Studying physical chromatin interactions in plants using Chromosome Conformation Capture (3C)

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

A high-quality genome assembly reveals adaptations underlying glossy, wax-coated leaves in the heat-tolerant wild raspberry Rubus leucanthus

With glossy, wax-coated leaves, Rubus leucanthus is one of the few heat-tolerant wild raspberry trees. To ascertain the underlying mechanism of heat tolerance, we generated a high-quality genome assembly with a genome size of 230.9 Mb and 24,918 protein-coding genes. Significantly expanded gene families were enriched in the flavonoid biosynthesis pathway and the circadian rhythm-plant pathway, enabling survival in subtropical areas by accumulating protective flavonoids and modifying photoperiodic responses. In contrast, plant-pathogen interaction and MAPK signaling involved in response to pathogens were significantly contracted. The well-known heat response elements (HSP70, HSP90, and HSFs) were reduced in R. leucanthus compared to two other heat-intolerant species, R. chingii and R. occidentalis, with transcriptome profiles further demonstrating their dispensable roles in heat stress response. At the same time, three significantly positively selected genes in the pathway of cuticular wax biosynthesis were identified, and may contribute to the glossy, wax-coated leaves of R. leucanthus. The thick, leathery, waxy leaves protect R. leucanthus against pathogens and herbivores, supported by the reduced R gene repertoire in R. leucanthus (355) compared to R. chingii (376) and R. occidentalis (449). Our study provides some insights into adaptive divergence between R. leucanthus and other raspberry species on heat tolerance.

Chromosome-level genome assembly of Helwingia omeiensis: the first genome in the family Helwingiaceae

Helwingia, a shrub of the monotypic cosmopolitan family Helwingiaceae, is distinguished by its inflorescence, in which flowers are borne on the midrib of the leaf-a trait not commonly observed in related plant families. Previous studies have investigated the development of this unusual structure using comparative anatomical methods. However, the scarcity of genomic data has hindered our understanding of the origins and evolutionary history of this uncommon trait at the molecular level. Here, we report the first high-quality genome of the family Helwingiaceae. Assembled using HiFi sequencing and Hi-C technologies, the genome of H. omeiensis is anchored to 19 chromosomes, with a total length of 2.75 Gb and a contig N50 length of 6.78 Mb. The BUSCO completeness score of the assembled genome was 98.2%. 53,951 genes were identified, of which 99.7% were annotated in at least one protein database. The high-quality reference genome of H. omeiensis provides an essential genetic resource and sheds light on the phylogeny and evolution of specific traits in the family Helwingiaceae.

Efficient purging of deleterious mutations contributes to the survival of a rare conifer

Cupressaceae is a conifer family rich in plants of horticultural importance, including Cupressus, Chamaecyparis, Juniperus, and Thuja, yet genomic surveys are lacking for this family. Cupressus gigantea, one of the many rare conifers that are threatened by climate change and anthropogenic habitat fragmentation, plays an ever-increasing role in ecotourism in Tibet. To infer how past climate change has shaped the population evolution of this species, we generated a de novo chromosome-scale genome (10.92 Gb) and compared the species' population history and genetic load with that of a widespread close relative, C. duclouxiana. Our demographic analyses, based on 83 re-sequenced individuals from multiple populations of the two species, revealed a sharp decline of population sizes during the first part of the Quaternary. However, populations of C. duclouxiana then started to recover, while C. gigantea populations continued to decrease until recently. The total genomic diversity of C. gigantea is smaller than that of C. duclouxiana, but contrary to expectations, C. gigantea has fewer highly and mildly deleterious mutations than C. duclouxiana, and simulations and statistical tests support purifying selection during prolonged inbreeding as the explanation. Our results highlight the evolutionary consequences of decreased population size on the genetic burden of a long-lived endangered conifer with large genome size and suggest that genetic purging deserves more attention in conservation management.

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Genome-wide identification analysis in wild-type Solanum pinnatisectum reveals some genes defending against Phytophthora infestans

Solanum pinnatisectum exhibits strong resistance to late blight caused by Phytophthora infestans but only an incomplete genome assembly based on short Illumina reads has been published. In this study, we generated the first chromosome-level draft genome for the wild-type potato species S. pinnatisectum in China using Oxford Nanopore technology sequencing and Hi-C technology. The high-quality assembled genome size is 664 Mb with a scaffold N50 value of 49.17 Mb, of which 65.87% was occupied by repetitive sequences, and predominant long terminal repeats (42.51% of the entire genome). The genome of S. pinnatisectum was predicted to contain 34,245 genes, of which 99.34% were functionally annotated. Moreover, 303 NBS-coding disease resistance (R) genes were predicted in the S. pinnatisectum genome to investigate the potential mechanisms of resistance to late blight disease. The high-quality chromosome-level reference genome of S. pinnatisectum is expected to provide potential valuable resources for intensively and effectively investigating molecular breeding and genetic research in the future.

A molecular mechanism for embryonic resetting of winter memory and restoration of winter annual growth habit in wheat

The staple food crop winter bread wheat (Triticum aestivum) acquires competence to flower in late spring after experiencing prolonged cold in temperate winter seasons, through the physiological process of vernalization. Prolonged cold exposure results in transcriptional repression of the floral repressor VERNALIZATION 2 (TaVRN2) and activates the expression of the potent floral promoter VERNALIZATION 1 (TaVRN1). Cold-induced TaVRN1 activation and TaVRN2 repression are maintained in post-cold vegetative growth and development, leading to an epigenetic 'memory of winter cold', enabling spring flowering. When and how the cold memory is reset in wheat is essentially unknown. Here we report that the cold-induced TaVRN1 activation is inherited by early embryos, but reset in subsequent embryo development, whereas TaVRN2 remains silenced through seed development, but is reactivated rapidly by light during seed germination. We further found that a chromatin reader mediates embryonic resetting of TaVRN1 and that chromatin modifications play an important role in the regulation of TaVRN1 expression and thus the floral transition, in response to developmental state and environmental cues. The findings define a two-step molecular mechanism for re-establishing vernalization requirement in common wheat, ensuring that each generation must experience winter cold to acquire competence to flower in spring.

A High-Quality Reference Genome Assembly of Prinsepia uniflora (Rosaceae)

This study introduces a meticulously constructed genome assembly at the chromosome level for the Rosaceae family species Prinsepia uniflora, a traditional Chinese medicinal herb. The final assembly encompasses 1272.71 megabases (Mb) distributed across 16 pseudochromosomes, boasting contig and super-scaffold N50 values of 2.77 and 79.32 Mb, respectively. Annotated within this genome is a substantial 875.99 Mb of repetitive sequences, with transposable elements occupying 777.28 Mb, constituting 61.07% of the entire genome. Our predictive efforts identified 49,261 protein-coding genes within the repeat-masked assembly, with 45,256 (91.87%) having functional annotations, 5127 (10.41%) demonstrating tandem duplication, and 2373 (4.82%) classified as transcription factor genes. Additionally, our investigation unveiled 3080 non-coding RNAs spanning 0.51 Mb of the genome sequences. According to our evolutionary study, P. uniflora underwent recent whole-genome duplication following its separation from Prunus salicina. The presented reference-level genome assembly and annotation for P. uniflora will significantly facilitate the in-depth exploration of genomic information pertaining to this species, offering substantial utility in comparative genomics and evolutionary analyses involving Rosaceae species.

Control of histone demethylation by nuclear-localized α-ketoglutarate dehydrogenase

Methylations on nucleosomal histones play fundamental roles in regulating eukaryotic transcription. Jumonji C domain-containing histone demethylases (JMJs) dynamically control the level of histone methylations. However, how JMJ activity is generally regulated is unknown. We found that the tricarboxylic acid cycle-associated enzyme α-ketoglutarate (α-KG) dehydrogenase (KGDH) entered the nucleus, where it interacted with various JMJs to regulate α-KG-dependent histone demethylations by JMJs, and thus controlled genome-wide gene expression in plants. We show that nuclear targeting is regulated by environmental signals and that KGDH is enriched at thousands of loci in Arabidopsis thaliana. Chromatin-bound KGDH catalyzes α-KG decarboxylation and thus may limit its local availability to KGDH-coupled JMJs, inhibiting histone demethylation. Thus, our results uncover a regulatory mechanism for histone demethylations by JMJs.

An upgraded method of high-throughput chromosome conformation capture (Hi-C 3.0) in cotton (Gossypium spp.)

High-throughput chromosome conformation capture (Hi-C) technology has been applied to explore the chromatin interactions and shed light on the biological functions of three-dimensional genomic features. However, it remains challenging to guarantee the high quality of Hi-C library in plants and hence the reliable capture of chromatin structures, especially loops, due to insufficient fragmentation and low efficiency of proximity ligations. To overcome these deficiencies, we optimized the parameters of the Hi-C protocol, principally the cross-linking agents and endonuclease fragmentation strategy. The double cross-linkers (FA+DSG) and double restriction enzymes (DpnII+DdeI) were utilized. Thus, a systematic in situ Hi-C protocol was designed using plant tissues embedded with comprehensive quality controls to monitor the library construction. This upgraded method, termed Hi-C 3.0, was applied to cotton leaves for trial. In comparison with the conventional Hi-C 2.0, Hi-C 3.0 can obtain more than 50% valid contacts at a given sequencing depth to improve the signal-to-noise ratio. Hi-C 3.0 can furthermore enhance the capturing of loops almost as twice as that of Hi-C 2.0. In addition, Hi-C 3.0 showed higher efficiency of compartment detection and identified compartmentalization more accurately. In general, Hi-C 3.0 contributes to the advancement of the Hi-C method in plants by promoting its capability on decoding the chromatin organization.

A high-quality chromosome-level Eutrema salsugineum genome, an extremophile plant model

BACKGROUND

Eutrema salsugineum (2n = 14), a halophyte in the family Brassicaceae, is an attractive model to study abiotic stress tolerance in plants. Two versions of E. salsugineum genomes that previously reported were based on relatively short reads; thus, the repetitive regions were difficult to characterize.

RESULTS

We report the sequencing and assembly of the E. salsugineum (Shandong accession) genome using long-read sequencing and chromosome conformation capture data. We generated Oxford Nanopore long reads at high depth (> 60X) of genome coverage with additional short reads for error correction. The new assembly has a total size of 295.5 Mb with 52.8% repetitive sequences, and the karyotype of E. salsugineum is consistent with the ancestral translocation Proto-Calepineae Karyotype structure in both order and orientation. Compared with previous assemblies, this assembly has higher contiguity, especially in the centromere region. Based on this new assembly, we predicted 25,399 protein-coding genes and identified the positively selected genes associated with salt and drought stress responses.

CONCLUSION

The new genome assembly will provide a valuable resource for future genomic studies and facilitate comparative genomic analysis with other plants.

Insights into the differentiation and adaptation within Circaeasteraceae from Circaeaster agrestis genome sequencing and resequencing

Circaeaster agrestis and Kingdonia uniflora are sister species that reproduce sexually and mainly asexually respectively, providing a good system for comparative genome evolution between taxa with different reproductive models. Comparative genome analyses revealed the two species have similar genome size, but C. agrestis encodes many more genes. The gene families specific to C. agrestis show significant enrichment of genes associated with defense response, while those gene families specific to K. uniflora are enriched in genes regulating root system development. Collinearity analyses revealed C. agrestis experienced two rounds of whole-genome duplication. Fst outlier test across 25 C. agrestis populations uncovered a close inter-relationship between abiotic stress and genetic variability. Genetic feature comparisons showed K. uniflora presents much higher genome heterozygosity, transposable element load, linkage disequilibrium degree, and π_N/π_S ratio. This study provides new insights into understanding the genetic differentiation and adaptation within ancient lineages characterized by multiple reproductive models.

Chromosome-Level Genome Assembly of Herpetospermum pedunculosum (Cucurbitaceae)

This study presents a chromosome-level reference genome assembly of a traditional Tibetan medicinal plant, Herpetospermum pedunculosum belonging to the Cucurbitaceae family. Following a combined PacBio high-fidelity sequencing and Hi-C analysis, a final H. pedunculosum genome assembly, 804.11 Mb in length was obtained, 90.45% of which was anchored into ten pseudochromosomes with a contig N50 of 24.39 Mb. In addition, 579.55 Mb repetitive sequences and 23,924 high-confidence protein-coding genes were annotated. Phylogenetic analysis revealed that H. pedunculosum was sister to a clade formed by cucumber, zucchini, and wax gourd. Further whole-genome duplication analysis revealed no recent polyploidization event in the H. pedunculosum genome. The high-quality H. pedunculosum genome presented here will be highly useful in investigating the molecular mechanisms underlying the biosynthesis of its active compounds and adaptation strategies to the extreme environment. It will also provide great insights into comparative genomic studies of Cucurbitaceae and flowering plants.

Karyotype evolution of the Asterids insights from the first genome sequences of the family Cornaceae

Cornaceae is a core representative family in Cornales, the earliest branching lineage in the Asterids on the life tree of angiosperms. This family includes the only genus Cornus, a group of ~55 species. These species occur widely in Northern Hemisphere and have been used as resources for horticultural ornaments, medicinal and industrial manufacturing. However, no any genome sequences are available for this family. Here, we reported a chromosome-level genome for Cornus controversa. This was generated using high-fidelity plus Hi-C sequencing, and totally ~771.80 Mb assembled sequences and 39,886 protein-coding genes were obtained. We provided evidence for a whole-genome duplication event (WGD) unique to C. controversa. The evolutionary features of this genome indicated that the expanded and unique genes might have contributed to response to stress, stimulus and defense. By using chromosome-level syntenic blocks shared between eight living genomes, we found high degrees of genomic diversification from the ancestral core-eudicot genome to the present-day genomes, suggesting an important role of WGD in genomic plasticity that leads to speciation and diversification. These results provide foundational insights on the evolutionary history of Cornaceae, as well as on the Asterids diversification.

Improved assembly and annotation of the sesame genome

Sesame (Sesamum indicum L.) is an important oilseed crop that produces abundant seed oil and has a pleasant flavor and high nutritional value. To date, several Illumina-based genome assemblies corresponding to different sesame genotypes have been published and widely used in genetic and genomic studies of sesame. However, these assemblies consistently showed low continuity with numerous gaps. Here, we reported a high-quality, reference-level sesame genome assembly by integrating PacBio high-fidelity sequencing and Hi-C technology. Our updated sesame assembly was 309.35 Mb in size with a high chromosome anchoring rate (97.54%) and contig N50 size (13.48 Mb), which were better than previously published genomes. We identified 163.38 Mb repetitive elements and 24,345 high-confidence protein-coding genes in the updated sesame assembly. Comparative genomic analysis showed that sesame shared an ancient whole-genome duplication event with two Lamiales species. A total of 2,782 genes were tandemly duplicated. We also identified several genes that were likely involved in fatty acid and triacylglycerol biosynthesis. Our improved sesame assembly and annotation will facilitate future genetic studies and genomics-assisted breeding of sesame.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

Complex chromosomal rearrangements induced by transposons in maize

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

Genomes and demographic histories of the endangered Bretschneidera sinensis (Akaniaceae)

BACKGROUND

Bretschneidera sinensis is an endangered relic tree species in the Akaniaceae family and is sporadically distributed in eastern Asia. As opposed to its current narrow and rare distribution, the fossil pollen of B. sinensis has been found to be frequent and widespread in the Northern Hemisphere during the Late Miocene. B. sinensis is also a typical mycorrhizal plant, and its annual seedlings exhibit high mortality rates in absence of mycorrhizal development. The chromosome-level high-quality genome of B. sinensis will help us to more deeply understand the survival and demographic histories of this relic species.

RESULTS

A total of 25.39 Gb HiFi reads and 109.17 Gb Hi-C reads were used to construct the chromosome-level genome of B. sinensis, which is 1.21 Gb in length with the contig N50 of 64.13 Mb and chromosome N50 of 146.54 Mb. The identified transposable elements account for 55.21% of the genome. A total of 45,839 protein-coding genes were predicted in B. sinensis. A lineage-specific whole-genome duplication was detected, and 7,283 lineage-specific expanded gene families with functions related to the specialized endotrophic mycorrhizal adaptation were identified. The historical effective population size (Ne) of B. sinensis was found to oscillate greatly in response to Quaternary climatic changes. The Ne of B. sinensis has decreased rapidly in the recent past, making its extant Ne extremely lower. Our additional evolutionary genomic analyses suggested that the developed mycorrhizal adaption might have been repeatedly disrupted by environmental changes caused by Quaternary climatic oscillations. The environmental changes and an already decreased population size during the Holocene may have led to the current rarity of B. sinensis.

CONCLUSION

This is a detailed report of the genome sequences for the family Akaniaceae distributed in evergreen forests in eastern Asia. Such a high-quality genomic resource may provide critical clues for comparative genomics studies of this family in the future.

The lncRNA MARS modulates the epigenetic reprogramming of the marneral cluster in response to ABA

Clustered organization of biosynthetic non-homologous genes is emerging as a characteristic feature of plant genomes. The co-regulation of clustered genes seems to largely depend on epigenetic reprogramming and three-dimensional chromatin conformation. In this study, we identified the long non-coding RNA (lncRNA) MARneral Silencing (MARS), localized inside the Arabidopsis marneral cluster, which controls the local epigenetic activation of its surrounding region in response to abscisic acid (ABA). MARS modulates the POLYCOMB REPRESSIVE COMPLEX 1 (PRC1) component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) binding throughout the cluster in a dose-dependent manner, determining H3K27me3 deposition and chromatin condensation. In response to ABA, MARS decoys LHP1 away from the cluster and promotes the formation of a chromatin loop bringing together the MARNERAL SYNTHASE 1 (MRN1) locus and a distal ABA-responsive enhancer. The enrichment of co-regulated lncRNAs in clustered metabolic genes in Arabidopsis suggests that the acquisition of novel non-coding transcriptional units may constitute an additional regulatory layer driving the evolution of biosynthetic pathways.

Chromosome-level genome assembly of a xerophytic plant, Haloxylon ammodendron

Haloxylon ammodendron is a xerophytic perennial shrub or small tree that has a high ecological value in anti-desertification due to its high tolerance to drought and salt stress. Here, we report a high-quality, chromosome-level genome assembly of H. ammodendron by integrating PacBio's high-fidelity sequencing and Hi-C technology. The assembled genome size was 685.4 Mb, of which 99.6% was assigned to nine pseudochromosomes with a contig N50 value of 23.6 Mb. Evolutionary analysis showed that both the recent substantial amplification of long terminal repeat retrotransposons and tandem gene duplication may have contributed to its genome size expansion and arid adaptation. An ample amount of low-GC genes was closely related to functions that may contribute to the desert adaptation of H. ammodendron. Gene family clustering together with gene expression analysis identified differentially expressed genes that may play important roles in the direct response of H. ammodendron to water-deficit stress. We also identified several genes possibly related to the degraded scaly leaves and well-developed root system of H. ammodendron. The reference-level genome assembly presented here will provide a valuable genomic resource for studying the genome evolution of xerophytic plants, as well as for further genetic breeding studies of H. ammodendron.

Stress Conditions Modulate the Chromatin Interactions Network in Arabidopsis

Stresses have been known to cause various responses like cellular physiology, gene regulation, and genome remodeling in the organism to cope and survive. Here, we assessed the impact of stress conditions on the chromatin-interactome network of Arabidopsis thaliana. We identified thousands of chromatin interactions in native as well as in salicylic acid treatment and high temperature conditions in a genome-wide fashion. Our analysis revealed the definite pattern of chromatin interactions and stress conditions could modulate the dynamics of chromatin interactions. We found the heterochromatic region of the genome actively involved in the chromatin interactions. We further observed that the establishment or loss of interactions in response to stress does not result in the global change in the expression profile of interacting genes; however, interacting regions (genes) containing motifs for known TFs showed either lower expression or no difference than non-interacting genes. The present study also revealed that interactions preferred among the same epigenetic state (ES) suggest interactions clustered the same ES together in the 3D space of the nucleus. Our analysis showed that stress conditions affect the dynamics of chromatin interactions among the chromatin loci and these interaction networks govern the folding principle of chromatin by bringing together similar epigenetic marks.

A chromosome-scale genome assembly for the holly (Ilex polyneura) provides insights into genomic adaptations to elevation in Southwest China

Southwest China is a plant diversity hotspot. The near-cosmopolitan genus Ilex (c. 664 spp., Aquifoliaceae) reaches its maximum diversity in this region, with many narrow-range and a few widespread species. Divergent selection on widespread species leads to local adaptation, with consequences for both conservation and utilization, but is counteracted by geneflow. Many Ilex species are utilized as teas, medicines, ornamentals, honey plants, and timber, but variation below the species level is largely uninvestigated. We therefore studied the widespread Ilex polyneura, which occupies most of the elevational range available and is cultivated for its decorative leafless branches with persistent red fruits. We assembled a chromosome-scale genome using approximately 100x whole genome long-read and short-read sequencing combined with Hi-C sequencing. The genome is approximately 727.1 Mb, with a contig N50 size of 5 124 369 bp and a scaffold N50 size of 36 593 620 bp, for which the BUSCO score was 97.6%, and 98.9% of the assembly was anchored to 20 pseudochromosomes. Out of 32 838 genes predicted, 96.9% were assigned functions. Two whole genome duplication events were identified. Using this genome as a reference, we conducted a population genomics study of 112 individuals from 21 populations across the elevation range using restriction site-associated DNA sequencing (RADseq). Most populations clustered into four clades separated by distance and elevation. Selective sweep analyses identified 34 candidate genes potentially under selection at different elevations, with functions related to responses to abiotic and biotic stresses. This first high-quality genome in the Aquifoliales will facilitate the further domestication of the genus.

Interrogating Global Chromatin Interaction Network by High-Throughput Chromosome Conformation Capture (Hi-C) in Plants

High-throughput chromosome conformation capture (Hi-C) enables the global quantification of chromatin interaction frequency in eukaryotic nuclei. This method is based on in situ Hi-C, in which chromatin is cross-linked with formaldehyde, then digested with restriction enzyme. Biotin-labeled nucleotide is incorporated before the spatially adjacent DNA ends are ligated, making it possible to enrich specifically the chimeric ligation products for deep sequencing. In this chapter, we describe a modified in situ Hi-C protocol for the global chromatin interaction analysis in plants.

Chromatin spatial organization of wild type and mutant peanuts reveals high-resolution genomic architecture and interaction alterations

BACKGROUND

Three-dimensional (3D) chromatin organization provides a critical foundation to investigate gene expression regulation and cellular homeostasis.

RESULTS

Here, we present the first 3D genome architecture maps in wild type and mutant allotetraploid peanut lines, which illustrate A/B compartments, topologically associated domains (TADs), and widespread chromatin interactions. Most peanut chromosomal arms (52.3%) have active regions (A compartments) with relatively high gene density and high transcriptional levels. About 2.0% of chromosomal regions switch from inactive to active (B-to-A) in the mutant line, harboring 58 differentially expressed genes enriched in flavonoid biosynthesis and circadian rhythm functions. The mutant peanut line shows a higher number of genome-wide cis-interactions than its wild-type. The present study reveals a new TAD in the mutant line that generates different chromatin loops and harbors a specific upstream AP2EREBP-binding motif which might upregulate the expression of the GA2ox gene and decrease active gibberellin (GA) content, presumably making the mutant plant dwarf.

CONCLUSIONS

Our findings will shed new light on the relationship between 3D chromatin architecture and transcriptional regulation in plants.

High-quality genome assembly of an important biodiesel plant, Euphorbia lathyris L

Caper spurge, Euphorbia lathyris L., is an important energy crop and medicinal crop. Here, we generated a high-quality, chromosome-level genome assembly of caper spurge using Oxford Nanopore sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼988.9 Mb in size, 99.8% of which could be grouped into 10 pseudochromosomes, with contig and scaffold N50 values of 32.6 and 95.7 Mb, respectively. A total of 651.4 Mb repetitive sequences and 36,342 protein-coding genes were predicted in the genome assembly. Comparative genomic analysis showed that caper spurge and castor bean clustered together. We found that no independent whole-genome duplication event had occurred in caper spurge after its split from castor bean, and recent substantial amplification of LTR retrotransposons (LTR-RTs) has contributed significantly to its genome expansion. Furthermore, based on gene homology searching, we identified a number of candidate genes involved in the biosynthesis of fatty acids and triacylglycerols. The reference genome presented here will be highly useful for the further study of the genetics, genomics, and breeding of this high-value crop, as well as for evolutionary studies of spurge family and angiosperms.

The spatial position effect: synthetic biology enters the era of 3D genomics

Microbial cell factories are critical to achieving green biomanufacturing. A position effect occurs when a synthetic gene circuit is expressed from different positions in the chassis strain genome. Here, we propose the concept of the 'spatial position effect,' which uses technologies in 3D genomics to reveal the spatial structure characteristics of the 3D genome of the chassis. On this basis, we propose to rationally design the integration sites of synthetic gene circuits, use reporter genes for preliminary screening, and integrate synthetic gene circuits into promising sites for further experiments. This approach can produce stable and efficient chassis strains for green biomanufacturing. The proposed spatial position effect brings synthetic biology into the era of 3D genomics.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

The histone variant H3.3 promotes the active chromatin state to repress flowering in Arabidopsis

The histone H3 family in animals and plants includes replicative H3 and nonreplicative H3.3 variants. H3.3 preferentially associates with active transcription, yet its function in development and transcription regulation remains elusive. The floral transition in Arabidopsis (Arabidopsis thaliana) involves complex chromatin regulation at a central flowering repressor FLOWERING LOCUS C (FLC). Here, we show that H3.3 upregulates FLC expression and promotes active histone modifications histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 36 trimethylation (H3K36me3) at the FLC locus. The FLC activator FRIGIDA (FRI) directly mediates H3.3 enrichment at FLC, leading to chromatin conformation changes and further induction of active histone modifications at FLC. Moreover, the antagonistic H3.3 and H2A.Z act in concert to activate FLC expression, likely by forming unstable nucleosomes ideal for transcription processing. We also show that H3.3 knockdown leads to H3K4me3 reduction at a subset of particularly short genes, suggesting the general role of H3.3 in promoting H3K4me3. The finding that H3.3 stably accumulates at FLC in the absence of H3K36me3 indicates that the H3.3 deposition may serve as a prerequisite for active histone modifications. Our results reveal the important function of H3.3 in mediating the active chromatin state for flowering repression.

The chromosome-level reference genome assembly for Dendrobium officinale and its utility of functional genomics research and molecular breeding study

Dendrobium officinale, an important medicinal plant of the genus Dendrobium in Orchidaceae family, has been used as traditional Chinese medicine (TCM) for nearly thousands of years. Here, we report the first chromosome-level reference genome of D. officinale, based on PacBio long-reads, Illumina short-reads and Hi-C data. The high-quality assembled genome is 1.23 Gb long, with contig N50 of 1.44 Mb. A total of 93.53% genome sequences were assembled into 19 pseudochromosomes with a super scaffold N50 of 63.07 Mb. Through comparative genomic analysis, we explored the expanded gene families of D. officinale, and also their impact on environmental adaptation and biosynthesis of secondary metabolites. We further performed detailed transcriptional analysis of D. officinale, and identified the candidate genes involved in the biosynthesis of three main active ingredients, including polysaccharides, alkaloids and flavonoids. In addition, the MODIFYING WALL LIGNIN-1 (MWL1) gene, which inferred from Genome-Wide Association Studies (GWAS) based on the resequencing date from D. officinale and five related species and their morphologic features, may contribute to the plant production (yield of stems) of D. officinale. Therefore, the high-quality reference genome reported in this study could benefits functional genomics research and molecular breeding of D. officinale.

Chromosome-level genome assembly of Sichuan pepper provides insights into apomixis, drought tolerance, and alkaloid biosynthesis

Sichuan pepper is a commonly used spice in Asian cuisine. Sanshools and wgx-50/gx-50 isolated from it have been shown to possess a wide spectrum of medicinal properties. Here we generated a chromosome-level genome assembly of one Sichuan pepper species Zanthoxylum armatum characterized by drought tolerance and apomixis. Analyses of functionally related genes suggested that increased gene copy number and expression level of drought-tolerant genes might play an important role in improving drought tolerance of Z. armatum. Moreover, a gene encoding an RWP-RK domain-containing protein was shown to contribute to apomixis in Z. armatum, which was further characterized by overexpression in Arabidopsis thaliana. Furthermore, based on gene homology searching and co-expression patterns of metabolite concentration and gene expressions, we identified a number of candidate genes involved in the biosynthesis of sanshools and wgx-50/gx-50. Taken together, our results yield valuable insights for understanding the evolution of apomixis, drought tolerance, and alkaloid biosynthesis in Z. armatum.

Chromosomal-Level Reference Genome of the Neotropical Tree Jacaranda mimosifolia D. Don

Jacaranda mimosifolia D. Don is a deciduous tree widely cultivated in the tropics and subtropics of the world. It is famous for its beautiful blue flowers and pinnate compound leaves. In addition, this tree has great potential in environmental monitoring, soil quality improvement, and medicinal applications. However, a genome resource for J. mimosifolia has not been reported to date. In this study, we constructed a chromosome-level genome assembly of J. mimosifolia using PacBio sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼707.32 Mb in size, 688.76 Mb (97.36%) of which could be grouped into 18 pseudochromosomes, with contig and scaffold N50 values of 16.77 and 39.98 Mb, respectively. A total of 30,507 protein-coding genes were predicted, 95.17% of which could be functionally annotated. Phylogenetic analysis among 12 plant species confirmed the close genetic relationship between J. mimosifolia and Handroanthus impetiginosus. Gene family clustering revealed 481 unique, 103 significantly expanded, and 16 significantly contracted gene families in the J. mimosifolia genome. This chromosome-level genome assembly of J. mimosifolia will provide a valuable genomic resource for elucidating the genetic bases of the morphological characteristics, adaption evolution, and active compounds biosynthesis of J. mimosifolia.

Chromatin architectural proteins regulate flowering time by precluding gene looping

Chromatin structure is critical for gene expression and many other cellular processes. In Arabidopsis thaliana, the floral repressor FLC adopts a self-loop chromatin structure via bridging of its flanking regions. This local gene loop is necessary for active FLC expression. However, the molecular mechanism underlying the formation of this class of gene loops is unknown. Here, we report the characterization of a group of linker histone-like proteins, named the GH1-HMGA family in Arabidopsis, which act as chromatin architecture modulators. We demonstrate that these family members redundantly promote the floral transition through the repression of FLC A genome-wide study revealed that this family preferentially binds to the 5' and 3' ends of gene bodies. The loss of this binding increases FLC expression by stabilizing the FLC 5' to 3' gene looping. Our study provides mechanistic insights into how a family of evolutionarily conserved proteins regulates the formation of local gene loops.

Inferring Single-Cell 3D Chromosomal Structures Based on the Lennard-Jones Potential

Reconstructing three-dimensional (3D) chromosomal structures based on single-cell Hi-C data is a challenging scientific problem due to the extreme sparseness of the single-cell Hi-C data. In this research, we used the Lennard-Jones potential to reconstruct both 500 kb and high-resolution 50 kb chromosomal structures based on single-cell Hi-C data. A chromosome was represented by a string of 500 kb or 50 kb DNA beads and put into a 3D cubic lattice for simulations. A 2D Gaussian function was used to impute the sparse single-cell Hi-C contact matrices. We designed a novel loss function based on the Lennard-Jones potential, in which the ε value, i.e., the well depth, was used to indicate how stable the binding of every pair of beads is. For the bead pairs that have single-cell Hi-C contacts and their neighboring bead pairs, the loss function assigns them stronger binding stability. The Metropolis-Hastings algorithm was used to try different locations for the DNA beads, and simulated annealing was used to optimize the loss function. We proved the correctness and validness of the reconstructed 3D structures by evaluating the models according to multiple criteria and comparing the models with 3D-FISH data.

JASMONATE-ZIM DOMAIN proteins engage Polycomb chromatin modifiers to modulate Jasmonate signaling in Arabidopsis

Jasmonate (JA) regulates various aspects of plant growth and development and stress responses, with prominent roles in male reproductive development and defenses against herbivores and necrotrophic pathogens. JASMONATE-ZIM DOMAIN (JAZ) proteins are key regulators in the JA signaling pathway and function to repress the expression of JA-responsive genes. Here, we show that JAZ proteins directly interact with several chromatin-associated Polycomb proteins to mediate repressive chromatin modifications at JA-responsive genes and, thus, their transcriptional repression in Arabidopsis. Genetic analyses revealed that the developmental defects, including anther and pollen abnormalities, resulting from loss or block of JA signaling were partially rescued by loss of Polycomb protein-mediated chromatin silencing (Polycomb repression). We further found that JAZ-mediated transcriptional repression during anther and pollen development requires Polycomb proteins at four key regulatory loci. Analysis of genome-wide occupancy of a Polycomb factor and transcriptome reprogramming in response to JA revealed that Polycomb repression is involved in the repression of various JA-responsive genes. Taken together, our study reveals an important chromatin-based mechanism for JAZ-mediated transcriptional repression and JA signaling in plants.

A chromosome-level genome assembly for the tertiary relict plant Tetracentron sinense oliv. (trochodendraceae)

Tetracentron sinense and Trochodendron aralioides are two Tertiary relict species of large trees in the family Trochodendraceae with narrow distributions on the mainland and islands of eastern Asia. They belong to the order Trochodendrales, which is one of the four early-diverged eudicot lineages. These two relict species provide a good system in which to examine genomic changes that occurred as they survived during repeated climatic oscillations in the Quaternary. We sequenced the genome of Te. sinense and compared it with that of Tr. aralioides. We found that Te. sinense has a smaller genome size (986.3 Mb) than that of Tr. aralioides (1610 Mb). Repetitive elements made the major contribution to the contrasting genome sizes in the two species, with most bursts of repeats occurring within the past four million years when the climate oscillated greatly. These species share two rounds of whole-genome duplications. The mainland species Te. sinense had a larger effective population size than the island species Tr. aralioides after the largest glaciation during the Quaternary climatic oscillation. However, soon after this recovery stage, the effective population sizes of both species continued to decrease, although the current effective population size of Te. sinense is still larger than that of Tr. aralioides. We recovered three distinctly diverged clades through resequencing the genomes of 50 individuals across the distributional range of Te. sinense in China. Our results provide an important genomic resource with which to examine early trait evolution in the core eudicots and assist efforts to conserve this relict tree species.

TOP1α, UPF1, and TTG2 regulate seed size in a parental dosage-dependent manner

Cues of maternal and paternal origins interact to control seed development, and the underlying molecular mechanisms are still far from clear. Here, we show that TOPOISOMERASE Iα (TOP1α), UP-FRAMESHIFT SUPPRESSOR 1 (UPF1), and TRANSPARENT TESTA GLABRA2 (TTG2) gametophytically, biparentally regulate seed size in Arabidopsis. TOP1α and UPF1 are mainly expressed in antipodal cells, and loss of their function leads to ectopic TTG2 expression in these female gametophytic cells. We further demonstrate that TOP1α and UPF1 directly repress TTG2 expression through affecting its chromatin status and determine its relative expression in antipodal cells versus sperm cells, which controls seed size in a dosage-dependent and parent-of-origin-dependent manner. The molecular interplay among these three genes explains their biparental gametophytic effect during diploidy and interploidy reciprocal crosses. Taken together, our findings reveal a molecular framework of parental interaction for seed size control.

Cues of maternal and paternal origin interact to control seed development, and the underlying molecular mechanisms are still far from clear. This study shows that in Arabidopsis, the relative dosage of the transcription factor TTG2 between antipodal cells and sperm cells at the beginning of seed development determines seed size under the control of TOP1α and UPF1.

Single-Cell Genomics and Epigenomics: Technologies and Applications in Plants

The development of genomics and epigenomics has allowed rapid advances in our understanding of plant biology. However, conventional bulk analysis dilutes cell-specific information by providing only average information, thereby limiting the resolution of genomic and functional genomic studies. Recent advances in single-cell sequencing technology concerning genomics and epigenomics open new avenues to dissect cell heterogeneity in multiple biological processes. Recent applications of these approaches to plants have provided exciting insights into diverse biological questions. We highlight the methodologies underlying the current techniques of single-cell genomics and epigenomics before covering their recent applications, potential significance, and future perspectives in plant biology.

Nuclear role for human Argonaute-1 as an estrogen-dependent transcription coactivator

In mammals, argonaute (AGO) proteins have been characterized for their roles in small RNA-mediated posttranscriptional and also in transcriptional gene silencing. Here, we report a different role for AGO1 in estradiol-triggered transcriptional activation in human cells. We show that in MCF-7 mammary gland cells, AGO1 associates with transcriptional enhancers of estrogen receptor α (ERα) and that this association is up-regulated by treating the cells with estrogen (E2), displaying a positive correlation with the activation of these enhancers. Moreover, we show that AGO1 interacts with ERα and that this interaction is also increased by E2 treatment, but occurs in the absence of RNA. We show that AGO1 acts positively as a coactivator in estradiol-triggered transcription regulation by promoting ERα binding to its enhancers. Consistently, AGO1 depletion decreases long-range contacts between ERα enhancers and their target promoters. Our results point to a role of AGO1 in transcriptional regulation in human cells that is independent from small RNA binding.

Heat stress-induced transposon activation correlates with 3D chromatin organization rearrangement in Arabidopsis

In higher eukaryotes, heterochromatin is mainly composed of transposable elements (TEs) silenced by epigenetic mechanisms. But, the silencing of certain heterochromatin-associated TEs is disrupted by heat stress. By comparing genome-wide high-resolution chromatin packing patterns under normal or heat conditions obtained through Hi-C analysis, we show here that heat stress causes global rearrangement of the 3D genome in Arabidopsis thaliana. Contacts between pericentromeric regions and distal chromosome arms, as well as proximal intra-chromosomal interactions along the chromosomes, are enhanced. However, interactions within pericentromeres and those between distal intra-chromosomal regions are decreased. Many inter-chromosomal interactions, including those within the KNOT, are also reduced. Furthermore, heat activation of TEs exhibits a high correlation with the reduction of chromosomal interactions involving pericentromeres, the KNOT, the knob, and the upstream and downstream flanking regions of the activated TEs. Together, our results provide insights into the relationship between TE activation and 3D genome reorganization.

Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement

Polyploidy is an evolutionary innovation for many animals and all flowering plants, but its impact on selection and domestication remains elusive. Here we analyze genome evolution and diversification for all five allopolyploid cotton species, including economically important Upland and Pima cottons. Although these polyploid genomes are conserved in gene content and synteny, they have diversified by subgenomic transposon exchanges that equilibrate genome size, evolutionary rate heterogeneities and positive selection between homoeologs within and among lineages. These differential evolutionary trajectories are accompanied by gene-family diversification and homoeolog expression divergence among polyploid lineages. Selection and domestication drive parallel gene expression similarities in fibers of two cultivated cottons, involving coexpression networks and N⁶-methyladenosine RNA modifications. Furthermore, polyploidy induces recombination suppression, which correlates with altered epigenetic landscapes and can be overcome by wild introgression. These genomic insights will empower efforts to manipulate genetic recombination and modify epigenetic landscapes and target genes for crop improvement.

Sequencing and genomic diversification of five allopolyploid cotton species provide insights into polyploid genome evolution and epigenetic landscapes for cotton improvement.

Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

Populus euphratica is well adapted to extreme desert environments and is an important model species for elucidating the mechanisms of abiotic stress resistance in trees. The current assembly of P. euphratica genome is highly fragmented with many gaps and errors, thereby impeding downstream applications. Here, we report an improved chromosome-level reference genome of P. euphratica (v2.0) using single-molecule sequencing and chromosome conformation capture (Hi-C) technologies. Relative to the previous reference genome, our assembly represents a nearly 60-fold improvement in contiguity, with a scaffold N50 size of 28.59 Mb. Using this genome, we have found that extensive expansion of Gypsy elements in P. euphratica led to its rapid increase in genome size compared to any other Salicaceae species studied to date, and potentially contributed to adaptive divergence driven by insertions near genes involved in stress tolerance. We also detected a wide range of unique structural rearrangements in P. euphratica, including 2,549 translocations, 454 inversions, 121 tandem and 14 segmental duplications. Several key genes likely to be involved in tolerance to abiotic stress were identified within these regions. This high-quality genome represents a valuable resource for poplar breeding and genetic improvement in the future, as well as comparative genomic analysis with other Salicaceae species.

A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: An Isatis genome

Isatis indigotica (2n = 14) is an important medicinal plant in China. Its dried leaves and roots (called Isatidis Folium and Isatidis Radix, respectively) are broadly used in traditional Chinese medicine for curing diseases caused by bacteria and viruses such as influenza and viral pneumonia. Various classes of compounds isolated from this species have been identified as effective ingredients. Previous studies based on transcriptomes revealed only a few candidate genes for the biosynthesis of these active compounds in this medicinal plant. Here, we report a high-quality chromosome-scale genome assembly of I. indigotica with a total size of 293.88 Mb and scaffold N50 = 36.16 Mb using single-molecule real-time long reads and high-throughput chromosome conformation capture techniques. We annotated 30,323 high-confidence protein-coding genes. Based on homolog searching and functional annotations, we identified many candidate genes involved in the biosynthesis of main active components such as indoles, terpenoids, and phenylpropanoids. In addition, we found that some key enzyme-coding gene families related to the biosynthesis of these components were expanded due to tandem duplications, which likely drove the production of these major active compounds and explained why I. indigotica has excellent antibacterial and antiviral activities. Our results highlighted the importance of genome sequencing in identifying candidate genes for metabolite synthesis in medicinal plants.

The tumor suppressive effect of long non-coding RNA FRMD6-AS2 in uteri corpus endometrial carcinoma

Uterine corpus endometrial cancer (UCEC) is one of the most common gynecological malignancies with increasing incidence and high morbidity and mortality. The currently acknowledged molecular mechanism of UCEC is still not adequate. Here, we reported that the expression of a novel long non-coding RNA (lncRNA) FRMD6-AS2 was reduced in UCEC compared to noncancerous endometrium tissues using the data from The Cancer Genome Atlas (TCGA) Project database. The gene ontology (GO) analysis on differential expressed targeted genes of FRMD6-AS2 in UCEC suggested that FRMD6-AS2 might impact with the function of actin-mediated cell movement and contraction. By over-expressing FRMD6-AS2 in UCEC cell lines, we observed that FRMD6-AS2 played a suppressive role in tumor growth, migration and invasion via activation of Hippo signaling pathway including FRMD6. Moreover, we also demonstrated that FRMD6-AS2 could interact with the 30 kb upstream beyond FRMD6 and facilitate the chromatin looping towards the promoter region of FRMD6 to enhance the expression of FRMD6. We concluded that lncRNA FRMD6-AS2 repressed UCEC, at least in part, by increasing FRMD6.

Back-spliced RNA from retrotransposon binds to centromere and regulates centromeric chromatin loops in maize

In most plants, centromeric DNA contains highly repetitive sequences, including tandem repeats and retrotransposons; however, the roles of these sequences in the structure and function of the centromere are unclear. Here, we found that multiple RNA sequences from centromeric retrotransposons (CRMs) were enriched in maize (Zea mays) centromeres, and back-spliced RNAs were generated from CRM1. We identified 3 types of CRM1-derived circular RNAs with the same back-splicing site based on the back-spliced sequences. These circular RNAs bound to the centromere through R-loops. Two R-loop sites inside a single circular RNA promoted the formation of chromatin loops in CRM1 regions. When RNA interference (RNAi) was used to target the back-splicing site of the circular CRM1 RNAs, the levels of R-loops and chromatin loops formed by these circular RNAs decreased, while the levels of R-loops produced by linear RNAs with similar binding sites increased. Linear RNAs with only one R-loop site could not promote chromatin loop formation. Higher levels of R-loops and lower levels of chromatin loops in the CRM1 regions of RNAi plants led to a reduced localization of the centromeric H3 variant (CENH3). Our work reveals centromeric chromatin organization by circular CRM1 RNAs via R-loops and chromatin loops, which suggested that CRM1 elements might help build a suitable chromatin environment during centromere evolution. These results highlight that R-loops are integral components of centromeric chromatin and proper centromere structure is essential for CENH3 localization.

Back-spliced RNA from retrotransposon binds to centromere and regulates centromeric chromatin loops in maize

In most plants, centromeric DNA contains highly repetitive sequences, including tandem repeats and retrotransposons; however, the roles of these sequences in the structure and function of the centromere are unclear. Here, we found that multiple RNA sequences from centromeric retrotransposons (CRMs) were enriched in maize (Zea mays) centromeres, and back-spliced RNAs were generated from CRM1. We identified 3 types of CRM1-derived circular RNAs with the same back-splicing site based on the back-spliced sequences. These circular RNAs bound to the centromere through R-loops. Two R-loop sites inside a single circular RNA promoted the formation of chromatin loops in CRM1 regions. When RNA interference (RNAi) was used to target the back-splicing site of the circular CRM1 RNAs, the levels of R-loops and chromatin loops formed by these circular RNAs decreased, while the levels of R-loops produced by linear RNAs with similar binding sites increased. Linear RNAs with only one R-loop site could not promote chromatin loop formation. Higher levels of R-loops and lower levels of chromatin loops in the CRM1 regions of RNAi plants led to a reduced localization of the centromeric H3 variant (CENH3). Our work reveals centromeric chromatin organization by circular CRM1 RNAs via R-loops and chromatin loops, which suggested that CRM1 elements might help build a suitable chromatin environment during centromere evolution. These results highlight that R-loops are integral components of centromeric chromatin and proper centromere structure is essential for CENH3 localization.

Molecular Basis of Natural Variation in Photoperiodic Flowering Responses

Plants exhibit different flowering behaviors in response to variable photoperiods across a wide geographical range. Here, we identify MYC3, a bHLH transcription factor, and its cis-element form the long-sought regulatory module responsible for cis-regulatory changes at the florigen gene FLOWERING LOCUS T (FT) that mediate natural variation in photoperiodic flowering responses in Arabidopsis. MYC3 is stabilized by DELLAs in the gibberellin pathway to suppress FT through binding the ACGGAT motif and antagonizing CONSTANS (CO) activation. Changing photoperiods modulate the relative abundance of MYC3 and CO, thus determining either of them as the predominant regulator for FT expression under different day lengths. Cis-regulatory changes in the MYC3 binding site at FT are associated with natural variation in day-length requirement for flowering in Arabidopsis accessions. Our findings reveal that environmental and developmental signals converge at MYC3 suppression of FT, an elementary event underlying natural variation in photoperiodic flowering responses.