3D genome architecture coordinates trans and cis regulation of differentially expressed ear and tassel genes in maize

Maintenance and dynamic reprogramming of chromatin organization during development

Controlled transcription of genes is critical for cell differentiation and development. Gene expression regulation therefore involves a multilayered control from nucleosome composition in histone variants and their post-translational modifications to higher-order folding of chromatin fibers and chromatin interactions in nuclear space. Recent technological advances have allowed gaining insight into these mechanisms, the interplay between local and higher-order chromatin organization, and the dynamic changes that occur during stress response and developmental transitions. In this review, we will discuss chromatin organization from the nucleosome to its three-dimensional structure in the nucleus, and consider how these different layers of organization are maintained during the cell cycle or rapidly reprogrammed during development.

BL-Hi-C reveals the 3D genome structure of Brassica crops with high sensitivity

High-throughput Chromatin Conformation Capture (Hi-C) technologies can be used to investigate the three-dimensional genomic structure of plants. However, the practical utility of these technologies is impeded by significant background noise, hindering their capability in detecting fine 3D genomic structures. In this study, we optimized the Bridge Linker Hi-C technology (BL-Hi-C) to comprehensively investigate the 3D chromatin landscape of Brassica rapa and Brassica oleracea. The Bouquet configuration of both B. rapa and B. oleracea was elucidated through the construction of a 3D genome simulation. The optimized BL-Hi-C exhibited lower background noise compared to conventional Hi-C methods. Taking this advantage, we used BL-Hi-C to identify FLC gene loops in Arabidopsis, B. rapa, and B. oleracea. We observed that gene loops of FLC2 exhibited conservation across Arabidopsis, B. rapa, and B. oleracea. While gene loops of syntenic FLCs exhibited conservation across B. rapa and B. oleracea, variations in gene loops were evident among multiple paralogs FLCs within the same species. Collectively, our findings highlight the high sensitivity of optimized BL-Hi-C as a powerful tool for investigating the fine 3D genomic organization.

Remodeling of the 3D chromatin architecture in the marine microalga Nannochloropsis oceanica during lipid accumulation

BACKGROUND

Genomic three-dimensional (3D) spatial organization plays a key role in shaping gene expression and associated chromatin modification, and it is highly sensitive to environmental stress conditions. In microalgae, exposure to nitrogen stress can drive lipid accumulation, yet the associated functional alterations in the spatial organization of the microalgal genome have yet to be effectively characterized.

RESULTS

Accordingly, the present study employed RNA-seq, Hi-C, and ChIP-seq approaches to explore the relationship between 3D chromosomal architecture and gene expression during lipid accumulation in the marine microalga Nannochloropsis oceanica in response to nitrogen deprivation (ND). These analyses revealed that ND resulted in various changes in chromosomal organization, including A/B compartment transitions, topologically associating domain (TAD) shifts, and the disruption of short-range interactions. Significantly higher levels of gene expression were evident in A compartments and TAD boundary regions relative to B compartments and TAD interior regions, consistent with observed histone modification enrichment in these areas. ND-induced differentially expressed genes (DEGs) were notably enriched in altered TAD-associated regions and regions exhibiting differential genomic contact. These DEGs were subjected to Gene Ontology (GO) term analyses that indicated they were enriched in the 'fatty acid metabolism', 'response to stress', 'carbon fixation' and 'photosynthesis' functional categories, in line with the ND treatment conditions used to conduct this study. These data indicate that Nannochloropsis cells exhibit a clear association between chromatin organization and transcriptional activity under nitrogen stress conditions. Pronounced and extensive histone modifications were evident in response to ND. Observed changes in chromatin architecture were linked to shifts in histone modifications and gene expression.

CONCLUSIONS

Overall, the reprogramming of many lipid metabolism-associated genes was evident under nitrogen stress conditions with respect to both histone modifications and chromosomal organization. Together these results revealed that higher-order chromatin architecture represents a new layer that can guide efforts to understand the transcriptional regulation of lipid metabolism in nitrogen-deprived microalgae.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Gene expression and expression quantitative trait loci analyses uncover natural variations underlying the improvement of important agronomic traits during modern maize breeding

Maize (Zea mays L.) is a major staple crop worldwide, and during modern maize breeding, cultivars with increased tolerance to high-density planting and higher yield per plant have contributed significantly to the increased yield per unit land area. Systematically identifying key agronomic traits and their associated genomic changes during modern maize breeding remains a significant challenge because of the complexity of genetic regulation and the interactions of the various agronomic traits, with most of them being controlled by numerous small-effect quantitative trait loci (QTLs). Here, we performed phenotypic and gene expression analyses for a set of 137 elite inbred lines of maize from different breeding eras in China. We found four yield-related traits are significantly improved during modern maize breeding. Through gene-clustering analyses, we identified four groups of expressed genes with distinct trends of expression pattern change across the historical breeding eras. In combination with weighted gene co-expression network analysis, we identified several candidate genes regulating various plant architecture- and yield-related agronomic traits, such as ZmARF16, ZmARF34, ZmTCP40, ZmPIN7, ZmPYL10, ZmJMJ10, ZmARF1, ZmSWEET15b, ZmGLN6 and Zm00001d019150. Further, by combining expression quantitative trait loci (eQTLs) analyses, correlation coefficient analyses and population genetics, we identified a set of candidate genes that might have been under selection and contributed to the genetic improvement of various agronomic traits during modern maize breeding, including a number of known key regulators of plant architecture, flowering time and yield-related traits, such as ZmPIF3.3, ZAG1, ZFL2 and ZmBES1. Lastly, we validated the functional variations in GL15, ZmPHYB2 and ZmPYL10 that influence kernel row number, flowering time, plant height and ear height, respectively. Our results demonstrates the effectiveness of our combined approaches for uncovering key candidate regulatory genes and functional variation underlying the improvement of important agronomic traits during modern maize breeding, and provide a valuable genetic resource for the molecular breeding of maize cultivars with tolerance for high-density planting.

Multi-omics analysis reveals spatiotemporal regulation and function of heteromorphic leaves in Populus

Despite the high economic and ecological importance of forests, our knowledge of the adaptive evolution of leaf traits remains very limited. Euphrates poplar (Populus euphratica), which has high tolerance to arid environment, has evolved four heteromorphic leaf forms, including narrow (linear and lanceolate) and broad (ovate and broad-ovate) leaves on different crowns. Here, we revealed the significant functional divergence of four P. euphratica heteromorphic leaves at physiological and cytological levels. Through global analysis of transcriptome and DNA methylation across tree and leaf developmental stages, we revealed that gene expression and DNA epigenetics differentially regulated key processes involving development and functional adaptation of heteromorphic leaves, such as hormone signaling pathways, cell division, and photosynthesis. Combined analysis of gene expression, methylation, ATAC-seq, and Hi-C-seq revealed longer interaction of 3D genome, hypomethylation, and open chromatin state upregulates IAA-related genes (such as PIN-FORMED1 and ANGUSTIFOLIA3) and promotes the occurrence of broad leaves while narrow leaves were associated with highly concentrated heterochromatin, hypermethylation, and upregulated abscisic acid pathway genes (such as Pyrabactin Resistance1-like10). Therefore, development of P. euphratica heteromorphic leaves along with functional divergence was regulated by differentially expressed genes, DNA methylation, chromatin accessibility, and 3D genome remodeling to adapt to the arid desert. This study advances our understanding of differential regulation on development and functional divergence of heteromorphic leaves in P. euphratica at the multi-omics level and provides a valuable resource for investigating the adaptive evolution of heteromorphic leaves in Populus.

Haplotype mapping of H3K27me3-associated chromatin interactions defines topological regulation of gene silencing in rice

Histone modification H3K27me3 is an important chromatin mark that plays vital roles in repressing expression of developmental genes. Here, we construct high-resolution 3D genome maps using long-read chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and characterize H3K27me3-associated chromatin interactions in an elite rice hybrid, Shanyou 63. We find that many H3K27me3-marked regions may function as silencer-like regulatory elements. The silencer-like elements can come into proximity with distal target genes via forming chromatin loops in 3D space of the nuclei, regulating gene silencing and plant traits. Natural and induced deletion of silencers upregulate expression of distal connected genes. Furthermore, we identify extensive allele-specific chromatin loops. We find that genetic variations alter allelic chromatin topology, thus modulating allelic gene imprinting in rice hybrids. In conclusion, the characterization of silencer-like regulatory elements and haplotype-resolved chromatin interaction maps provide insights into the understanding of molecular mechanisms underlying allelic gene silencing and plant trait controlling.

Transcriptome and DNA methylome divergence of inflorescence development between two ecotypes in Panicum hallii

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the two major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMRs-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Toward cis-regulation in soybean: a 3D genome scope

In eukaryotic cells, 3D genome plays an important role in the regulation of gene spatiotemporal expression, which is essential for the biological and developmental processes in a life cycle. In the past decade, the development of high-throughput technologies greatly enhances our ability to map the 3D genome organization, identifies multiple 3D genome structures, and investigates the functional role of 3D genome organization in gene regulation, which facilitates our understandings of cis-regulatory landscape and biological development. Comparing with the comprehensive analyses of 3D genome in mammals and model plants, the progress in soybean is much less. Future development and application of tools to precisely manipulate 3D genome structure at different levels will significantly strengthen the functional genome study and molecular breeding in soybean. Here, we review the recent progresses in 3D genome study and discuss future directions, which may help to improve soybean 3D functional genome study and molecular breeding.

Reorganization of three-dimensional chromatin architecture in Medicago truncatula under phosphorus deficiency

Emerging evidence reveals that the three-dimensional (3D) chromatin architecture plays a key regulatory role in various biological processes of plants. However, information on the 3D chromatin architecture of the legume model plant Medicago truncatula and its potential roles in the regulation of response to mineral nutrient deficiency are very limited. Using high-resolution chromosome conformation capture sequencing, we identified the 3D genome structure of M. truncatula in terms of A/B compartments, topologically associated domains (TADs) and chromatin loops. The gene density, expressional level, and active histone modification were higher in A compartments than in B compartments. Moreover, we analysed the 3D chromatin architecture reorganization in response to phosphorus (P) deficiency. The intra-chromosomal cis-interaction proportion was increased by P deficiency, and a total of 748 A/B compartment switch regions were detected. In these regions, density changes in H3K4me3 and H3K27ac modifications were associated with expression of P deficiency-responsive genes involved in root system architecture and hormonal responses. Furthermore, these genes enhanced P uptake and mobilization by increasing root surface area and strengthening signal transduction under P deficiency. These findings advance our understanding of the potential roles of 3D chromatin architecture in responses of plants in general, and in particular in M. truncatula, to P deficiency.

Differential chromatin binding preference is the result of the neo-functionalization of the TB1 clade of TCP transcription factors in grasses

The understanding of neo-functionalization of plant transcription factors (TFs) after gene duplication has been extensively focused on changes in protein-protein interactions, the expression pattern of TFs, or the variation of cis-elements bound by TFs. Yet, the main molecular role of a TF, that is, its specific chromatin binding for the direct regulation of target gene expression, continues to be mostly overlooked. Here, we studied the TB1 clade of the TEOSINTE BRANCHED 1, CYCLOIDEA, PROLIFERATING CELL FACTORS (TCP) TF family within the grasses (Poaceae). We identified an Asp/Gly amino acid replacement within the TCP domain, originated within a paralog TIG1 clade exclusive for grasses. The heterologous expression of Zea mays TB1 and its two paralogs BAD1 and TIG1 in Arabidopsis mutant plants lacking the TB1 ortholog BRC1 revealed distinct functions in plant development. Notably, the Gly acquired in the TIG1 clade does not impair TF homodimerization and heterodimerization, while it modulates chromatin binding preferences. We found that in vivo TF recognition of target promoters depends on this Asp/Gly mutation and directly impacts downstream gene expression and subsequent plant development. These results provided new insights into how natural selection fine-tunes gene expression regulation after duplication of TFs to define plant architecture.

De novo genome assembly and analyses of 12 founder inbred lines provide insights into maize heterosis

Hybrid maize displays superior heterosis and contributes over 30% of total worldwide cereal production. However, the molecular mechanisms of heterosis remain obscure. Here we show that structural variants (SVs) between the parental lines have a predominant role underpinning maize heterosis. De novo assembly and analyses of 12 maize founder inbred lines (FILs) reveal abundant genetic variations among these FILs and, through expression quantitative trait loci and association analyses, we identify several SVs contributing to genomic and phenotypic differentiations of various heterotic groups. Using a set of 91 diallel-cross F₁ hybrids, we found strong positive correlations between better-parent heterosis of the F₁ hybrids and the numbers of SVs between the parental lines, providing concrete genomic support for a prevalent role of genetic complementation underlying heterosis. Further, we document evidence that SVs in both ZAR1 and ZmACO2 contribute to yield heterosis in an overdominance fashion. Our results should promote genomics-based breeding of hybrid maize.

Cis-Regulation by NACs: A Promising Frontier in Wheat Crop Improvement

Crop traits are controlled by multiple genes; however, the complex spatio-temporal transcriptional behavior of genes cannot be fully understood without comprehending the role of transcription factors (TFs) and the underlying mechanisms of the binding interactions of their cis-regulatory elements. NAC belongs to one of the largest families of plant-specific TFs and has been associated with the regulation of many traits. This review provides insight into the cis-regulation of genes by wheat NACs (TaNACs) for the improvement in yield-related traits, including phytohormonal homeostasis, leaf senescence, seed traits improvement, root modulation, and biotic and abiotic stresses in wheat and other cereals. We also discussed the current potential, knowledge gaps, and prospects of TaNACs.

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.

The chromatin accessibility landscape of pistils and anthers in rice

Transcription activation is tightly associated with the openness of chromatin and allows direct contact between transcriptional regulators and their targeted DNA for gene expression. However, there are limited studies on the annotation of open chromatin regions (OCRs) in rice (Oryza sativa), especially those in reproductive organs. Here, we characterized OCRs in rice pistils and anthers with an assay for transposase-accessible chromatin using sequencing. Despite a large overlap, we found more OCRs in pistils than in anthers. These OCRs were enriched in gene transcription start sites (TSSs) and showed tight associations with gene expression. Transcription factor (TF) binding motifs were enriched at these OCRs as validated by TF chromatin immunoprecipitation followed by sequencing. Pistil-specific OCRs provided potential regulatory networks by binding directly to the targets, indicating that pistil-specific OCRs may be indicators of cis-regulatory elements in regulating pistil development, which are absent in anthers. We also found that open chromatin of pistils and anthers responded differently to low temperature (LT). These data offer a comprehensive overview of OCRs regulating reproductive organ development and LT responses in rice.

Comprehensive transcriptional variability analysis reveals gene networks regulating seed oil content of Brassica napus

BACKGROUND

Regulation of gene expression plays an essential role in controlling the phenotypes of plants. Brassica napus (B. napus) is an important source for the vegetable oil in the world, and the seed oil content is an important trait of B. napus.

RESULTS

We perform a comprehensive analysis of the transcriptional variability in the seeds of B. napus at two developmental stages, 20 and 40 days after flowering (DAF). We detect 53,759 and 53,550 independent expression quantitative trait loci (eQTLs) for 79,605 and 76,713 expressed genes at 20 and 40 DAF, respectively. Among them, the local eQTLs are mapped to the adjacent genes more frequently. The adjacent gene pairs are regulated by local eQTLs with the same open chromatin state and show a stronger mode of expression piggybacking. Inter-subgenomic analysis indicates that there is a feedback regulation for the homoeologous gene pairs to maintain partial expression dosage. We also identify 141 eQTL hotspots and find that hotspot87-88 co-localizes with a QTL for the seed oil content. To further resolve the regulatory network of this eQTL hotspot, we construct the XGBoost model using 856 RNA-seq datasets and the Basenji model using 59 ATAC-seq datasets. Using these two models, we predict the mechanisms affecting the seed oil content regulated by hotspot87-88 and experimentally validate that the transcription factors, NAC13 and SCL31, positively regulate the seed oil content.

CONCLUSIONS

We comprehensively characterize the gene regulatory features in the seeds of B. napus and reveal the gene networks regulating the seed oil content of B. napus.

A cost-effective tsCUT&Tag method for profiling transcription factor binding landscape

Knowledge of the transcription factor binding landscape (TFBL) is necessary to analyze gene regulatory networks for important agronomic traits. However, a low-cost and high-throughput in vivo chromatin profiling method is still lacking in plants. Here, we developed a transient and simplified cleavage under targets and tagmentation (tsCUT&Tag) that combines transient expression of transcription factor proteins in protoplasts with a simplified CUT&Tag without nucleus extraction. Our tsCUT&Tag method provided higher data quality and signal resolution with lower sequencing depth compared with traditional ChIP-seq. Furthermore, we developed a strategy combining tsCUT&Tag with machine learning, which has great potential for profiling the TFBL across plant development.

Plant synthetic epigenomic engineering for crop improvement

Efforts have been directed to redesign crops with increased yield, stress adaptability, and nutritional value through synthetic biology-the application of engineering principles to biology. A recent expansion in our understanding of how epigenetic mechanisms regulate plant development and stress responses has unveiled a new set of resources that can be harnessed to develop improved crops, thus heralding the promise of "synthetic epigenetics." In this review, we summarize the latest advances in epigenetic regulation and highlight how innovative sequencing techniques, epigenetic editing, and deep learning-driven predictive tools can rapidly extend these insights. We also proposed the future directions of synthetic epigenetics for the development of engineered smart crops that can actively monitor and respond to internal and external cues throughout their life cycles.

Characterization of the MADS-Box Gene Family in Akebia trifoliata and Their Evolutionary Events in Angiosperms

As the largest clade of modern plants, flower plants have evolved a wide variety of flowers and fruits. MADS-box genes play key roles in regulating plant morphogenesis, while basal eudicots have an evolutionarily important position of acting as an evolutionary bridge between basal angiosperms and core eudicots. Akebia trifoliata is an important member of the basal eudicot group. To study the early evolution of angiosperms, we identified and characterized the MADS-Box gene family on the whole-genome level of A. trifoliata. There were 47 MADS-box genes (13 type I and 34 type II genes) in the A. trifoliata genome; type I genes had a greater gene length and coefficient of variation and a smaller exon number than type II genes. A total of 27 (57.4%) experienced whole or segmental genome duplication and purifying selection. A transcriptome analysis suggested that three and eight genes were involved in whole fruit and seed development, respectively. The diversification and phylogenetic analysis of 1479 type II MADS-box genes of 22 angiosperm species provided some clues indicating that a γ whole genome triplication event of eudicots possibility experienced a two-step process. These results are valuable for improving A. trifoliata fruit traits and theoretically elucidating evolutionary processes of angiosperms, especially eudicots.

The Heterogeneity in the Landscape of Gene Dominance in Maize is Accompanied by Unique Chromatin Environments

Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Transcription factors (TFs) play a critical role in determining cell fate decisions by integrating developmental and environmental signals through binding to specific cis-regulatory modules and regulating spatio-temporal specificity of gene expression patterns. Precise identification of functional TF binding sites in time and space not only will revolutionize our understanding of regulatory networks governing cell fate decisions but is also instrumental to uncover how genetic variations cause morphological diversity or disease. In this review, we discuss recent advances in mapping TF binding sites and characterizing the various parameters underlying the complexity of binding site recognition by TFs.

Integrated Analysis of Hi-C and RNA-Seq Reveals the Molecular Mechanism of Autopolyploid Growth Advantages in Pak Choi (Brassica rapa ssp. chinensis)

Polyploids generated by the replication of a single genome (autopolyploid) or synthesis of two or more distinct genomes (allopolyploid) usually show significant advantages over their diploid progenitors in biological characteristics, including growth and development, nutrient accumulation, and plant resistance. Whereas, the impacts of genomic replication on transcription regulation and chromatin structure in pak choi have not been explored fully. In this study, we observed the transcriptional and genomic structural alterations between diploid B. rapa (AA) and artificial autotetraploid B. rapa (AAAA) using RNA-seq and Hi-C. RNA-seq revealed 1,786 differentially expressed genes (DEGs) between the diploids and autotetraploids, including 717 down-regulated and 1,069 up-regulated genes in autotetraploids. Of all the 1,786 DEGs, 23 DEGs (10 down-regulated DEGs in autotetraploids) were involved in Compartment A-B shifts, while 28 DEGs (20 up-regulated DEGs in autotetraploids) participated in Compartment B-A shifts. Moreover, there were 15 DEGs in activated topologically associating domains (TADs) (9 up-regulated DEGs in diploids) and 80 DEGs in repressed TADs (49 down-regulated DEGs in diploids). Subsequently, eight DEGs with genomic structural variants were selected as potential candidate genes, including four DEGs involved in photosynthesis (BraA01003143, BraA09002798, BraA04002224, and BraA08000594), three DEGs related to chloroplast (BraA05002974, BraA05001662, and BraA04001148), and one DEG associated with disease resistance (BraA09004451), which all showed high expression in autotetraploids. Overall, our results demonstrated that integrative RNA-seq and Hi-C analysis can identify related genes to phenotypic traits and also provided new insights into the molecular mechanism of the growth advantage of polyploids.

Spatial Features and Functional Implications of Plant 3D Genome Organization

The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.

Fine-tuning shoot meristem size to feed the world

In order to maintain food security for the world's growing population, crop yields need to be significantly improved. Domestication and crop improvement involve modification of traits such as fruit size and seed number to optimize productivity. Although these traits are selected at the mature stage, they are determined during the development of shoot meristem, a tissue that forms successive meristems and reproductive organs that make edible fruits or seeds. Therefore, the architecture of reproductive organs and yield-related traits are determined during the maturation of shoot meristem. Here, we highlight recent progress in understanding how shoot meristem size affects yield-related traits and outline the strategies to fine-tune meristem regulatory genes to meet the demands of a growing population and promote sustainable agriculture.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

Dynamic 3D genome architecture of cotton fiber reveals subgenome-coordinated chromatin topology for 4-staged single-cell differentiation

Background

Despite remarkable advances in our knowledge of epigenetically mediated transcriptional programming of cell differentiation in plants, little is known about chromatin topology and its functional implications in this process.

Results

To interrogate its significance, we establish the dynamic three-dimensional (3D) genome architecture of the allotetraploid cotton fiber, representing a typical single cell undergoing staged development in plants. We show that the subgenome-relayed switching of the chromatin compartment from active to inactive is coupled with the silencing of developmentally repressed genes, pinpointing subgenome-coordinated contribution to fiber development. We identify 10,571 topologically associating domain-like (TAD-like) structures, of which 25.6% are specifically organized in different stages and 75.23% are subject to partition or fusion between two subgenomes. Notably, dissolution of intricate TAD-like structure cliques showing long-range interactions represents a prominent characteristic at the later developmental stage. Dynamic chromatin loops are found to mediate the rewiring of gene regulatory networks that exhibit a significant difference between the two subgenomes, implicating expression bias of homologous genes.

Conclusions

This study sheds light on the spatial-temporal asymmetric chromatin structures of two subgenomes in the cotton fiber and offers a new insight into the regulatory orchestration of cell differentiation in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02616-y.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Chromatin and regulatory differentiation between bundle sheath and mesophyll cells in maize

C₄ plants partition photosynthesis enzymes between the bundle sheath (BS) and the mesophyll (M) cells for the better delivery of CO₂ to RuBisCO and to reduce photorespiration. To better understand how C₄ photosynthesis is regulated at the transcriptional level, we performed RNA-seq, ATAC-seq, ChIP-seq and Bisulfite-seq (BS-seq) on BS and M cells isolated from maize leaves. By integrating differentially expressed genes with chromatin features, we found that chromatin accessibility coordinates with epigenetic features, especially H3K27me3 modification and CHH methylation, to regulate cell type-preferentially enriched gene expression. Not only the chromatin-accessible regions (ACRs) proximal to the genes (pACRs) but also the distal ACRs (dACRs) are determinants of cell type-preferentially enriched expression. We further identified cell type-preferentially enriched motifs, e.g. AAAG for BS cells and TGACC/T for M cells, and determined their corresponding transcription factors: DOFs and WRKYs. The complex interaction between cis and trans factors in the preferential expression of C₄ genes was also observed. Interestingly, cell type-preferentially enriched gene expression can be fine-tuned by the coordination of multiple chromatin features. Such coordination may be critical in ensuring the cell type-specific function of key C₄ genes. Based on the observed cell type-preferentially enriched expression pattern and coordinated chromatin features, we predicted a set of functionally unknown genes, e.g. Zm00001d042050 and Zm00001d040659, to be potential key C₄ genes. Our findings provide deep insight into the architectures associated with C₄ gene expression and could serve as a valuable resource to further identify the regulatory mechanisms present in C₄ species.

Single-cell analysis of cis-regulatory elements

Plant tissues and organs are composed of functionally discrete cell types that are all defined by the same genome sequence. Cell-type variation in part arises from differential accessibility of cis-regulatory elements that encode the blueprints for transcriptional programs underlying cell identity and function. Owing to technical limitations, the role of cis-regulatory elements in cell identity maintenance, differentiation, and functional specialization has remained relatively unexplored in plant systems. Single-cell profiling has emerged as a powerful tool to circumvent these past obstacles by enabling unbiased charting of transcriptional and cis-regulatory states at the resolution of individual cells. Here, we review state-of-the-art single-cell approaches and analytical frameworks that have paved the way for establishing the link between cellular phenotypic variation and cis-regulatory mechanisms in plants.

Open chromatin interaction maps reveal functional regulatory elements and chromatin architecture variations during wheat evolution

BACKGROUND

Bread wheat (Triticum aestivum) is an allohexaploid that is generated by two subsequent allopolyploidization events. The large genome size (16 Gb) and polyploid complexity impede our understanding of how regulatory elements and their interactions shape chromatin structure and gene expression in wheat. The open chromatin enrichment and network Hi-C (OCEAN-C) is a powerful antibody-independent method to detect chromatin interactions between open chromatin regions throughout the genome.

RESULTS

Here we generate open chromatin interaction maps for hexaploid wheat and its tetraploid and diploid relatives using OCEAN-C. The anchors of chromatin loops show high chromatin accessibility and are concomitant with several active histone modifications, with 67% of them interacting with multiple loci. Binding motifs of various transcription factors are significantly enriched in the hubs of open chromatin interactions (HOCIs). The genes linked by HOCIs represent higher expression level and lower coefficient expression variance than the genes linked by other loops, which suggests HOCIs may coordinate co-expression of linked genes. Thousands of interchromosomal loops are identified, while limited interchromosomal loops (0.4%) are identified between homoeologous genes in hexaploid wheat. Moreover, we find structure variations contribute to chromatin interaction divergence of homoeologs and chromatin topology changes between different wheat species. The genes with discrepant chromatin interactions show expression alteration in hexaploid wheat compared with its tetraploid and diploid relatives.

CONCLUSIONS

Our results reveal open chromatin interactions in different wheat species, which provide new insights into the role of open chromatin interactions in gene expression during the evolution of polyploid wheat.

Comparison of dynamic 3D chromatin architecture uncovers heterosis for leaf size in Brassica napus

INTRODUCTION

Heterosis is the major event driving plant development and promoting crop breeding, but the molecular bases for this phenomenon remain elusive.

OBJECTIVES

We aim to explore the effect of three-dimensional (3D) chromatin architecture on the underlying mechanism of heterosis.

METHODS

Here, we constructed the North Carolina II (NC-II) population to select superior and inferior heterosis sets by comparing mid-parent heterosis (MPH) in Brassica napus. To decipher the impact of 3D chromatin architecture on the underlying mechanism of heterosis, we combined genetics, transcriptomics and 3D genomics approaches.

RESULTS

We suggest that F1 hybrids with superior heterosis tend to contain more transcriptionally active A compartments compared with F1 hybrids with inferior heterosis, and approximately 19-21% compartment significantly altered in the F1 hybrids relative to the parental lines. Further analyses show that chromatin compartments correlate with genetic variance among parents, which may form the basis for differentially active chromatin compartments. Having more A compartments in F1 hybrids confers a more accessible chromatin circumstance, which promotes a higher proportion of highly expressed ELD (expression level dominance) genes in superior heterosis F1 hybrids (46-64%) compared with inferior heterosis F1 hybrids (22-31%). Moreover, genes related to hormones which affect plant growth, are more up-regulated with changes of 3D genome architecture, and we validate that increased hormone content contributes to cell proliferation and expansion by influencing the key genes of cell cycle thereby promoting leaf size.

CONCLUSION

Dynamic 3D chromatin architecture correlates with genetic variance among parents and contributes to heterosis in Brassica napus.

Mapping open chromatin by ATAC-seq in bread wheat

Gene transcription is largely regulated by cis-regulatory elements. Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is an emerging technology that can accurately map cis-regulatory elements in animals and plants. However, the presence of cell walls and chloroplasts in plants hinders the extraction of high-quality nuclei, thereby affects the quality of ATAC-seq data. Meanwhile, it is tricky to perform ATAC-seq with different tissue types, especially for those with limited size and amount. Moreover, with rapid growth of ATAC-seq datasets from plants, powerful and easy-to-use data analysis pipelines for ATAC-seq, especially for wheat is lacking. Here, we provided an all-in-one solution for mapping open chromatin in wheat including both experimental and data analysis procedure. We efficiently obtained nuclei with less cell debris from various wheat tissues. High-quality ATAC-seq data from young spike and ovary, which are hard to harvest were generated. We determined that the saturation sequencing depth of wheat ATAC-seq is about 16 Gb. Particularly, we developed a powerful and easy-to-use online pipeline to analyze the wheat ATAC-seq data and this pipeline can be easily extended to other plant species. The method developed here will facilitate plant regulatory genome study not only for wheat but also for other plant species.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

Plant 3D Chromatin Organization: Important Insights from Chromosome Conformation Capture Analyses of the Last 10 Years

Over the past few decades, eukaryotic linear genomes and epigenomes have been widely and extensively studied for understanding gene expression regulation. More recently, the three-dimensional (3D) chromatin organization was found to be important for determining genome functionality, finely tuning physiological processes for appropriate cellular responses. With the development of visualization techniques and chromatin conformation capture (3C)-based techniques, increasing evidence indicates that chromosomal architecture characteristics and chromatin domains with different epigenetic modifications in the nucleus are correlated with transcriptional activities. Subsequent studies have further explored the intricate interplay between 3D genome organization and the function of interacting regions. In this review, we summarize spatial distribution patterns of chromatin, including chromatin positioning, configurations and domains, with a particular focus on the effect of a unique form of interaction between varieties of factors that shape the 3D genome conformation in plants. We further discuss the methods, advantages and limitations of various 3C-based techniques, highlighting the applications of these technologies in plants to identify chromatin domains, and address their dynamic changes and functional implications in evolution, and adaptation to development and changing environmental conditions. Moreover, the future implications and emerging research directions of 3D genome organization are discussed.

The arches and spandrels of maize domestication, adaptation, and improvement

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

Modeling the 3D genome of plants

Chromosomes are the carriers of inheritable traits and define cell function and development. This is not only based on the linear DNA sequence of chromosomes but also on the additional molecular information they are associated with, including the transcription machinery, histone modifications, and their three-dimensional folding. The synergistic application of experimental approaches and computer simulations has helped to unveil how these organizational layers of the genome interplay in various organisms. However, such multidisciplinary approaches are still rarely explored in the plant kingdom. Here, we provide an overview of our current knowledge on plant 3D genome organization and review recent efforts to integrate cutting-edge experiments from microscopy and next-generation sequencing approaches with theoretical models. Building on these recent approaches, we propose possible avenues to extend the application of theoretical modeling in the characterization of the 3D genome organization in plants.

An ethylene biosynthesis enzyme controls quantitative variation in maize ear length and kernel yield

Maize ear size and kernel number differ among lines, however, little is known about the molecular basis of ear length and its impact on kernel number. Here, we characterize a quantitative trait locus, qEL7, to identify a maize gene controlling ear length, flower number and fertility. qEL7 encodes 1-aminocyclopropane-1- carboxylate oxidase2 (ACO2), a gene that functions in the final step of ethylene biosynthesis and is expressed in specific domains in developing inflorescences. Confirmation of qEL7 by gene editing of ZmACO2 leads to a reduction in ethylene production in developing ears, and promotes meristem and flower development, resulting in a ~13.4% increase in grain yield per ear in hybrids lines. Our findings suggest that ethylene serves as a key signal in inflorescence development, affecting spikelet number, floral fertility, ear length and kernel number, and also provide a tool to improve grain productivity by optimizing ethylene levels in maize or in other cereals.

Considerable genetic variation exists in maize ear size and kernel number. Here the authors show that variation in a gene encoding an ethylene biosynthetic enzyme impacts ear length, flower fertility and kernel yield suggesting an important role for ethylene signaling during inflorescence development.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

Plant 3D genomics: the exploration and application of chromatin organization

Eukaryotic genomes are highly folded for packing into higher-order chromatin structures in the nucleus. With the emergence of state-of-the-art chromosome conformation capture methods and microscopic imaging techniques, the spatial organization of chromatin and its functional implications have been interrogated. Our knowledge of 3D chromatin organization in plants has improved dramatically in the past few years, building on the early advances in animal systems. Here, we review recent advances in 3D genome mapping approaches, our understanding of the sophisticated organization of spatial structures, and the application of 3D genomic principles in plants. We also discuss directions for future developments in 3D genomics in plants.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Genomic Features of Open Chromatin Regions (OCRs) in Wild Soybean and Their Effects on Gene Expressions

Transcription activation is tightly associated with the openness of chromatin, which allows direct contact between transcriptional regulators, such as transcription factors, and their targeted DNA for downstream gene activation. However, the annotation of open chromatin regions (OCRs) in the wild soybean (Glycine soja) genome is limited. We performed assay for transposase-accessible chromatin using sequencing (ATAC-seq) and successfully identified 22,333 OCRs in the leaf of W05 (a wild soybean accession). These OCRs were enriched in gene transcription start sites (TSS) and were positively correlated with downstream gene expression. Several known transcription factor (TF)-binding motifs were also enriched at the OCRs. A potential regulatory network was constructed using these transcription factors and the OCR-marked genes. Furthermore, by overlapping the OCR distribution with those of histone modifications from chromatin immunoprecipitation followed by sequencing (ChIP-seq), we found that the distribution of the activation histone mark, H3K4me3, but not that of the repressive H3K27me3 mark, was closely associated with OCRs for gene activation. Several putative enhancer-like distal OCRs were also found to overlap with LincRNA-encoding loci. Moreover, our data suggest that homologous OCRs could potentially influence homologous gene expression. Hence, the duplication of OCRs might be essential for plant genome architecture as well as for regulating gene expression.

Rapid and Low-Input Profiling of Histone Marks in Plants Using Nucleus CUT&Tag

Characterizing genome-wide histone posttranscriptional modifications and transcriptional factor occupancy is crucial for deciphering their biological functions. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a powerful method for genome-wide profiling of histone modifications and transcriptional factor-binding sites. However, the current ChIP-seq experimental procedure in plants requires significant material and several days for completion. CUT&Tag is an alternative method of ChIP-seq for low-sample and single-cell epigenomic profiling using protein A-Tn5 transposase fusion proteins (PAT). In this study, we developed a nucleus CUT&Tag (nCUT&Tag) protocol based on the live-cell CUT&Tag technology. Our results indicate that nCUT&Tag could be used for histone modifications profiling in both monocot rice and dicot rapeseed using crosslinked or fresh tissues. In addition, both active and repressive histone marks such as H3K4me3 and H3K9me2 can be identified using our nCUT&Tag. More importantly, all the steps in nCUT&Tag can be finished in only 1 day, and the assay can be performed with as little as 0.01 g of plant tissue as starting materials. Therefore, our results demonstrate that nCUT&Tag is an efficient alternative strategy for plant epigenomic studies.

Improving architectural traits of maize inflorescences

The domestication and improvement of maize resulted in radical changes in shoot architecture relative to its wild progenitor teosinte. In particular, critical modifications involved a reduction of branching and an increase in inflorescence size to meet the needs for human consumption and modern agricultural practices. Maize is a major contributor to global agricultural production by providing large and inexpensive quantities of food, animal feed, and ethanol. Maize is also a classic system for studying the genetic regulation of inflorescence formation and its enlarged female inflorescences directly influence seed production and yield. Studies on the molecular and genetic networks regulating meristem proliferation and maintenance, including receptor-ligand interactions, transcription factor regulation, and hormonal control, provide important insights into maize inflorescence development and reveal potential avenues for the targeted modification of specific architectural traits. In this review, we summarize recent findings on the molecular mechanisms controlling inflorescence formation and discuss how this knowledge can be applied to improve maize productivity in the face of present and future environmental challenges.

Enhancing grain-yield-related traits by CRISPR-Cas9 promoter editing of maize CLE genes

Several yield-related traits selected during crop domestication and improvement^1,2 are associated with increases in meristem size³, which is controlled by CLE peptide signals in the CLAVATA-WUSCHEL pathway^4-13. Here, we engineered quantitative variation for yield-related traits in maize by making weak promoter alleles of CLE genes, and a null allele of a newly identified partially redundant compensating CLE gene, using CRISPR-Cas9 genome editing. These strategies increased multiple maize grain-yield-related traits, supporting the enormous potential for genomic editing in crop enhancement.

Insights into the evo-devo of plant reproduction using next-generation sequencing approaches

The development of plant model organisms has traditionally been analyzed using resource-heavy, tailored applications that are not easily transferable to distantly related non-model taxa. Thus, our understanding of plant development has been limited to a subset of traits, and evolutionary studies conducted most effectively either across very wide [e.g. Arabidopsis thaliana and Oryza sativa (rice)] or narrow (i.e. population level) phylogenetic distances. As plant biologists seek to capitalize on natural diversity for crop improvement, enhance ecosystem functioning, and better understand plant responses to climate change, high-throughput and broadly applicable forms of existing molecular biology assays are becoming an invaluable resource. Next-generation sequencing (NGS) is increasingly becoming a powerful tool in evolutionary developmental biology (evo-devo) studies, particularly through its application to understanding trait evolution at different levels of gene regulation. Here, I review some of the most common and emerging NGS-based methods, using exemplar studies in reproductive plant evo-devo to illustrate their potential.

Predicting Genome Architecture: Challenges and Solutions

Genome architecture plays a pivotal role in gene regulation. The use of high-throughput methods for chromatin profiling and 3-D interaction mapping provide rich experimental data sets describing genome organization and dynamics. These data challenge development of new models and algorithms connecting genome architecture with epigenetic marks. In this review, we describe how chromatin architecture could be reconstructed from epigenetic data using biophysical or statistical approaches. We discuss the applicability and limitations of these methods for understanding the mechanisms of chromatin organization. We also highlight the emergence of new predictive approaches for scoring effects of structural variations in human cells.