Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis -regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Crossover hotspots are strikingly distinct among rice subspecies, despite their relatively recent divergence, and associate with DNA transposons and distinct pa

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana . We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana . Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species . We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131–157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

The highly transformable maize inbred line LH244 represents an attractive model for gene discovery and genome engineering. However, the lack of a high‐quality genome assembly has limited its utility in functional genomics research. Here, we present a 2.29 Gb near‐complete assembly of the LH244 maize genome, with an overall base accuracy of 99.998%. Except for five gaps associated with super‐long thymine–adenine–guanine (TAG) repeat arrays, all the genome sequences were assembled from telomere to telomere (T2T). Comparative analysis revealed high genetic similarity between LH244 and B73, including 80.06% genome‐wide synteny and 90.92% of genes nearly identical. The LH244 genome was also compared with the complete Mo17 genome and revealed extensive intraspecific genomic variations. A total of 14 megabase‐scale structural variations (SVs) were identified, including a 3.15 Mb insertion, harboring 95 genes, within the 45S rDNA array of LH244 but not in the Mo17 genome. In addition, there were five knob arrays, with an average size of 21.76 Mb and the longest of 38.70 Mb, only existing in the LH244 genome. Despite the substantial variation in knob abundance, knob‐6S and knob‐8L were highly conserved between LH244 and Mo17, showing strong synteny and sequence identity, as well as consistent insertion patterns of genes and transposable elements (TEs). Overall, our study provides a near‐complete reference genome of an important transformable maize germplasm, which will serve as a much‐needed resource for functional genomics and genome editing of maize.