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.

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.

Soybean production is severely hampered by saline-alkaline stress caused by saline-alkalization. Plants have aldehydrogenase (ALDH) family members that convert reactive aldehydes to carboxylic acids to remove active aldehyde molecules. However, little is known about the increased saline-alkali tolerance caused by the ALDH function in soybean. Here, we introduced a previously identified ALDH coding gene AhALDH3H1 from Arachis hypogaea into the soybean genome to investigate its critical role in response to saline-alkali stress. Transgenic soybean with increased aldehyde dehydrogenase activity showed significant tolerance to saline-alkali stress. It reduced malondialdehyde (MDA) content compared to its receptor, suggesting that over-expression of AhALDH3H1 accelerated soybean tolerance to saline-alkali stress by increasing aldehyde dehydrogenase activity, which is responsible for scavenging toxic MDA. To further analyze the inner mechanisms that allow transgenic plants to tolerate saline-alkali stress, we sequenced the transcriptome and metabolome of P3 (wild type, WT) and transgenic lines which were separately treated with water and a saline-alkali solution. When subjected to saline-alkali stress, the integrated analysis of the transcriptome and metabolome suggested that several genes related to cell wall structure crucial for preserving cell wall extensibility and plasticity were largely responsible for restoring homeostasis within the transgenic cells compared to WT. Metabolites, including both necessary ingredients for cell wall genesis and harmful production produced during the saline-alkali stress response, could be transported efficiently with the help of the ABC transporter, reducing the negative effects of saline-alkali stress. These findings suggest that introducing AhALDH3H1 increases transgenic soybean tolerance to saline-alkali stress may through cell wall structure maintenance and metabolites transport.

The soybean pod borer, Leguminivora glycinivorella Matsumura (Lepidoptera: Tortricidae), is an economically significant soybean pest in northeastern Asia. Serine proteases are crucial enzymes responsible for protein digestion in herbivorous lepidopterans. In this study, a gene (Spbtry1) encoding a soybean pod borer serine protease was cloned from the organism's midgut. The Spbtry1 open reading frame encoded a 269 amino acid protein with a predicted molecular mass of approximately 29 kDa. Alignment of Spbtry1 with trypsins and chymotrypsins from other insects revealed a high degree of conservation in the putative catalytic domain region. Analysis by reverse transcriptase polymerase chain reaction indicated that Spbtry1 was specifically expressed in the midgut, its transcript existed constitutively in the larval stage, and its expression was highest in the 3rd instar larval stage. RNA interference indicated that Spbtry1 expression levels decreased on diets containing Spbtry1 double-stranded RNA (ds... Resumen El barrenador de la vaina de la soja, Leguminivora glycinivorella Matsumura (Lepidoptera: Tortricidae), es una plaga económicamente significativa de la soja en el noreste de Asia. Las serinas proteasas son enzimas cruciales responsables de la digestión de proteínas en los lepidópteros herbívoros. En este estudio, se clonó un gen (Spbtry1) que codifica una serina proteasa del intestino medio del organismo del barrenador de vaina de la soja. El marco de lectura abierto Spbtry1 codificó una proteína de 269 aminoácidos con una masa molecular predicha de aproximadamente 29 kDa. La alineación de Spbtry1 con tripsinas y quimotripsinas de otros insectos reveló un alto grado de conservación en la región del dominio catalítico putativo. El análisis por la reacción en cadena de la polimerasa de la transcriptasa inversa indicó que Spbtry1 se expresó específicamente en el intestino medio, su transcripción existió constitutivamente en la fase larvaria, y su expresión fue más alta en la fase larvaria del 3er ins...