Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato.

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening," is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200 bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene and the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

Cold-induced deposition of bivalent H3K4me3-H3K27me3 modification and nucleosome depletion in Arabidopsis.

Epigenetic regulation of gene expression plays a crucial role in plant development and environmental adaptation. The H3K4me3 and H3K27me3 have not only been discovered in the regulation of gene expression in multiple biological processes but also in responses to abiotic stresses in plants. However, evidence for the presence of both H3K4me3 and H3K27me3 on the same nucleosome is sporadic. Cold-induced deposition of bivalent H3K4me3-H3K27me3 modifications and nucleosome depletion over a considerable number of active genes is documented in potato tubers and provides clues on an additional role of the bivalent modifications. Limited by the available information of genes encoding PcG/TrxG proteins as well as their corresponding mutants in potatoes, the molecular mechanism underlying the cold-induced deposition of the bivalent mark remains elusive. In this study, we found a similar deposition of the bivalent H3K4me3-H3K27me3 mark over 2129 active genes in cold-treated Arabidopsis Col-0 seedlings. The expression levels of the bivalent mark-associated genes tend to be independent of bivalent modification levels. However, these genes were associated with greater chromatin accessibility, presumably to provide a distinct chromatin environment for gene expression. In mutants clf28 and lhp1, failure to deposit H3K27me3 in active genes upon cold treatment implies that the CLF is potentially involved in cold-induced deposition of H3K27me3, with assistance from LHP1. Failure to deposit H3K4me3 during cold treatment in atx1-2 suggests a regulatory role of ATX1 in the deposition of H3K4me3. In addition, we observed a cold-induced global reduction in nucleosome occupancy, which is potentially mediated by LHP1 in an H3K27me3-dependent manner.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii.

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Genomic rearrangements and evolutionary changes in 3D chromatin topologies in the cotton tribe (Gossypieae).

BACKGROUND: Analysis of the relationship between chromosomal structural variation (synteny breaks) and 3D-chromatin architectural changes among closely related species has the potential to reveal causes and correlates between chromosomal change and chromatin remodeling. Of note, contrary to extensive studies in animal species, the pace and pattern of chromatin architectural changes following the speciation of plants remain unexplored; moreover, there is little exploration of the occurrence of synteny breaks in the context of multiple genome topological hierarchies within the same model species. RESULTS: Here we used Hi-C and epigenomic analyses to characterize and compare the profiles of hierarchical chromatin architectural features in representative species of the cotton tribe (Gossypieae), including Gossypium arboreum, Gossypium raimondii, and Gossypioides kirkii, which differ with respect to chromosome rearrangements. We found that (i) overall chromatin architectural territories were preserved in Gossypioides and Gossypium, which was reflected in their similar intra-chromosomal contact patterns and spatial chromosomal distributions; (ii) the non-random preferential occurrence of synteny breaks in A compartment significantly associate with the B-to-A compartment switch in syntenic blocks flanking synteny breaks; (iii) synteny changes co-localize with open-chromatin boundaries of topologically associating domains, while TAD stabilization has a greater influence on regulating orthologous expression divergence than do rearrangements; and (iv) rearranged chromosome segments largely maintain ancestral in-cis interactions. CONCLUSIONS: Our findings provide insights into the non-random occurrence of epigenomic remodeling relative to the genomic landscape and its evolutionary and functional connections to alterations of hierarchical chromatin architecture, on a known evolutionary timescale.

Development of homozygous tetraploid potato and whole genome doubling-induced the enrichment of H3K27ac and potentially enhanced resistance to cold-induced sweetening in tubers.

Polyploid plants typically display advantages on some agronomically important traits over their diploid counterparts. Extensive studies have shown genetic, transcriptomic, and epigenetic dynamics upon polyploidization in multiple plant species. However, few studies have unveiled those alternations imposed only by ploidy level, without any interference from heterozygosity. Cultivated potato is highly heterozygous. Thus, in this study, we developed two homozygous autotetraploid lines and one homozygous diploid line in parallel from a homozygous diploid potato. We confirmed their ploidy levels using chloroplast counting and karyotyping. Oligo-FISH and genome re-sequencing validated that these potato lines are nearly homozygous. We investigated variations in phenotypes, transcription, and histone modifications between two ploidies. Both autotetraploid lines produced larger but fewer tubers than the diploid line. Interestingly, each autotetraploid line displayed ploidy-related differential expression for various genes. We also discovered a genome-wide enrichment of H3K27ac in genic regions upon whole-genome doubling (WGD). However, such enrichment was not associated with the differential gene expression between two ploidies. The tetraploid lines may exhibit better resistance to cold-induced sweetening (CIS) than the diploid line in tubers, potentially regulated through the expression of CIS-related key genes, which seems to be associated with the levels of H3K4me3 in cold-stored tubers. These findings will help to understand the impacts of autotetraploidization on dynamics of phenotypes, transcription, and histone modifications, as well as on CIS-related genes in response to cold storage.

Chromatin architectural alterations due to null mutation of a major CG methylase in rice.

Associations between 3D chromatin architectures and epigenetic modifications have been characterized in animals. However, any impact of DNA methylation on chromatin architecture in plants is understudied, which is confined to Arabidopsis thaliana. Because plant species differ in genome size, composition, and overall chromatin packing, it is unclear to what extent findings from A. thaliana hold in other species. Moreover, the incomplete chromatin architectural profiles and the low-resolution high-throughput chromosome conformation capture (Hi-C) data from A. thaliana have hampered characterizing its subtle chromatin structures and their associations with DNA methylation. We constructed a high-resolution Hi-C interaction map for the null OsMET1-2 (the major CG methyltransferase in rice) mutant (osmet1-2) and isogenic wild-type rice (WT). Chromatin structural changes occurred in osmet1-2, including intra-/inter-chromosomal interactions, compartment transition, and topologically associated domains (TAD) variations. Our findings provide novel insights into the potential function of DNA methylation in TAD formation in rice and confirmed DNA methylation plays similar essential roles in chromatin packing in A. thaliana and rice.

Quasi-Consensus of Time-Varying Multi-Agent Systems with External Inputs under Deception Attacks.

The quasi-consensus of a class of nonlinear time-varying multi-agent systems suffering from both external inputs and deception attacks is studied in this paper. This is different from a time-varying matrix, which is assumed to be bounded; further reasonable assumptions are supposed. In addition, impulsive deception attacks modeled with Bernoulli variables are considered. Sufficient conditions to achieve quasi-consensus are given, and the upper bounds of the error state related to the deception attacks is derived. Finally, a numerical simulation example is provided to show the validity of the obtained results.

FtMYB18 acts as a negative regulator of anthocyanin/proanthocyanidin biosynthesis in Tartary buckwheat.

KEY MESSAGE: FtMYB18 plays a role in the repression of anthocyanins and proanthocyanidins accumulation by strongly down-regulating the CHS and DFR genes in Tartary buckwheat, and the C5 motif plays an important role in this process. Anthocyanins and proanthocyanidins (PAs) are important flavonoids in Tartary buckwheat (Fagopyrum tataricum Gaertn.), which provides various vibrant color and stronge abiotic stress resistance. Their synthesis is generally regulated by MYB transcription factors at transcription level. However, the negative regulations of MYB and their effects on flavonol metabolism are poorly understood. A SG4-like MYB subfamily TF, FtMYB18, containing C5 motif was identified from Tartary buckwheat. The expression of FtMYB18 was not only showed a negative correlation with anthocyanins and PAs content but also strongly respond to MeJA and ABA. As far as the transgenic lines with FtMYB18 overexpression, anthocyanins and PAs accumulations were decreased through down-regulating expression levels of NtCHS and NtDFR in tobacco, AtDFR and AtTT12 in Arabidopsis, FtCHS, FtDFR and FtANS in Tartary buckwheat hairy roots, respectively. However, FtMYB18 showed no effect on the FLS gene expression and the metabolites content in flavonol synthesis branch. The further molecular interaction analysis indicated FtMYB18 could mediate the inhibition of anthocyanins and PAs synthesis by forming MBW transcriptional complex with FtTT8 and FtTTG1, or MYB-JAZ complex with FtJAZ1/-3/-4/-7. Importantly, in FtMYB18 mutant lines with C5 motif deletion (FtMYB18-C), both of anthocyanins and PAs accumulations had recovered to the similar level as that in wild type, which was attributed to the weakened MBW complex activity or the deficient molecular interaction between FtMYB18ΔC5 with FtJAZ3/-4. The results showed that FtMYB18 could suppress anthocyanins and PAs synthesis at transcription level through the specific interaction of C5 motif with other proteins in Tartary buckwheat.

Cold stress induces enhanced chromatin accessibility and bivalent histone modifications H3K4me3 and H3K27me3 of active genes in potato.

BACKGROUND: Cold stress can greatly affect plant growth and development. Plants have developed special systems to respond to and tolerate cold stress. While plant scientists have discovered numerous genes involved in responses to cold stress, few studies have been dedicated to investigation of genome-wide chromatin dynamics induced by cold or other abiotic stresses. RESULTS: Genomic regions containing active cis-regulatory DNA elements can be identified as DNase I hypersensitive sites (DHSs). We develop high-resolution DHS maps in potato (Solanum tuberosum) using chromatin isolated from tubers stored under room (22 °C) and cold (4 °C) conditions. We find that cold stress induces a large number of DHSs enriched in genic regions which are frequently associated with differential gene expression in response to temperature variation. Surprisingly, active genes show enhanced chromatin accessibility upon cold stress. A large number of active genes in cold-stored tubers are associated with the bivalent H3K4me3-H3K27me3 mark in gene body regions. Interestingly, upregulated genes associated with the bivalent mark are involved in stress response, whereas downregulated genes with the bivalent mark are involved in developmental processes. In addition, we observe that the bivalent mark-associated genes are more accessible than others upon cold stress. CONCLUSIONS: Collectively, our results suggest that cold stress induces enhanced chromatin accessibility and bivalent histone modifications of active genes. We hypothesize that in cold-stored tubers, the bivalent H3K4me3-H3K27me3 mark represents a distinct chromatin environment with greater accessibility, which may facilitate the access of regulatory proteins required for gene upregulation or downregulation in response to cold stress.

Historical Meiotic Crossover Hotspots Fueled Patterns of Evolutionary Divergence in Rice.

Recombination plays an integral role in the creation of novel genetic variation in sexually reproducing species. Despite this important role, the determinants and evolution of crossover hotspots have remained poorly understood in plants. Here, we present a comparative analysis of two rice (Oryza sativa) historical recombination maps from two subspecies (indica and japonica) using 150 resequenced genomes. Fine-scale recombination rates and crossover hotspots were validated by comparison with a consensus genetic map and empirically derived crossovers, respectively. Strikingly, nearly 80% of crossover hotspots were unique to each subspecies, despite their relatively recent divergence and broad-scale correlated recombination rates. Crossover hotspots were enriched with Stowaway and P instability factor (PIF)/Harbinger transposons and overlapped accessible chromatin regions. Increased nucleotide diversity and signatures of population differentiation augmented by Stowaway and PIF/Harbinger transposons were prevalent at subspecies-specific crossover hotspots. Motifs derived from lineage-specific indica and japonica crossover hotspots were nearly identical in the two subspecies, implicating a core set of crossover motifs in rice. Finally, Stowaway and PIF/Harbinger transposons were associated with stabilized G/C bias within highly active hotspots, suggesting that hotspot activity can be fueled by de novo variation. These results provide evolutionary insight into historical crossover hotspots as potentially powerful drivers of sequence and subspecies evolution in plants.

Amplification and adaptation of centromeric repeats in polyploid switchgrass species.

Centromeres in most higher eukaryotes are composed of long arrays of satellite repeats from a single satellite repeat family. Why centromeres are dominated by a single satellite repeat and how the satellite repeats originate and evolve are among the most intriguing and long-standing questions in centromere biology. We identified eight satellite repeats in the centromeres of tetraploid switchgrass (Panicum virgatum). Seven repeats showed characteristics associated with classical centromeric repeats with monomeric lengths ranging from 166 to 187 bp. Interestingly, these repeats share an 80-bp DNA motif. We demonstrate that this 80-bp motif may dictate translational and rotational phasing of the centromeric repeats with the cenH3 nucleosomes. The sequence of the last centromeric repeat, Pv156, is identical to the 5S ribosomal RNA genes. We demonstrate that a 5S ribosomal RNA gene array was recruited to be the functional centromere for one of the switchgrass chromosomes. Our findings reveal that certain types of satellite repeats, which are associated with unique sequence features and are composed of monomers in mono-nucleosomal length, are favorable for centromeres. Centromeric repeats may undergo dynamic amplification and adaptation before the centromeres in the same species become dominated by the best adapted satellite repeat.

Meiotic crossovers are associated with open chromatin and enriched with Stowaway transposons in potato.

BACKGROUND: Meiotic recombination is the foundation for genetic variation in natural and artificial populations of eukaryotes. Although genetic maps have been developed for numerous plant species since the late 1980s, few of these maps have provided the necessary resolution needed to investigate the genomic and epigenomic features underlying meiotic crossovers. RESULTS: Using a whole genome sequencing-based approach, we developed two high-density reference-based haplotype maps using diploid potato clones as parents. The vast majority (81%) of meiotic crossovers were mapped to less than 5 kb. The fine-scale accuracy of crossover detection was validated by Sanger sequencing for a subset of ten crossover events. We demonstrate that crossovers reside in genomic regions of "open chromatin", which were identified based on hypersensitivity to DNase I digestion and association with H3K4me3-modified nucleosomes. The genomic regions spanning crossovers were significantly enriched with the Stowaway family of miniature inverted-repeat transposable elements (MITEs). The occupancy of Stowaway elements in gene promoters is concomitant with an increase in recombination rate. A generalized linear model identified the presence of Stowaway elements as the third most important genomic or chromatin feature behind genes and open chromatin for predicting crossover formation over 10-kb windows. CONCLUSIONS: Collectively, our results suggest that meiotic crossovers in potato are largely determined by the local chromatin status, marked by accessible chromatin, H3K4me3-modified nucleosomes, and the presence of Stowaway transposons.

Recurrent establishment of de novo centromeres in the pericentromeric region of maize chromosome 3.

Centromeres can arise de novo from non-centromeric regions, which are often called "neocentromeres." Neocentromere formation provides the best evidence for the concept that centromere function is not determined by the underlying DNA sequences, but controlled by poorly understood epigenetic mechanisms. Numerous neocentromeres have been reported in several plant and animal species. However, it has been elusive how and why a specific chromosomal region is chosen to be a new centromere during the neocentromere activation events. We report recurrent establishment of neocentromeres in a pericentromeric region of chromosome 3 in maize (Zea mays). This latent region is located in the short arm and is only 2 Mb away from the centromere (Cen3) of chromosome 3. At least three independent neocentromere activation events, which were likely induced by different mechanisms, occurred within this latent region. We mapped the binding domains of CENH3, the centromere-specific H3 histone variant, of the three neocentromeres and analyzed the genomic and epigenomic features associated with Cen3, the de novo centromeres and an inactivated centromere derived from an ancestral chromosome. Our results indicate that lack of genes and transcription and a relatively high level of DNA methylation in this pericentromeric region may provide a favorable chromatin environment for neocentromere activation.

Isolation and Proteomics Analysis of Barley Centromeric Chromatin Using PICh.

Identification of proteins that are directly or indirectly associated with a specific DNA sequence is often an important goal in molecular biology research. Proteomics of isolated chromatin fragments (PICh) is a technique used to isolate chromatin that contains homologous DNA sequence to a specific nucleic acid probe. All proteins directly and indirectly associated with the DNA sequences that hybridize to the probe are then identified by proteomics.1 We used the PICh technique to isolate chromatin associated with the centromeres of barley (Hordeum vulgare) by using a 2'-deoxy-2'fluoro-ribonucleotides (2'-F RNA) probe that is homologous to the AGGGAG satellite DNA specific to barley centromeres. Proteins associated with the barley centromeric chromatin were then isolated and identified by mass spectrometry. Both alpha-cenH3 and beta-cenH3, the two centromeric histone H3 variants associated with barley centromeres, were positively identified. Interestingly, several different H2A and H2B variants were recovered in the PIChed chromatin. The limitations and future potential of PICh in plant chromatin research are discussed.

Genome Reduction Uncovers a Large Dispensable Genome and Adaptive Role for Copy Number Variation in Asexually Propagated Solanum tuberosum.

Clonally reproducing plants have the potential to bear a significantly greater mutational load than sexually reproducing species. To investigate this possibility, we examined the breadth of genome-wide structural variation in a panel of monoploid/doubled monoploid clones generated from native populations of diploid potato (Solanum tuberosum), a highly heterozygous asexually propagated plant. As rare instances of purely homozygous clones, they provided an ideal set for determining the degree of structural variation tolerated by this species and deriving its minimal gene complement. Extensive copy number variation (CNV) was uncovered, impacting 219.8 Mb (30.2%) of the potato genome with nearly 30% of genes subject to at least partial duplication or deletion, revealing the highly heterogeneous nature of the potato genome. Dispensable genes (>7000) were associated with limited transcription and/or a recent evolutionary history, with lower deletion frequency observed in genes conserved across angiosperms. Association of CNV with plant adaptation was highlighted by enrichment in gene clusters encoding functions for environmental stress response, with gene duplication playing a part in species-specific expansions of stress-related gene families. This study revealed unique impacts of CNV in a species with asexual reproductive habits and how CNV may drive adaption through evolution of key stress pathways.

Silencing of vacuolar invertase and asparagine synthetase genes and its impact on acrylamide formation of fried potato products.

Acrylamide is produced in a wide variety of carbohydrate-rich foods during high-temperature cooking. Dietary acrylamide is a suspected human carcinogen, and health concerns related to dietary acrylamide have been raised worldwide. French fries and potato chips contribute a significant proportion to the average daily intake of acrylamide, especially in developed countries. One way to mitigate health concerns related to acrylamide is to develop potato cultivars that have reduced contents of the acrylamide precursors asparagine, glucose and fructose in tubers. We generated a large number of silencing lines of potato cultivar Russet Burbank by targeting the vacuolar invertase gene VInv and the asparagine synthetase genes StAS1 and StAS2 with a single RNA interference construct. The transcription levels of these three genes were correlated with reducing sugar (glucose and fructose) and asparagine content in tubers. Fried potato products from the best VInv/StAS1/StAS2-triple silencing lines contained only one-fifteenth of the acrylamide content of the controls. Interestingly, the extent of acrylamide reduction of the best triple silencing lines was similar to that of the best VInv-single silencing lines developed previously from the same potato cultivar Russet Burbank. These results show that an acrylamide mitigation strategy focused on developing potato cultivars with low reducing sugars is likely to be an effective and sufficient approach for minimizing the acrylamide-forming potential of French fry processing potatoes.

Three potato centromeres are associated with distinct haplotypes with or without megabase-sized satellite repeat arrays.

We report discoveries of different haplotypes associated with the centromeres of three potato chromosomes, including haplotypes composed of long arrays of satellite repeats and haplotypes lacking the same repeats. These results are in favor of the hypothesis that satellite repeat-based centromeres may originate from neocentromeres that lack repeats.

High-resolution mapping of open chromatin in the rice genome.

Gene expression is controlled by the complex interaction of transcription factors binding to promoters and other regulatory DNA elements. One common characteristic of the genomic regions associated with regulatory proteins is a pronounced sensitivity to DNase I digestion. We generated genome-wide high-resolution maps of DNase I hypersensitive (DH) sites from both seedling and callus tissues of rice (Oryza sativa). Approximately 25% of the DH sites from both tissues were found in putative promoters, indicating that the vast majority of the gene regulatory elements in rice are not located in promoter regions. We found 58% more DH sites in the callus than in the seedling. For DH sites detected in both the seedling and callus, 31% displayed significantly different levels of DNase I sensitivity within the two tissues. Genes that are differentially expressed in the seedling and callus were frequently associated with DH sites in both tissues. The DNA sequences contained within the DH sites were hypomethylated, consistent with what is known about active gene regulatory elements. Interestingly, tissue-specific DH sites located in the promoters showed a higher level of DNA methylation than the average DNA methylation level of all the DH sites located in the promoters. A distinct elevation of H3K27me3 was associated with intergenic DH sites. These results suggest that epigenetic modifications play a role in the dynamic changes of the numbers and DNase I sensitivity of DH sites during development.

Molecular characterization of a wheat -Thinopyrum ponticum partial amphiploid and its derived substitution line for resistance to stripe rust.

Stripe rust (caused by Puccinia striiformis) occurs annually in most wheat-growing areas of the world. Thinopyrum ponticum has provided novel rust resistance genes to protect wheat from this fungal disease. Wheat - Th. ponticum partial amphiploid line 7430 and a substitution line X005 developed from crosses between wheat and 7430 were resistant to stripe rust isolates from China. Genomic in situ hybridization (GISH) analysis using Pseudoroegneria spicata genomic DNA as a probe demonstrated that the partial amphiploid line 7430 contained ten J(s) and six J genome chromosomes, and line X005 had a pair of J(s)-chromosomes. Giemsa-C banding further revealed that both lines 7430 and X005 were absent of wheat chromosomes 6B. The EST based PCR confirmed that the introduced J(s) chromosomes belonging to linkage group 6, indicating that line X005 was a 6J(s)/6B substitution line. Both resistance observation and sequence characterized amplified region (SCAR) markers displayed that the introduced chromosomes 6J(s) were responsible for the stripe rust resistances. Therefore, lines 7430 and X005 can be used as a donor in wheat breeding for stripe rust resistance.