Identification of AtALKBH1A and AtALKBH1D as DNA N6-adenine demethylases in Arabidopsis thaliana

DNA N6-methyladenine (6 mA) has recently been discovered as a novel DNA modification in animals and plants. In mammals, AlkB homolog 1 (ALKBH1) has been identified as a DNA 6 mA demethylase. ALKBH1 tightly controls the DNA 6 mA methylation level of mammalian genomes and plays important role in regulating gene expression. DNA 6 mA methylation has also been reported to exist in plant genomes, however, the plant DNA 6 mA demethylases and their function remain largely unknown. Here we identify homologs of ALKBH1 as DNA 6 mA demethylases in Arabidopsis. We discover that there are four homologs of ALKBH1, AtALKBH1A, AtALKBH1B, AtALKBH1C and AtALKBH1D, in Arabidopsis. In vitro enzymatic activity studies reveal that AtALKBH1A and 1D can efficiently erase DNA 6 mA methylation. Loss of function of AtALKBH1A and AtALKBH1D causes elevated DNA 6 mA methylation levels in vivo. atalkbh1a/1d mutant displays delayed seed gemination. Based on our RNA-seq data, we find some regulators of seed gemination are dysregulated in atalkbh1a/1d, and the dysregulation is correlated with changes of DNA 6 mA methylation levels. This study identifies plant DNA 6 mA demethylases and reports the function of DNA 6 mA methylation in regulating seed germination.

Population epigenetics: DNA methylation in the plant omics era

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

Molecular basis of chromatin remodelling by DDM1 involved in plant DNA methylation

Eukaryotic gene regulation occurs at the chromatin level, which requires changing the chromatin structure by a group of ATP-dependent DNA translocases-namely, the chromatin remodellers1. In plants, chromatin remodellers function in various biological processes and possess both conserved and plant-specific components2-5. DECREASE IN DNA METHYLATION 1 (DDM1) is a plant chromatin remodeller that plays a key role in the maintenance DNA methylation6-11. Here we determined the structures of Arabidopsis DDM1 in complex with nucleosome in ADP-BeFx-bound, ADP-bound and nucleotide-free conformations. We show that DDM1 specifically recognizes the H4 tail and nucleosomal DNA. The conformational differences between ADP-BeFx-bound, ADP-bound and nucleotide-free DDM1 suggest a chromatin remodelling cycle coupled to ATP binding, hydrolysis and ADP release. This, in turn, triggers conformational changes in the DDM1-bound nucleosomal DNA, which alters the nucleosome structure and promotes DNA sliding. Together, our data reveal the molecular basis of chromatin remodelling by DDM1.

Molecular mechanisms of the RNA polymerases in plant RNA-directed DNA methylation

In plants, two atypical DNA-dependent RNA polymerases, RNA polymerase IV (Pol IV) and Pol V, and an RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) together produce noncoding RNAs (ncRNAs) to guide the plant-specific RNA-directed DNA methylation (RdDM). Although both Pol IV and Pol V have evolved from the canonical Pol II, they have adapted to different roles in RdDM. The mechanisms of their adaptation are key to understanding plant DNA methylation and the divergent evolution of polymerases. In this review, we summarize insights that have emerged from recent structural studies of Pol IV, Pol V, and RDR2 and discuss their structural features critical for efficient ncRNA production in RdDM.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

Base Excision DNA Repair in Plants: Arabidopsis and Beyond

Base excision DNA repair (BER) is a key pathway safeguarding the genome of all living organisms from damage caused by both intrinsic and environmental factors. Most present knowledge about BER comes from studies of human cells, E. coli, and yeast. Plants may be under an even heavier DNA damage threat from abiotic stress, reactive oxygen species leaking from the photosynthetic system, and reactive secondary metabolites. In general, BER in plant species is similar to that in humans and model organisms, but several important details are specific to plants. Here, we review the current state of knowledge about BER in plants, with special attention paid to its unique features, such as the existence of active epigenetic demethylation based on the BER machinery, the unexplained diversity of alkylation damage repair enzymes, and the differences in the processing of abasic sites that appear either spontaneously or are generated as BER intermediates. Understanding the biochemistry of plant DNA repair, especially in species other than the Arabidopsis model, is important for future efforts to develop new crop varieties.

Molecular basis of the plant ROS1-mediated active DNA demethylation

Active DNA demethylation plays a crucial role in eukaryotic gene imprinting and antagonizing DNA methylation. The plant-specific REPRESSOR OF SILENCING 1/DEMETER (ROS1/DME) family of enzymes directly excise 5-methyl-cytosine (5mC), representing an efficient DNA demethylation pathway distinct from that of animals. Here, we report the cryo-electron microscopy structures of an Arabidopsis ROS1 catalytic fragment in complex with substrate DNA, mismatch DNA and reaction intermediate, respectively. The substrate 5mC is flipped-out from the DNA duplex and subsequently recognized by the ROS1 base-binding pocket through hydrophobic and hydrogen-bonding interactions towards the 5-methyl group and Watson-Crick edge respectively, while the different protonation states of the bases determine the substrate preference for 5mC over T:G mismatch. Together with the structure of the reaction intermediate complex, our structural and biochemical studies revealed the molecular basis for substrate specificity, as well as the reaction mechanism underlying 5mC demethylation by the ROS1/DME family of plant-specific DNA demethylases.

Factor of DNA methylation 1 affects woodland strawberry plant stature and organ size via DNA methylation

RNA-directed DNA methylation (RdDM) is an epigenetic process that directs silencing to specific genomic regions and loci. The biological functions of RdDM are not well studied in horticultural plants. Here, we isolated the ethyl methane-sulfonate-induced mutant reduced organ size (ros) producing small leaves, flowers, and fruits in woodland strawberry (Fragaria vesca) due to reduced cell numbers compared with that in the wild-type (WT). The candidate mutation causes a premature stop codon in FvH4_6g28780, which shares high similarity to Arabidopsis (Arabidopsis thaliana) Factor of DNA Methylation1 (FDM1) encoding an RdDM pathway component and was named FveFDM1. Consistently, the fvefdm1CR mutants generated by CRISPR/Cas9 also produced smaller organs. Overexpressing FveFDM1 in an Arabidopsis fdm1-1 fdm2-1 double mutant restored DNA methylation at the RdDM target loci. FveFDM1 acts in a protein complex with its homolog Involved in De Novo 2 (FveIDN2). Furthermore, whole-genome bisulfite sequencing revealed that DNA methylation, especially in the CHH context, was remarkably reduced throughout the genome in fvefdm1. Common and specific differentially expressed genes were identified in different tissues of fvefdm1 compared to in WT tissues. DNA methylation and expression levels of several gibberellic acid (GA) biosynthesis and cell cycle genes were validated. Moreover, the contents of GA and auxin were substantially reduced in the young leaves of fvefdm1 compared to in the WT. However, exogenous application of GA and auxin could not recover the organ size of fvefdm1. In addition, expression levels of FveFDM1, FveIDN2, Nuclear RNA Polymerase D1 (FveNRPD1), Domains Rearranged Methylase 2 (FveDRM2), and cell cycle genes were greatly induced by GA treatment. Overall, our work demonstrated the critical roles of FveFDM1 in plant growth and development via RdDM-mediated DNA methylation in horticultural crops.

A diRNA-protein scaffold module mediates SMC5/6 recruitment in plant DNA repair

In eukaryotes, the STRUCTURAL MAINTENANCE OF CHROMOSOME 5/6 (SMC5/6) complex is critical to maintaining chromosomal structures around double-strand breaks (DSBs) in DNA damage repair. However, the recruitment mechanism of this conserved complex at DSBs remains unclear. In this study, using Arabidopsis thaliana as a model, we found that SMC5/6 localization at DSBs is dependent on the protein scaffold containing INVOLVED IN DE NOVO 2 (IDN2), CELL DIVISION CYCLE 5 (CDC5), and ALTERATION/DEFICIENCY IN ACTIVATION 2B (ADA2b), whose recruitment is further mediated by DNA-damage-induced RNAs (diRNAs) generated from DNA regions around DSBs. The physical interactions of protein components including SMC5-ADA2b, ADA2b-CDC5, and CDC5-IDN2 result in formation of the protein scaffold. Further analysis indicated that the DSB localization of IDN2 requires its RNA-binding activity and ARGONAUTE 2 (AGO2), indicating a role for the AGO2-diRNA complex in this process. Given that most of the components in the scaffold are conserved, the mechanism presented here, which connects SMC5/6 recruitment and small RNAs, will improve our understanding of DNA repair mechanisms in eukaryotes.

The histone H3.1 variant regulates TONSOKU-mediated DNA repair during replication

The tail of replication-dependent histone H3.1 varies from that of replication-independent H3.3 at the amino acid located at position 31 in plants and animals, but no function has been assigned to this residue to demonstrate a unique and conserved role for H3.1 during replication. We found that TONSOKU (TSK/TONSL), which rescues broken replication forks, specifically interacts with H3.1 via recognition of alanine 31 by its tetratricopeptide repeat domain. Our results indicate that genomic instability in the absence of ATXR5/ATXR6-catalyzed histone H3 lysine 27 monomethylation in plants depends on H3.1, TSK, and DNA polymerase theta (Pol θ). This work reveals an H3.1-specific function during replication and a common strategy used in multicellular eukaryotes for regulating post-replicative chromatin maturation and TSK, which relies on histone monomethyltransferases and reading of the H3.1 variant.

The Arabidopsis ATR-SOG1 signaling module regulates pleiotropic developmental adjustments in response to 3'-blocked DNA repair intermediates

Base excision repair and active DNA demethylation produce repair intermediates with DNA molecules blocked at the 3'-OH end by an aldehyde or phosphate group. However, both the physiological consequences of these accumulated single-strand DNAs break with 3'-blocked ends (DNA 3'-blocks) and the signaling pathways responding to unrepaired DNA 3'-blocks remain unclear in plants. Here, we investigated the effects of DNA 3'-blocks on plant development using the zinc finger DNA 3'-phosphoesterase (zdp) AP endonuclease2 (ape2) double mutant, in which 3'-blocking residues are poorly repaired. The accumulation of DNA 3'-blocked triggered diverse developmental defects that were dependent on the ATM and RAD3-related (ATR)-suppressor of gamma response 1 (SOG1) signaling module. SOG1 mutation rescued the developmental defects of zdp ape2 leaves by preventing cell endoreplication and promoting cell proliferation. However, SOG1 mutation caused intensive meristematic cell death in the radicle of zdp ape2 following germination, resulting in rapid termination of radicle growth. Notably, mutating FORMAMIDOPYRIMIDINE DNA GLYCOSYLASE (FPG) in zdp ape2 sog1 partially recovered its radicle growth, demonstrating that DNA 3'-blocks generated by FPG caused the meristematic defects. Surprisingly, despite lacking a functional radicle, zdp ape2 sog1 mutants compensated the lack of root growth by generating anchor roots having low levels of DNA damage response. Our results reveal dual roles of SOG1 in regulating root establishment when seeds germinate with excess DNA 3'-blocks.

Optimization of T-DNA configuration with UBIQUITIN10 promoters and tRNA-sgRNA complexes promotes highly efficient genome editing in allotetraploid tobacco

KEY MESSAGE

Combination of UBIQUITIN10 promoter-directed CAS9 and tRNA-gRNA complexes in gene-editing assay induces 80% mutant phenotype with a knockout of the four allelic copies in the T0 generation of allotetraploid tobaccos. While gene-editing methodologies, such as CRISPR-Cas9, have been developed and successfully used in many plant species, their use remains challenging, because they most often rely on stable or transient transgene expression. Regrettably, in all plant species, transformation causes epigenetic effects such as gene silencing and variable transgene expression. Here, UBIQUITIN10 promoters from several plant species were characterized and showed their capacity to direct high levels of transgene expression in transient and stable transformation assays, which in turn was used to improve the selection process of regenerated transformants. Furthermore, we compared various sgRNAs delivery systems and showed that the combination of UBIQUITIN10 promoters and tRNA-sgRNA complexes produced 80% mutant phenotype with a complete knockout of the four allelic copies, while the remaining 20% exhibited weaker phenotype, which suggested partial allelic knockout, in the T0 generation of the allotetraploid Nicotiana tabacum. These data provide valuable information to optimize future designs of gene editing constructs for plant research and crop improvement and open the way for valuable gene editing projects in non-model Solanaceae species.

The transcription factor ORA59 exhibits dual DNA binding specificity that differentially regulates ethylene- and jasmonic acid-induced genes in plant immunity

Jasmonic acid (JA) and ethylene (ET) signaling modulate plant defense against necrotrophic pathogens in a synergistic and interdependent manner, while JA and ET also have independent roles in certain processes, e.g. in responses to wounding and flooding, respectively. These hormone pathways lead to transcriptional reprogramming, which is a major part of plant immunity and requires the roles of transcription factors. ET response factors are responsible for the transcriptional regulation of JA/ET-responsive defense genes, of which ORA59 functions as a key regulator of this process and has been implicated in the JA-ET crosstalk. We previously demonstrated that Arabidopsis (Arabidopsis thaliana) GDSL LIPASE 1 (GLIP1) depends on ET for gene expression and pathogen resistance. Here, promoter analysis of GLIP1 revealed ERELEE4 as the critical cis-element for ET-responsive GLIP1 expression. In a yeast one-hybrid screening, ORA59 was isolated as a specific transcription factor that binds to the ERELEE4 element, in addition to the well-characterized GCC box. We found that ORA59 regulates JA/ET-responsive genes through direct binding to these elements in gene promoters. Notably, ORA59 exhibited a differential preference for GCC box and ERELEE4, depending on whether ORA59 activation is achieved by JA and ET, respectively. JA and ET induced ORA59 phosphorylation, which was required for both activity and specificity of ORA59. Furthermore, RNA-seq and virus-induced gene silencing analyses led to the identification of ORA59 target genes of distinct functional categories in JA and ET pathways. Our results provide insights into how ORA59 can generate specific patterns of gene expression dynamics through JA and ET hormone pathways.

Changes in Epigenetic Patterns Related to DNA Replication in Vicia faba Root Meristem Cells under Cadmium-Induced Stress Conditions

Experiments on Vicia faba root meristem cells exposed to 150 µM cadmium chloride (CdCl2) were undertaken to analyse epigenetic changes, mainly with respect to DNA replication stress. Histone modifications examined by means of immunofluorescence labeling included: (1) acetylation of histone H3 on lysine 56 (H3K56Ac), involved in transcription, S phase, and response to DNA damage during DNA biosynthesis; (2) dimethylation of histone H3 on lysine 79 (H3K79Me2), correlated with the replication initiation; (3) phosphorylation of histone H3 on threonine 45 (H3T45Ph), engaged in DNA synthesis and apoptosis. Moreover, immunostaining using specific antibodies against 5-MetC-modified DNA was used to determine the level of DNA methylation. A significant decrease in the level of H3K79Me2, noted in all phases of the CdCl2-treated interphase cell nuclei, was found to correspond with: (1) an increase in the mean number of intranuclear foci of H3K56Ac histones (observed mainly in S-phase), (2) a plethora of nuclear and nucleolar labeling patterns (combined with a general decrease in H3T45Ph), and (3) a decrease in DNA methylation. All these changes correlate well with a general viewpoint that DNA modifications and post-translational histone modifications play an important role in gene expression and plant development under cadmium-induced stress conditions.

Mechanistic and genetic basis of single-strand templated repair at Cas12a-induced DNA breaks in Chlamydomonas reinhardtii

Single-stranded oligodeoxynucleotides (ssODNs) are widely used as DNA repair templates in CRISPR/Cas precision genome editing. However, the underlying mechanisms of single-strand templated DNA repair (SSTR) are inadequately understood, constraining rational improvements to precision editing. Here we study SSTR at CRISPR/Cas12a-induced DNA double-strand breaks (DSBs) in the eukaryotic model green microalga Chlamydomonas reinhardtii. We demonstrate that ssODNs physically incorporate into the genome during SSTR at Cas12a-induced DSBs. This process is genetically independent of the Rad51-dependent homologous recombination and Fanconi anemia pathways, is strongly antagonized by non-homologous end-joining, and is mediated almost entirely by the alternative end-joining enzyme polymerase θ. These findings suggest differences in SSTR between C. reinhardtii and animals. Our work illustrates the promising potentially of C. reinhardtii as a model organism for studying nuclear DNA repair.

Kleisin NSE4 of the SMC5/6 complex is necessary for DNA double strand break repair, but not for recovery from DNA damage in Physcomitrella (Physcomitrium patens)

Kleisin NSE4 and circular form of SMC5/6 is indispensable for DSB repair and necessary for gene targeting but is not enough for recovery of cells from DNA damage in Physcomitrella. Structural maintenance of chromosomes (SMC) complexes are involved in cohesion, condensation and maintenance of genome stability. Based on the sensitivity of mutants to genotoxic stress the SMC5/6 complex is thought to play a prominent role in DNA stabilization during repair by tethering DNA at the site of lesion by a heteroduplex of SMC5 and SMC6 encircled with non-SMC components NSE1, NSE3 and kleisin NSE4. In this study, we tested how formation of the SMC5/6 circular structure affects mutant sensitivity to DNA damage, kinetics of DSB repair and gene targeting. In the moss Physcomitrella (Physcomitrium patens), SMC6 and NSE4 are essential single copy genes and this is why we used blocking of transcription to reveal their mutated phenotype. Even slight reduction of transcript levels by dCas9 binding was enough to obtain stable lines with severe DSB repair defects and specific bleomycin sensitivity. We show that survival after bleomycin or MMS treatment fully depends on active SMC6, whereas attenuation of NSE4 has little or negligible effect. We conclude that circularization of SMC5/6 provided by the kleisin NSE4 is indispensable for the DSB repair, nevertheless there are other functions associated with the SMC5/6 complex, which are critical to survive DNA damage.

Genome-wide detection of cytosine methylations in plant from Nanopore data using deep learning

In plants, cytosine DNA methylations (5mCs) can happen in three sequence contexts as CpG, CHG, and CHH (where H = A, C, or T), which play different roles in the regulation of biological processes. Although long Nanopore reads are advantageous in the detection of 5mCs comparing to short-read bisulfite sequencing, existing methods can only detect 5mCs in the CpG context, which limits their application in plants. Here, we develop DeepSignal-plant, a deep learning tool to detect genome-wide 5mCs of all three contexts in plants from Nanopore reads. We sequence Arabidopsis thaliana and Oryza sativa using both Nanopore and bisulfite sequencing. We develop a denoising process for training models, which enables DeepSignal-plant to achieve high correlations with bisulfite sequencing for 5mC detection in all three contexts. Furthermore, DeepSignal-plant can profile more 5mC sites, which will help to provide a more complete understanding of epigenetic mechanisms of different biological processes.

HvLEAFY controls the early stages of floral organ specification and inhibits the formation of multiple ovaries in barley

Transition to the reproductive phase, inflorescence formation and flower development are crucial elements that ensure maximum reproductive success in a plant's life cycle. To understand the regulatory mechanisms underlying correct flower development in barley (Hordeum vulgare), we characterized the multiovary 5 (mov5.o) mutant. This mutant develops abnormal flowers that exhibit mosaic floral organs typified by multiple carpels at the total or partial expense of stamens. Genetic mapping positioned mov5 on the long arm of chromosome 2H, incorporating a region that encodes HvLFY, the barley orthologue of LEAFY from Arabidopsis. Sequencing revealed that, in mov5.o plants, HvLFY contains a single amino acid substitution in a highly conserved proline residue. CRISPR-mediated knockout of HvLFY replicated the mov5.o phenotype, suggesting that HvLFY_mov5 represents a loss of function allele. In heterologous assays, the HvLFY_mov5 polymorphism influenced protein-protein interactions and affinity for a putative binding site in the promoter of HvMADS58, a C-class MADS-box gene. Moreover, molecular analysis indicated that HvLFY interacts with HvUFO and regulates the expression of floral homeotic genes including HvMADS2, HvMADS4 and HvMADS16. Other distinct changes in expression differ from those reported in the rice LFY mutants apo2/rfl, suggesting that LFY function in the grasses is modulated in a species-specific manner. This pathway provides a key entry point for the study of LFY function and multiple ovary formation in barley, as well as cereal species in general.

Development and evaluation of a novel loop mediated isothermal amplification coupled with TaqMan probe assay for detection of genetically modified organism with NOS terminator

In this study, we aim to develop a novel loop mediated isothermal amplification (LAMP) coupled with TaqMan (LAMP-TaqMan) method for quick qualitative detection of genetically modified organism (GMOs). We designed four LAMP primers and one TaqMan probe for the LAMP-TaqMan detection method to detect the nopaline synthase gene (NOS) terminator in GMOs. This assay enabled the amplification of DNA within ~20 min at a constant temperature of 65 °C. This assay detected as few as five copies of target sequences, which had a high specificity similar to the TaqMan qPCR method. Furthermore, the LAMP-TaqMan detection method was successfully used to amplify and detect DNA from food samples of the major crops (soybean, maize, rice, etc.). In summary, a novel LAMP-TaqMan assay has been developed, which has the similar sensitivity but takes less time than the TaqMan qPCR method. This method offers a novel approach for rapid detection of GMOs in foods.

DNA demethylases remodel DNA methylation in rice gametes and zygote and are required for reproduction

Fertilization constitutes a critical step in the plant life cycle during which the gamete genomes undergo chromatin dynamics in preparation for embryogenesis. In mammals, parental chromatin is extensively reprogrammed through the global erasure of DNA methylation. However, in flowering plants it remains unclear whether and how DNA methylation is remodeled in gametes and after fertilization in the zygote. In this study, we characterize DNA methylation patterns and investigate the function of DNA glycosylases in rice eggs, sperm, and unicellular zygotes and during embryogenesis. We found that DNA methylation is locally reconfigured after fertilization and is intensified during embryogenesis. Genetic, epigenomic, and transcriptomic analysis revealed that three rice DNA glycosylases, DNG702, DNG701, and DNG704, demethylate DNA at distinct genomic regions in the gametes and the zygote, and are required for zygotic gene expression and development. Collectively, these results indicate that active DNA demethylation takes place in the gametes and the zygote to locally remodel DNA methylation, which is critical for egg and zygote gene expression and reproduction in rice.

Connections between the Cell Cycle and the DNA Damage Response in Plants

Due to their sessile lifestyle, plants are especially exposed to various stresses, including genotoxic stress, which results in altered genome integrity. Upon the detection of DNA damage, distinct cellular responses lead to cell cycle arrest and the induction of DNA repair mechanisms. Interestingly, it has been shown that some cell cycle regulators are not only required for meristem activity and plant development but are also key to cope with the occurrence of DNA lesions. In this review, we first summarize some important regulatory steps of the plant cell cycle and present a brief overview of the DNA damage response (DDR) mechanisms. Then, the role played by some cell cycle regulators at the interface between the cell cycle and DNA damage responses is discussed more specifically.

The plastome sequence of Bactris gasipaes and evolutionary analysis in tribe Cocoseae (Arecaceae)

The family Arecaceae is distributed throughout tropical and subtropical regions of the world. Among the five subfamilies, Arecoideae is the most species-rich and still contains some ambiguous inter-generic relationships, such as those within subtribes Attaleinae and Bactridineae. The hypervariable regions of plastid genomes (plastomes) are interesting tools to clarify unresolved phylogenetic relationships. We sequenced and characterized the plastome of Bactris gasipaes (Bactridinae) and compared it with eight species from the three Cocoseae sub-tribes (Attaleinae, Bactridinae, and Elaeidinae) to perform comparative analysis and to identify hypervariable regions. The Bactris gasipaes plastome has 156,646 bp, with 113 unique genes. Among them, four genes have an alternative start codon (cemA, rps19, rpl2, and ndhD). Plastomes are highly conserved within tribe Cocoseae: 97.3% identity, length variation of ~2 kb, and a single ~4.5 kb inversion in Astrocaryum plastomes. The LSC/IR and IR/SSC junctions vary among the subtribes: in Bactridinae and Elaeidinae the rps19 gene is completely contained in the IR region; in the subtribe Attaleinae the rps19 gene is only partially contained in the IRs. The hypervariable regions selected according to sequence variation (SV%) and frequency of parsimony informative sites (PIS%) revealed plastome regions with great potential for molecular analysis. The ten regions with greatest SV% showed higher variation than the plastid molecular markers commonly used for phylogenetic analysis in palms. The phylogenetic trees based on the plastomes and the hypervariable regions (SV%) datasets had well-resolved relationships, with consistent topologies within tribe Cocoseae, and confirm the monophyly of the subtribes Bactridinae and Attaleinae.

RNase H1C collaborates with ssDNA binding proteins WHY1/3 and recombinase RecA1 to fulfill the DNA damage repair in Arabidopsis chloroplasts

Proper repair of damaged DNA is crucial for genetic integrity and organismal survival. As semi-autonomous organelles, plastids have their own genomes whose integrity must be preserved. Several factors have been shown to participate in plastid DNA damage repair; however, the underlying mechanism remains unclear. Here, we elucidate a mechanism of homologous recombination (HR) repair in chloroplasts that involves R-loops. We find that the recombinase RecA1 forms filaments in chloroplasts during HR repair, but aggregates as puncta when RNA:DNA hybrids accumulate. ssDNA-binding proteins WHY1/3 and chloroplast RNase H1 AtRNH1C are recruited to the same genomic sites to promote HR repair. Depletion of AtRNH1C or WHY1/3 significantly suppresses the binding of RNA polymerase to the damaged DNA, thus reducing HR repair and modulating microhomology-mediated double-strand break repair. Furthermore, we show that DNA polymerase IB works with AtRNH1C genetically to complete the DNA damage repair process. This study reveals the positive role of R-loops in facilitating the activities of WHY1/3 and RecA1, which in turn secures HR repair and organellar development.

Redox feedback regulation of ANAC089 signaling alters seed germination and stress response

The interplay between the phytohormone abscisic acid (ABA) and the gasotransmitter nitric oxide (NO) regulates seed germination and post-germinative seedling growth. We show that GAP1 (germination in ABA and cPTIO 1) encodes the transcription factor ANAC089 with a critical membrane-bound domain and extranuclear localization. ANAC089 mutants lacking the membrane-tethered domain display insensitivity to ABA, salt, and osmotic and cold stresses, revealing a repressor function. Whole-genome transcriptional profiling and DNA-binding specificity reveals that ANAC089 regulates ABA- and redox-related genes. ANAC089 truncated mutants exhibit higher NO and lower ROS and ABA endogenous levels, alongside an altered thiol and disulfide homeostasis. Consistently, translocation of ANAC089 to the nucleus is directed by changes in cellular redox status after treatments with NO scavengers and redox-related compounds. Our results reveal ANAC089 to be a master regulator modulating redox homeostasis and NO levels, able to repress ABA synthesis and signaling during Arabidopsis seed germination and abiotic stress.

SOG1 transcription factor promotes the onset of endoreduplication under salinity stress in Arabidopsis

As like in mammalian system, the DNA damage responsive cell cycle checkpoint functions play crucial role for maintenance of genome stability in plants through repairing of damages in DNA and induction of programmed cell death or endoreduplication by extensive regulation of progression of cell cycle. ATM and ATR (ATAXIA-TELANGIECTASIA-MUTATED and -RAD3-RELATED) function as sensor kinases and play key role in the transmission of DNA damage signals to the downstream components of cell cycle regulatory network. The plant-specific NAC domain family transcription factor SOG1 (SUPPRESSOR OF GAMMA RESPONSE 1) plays crucial role in transducing signals from both ATM and ATR in presence of double strand breaks (DSBs) in the genome and found to play crucial role in the regulation of key genes involved in cell cycle progression, DNA damage repair, endoreduplication and programmed cell death. Here we report that Arabidopsis exposed to high salinity shows generation of oxidative stress induced DSBs along with the concomitant induction of endoreduplication, displaying increased cell size and DNA ploidy level without any change in chromosome number. These responses were significantly prominent in SOG1 overexpression line than wild-type Arabidopsis, while sog1 mutant lines showed much compromised induction of endoreduplication under salinity stress. We have found that both ATM-SOG1 and ATR-SOG1 pathways are involved in the salinity mediated induction of endoreduplication. SOG1was found to promote G2-M phase arrest in Arabidopsis under salinity stress by downregulating the expression of the key cell cycle regulators, including CDKB1;1, CDKB2;1, and CYCB1;1, while upregulating the expression of WEE1 kinase, CCS52A and E2Fa, which act as important regulators for induction of endoreduplication. Our results suggest that Arabidopsis undergoes endoreduplicative cycle in response to salinity induced DSBs, showcasing an adaptive response in plants under salinity stress.

Maize Endochitinase Expression in Response to Fall Armyworm Herbivory

A large percentage of crop loss is due to insect damage, especially caterpillar damage. Plant chitinases are considered excellent candidates to combat these insects since they can degrade chitin in peritrophic matrix (PM), an important protective structure in caterpillar midgut. Compared to chemical insecticides, chitinases could improve host plant resistance and be both economically and environmentally advantageous. The focus of this research was to find chitinase candidates that could improve plant resistance by effectively limiting caterpillar damage. Five classes of endochitinase (I-V) genes were characterized in the maize genome, and we isolated and cloned four chitinase genes (chitinase A, chitinase B, chitinase I, and PRm3) present in two maize (Zea mays L.) inbred lines Mp708 and Tx601, with different levels of resistance to caterpillar pests. We also investigated the expression of these maize chitinases in response to fall armyworm (Spodoptera frugiperda, FAW) attack. The results indicated that both chitinase transcript abundance and enzymatic activity increased in response to FAW feeding and mechanical wounding. Furthermore, chitinases retained activity inside the caterpillar midgut and enzymatic activity was detected in the food bolus and frass. When examined under scanning electron microscopy, PMs from Tx601-fed caterpillars showed structural damage when compared to diet controls. Analysis of chitinase transcript abundance after caterpillar feeding and proteomic analysis of maize leaf trichomes in the two inbreds implicated chitinase PRm3 found in Tx601 as a potential insecticidal protein.

A Dual-Color Tyr-FISH Method for Visualizing Genes/Markers on Plant Chromosomes to Create Integrated Genetic and Cytogenetic Maps

In situ imaging of molecular markers on a physical chromosome is an indispensable tool for refining genetic maps and validation genome assembly at the chromosomal level. Despite the tremendous progress in genome sequencing, the plant genome assembly at the chromosome level remains a challenge. Recently developed optical and Hi-C mapping are aimed at assistance in genome assembly. For high confidence in the genome assembly at chromosome level, more independent approaches are required. The present study is aimed at refining an ultrasensitive Tyr-FISH technique and developing a reliable and simple method of in situ mapping of a short unique DNA sequences on plant chromosomes. We have carefully analyzed the critical steps of the Tyr-FISH to find out the reasons behind the flaws of this technique. The accurate visualization of markers/genes appeared to be significantly dependent on the means of chromosome slide preparation, probe design and labeling, and high stringency washing. Appropriate adjustment of these steps allowed us to detect a short DNA sequence of 1.6 Kb with a frequency of 51.6%. Based on our results, we developed a more reliable and simple protocol for dual-color Tyr-FISH visualization of unique short DNA sequences on plant chromosomes. This new protocol can allow for more accurate determination of the physical distance between markers and can be applied for faster integration of genetic and cytogenetic maps.

Physiological and genetic characterization of heat stress effects in a common bean RIL population

Heat stress is a major abiotic stress factor reducing crop productivity and climate change models predict increasing temperatures in many production regions. Common bean (Phaseolus vulgaris L.) is an important crop for food security in the tropics and heat stress is expected to cause increasing yield losses. To study physiological responses and to characterize the genetics of heat stress tolerance, we evaluated the recombinant inbred line (RIL) population IJR (Indeterminate Jamaica Red) x AFR298 of the Andean gene pool. Heat stress (HS) conditions in the field affected many traits across the reproductive phase. High nighttime temperatures appeared to have larger effects than maximum daytime temperatures. Yield was reduced compared to non-stress conditions by 37% and 26% in 2016 and 2017 seasons, respectively. The image analysis tool HYRBEAN was developed to evaluate pollen viability (PolVia). A significant reduction of PolVia was observed in HS and higher viability was correlated with yield only under stress conditions. In susceptible lines the reproductive phase was extended and defects in the initiation of seed, seed fill and seed formation were identified reducing grain quality. Higher yields under HS were correlated with early flowering, high pollen viability and effective seed filling. Quantitative trait loci (QTL) analysis revealed a QTL for both pod harvest index and PolVia on chromosome Pv05, for which the more heat tolerant parent IJR contributed the positive allele. Also, on chromosome Pv08 a QTL from IJR improved PolVia and the yield component pods per plant. HS affected several traits during the whole reproductive development, from floral induction to grain quality traits, indicating a general heat perception affecting many reproductive processes. Identification of tolerant germplasm, indicator traits for heat tolerance and molecular tools will help to breed heat tolerant varieties to face future climate change effects.

EIN2-directed histone acetylation requires EIN3-mediated positive feedback regulation in response to ethylene

Ethylene is an important phytohormone with pleotropic roles in plant growth, development, and stress responses. ETHYLENE INSENSITIVE2 (EIN2) mediates the transduction of the ethylene signal from the endoplasmic reticulum membrane to the nucleus, where its C-terminus (EIN2-C) regulates histone acetylation to mediate transcriptional regulation by EIN3. However, no direct interaction between EIN2-C and EIN3 has been detected. To determine how EIN2-C and EIN3 act together, we followed a synthetic approach and engineered a chimeric EIN2-C with EIN3 DNA-binding activity but lacking its transactivation activity (EIN2C-EIN3DB). The overexpression of EIN2C-EIN3DB in either wild-type or in the ethylene-insensitive mutant ein3-1 eil1-1 led to a partial constitutive ethylene response. Chromatin immunoprecipitation sequencing showed that EIN2C-EIN3DB has DNA-binding activity, indicating that EIN3DB is functional in EIN2C-EIN3DB. Furthermore, native EIN3 protein levels determine EIN2C-EIN3DB binding activity and binding targets in a positive feedback loop by interacting with EIN2C-EIN3DB to form a heterodimer. Additionally, although EIN3 does not direct affect histone acetylation levels in the absence of EIN2, it is required for the ethylene-induced elevation of H3K14Ac and H3K23Ac in the presence of EIN2. Together, we reveal efficient and specific DNA-binding by dimerized EIN3 in the presence of ethylene to mediate positive feedback regulation, which is required for EIN2-directed elevation of histone acetylation to integrate into an EIN3-dependent transcriptional activation.

Identification of regulatory factors promoting embryogenic callus formation in barley through transcriptome analysis

BACKGROUND

Barley is known to be recalcitrant to tissue culture, which hinders genetic transformation and its biotechnological application. To date, the ideal explant for transformation remains limited to immature embryos; the mechanism underlying embryonic callus formation is elusive.

RESULTS

This study aimed to uncover the different transcription regulation pathways between calli formed from immature (IME) and mature (ME) embryos through transcriptome sequencing. We showed that incubation of embryos in an auxin-rich medium caused dramatic changes in gene expression profiles within 48 h. Overall, 9330 and 11,318 differentially expressed genes (DEGs) were found in the IME and ME systems, respectively. 3880 DEGs were found to be specific to IME_0h/IME_48h, and protein phosphorylation, regulation of transcription, and oxidative-reduction processes were the most common gene ontology categories of this group. Twenty-three IAA, fourteen ARF, eight SAUR, three YUC, and four PIN genes were found to be differentially expressed during callus formation. The effect of callus-inducing medium (CIM) on IAA genes was broader in the IME system than in the ME system, indicating that auxin response participates in regulating cell reprogramming during callus formation. BBM, LEC1, and PLT2 exhibited a significant increase in expression levels in the IME system but were not activated in the ME system. WUS showed a more substantial growth trend in the IME system than in the ME system, suggesting that these embryonic, shoot, and root meristem genes play crucial roles in determining the acquisition of competency. Moreover, epigenetic regulators, including SUVH3A, SUVH2A, and HDA19B/703, exhibited differential expression patterns between the two induction systems, indicating that epigenetic reprogramming might contribute to gene expression activation/suppression in this process. Furthermore, we examined the effect of ectopic expression of HvBBM and HvWUS on Agrobacterium-mediated barley transformation. The transformation efficiency in the group expressing the PLTPpro:HvBBM + Axig1pro:HvWUS construct was increased by three times that in the control (empty vector) because of enhanced plant regeneration capacity.

CONCLUSIONS

We identified some regulatory factors that might contribute to the differential responses of the two explants to callus induction and provide a promising strategy to improve transformation efficiency in barley.

N uptake, assimilation and isotopic fractioning control δ 15N dynamics in plant DNA: A heavy labelling experiment on Brassica napus L

In last decades, a large body of evidence clarified nitrogen isotope composition (δ¹⁵N) patterns in plant leaves, roots and metabolites, showing isotopic fractionation along N uptake and assimilation pathways, in relation to N source and use efficiency, also suggesting ¹⁵N depletion in plant DNA. Here we present a manipulative experiment on Brassica napus var. oleracea, where we monitored δ ¹⁵N of purified, lyophilized DNA and source leaf and root materials, over a 60-days growth period starting at d 60 after germination, in plants initially supplied with a heavy labelled (δ ¹⁵N_Air-N2 = 2100 mUr) ammonium nitrate solution covering nutrient requirements for the whole observation period (470 mg N per plant) and controlling for the labelled N species (ṄH₄, ṄO₃ and both). Dynamics of Isotopic Ratio Mass Spectrometry (IRMS) data for the three treatments showed that: (1) leaf and root δ ¹⁵N dynamics strictly depend on the labelled chemical species, with ṄH₄, ṄO₃ and ṄH₄ṄO₃ plants initially showing higher, lower and intermediate values, respectively, then converging due to the progressive NH₄⁺ depletion from the nutrient solution; (2) in ṄH₄ṄO₃, where δ¹⁵N was not affected by the labelled chemical species, we did not observe isotopic fractionation associated to inorganic N uptake; (3) δ¹⁵N values in roots compared to leaves did not fully support patterns predicted by differences in assimilation rates of NH₄⁺ and NO₃^-; (4) DNA is depleted in ¹⁵N compared to the total N pools of roots and leaves, likely due to enzymatic discrimination during purine biosynthesis. In conclusion, while our experimental setup did not allow to assess the fractionation coefficient (ε) associated to DNA bases biosynthesis, this is the first study specifically reporting on dynamics of specific plant molecular pools such as nucleic acids over a long observation period with a heavy labelling technique.

The Effects of Cold Atmospheric Pressure Plasma on Germination Parameters, Enzyme Activities and Induction of DNA Damage in Barley

Climate change, environmental pollution and pathogen resistance to available chemical agents are part of the problems that the food industry has to face in order to ensure healthy food for people and livestock. One of the promising solutions to these problems is the use of cold atmospheric pressure plasma (CAPP). Plasma is suitable for efficient surface decontamination of seeds and food products, germination enhancement and obtaining higher yields in agricultural production. However, the plasma effects vary due to plasma source, treatment conditions and seed type. In our study, we tried to find the proper conditions for treatment of barley grains by diffuse coplanar surface barrier discharge, in which positive effects of CAPP, such as enhanced germination or decontamination effects, would be maximized and harmful effects, such as oxidation and genotoxic potential, minimized. Besides germination parameters, we evaluated DNA damage and activities of various germination and antioxidant enzymes in barley seedlings. Plasma exposure resulted in changes in germination parameters and enzyme activities. Longer exposures had also genotoxic effects. As such, our findings indicate that appropriate plasma exposure conditions need to be carefully optimized in order to preserve germination, oxidation balance and genome stability, should CAPP be used in agricultural practice.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq

We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

Nucleolar rDNA folds into condensed foci with a specific combination of epigenetic marks

Arabidopsis thaliana 45S ribosomal genes (rDNA) are located in tandem arrays called nucleolus organizing regions on the termini of chromosomes 2 and 4 (NOR2 and NOR4) and encode rRNA, a crucial structural element of the ribosome. The current model of rDNA organization suggests that inactive rRNA genes accumulate in the condensed chromocenters in the nucleus and at the nucleolar periphery, while the nucleolus delineates active genes. We challenge the perspective that all intranucleolar rDNA is active by showing that a subset of nucleolar rDNA assembles into condensed foci marked by H3.1 and H3.3 histones that also contain the repressive H3K9me2 histone mark. By using plant lines containing a low number of rDNA copies, we further found that the condensed foci relate to the folding of rDNA, which appears to be a common mechanism of rDNA regulation inside the nucleolus. The H3K9me2 histone mark found in condensed foci represents a typical modification of bulk inactive rDNA, as we show by genome-wide approaches, similar to the H2A.W histone variant. The euchromatin histone marks H3K27me3 and H3K4me3, in contrast, do not colocalize with nucleolar foci and their overall levels in the nucleolus are very low. We further demonstrate that the rDNA promoter is an important regulatory region of the rDNA, where the distribution of histone variants and histone modifications are modulated in response to rDNA activity.

Looking for Hidden Enemies of Metabarcoding: Species Composition, Habitat and Management Can Strongly Influence DNA Extraction while Examining Grassland Communities

Despite the raising preoccupation, the critical question of how the plant community is composed belowground still remains unresolved, particularly for the conservation priority types of vegetation. The usefulness of metabarcoding analysis of the belowground parts of the plant community is subjected to a considerable bias, that often impedes detection of all species in a sample due to insufficient DNA quality or quantity. In the presented study we have attempted to find environmental factors that determine the amount and quality of DNA extracted from total plant tissue from above- and belowground samples (1000 and 10,000 cm2). We analyzed the influence of land use intensity, soil properties, species composition, and season on DNA extraction. The most important factors for DNA quality were vegetation type, soil conductometry (EC), and soil pH for the belowground samples. The species that significantly decreased the DNA quality were Calamagrostis epigejos, Coronilla varia, and Holcus lanatus. For the aboveground part of the vegetation, the season, management intensity, and certain species-with the most prominent being Centaurea rhenana and Cirsium canum-have the highest influence. Additionally, we found that sample size, soil granulation, MgO, organic C, K2O, and total soil N content are important for DNA extraction effectiveness. Both low EC and pH reduce significantly the yield and quality of DNA. Identifying the potential inhibitors of DNA isolation and predicting difficulties of sampling the vegetation plots for metabarcoding analysis will help to optimize the universal, low-cost multi-stage DNA extraction procedure in molecular ecology studies.

Deep6mA: A deep learning framework for exploring similar patterns in DNA N6-methyladenine sites across different species

N6-methyladenine (6mA) is an important DNA modification form associated with a wide range of biological processes. Identifying accurately 6mA sites on a genomic scale is crucial for under-standing of 6mA’s biological functions. However, the existing experimental techniques for detecting 6mA sites are cost-ineffective, which implies the great need of developing new computational methods for this problem. In this paper, we developed, without requiring any prior knowledge of 6mA and manually crafted sequence features, a deep learning framework named Deep6mA to identify DNA 6mA sites, and its performance is superior to other DNA 6mA prediction tools. Specifically, the 5-fold cross-validation on a benchmark dataset of rice gives the sensitivity and specificity of Deep6mA as 92.96% and 95.06%, respectively, and the overall prediction accuracy is 94%. Importantly, we find that the sequences with 6mA sites share similar patterns across different species. The model trained with rice data predicts well the 6mA sites of other three species: Arabidopsis thaliana, Fragaria vesca and Rosa chinensis with a prediction accuracy over 90%. In addition, we find that (1) 6mA tends to occur at GAGG motifs, which means the sequence near the 6mA site may be conservative; (2) 6mA is enriched in the TATA box of the promoter, which may be the main source of its regulating downstream gene expression.

DNA N6 methyladenine (6mA) is a newly recognized methylation modification in eukaryotes. It exists widely and conservatively in organisms, and its modification level changes dynamically in the whole life cycle. This study proposes an algorithm based on a deep learning framework including LSTM and CNN to predict 6mA sites. The results showed that our method could accurately predict the 6mA sites in different species, which means DNA sub-sequences containing 6mA sites among species have certain conservation. Importantly, we found that 6mA methylation in most different species is more likely to occur on the GAGG motif. In addition, we also found that 6mA is rich in the promoter’s TATA box, which may be a mechanism of regulating downstream gene expression.

Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae)

As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.

Distribution and degradation of DNA from non-genetically and genetically modified soybean (Roundup Ready): Impact of soybean protein concentrate and soybean protein isolate preparation

To improve genetically modified product labelling legislation and promote the development of genetically modified foods, the mass variations of genomic DNA and length distributions of DNA fragments in non-genetically and genetically modified soybean (Roundup Ready) and the variations in transgenic contents during soybean protein concentrate (SPC) and soybean protein isolate (SPI) preparation were monitored. The material masses throughout the process conformed to the law of mass conservation, and amounts of DNA recovered decreased with SPC and SPI preparation. The successive steps of ethanol extraction led to a decrease in the size distribution of the recovered DNA. For the LECTIN and CP4 EPSPS targets investigated, longer fragments (>800 bp) were more affected than smaller fragments (<200 bp), and both targets degraded slowly upon degradation into small fragments. DNA distribution and degradation thereby affected GMO quantification. After preparation procedures, the transgenic contents of SPC and SPI products were higher than that of raw soybean.

Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia

Subgenus Perovskia of the extended genus of Salvia comprises several Central Asian medicinal and aromatic species, of which S. yangii and S. abrotanoides are the most widespread. These plants are cultivated in Europe as robust ornamentals, and several cultivars are available. However, their medicinal potential remains underutilized because of limited information about their phytochemical and genetic diversity. Thus, we combined an ultra-high performance liquid chromatography quadrupole time of flight mass spectrometry (UHPLC-QTOF-MS) based metabolomics with DNA barcoding approach based on trnH-psbA and ITS2 barcodes to clarify the relationships between these two taxa. Metabolomic analysis demonstrated that aerial parts are more similar than roots and none of the major compounds stand out as distinct. Sugiol in S. yangii leaves and carnosic acid quinone in S. abrotanoides were mostly responsible for their chemical differentiation, whereas in roots the distinction was supported by the presence of five norditerpenoids in S. yangii and two flavonoids and one norditerpenoid in S. abrotanoides. To verify the metabolomics-based differentiation, we performed DNA authentication that revealed S. yangii and S. abrotanoides to be very closely related but separate species. We demonstrated that DNA barcoding coupled with parallel LC-MS profiling constitutes a powerful tool in identification of taxonomically close Salvia species.

New Insights into Plant Extracellular DNA. A Study in Soybean Root Extracellular Trap

exDNA is found in various organisms, including plants. However, plant exDNA has thus far received little attention related to its origin and role in the RET (root extracellular trap). In this study, we performed the first high-throughput genomic sequencing of plant exDNA from a Fabaceae with worldwide interest: soybean (Glycine max (L.) Merr.). The origin of this exDNA was first investigated in control condition, and the results show high-coverage on organelles (mitochondria/plastid) DNA relative to nuclear DNA, as well as a mix of coding and non-coding sequences. In the second part of this study, we investigated if exDNA release was modified during an elicitation with PEP-13 (a peptide elicitor from oomycete genus Phytophthora). Our results show that treatment of roots with PEP-13 does not affect the composition of exDNA.

Unusual DNA-binding properties of the Arabidopsis thaliana WRKY50 transcription factor at target gene promoters

WRKY50 from A. thaliana requires WT-boxes at target gene promoters for activation and binding. Based on the genome-wide prediction of WRKY50 target genes and the similarity of a WRKY50 binding site to WT-boxes in microbe-associated molecular pattern (MAMP)-responsive cis-regulatory modules (CRM), four WT-box containing CRMs from the promoter region of three WRKY50 target genes were investigated for their interaction with WRKY50. These target genes are DJ1E, WRKY30 and ATBBE4. Two of the four CRMs, one from DJ1E and one from WRKY30, were able to activate reporter gene expression in the presence of WRKY50. Activation requires the WT-boxes GGACTTTT, GGACTTTG from DJ1E and GGACTTTC from WRKY30. WRKY50 does not activate a second CRM from WRKY30 and the CRM from ATBBE4, both containing the WT-box TGACTTTT. In vitro gel-shift assays demonstrate WT-box-specific binding of the WRKY50 DNA-binding domain to all four CRMs. This work shows a high flexibility of WRKY50 binding site recognition beyond the classic W-box TTGACC/T.

Profiling Protein-DNA Interactions by Chromatin Immunoprecipitation in Arabidopsis

In plant cells, transcription factors play an important role in the regulation of gene expression, which eventually leads to the formation of complex phenotypes. Although chromatin immunoprecipitation (ChIP) involves a lengthy process that requires up to 4 days to complete, it is a powerful technique to investigate the interactions between transcription factors and their target sequences in vivo. Here, we describe a detailed ChIP protocol, focusing on ChIP-qPCR, from material collection to data analyses. Moreover, we explain multiple checkpoints for the quality control of ChIP-qPCR data to ensure the success of this protocol. As this protocol is robust, it can be adapted to other plant materials and plant species, and it can be used for genome-wide profiling experiments, including ChIP-chip and ChIP-seq analyses. We believe that our ChIP-qPCR protocol facilitates research on the interactions between plant transcription factors and their target sequences in vivo.

A moso bamboo transcription factor, Phehdz1, positively regulates the drought stress response of transgenic rice

78 HD-Zip family genes in Phyllostachys edulis were analyzed. Overexpression of Phehdz1 can improve the drought tolerance of transgenic rice and affect its secondary metabolism. Many studies suggested homeodomain-leucine zipper (HD-Zip) transcription factors are important regulators of plant growth and development, signal transduction, and responses to environmental stresses. In this study, 78 moso bamboo (Phyllostachys edulis) HD-Zip genes were investigated and classified into four subfamilies (HD-Zip I-IV). Additionally, Phehdz1 (HD-Zip I gene) was isolated and confirmed to be highly expressed in the roots. A quantitative real-time PCR analysis indicated Phehdz1 expression was significantly induced by drought, high salinity, and abscisic acid (ABA). A transient expression assay proved that Phehdz1 was localized in the nucleus of tobacco cells. Moreover, it could bind to the core region encoded by the H-box sequence (CAATAATTG) in yeast. In response to mannitol treatments, the Phehdz1-overexpressing transgenic rice had a higher germination rate and longer shoots than the wild-type controls. Moreover, Phehdz1-overexpressing rice plants had a higher survival rate as well as higher relative water and proline contents, but a lower malondialdehyde content, than the WT plants after a 30% polyethylene glycol 6000 treatment. Accordingly, the overexpression of Phehdz1 enhances the drought tolerance of transgenic rice. Many of the differentially expressed genes identified by a transcriptome analysis are involved in MAPK signal transduction and the biosynthesis of secondary metabolites. Thus, the overexpression of Phehdz1 enhances the drought stress tolerance of transgenic rice, while also potentially modulating the expression of metabolism-related genes.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Structural determinants for NF-Y subunit organization and NF-Y/DNA association in plants

NF-Y transcription factor comprises three subunits: NF-YA, NF-YB and NF-YC. NF-YB and NF-YC dimerize through their histone fold domain (HFD), which can bind DNA in a non-sequence-specific fashion while serving as a scaffold for NF-YA trimerization. Upon trimerization, NF-YA specifically recognizes the CCAAT box sequence on promoters and enhancers. In plants, each NF-Y subunit is encoded by several genes giving rise to hundreds of potential heterotrimeric combinations. In addition, plant NF-YBs and NF-YCs interact with other protein partners to recognize a plethora of genomic motifs, as the CCT protein family that binds CORE sites. The NF-Y subunit organization and its DNA-binding properties, together with the NF-Y HFD capacity to adapt different protein modules, represent plant-specific features that play a key role in development, growth and reproduction. Despite their relevance, these features are still poorly understood at the molecular level. Here, we present the structures of Arabidopsis and rice NF-YB/NF-YC dimers, and of an Arabidopsis NF-Y trimer in complex with the FT CCAAT box, together with biochemical data on NF-Y mutants. The dimeric structures identify the key residues for NF-Y HFD stabilization. The NF-Y/DNA structure and the mutation experiments shed light on HFD trimerization interface properties and the NF-YA sequence appetite for the bases flanking the CCAAT motif. These data explain the logic of plant NF-Y gene expansion: the trimerization adaptability and the flexible DNA-binding rules serve the scopes of accommodating the large number of NF-YAs, CCTs and possibly other NF-Y HFD binding partners and a diverse audience of genomic motifs.

Characterization of the N- terminal domain of Mre11 protein from rice (OsMre11) Oryza sativa

Chromosomal breaks occur in the genome of all living organisms upon exposure to ionizing radiation, xenobiotics and as intermediates during normal cell cycle progression. Most of the information on DNA repair process has emerged from bacteria, human, mice, and yeast while information on plant DNA repair genes and proteins is limited. Among other DNA repair proteins, MRE11 forms the core of the MRN (Mre11-Rad50-Nbs1) complex and is the first responder to double strand breaks (DSBs), promotes repair either by Non-Homologous End Joining (NHEJ) or Homologous Recombination (HR). Till date, MRE11 has not been biochemically characterized from plant systems. Here, we report the in vitro biochemical activities of Oryza sativa MRE11. We cloned and purified the N- terminal region of OsMre11, which represents both the nuclease and DNA binding domains. The N- terminal end of OsMre11-N protein (∼55.0 kDa) showed binding activity with dsDNA, ssDNA and G-quadruplex DNA. Tryptophan fluorescence analysis also showed that OsMre11-N protein binds to ssDNA, dsDNA and G4 DNA in a protein concentration dependant manner. Additionally, OsMre11 protein showed exonuclease activity only in the presence of Mn²⁺. A protein concentration dependant endonuclease activity also was observed and was enhanced in the presence of Mn²⁺, Mg²⁺ and Ca²⁺. Put together, OsMre11 has properties similar to its counterparts in yeast and humans and may play an important role in cellular response to DNA damage in plants, especially rice.

Nonhost Disease Resistance in Pea: Chitosan's Suggested Role in DNA Minor Groove Actions Relative to Phytoalexin-Eliciting Anti-Cancer Compounds

A stable intense resistance called "nonhost resistance" generates a complete multiple-gene resistance against plant pathogenic species that are not pathogens of pea such as the bean pathogen, Fusarium solani f. sp. phaseoli (Fsph). Chitosan is a natural nonhost resistance response gene activator of defense responses in peas. Chitosan may share with cancer-treatment compounds, netropsin and some anti-cancer drugs, a DNA minor groove target in plant host tissue. The chitosan heptamer and netropsin have the appropriate size and charge to reside in the DNA minor groove. The localization of a percentage of administered radio-labeled chitosan in the nucleus of plant tissue in vivo indicates its potential to transport to site(s) within the nuclear chromatin (1,2). Other minor groove-localizing compounds administered to pea tissue activate the same secondary plant pathway that terminates in the production of the anti-fungal isoflavonoid, pisatin an indicator of the generated resistance response. Some DNA minor groove compounds also induce defense genes designated as "pathogenesis-related" (PR) genes. Hypothetically, DNA targeting components alter host DNA in a manner enabling the transcription of defense genes previously silenced or minimally expressed. Defense-response-elicitors can directly (a) target host DNA at the site of transcription or (b) act by a series of cascading events beginning at the cell membrane and indirectly influence transcription. A single defense response, pisatin induction, induced by chitosan and compounds with known DNA minor groove attachment potential was followed herein. A hypothesis is formulated suggesting that this DNA target may be accountable for a portion of the defense response generated in nonhost resistance.