Promoter-proximal introns in Arabidopsis thaliana are enriched in dispersed signals that elevate gene expression

Cas12a-mediated gene targeting by sequential transformation strategy in Arabidopsis thaliana

Gene targeting (GT) allows precise manipulation of genome sequences, such as knock-ins and sequence substitutions, but GT in seed plants remains a challenging task. Engineered sequence-specific nucleases (SSNs) are known to facilitate GT via homology-directed repair (HDR) in organisms. Here, we demonstrate that Cas12a and a temperature-tolerant Cas12a variant (ttCas12a) can efficiently establish precise and heritable GT at two loci in Arabidopsis thaliana (Arabidopsis) through a sequential transformation strategy. As a result, ttCas12a showed higher GT efficiency than unmodified Cas12a. In addition, the efficiency of transcriptional and translational enhancers for GT via sequential transformation strategy was also investigated. These enhancers and their combinations were expected to show an increase in GT efficiency in the sequential transformation strategy, similar to previous reports of all-in-one strategies, but only a maximum twofold increase was observed. These results indicate that the frequency of double strand breaks (DSBs) at the target site is one of the most important factors determining the efficiency of genetic GT in plants. On the other hand, a higher frequency of DSBs does not always lead to higher efficiency of GT, suggesting that some additional factors are required for GT via HDR. Therefore, the increase in DSB can no longer be expected to improve GT efficiency, and a new strategy needs to be established in the future. This research opens up a wide range of applications for precise and heritable GT technology in plants.

The first intron of ARF7 is required for expression in root tips

Auxin regulates plant growth and development through the transcription factors of the AUXIN RESPONSE FACTOR (ARF) gene family. ARF7 is one of five activators that bind DNA and elicit downstream transcriptional responses. In roots, ARF7 regulates growth, gravitropism and redundantly with ARF19, lateral root organogenesis. In this study we analyzed ARF7 cis-regulation, using different non-coding sequences of the ARF7 locus to drive GFP. We show that constructs containing the first intron led to increased signal in the root tip. Although bioinformatics analyses predicted several transcription factor binding sites in the first intron, we were unable to significantly alter expression of GFP in the root by mutating these. We instead observed the intronic sequences needed to be present within the transcribed sequences to drive expression in the root meristem. These data support a mechanism by which intron-mediated enhancement regulates the tissue specific expression of ARF7 in the root meristem.

Building a pipeline to identify and engineer constitutive and repressible promoters

To support the increasingly complex circuits needed for plant synthetic biology applications, additional constitutive promoters are essential. Reusing promoter parts can lead to difficulty in cloning, increased heterogeneity between transformants, transgene silencing and trait instability. We have developed a pipeline to identify genes that have stable expression across a wide range of Arabidopsis tissues at different developmental stages and have identified a number of promoters that are well expressed in both transient (Nicotiana benthamiana) and stable (Arabidopsis) transformation assays. We have also introduced two genome-orthogonal gRNA target sites in a subset of the screened promoters, converting them into NOR logic gates. The work here establishes a pipeline to screen for additional constitutive promoters and can form the basis of constructing more complex information processing circuits in the future.

The intragenic cis-elements mediate temperature response of RrKSN

Gene regulation via intragenic sequences is becoming more recognized in many eukaryotes. However, the intragenic sequences mediated gene expressions in response to environmental stimuli have been largely uncharacterized. Here, we showed that the first intron of RrKSN from the Rosa rugosa cultivar 'Purple branch' had a positive effect on RrKSN expression, and the effect depends on its position and orientation. Further analyses revealed that the four adjacent cis-elements (T)CGATT/AATCG(A) within the first intron were critical for the positive regulation, and the RrKSN promotion was significantly suppressed with mutations of these elements. These cis-elements were further evidenced as binding sites for RrARR1, the homologous of Arabidopsis type-B ARABIDOPSIS RESPONSE REGULATOR 1 (ARR1) transcription factor. The first intron-mediated RrKSN expression was enhanced with over-expressing of RrARR1, but abolished with RrARR1 silencing in rose seedlings. Moreover, the expression difference of RrKSN between 16°C and 28°C was eliminated along with RrARR1-silencing. Taken together, these results suggested both RrARR1 and its binding elements are required for the first intron-mediated RrKSN expression in response to varying temperatures. Therefore, our results reveal a unique intragenic regulation mechanism of gene expression by which plants perceive the signal of ambient temperature in rose.

Comprehensive analyses of microtubule-associated protein MAP65 family genes in Cucurbitaceae and CsaMAP65s expression profiles in cucumber

MAP65 is a microtubule-binding protein family in plants and plays crucial roles in regulating cell growth and development, intercellular communication, and plant responses to various environmental stresses. However, MAP65s in Cucurbitaceae are still less understood. In this study, a total of 40 MAP65s were identified from six Cucurbitaceae species (Cucumis sativus L., Citrullus lanatus, Cucumis melo L., Cucurbita moschata, Lagenaria siceraria, and Benincasa hispida) and classified into five groups by phylogenetic analysis according to gene structures and conserved domains. A conserved domain (MAP65_ASE1) was found in all MAP65 proteins. In cucumber, we isolated six CsaMAP65s with different expression patterns in tissues including root, stem, leaf, female flower, male flower, and fruit. Subcellular localizations of CsaMAP65s verified that all CsaMAP65s were localized in microtubule and microfilament. Analyses of the promoter regions of CsaMAP65s have screened different cis-acting regulatory elements involved in growth and development and responses to hormone and stresses. In addition, CsaMAP65-5 in leaves was significantly upregulated by salt stress, and this promotion effect was higher in cucumber cultivars with salt tolerant than that without salt tolerant. CsaMAP65-1 in leaves was significantly upregulated by cold stress, and this promotion was higher in cold-tolerant cultivar than intolerant cultivar. With the genome-wide characterization and phylogenetic analysis of Cucurbitaceae MAP65s, and the expression profile of CsaMAP65s in cucumber, this study laid a foundation for further study on MAP65 functions in developmental processes and responses to abiotic stress in Cucurbitaceae species.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

The position of the longest intron is related to biological functions in some human genes

The evidence that introns can influence different levels of transfer of genetic information between DNA and the final product is increasing. Longer first introns were found to be a general property of eukaryotic gene structure and shown to contain a higher fraction of conserved sequence and different functional elements. Our work brings more precise information about the position of the longest introns in human protein-coding genes and possible connection with biological function and gene expression. According to our results, the position of the longest intron can be localized to the first third of introns in 64%, the second third in 19%, and the third in 17%, with notable peaks at the middle and last introns of approximately 5% and 6%, respectively. The median lengths of the longest introns decrease with increasing distance from the start of the gene from approximately 15,000 to 5,000 bp. We have shown that the position of the longest intron is in some cases linked to the biological function of the given gene. For example, DNA repair genes have the longest intron more often in the second or third. In the distribution of gene expression according to the position of the longest intron, tissue-specific profiles can be traced with the highest expression usually at the absolute positions of intron 1 and 2. In this work, we present arguments supporting the hypothesis that the position of the longest intron in a gene is another biological factor modulating the transmission of genetic information. The position of the longest intron is related to biological functions in some human genes.

Splicing-mediated activation of SHAGGY-like kinases underpinning carbon partitioning in Arabidopsis seeds

Glycogen synthase kinase 3 (GSK3) family members serve as signaling hubs for plant development and stress responses, yet the underlying mechanism of their transcriptional regulation remains a long-standing mystery. Here we show that the transcription of SHAGGY-like kinase 11/12 (SK11/12), two members of the GSK3 gene family, is promoted by the splicing factor SmD1b, which is essential for distributing carbon sources into storage and protective components in Arabidopsis seeds. The chromatin recruitment of SmD1b at the SK11/12 loci promotes their transcription associated with co-transcriptional splicing of the first introns in the 5'-untranslated region of SK11/12. The loss of SmD1b function generates transcripts with unspliced introns that create disruptive R-loops to hamper the transcriptional elongation of SK11/12, in addition to compromising the recruitment of RNA polymerase II to the SK11/12 genomic regions. These effects imposed by SmD1b determine the transcription of SK11/12 to confer a key switch of carbon flow among metabolic pathways in zygotic and maternal tissues in seeds.

Orphan gene in Littorina: An unexpected role of symbionts in the host evolution

Mechanisms of reproductive isolation between closely related sympatric species are of high evolutionary significance as they may function as initial drivers of speciation and protect species integrity afterwards. Proteins involved in the establishment of reproductive barriers often evolve fast and may be key players in cessation of gene flow between the incipient species. The five Atlantic Littorina (Neritrema) species represent a notable example of recent radiation. The geographic ranges of these young species largely overlap and the mechanisms of reproductive isolation are poorly understood. In this study, we performed a detailed analysis of the reproductive protein LOSP, previously identified in Littorina. We showed that this protein is evolutionary young and taxonomically restricted to the genus Littorina. It has high sequence variation both within and between Littorina species, which is compatible with its presumable role in the reproductive isolation. The strongest differences in the LOSP structure were detected between Littorina subgenera with distinctive repetitive motifs present exclusively in the Neritrema species, but not in L. littorea. Moreover, the sequence of these repetitive structural elements demonstrates a high homology with genetic elements of bacteria, identified as components of Littorina associated microbiomes. We suggest that these elements were acquired from a symbiotic bacterial donor via horizontal genetic transfer (HGT), which is indirectly confirmed by the presence of multiple transposable elements in the LOSP flanking and intronic regions. Furthermore, we hypothesize that this HGT-driven evolutionary innovation promoted LOSP function in reproductive isolation, which might be one of the factors determining the intensive cladogenesis in the Littorina (Neritrema) lineage in contrast to the anagenesis in the L. littorea clade.

Interplay between gene expression and gene architecture as a consequence of gene and genome duplications: evidence from metabolic genes of Arabidopsis thaliana

Gene and genome duplications have been widespread during the evolution of flowering plant which resulted in the increment of biological complexity as well as creation of plasticity of a genome helping the species to adapt to changing environments. Duplicated genes with higher evolutionary rates can act as a mechanism of generating novel functions in secondary metabolism. In this study, we explored duplication as a potential factor governing the expression heterogeneity and gene architecture of Primary Metabolic Genes (PMGs) and Secondary Metabolic Genes (SMGs) of Arabidopsis thaliana. It is remarkable that different types of duplication processes controlled gene expression and tissue specificity differently in PMGs and SMGs. A complex relationship exists between gene architecture and expression patterns of primary and secondary metabolic genes. Our study reflects, expression heterogeneity and gene structure variation of primary and secondary metabolism in Arabidopsis thaliana are partly results of duplication events of different origins. Our study suggests that duplication has differential effect on PMGs and SMGs regarding expression pattern by controlling gene structure, epigenetic modifications, multifunctionality and subcellular compartmentalization. This study provides an insight into the evolution of metabolism in plants in the light of gene and genome scale duplication.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01188-2.

Genome-Wide Identification and Expression Profiling of Potassium Transport-Related Genes in Vigna radiata under Abiotic Stresses

Potassium (K+) is one of the most important cations that plays a significant role in plants and constitutes up to 10% of plants' dry weight. Plants exhibit complex systems of transporters and channels for the distribution of K+ from soil to numerous parts of plants. In this study, we have identified 39 genes encoding putative K+ transport-related genes in Vigna radiata. Chromosomal mapping of these genes indicated an uneven distribution across eight out of 11 chromosomes. Comparative phylogenetic analysis of different plant species, i.e., V. radiata, Glycine max, Cicer arietinum, Oryza sativa, and Arabidopsis thaliana, showed their strong conservation in different plant species. Evolutionary analysis of these genes suggests that gene duplication is a major route of expansion for this family in V. radiata. Comprehensive promoter analysis identified several abiotic stresses related to cis-elements in the promoter regions of these genes, suggesting their role in abiotic stress tolerance. Our additional analyses indicated that abiotic stresses adversely affected the chlorophyll concentration, carotenoids, catalase, total soluble protein concentration, and the activities of superoxide and peroxidase in V. radiata. It also disturbs the ionic balance by decreasing the uptake of K+ content and increasing the uptake of Na+. Expression analysis from high-throughput sequencing data and quantitative real-time PCR experiments revealed that several K+ transport genes were expressed in different tissues (seed, flower, and pod) and in abiotic stress-responsive manners. A highly significant variation of expression was observed for VrHKT (1.1 and 1.2), VrKAT (1 and 2) VrAKT1.1, VrAKT2, VrSKOR, VrKEA5, VrTPK3, and VrKUP/HAK/KT (4, 5, and 8.1) in response to drought, heat or salinity stress. It reflected their potential roles in plant growth, development, or stress adaptations. The present study gives an in-depth understanding of K+ transport system genes in V. radiata and will serve as a basis for a functional analysis of these genes.

Genome-Scale Computational Identification and Characterization of UTR Introns in Atalantia buxifolia

Accumulated evidence has shown that CDS introns (CIs) play important roles in regulating gene expression. However, research on UTR introns (UIs) is limited. In this study, UIs (including 5′UTR and 3′UTR introns (5UIs and 3UIs)) were identified from the Atalantia buxifolia genome. The length and nucleotide distribution characteristics of both 5UIs and 3UIs and the distributions of cis-acting elements and transcription factor binding sites (TFBSs) in 5UIs were investigated. Moreover, PageMan enrichment analysis was applied to show the possible roles of transcripts containing UIs (UI-Ts). In total, 1077 5UIs and 866 3UIs were identified from 897 5UI-Ts and 670 3UI-Ts, respectively. Among them, 765 (85.28%) 5UI-Ts and 527 (78.66%) 3UI-Ts contained only one UI, and 94 (6.38%) UI-Ts contained both 5UI and 3UI. The UI density was lower than that of CDS introns, but their mean and median intron sizes were ~2 times those of the CDS introns. The A. buxifolia 5UIs were rich in gene-expression-enhancement-related elements and contained many TFBSs for BBR-BPC, MIKC_MADS, AP2 and Dof TFs, indicating that 5UIs play a role in regulating or enhancing the expression of downstream genes. Enrichment analysis revealed that UI-Ts involved in ‘not assigned’ and ‘RNA’ pathways were significantly enriched. Noteworthily, 119 (85.61%) of the 3UI-Ts were genes encoding pentatricopeptide (PPR) repeat-containing proteins. These results will be helpful for the future study of the regulatory roles of UIs in A. buxifolia.

Ectopic expression of VRT-A2 underlies the origin of Triticum polonicum and Triticum petropavlovskyi with long outer glumes and grains

Polish wheat (Triticum polonicum) is a unique tetraploid wheat species characterized by an elongated outer glume. The genetic control of the long-glume trait by a single semi-dominant locus, P1 (from Polish wheat), was established more than 100 years ago, but the underlying causal gene and molecular nature remain elusive. Here, we report the isolation of VRT-A2, encoding an SVP-clade MADS-box transcription factor, as the P1 candidate gene. Genetic evidence suggests that in T. polonicum, a naturally occurring sequence rearrangement in the intron-1 region of VRT-A2 leads to ectopic expression of VRT-A2 in floral organs where the long-glume phenotype appears. Interestingly, we found that the intron-1 region is a key ON/OFF molecular switch for VRT-A2 expression, not only because it recruits transcriptional repressors, but also because it confers intron-mediated transcriptional enhancement. Genotypic analyses using wheat accessions indicated that the P1 locus is likely derived from a single natural mutation in tetraploid wheat, which was subsequently inherited by hexaploid T. petropavlovskyi. Taken together, our findings highlight the promoter-proximal intron variation as a molecular basis for phenotypic differentiation, and thus species formation in Triticum plants.

Cooperative roles of introns 1 and 2 of tobacco resistance gene N in enhanced N transcript expression and antiviral defense responses

The tobacco virus resistance gene N contains four introns. Transient expression of transcripts from an N transgene containing these introns and driven by the native promoter in the presence of the elicitor of tobacco mosaic virus resulted in its increased expression. The requirement of the native promoter, the elicitor, or the individual introns for enhanced expression of N has not been fully studied. Here, we determined that 35S promoter-driven N transcript expression could be enhanced in the presence of the four introns regardless of the co-expression of the virus elicitor in tobacco. Function analyses using a series of N transgenes with different combination of introns revealed that the presence of intron 1 more so than intron 2 allowed higher accumulation of premature and mature N transcripts; however, both introns were important for not only enhanced gene expression but also for induction of cell death in tobacco and induced local resistance to spread of virus in Nicotiana benthamiana. Our findings indicate that introns 1 and 2 cooperatively contribute to N expression and virus resistance.

Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana

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.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

Identification of cis-regulatory motifs in first introns and the prediction of intron-mediated enhancement of gene expression in Arabidopsis thaliana

BACKGROUND

Intron mediated enhancement (IME) is the potential of introns to enhance the expression of its respective gene. This essential function of introns has been observed in a wide range of species, including fungi, plants, and animals. However, the mechanisms underlying the enhancement are as of yet poorly understood. The goal of this study was to identify potential IME-related sequence motifs and genomic features in first introns of genes in Arabidopsis thaliana.

RESULTS

Based on the rationale that functional sequence motifs are evolutionarily conserved, we exploited the deep sequencing information available for Arabidopsis thaliana, covering more than one thousand Arabidopsis accessions, and identified 81 candidate hexamer motifs with increased conservation across all accessions that also exhibit positional occurrence preferences. Of those, 71 were found associated with increased correlation of gene expression of genes harboring them, suggesting a cis-regulatory role. Filtering further for effect on gene expression correlation yielded a set of 16 hexamer motifs, corresponding to five consensus motifs. While all five motifs represent new motif definitions, two are similar to the two previously reported IME-motifs, whereas three are altogether novel. Both consensus and hexamer motifs were found associated with higher expression of alleles harboring them as compared to alleles containing mutated motif variants as found in naturally occurring Arabidopsis accessions. To identify additional IME-related genomic features, Random Forest models were trained for the classification of gene expression level based on an array of sequence-related features. The results indicate that introns contain information with regard to gene expression level and suggest sequence-compositional features as most informative, while position-related features, thought to be of central importance before, were found with lower than expected relevance.

CONCLUSIONS

Exploiting deep sequencing and broad gene expression information and on a genome-wide scale, this study confirmed the regulatory role on first-introns, characterized their intra-species conservation, and identified a set of novel sequence motifs located in first introns of genes in the genome of the plant Arabidopsis thaliana that may play a role in inducing high and correlated gene expression of the genes harboring them.

The Arabidopsis condensin CAP-D subunits arrange interphase chromatin

Condensins are best known for their role in shaping chromosomes. Other functions such as organizing interphase chromatin and transcriptional control have been reported in yeasts and animals, but little is known about their function in plants. To elucidate the specific composition of condensin complexes and the expression of CAP-D2 (condensin I) and CAP-D3 (condensin II), we performed biochemical analyses in Arabidopsis. The role of CAP-D3 in interphase chromatin organization and function was evaluated using cytogenetic and transcriptome analysis in cap-d3 T-DNA insertion mutants. CAP-D2 and CAP-D3 are highly expressed in mitotically active tissues. In silico and pull-down experiments indicate that both CAP-D proteins interact with the other condensin I and II subunits. In cap-d3 mutants, an association of heterochromatic sequences occurs, but the nuclear size and the general histone and DNA methylation patterns remain unchanged. Also, CAP-D3 influences the expression of genes affecting the response to water, chemicals, and stress. The expression and composition of the condensin complexes in Arabidopsis are similar to those in other higher eukaryotes. We propose a model for the CAP-D3 function during interphase in which CAP-D3 localizes in euchromatin loops to stiffen them and consequently separates centromeric regions and 45S rDNA repeats.

Functional Characterization of GhACX3 Gene Reveals Its Significant Role in Enhancing Drought and Salt Stress Tolerance in Cotton

The acyl-coenzyme A oxidase 3 (ACX3) gene involved in the β-oxidation pathway plays a critical role in plant growth and development as well as stress response. Earlier on, studies focused primarily on the role of β-oxidation limited to fatty acid breakdown. However, ACX3 peroxisomal β-oxidation pathways result in a downstream cascade of events that act as a transduction of biochemical and physiological responses to stress. A role that is yet to be studied extensively. In this study, we identified 20, 18, 22, 23, 20, 11, and 9 proteins in Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum, G. darwinii, G. arboretum, and G. raimondii genomes, respectively. The tetraploid cotton genome had protein ranging between 18 and 22, while diploids had between 9 and 11. After analyzing the gene family evolution or selection pressure, we found that this gene family undergoes purely segmental duplication both in diploids and tetraploids. W-Box (WRKY-binding site), ABRE, CAAT-Box, TATA-box, MYB, MBS, LTR, TGACG, and CGTCA-motif are abiotic stress cis-regulatory elements identified in this gene family. All these are the binding sites for abiotic stress transcription factors, indicating that this gene is essential. Genes found in G. hirsutum showed a clear response to drought and salinity stress, with higher expression under drought and salt stress, particularly in the leaf and root, according to expression analysis. We selected Gh_DO1GO186, one of the highly expressed genes, for functional characterization. We functionally characterized the GhACX3 gene through overexpression and virus-induced gene silencing (VIGS). Overexpression of this gene enhanced tolerance under stress, which was exhibited by the germination assay. The overexpressed seed growth rate was faster relative to control under drought and salt stress conditions. The survival rate was also higher in overexpressed plants relative to control plants under stress. In contrast, the silencing of the GhACX3 gene in cotton plants resulted in plants showing the stress susceptibility phenotype and reduced root length compared to control. Biochemical analysis also demonstrated that GhACX3-silenced plants experienced oxidative stress while the overexpressed plants did not. This study has revealed the importance of the ACX3 family during stress tolerance and can breed stress-resilient cultivar.

Analysis of CACTA transposase genes unveils the mechanism of intron loss and distinct small RNA silencing pathways underlying divergent evolution of Brassica genomes

In comparison with retrotransposons, DNA transposons make up a smaller proportion of most plant genomes. However, these elements are often proximal to genes to affect gene expression depending on the activity of the transposons, which is largely reflected by the activity of the transposase genes. Here, we show that three AT-rich introns were retained in the TNP2-like transposase genes of the Bot1 (Brassica oleracea transposon 1) CACTA transposable elements in Brassica oleracea, but were lost in the majority of the Bot1 elements in Brassica rapa. A recent burst of transposition of Bot1 was observed in B. oleracea, but not in B. rapa. This burst of transposition is likely related to the activity of the TNP2-like transposase genes as the expression values of the transposase genes were higher in B. oleracea than in B. rapa. In addition, distinct populations of small RNAs (21, 22 and 24 nt) were detected from the Bot1 elements in B. oleracea, but the vast majority of the small RNAs from the Bot1 elements in B. rapa are 24 nt in length. We hypothesize that the different activity of the TNP2-like transposase genes is likely associated with the three introns, and intron loss is likely reverse transcriptase mediated. Furthermore, we propose that the Bot1 family is currently undergoing silencing in B. oleracea, but has already been silenced in B. rapa. Taken together, our data provide new insights into the differentiation of transposons and their role in the asymmetric evolution of these two closely related Brassica species.

Multidimensional analysis of actin depolymerising factor family in pigeon pea under different environmental stress revealed specific response genes in each subgroup

Actin depolymerising factor (ADF) is an actin binding protein that is ubiquitous in animal and plant cells. It plays an important role in plant growth and development, as well as resistance to biotic and abiotic stress. The research of plant ADF family has been restricted to Arabidopsis thaliana (L.) Heynh. and some herb crops, but no woody cash crops have been reported to date. All members of the Cajanus cajan (L.) Millsp. ADF (CcADF) family were identified from the pigeon pea genome, and distributed among the four subfamilies by phylogenetic analysis. CcADFs were relatively conservative in gene structure evolution, protein structure and functional expression, and different CcADFs showed specific expression patterns under different treatments. The expression characteristics of several key CcADFs were revealed by analysing the stress response pattern of CcADFs and the time series RNA-seq of aluminium stress. Among them, CcADF9 in the first subgroup specifically responded to aluminium stress in the roots; CcADF3 in the second subgroup intensively responded to fungal infection in the leaves; and CcADF2 in the fourth subgroup positively responded to various stress treatments in different tissues. This study extended the relationship between plant ADF family and aluminium tolerance, as well as adding to the understanding of CcADF family in woody crops.

Introns mediate post-transcriptional enhancement of nuclear gene expression in the green microalga Chlamydomonas reinhardtii

Efficient nuclear transgene expression in the green microalga Chlamydomonas reinhardtii is generally hindered by low transcription rates. Introns can increase transcript abundance by a process called Intron-Mediated Enhancement (IME) in this alga and has been broadly observed in other eukaryotes. However, the mechanisms of IME in microalgae are poorly understood. Here, we identified 33 native introns from highly expressed genes in C. reinhardtii selected from transcriptome studies as well as 13 non-native introns. We investigated their IME capacities and probed the mechanism of action by modification of splice sites, internal sequence motifs, and position within transgenes. Several introns were found to elicit strong IME and found to be broadly applicable in different expression constructs. We determined that IME in C. reinhardtii exclusively occurs from introns within transcribed ORFs regardless of the promoter and is not induced by traditional enhancers of transcription. Our results elucidate some mechanistic details of IME in C. reinhardtii, which are similar to those observed in higher plants yet underly distinctly different induction processes. Our findings narrow the focus of targets responsible for algal IME and provides evidence that introns are underestimated regulators of C. reinhardtii nuclear gene expression.

Although many genetic tools and basic transformation strategies exist for the model microalga Chlamydomonas reinhardtii, high-level genetic engineering with this organism is hindered by its inherent recalcitrance to foreign gene expression and limited knowledge of responsible expression regulators. In this work, we characterized the dynamics of 33 endogenous and 13 non-native introns and their effect on gene expression as artificial insertions into codon optimized transgenes. We found that introns from different origins have the capacity to increase gene expression rates. Intron-mediated enhancement was observed exclusively when these elements were placed in transcripts but not outside of transcribed mRNA regions. Insertion of different endogenous introns into coding sequences was found to positively affect expression rates through a synergy of additive transcription enhancement and exon length reduction, similar to those natively found in the C. reinhardtii genome. Our results indicate that intensive mRNA processing plays an underestimated role in the regulation of native gene expression in C. reinhardtii. In addition to internal sequence motifs, the location of artificially introduced introns greatly affected transgene expression levels. This work is highly valuable to the greater microalgal and synthetic biology research communities and contributes to broadening our understanding of eukaryotic intron-mediated enhancement.

First Come, First Served: Sui Generis Features of the First Intron

Most of the transcribed genes in eukaryotic cells are interrupted by intervening sequences called introns that are co-transcriptionally removed from nascent messenger RNA through the process of splicing. In Arabidopsis, 79% of genes contain introns and more than 60% of intron-containing genes undergo alternative splicing (AS), which ostensibly is considered to increase protein diversity as one of the intrinsic mechanisms for fitness to the varying environment or the internal developmental program. In addition, recent findings have prevailed in terms of overlooked intron functions. Here, we review recent progress in the underlying mechanisms of intron function, in particular by focusing on unique features of the first intron that is located in close proximity to the transcription start site. The distinct deposition of epigenetic marks and nucleosome density on the first intronic DNA sequence, the impact of the first intron on determining the transcription start site and elongation of its own expression (called intron-mediated enhancement, IME), translation control in 5′-UTR, and the new mechanism of the trans-acting function of the first intron in regulating gene expression at the post-transcriptional level are summarized.

Conserved exitrons of FLAGELLIN-SENSING 2 (FLS2) across dicot plants and their functions

The alternative splicing of pattern recognition receptor genes regulates immune signalling in mammals, but in plants its role is still unknown. Here, we detected alternatively spliced introns (exitrons) in the first annotated exons of FLAGELLIN-SENSING 2 (FLS2) genes in all the examined dicot plants across nine families. The 5' splice site (SS) regions were conserved and with rare synonymous substitutions. Point mutations and gene swaps indicated that the position and efficiency of exitron splicing primarily depended on the nucleotide sequences of FLS2 genes. Single-nucleotide mutations in the invariable codon carrying 5' SS dramatically altered the accumulation of poplar and tomato FLS2 transcripts, indicating the 5'-proximal exitrons of FLS2 function as stimulatory introns on gene expression. The 3' SSs of exitrons are diverse and can be changed by 1-2 nucleotide mutations in Salicaceae FLS2. The alternative transcripts (ATs) of poplar and tobacco FLS2, which encode small secreted proteins, were specifically induced by flg22, and one such AT from tobacco FLS2 suppressed flg22-induced response. Our results indicated that the exitrons of FLS2 genes regulate the accumulation of transcripts by an intron mediated enhancement (IME) mechanism and some ATs have the potential to encode suppressors for FLS2 pathway.

The transcription factor HAND2 up-regulates transcription of the IL15 gene in human endometrial stromal cells

Cyclic changes of the human endometrium, such as proliferation, secretion, and decidualization, occur during regular menstrual cycles. Heart- and neural crest derivatives-expressed transcript 2 (HAND2) is a key transcription factor in progestin-induced decidualization of human endometrial stromal cells (ESCs). It has been suggested that HAND2 regulates interleukin 15 (IL15), a key immune factor required for the activation and survival of uterine natural killer (uNK) cells. Activated uNK cells can promote spiral artery remodeling and secrete cytokines to induce immunotolerance. To date, no studies have evaluated the transcription factors that regulate IL15 expression in human ESCs. In the present study, we examined whether HAND2 controls IL15 transcriptional regulation in human ESCs. Quantitative RT-PCR and histological analyses revealed that HAND2 and IL15 levels increase considerably in the secretory phase of human endometrium tissues. Results from ChIP-quantitative PCR suggested that HAND2 binds to a putative HAND2 motif, which we identified in the upstream region of the human IL15 gene through in silico analysis. Using a luciferase reporter assay, we found that the upstream region of the human IL15 gene up-regulates reporter gene activities in response to estradiol and a progestin representative (medroxyprogesterone) in ESCs. The upstream region of the human IL15 gene also exhibited increasing responsiveness to transfection with a HAND2 expression vector. Of note, deletion and substitution variants of the putative HAND2 motif in the upstream region of IL15 did not respond to HAND2 transfection. These findings confirm that HAND2 directly up-regulates human IL15 transcription in ESCs.

Native elongation transcript sequencing reveals temperature dependent dynamics of nascent RNAPII transcription in Arabidopsis

Temperature profoundly affects the kinetics of biochemical reactions, yet how large molecular complexes such as the transcription machinery accommodate changing temperatures to maintain cellular function is poorly understood. Here, we developed plant native elongating transcripts sequencing (plaNET-seq) to profile genome-wide nascent RNA polymerase II (RNAPII) transcription during the cold-response of Arabidopsis thaliana with single-nucleotide resolution. Combined with temporal resolution, these data revealed transient genome-wide reprogramming of nascent RNAPII transcription during cold, including characteristics of RNAPII elongation and thousands of non-coding transcripts connected to gene expression. Our results suggest a role for promoter-proximal RNAPII stalling in predisposing genes for transcriptional activation during plant-environment interactions. At gene 3'-ends, cold initially facilitated transcriptional termination by limiting the distance of read-through transcription. Within gene bodies, cold reduced the kinetics of co-transcriptional splicing leading to increased intragenic stalling. Our data resolved multiple distinct mechanisms by which temperature transiently altered the dynamics of nascent RNAPII transcription and associated RNA processing, illustrating potential biotechnological solutions and future focus areas to promote food security in the context of a changing climate.

Putative cis-Regulatory Elements Predict Iron Deficiency Responses in Arabidopsis Roots

Plant iron deficiency (-Fe) activates a complex regulatory network that coordinates root Fe uptake and distribution to sink tissues. In Arabidopsis (Arabidopsis thaliana), FER-LIKE FE DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT), a basic helix-loop-helix (bHLH) transcription factor (TF), regulates root Fe acquisition genes. Many other -Fe-induced genes are FIT independent, and instead regulated by other bHLH TFs and by yet unknown TFs. The cis-regulatory code, that is, the cis-regulatory elements (CREs) and their combinations that regulate plant -Fe-responses, remains largely elusive. Using Arabidopsis root transcriptome data and coexpression clustering, we identified over 100 putative CREs (pCREs) that predicted -Fe-induced gene expression in computational models. To assess pCRE properties and possible functions, we used large-scale in vitro TF binding data, positional bias, and evolutionary conservation. As one example, our approach uncovered pCREs resembling IDE1 (iron deficiency-responsive element 1), a known grass -Fe response CRE. Arabidopsis IDE1-likes were associated with FIT-dependent gene expression, more specifically with biosynthesis of Fe-chelating compounds. Thus, IDE1 seems to be conserved in grass and nongrass species. Our pCREs matched among others in vitro binding sites of B3, NAC, bZIP, and TCP TFs, which might be regulators of -Fe responses. Altogether, our findings provide a comprehensive source of cis-regulatory information for -Fe-responsive genes that advance our mechanistic understanding and inform future efforts in engineering plants with more efficient Fe uptake or transport systems.

Shoot sodium exclusion in salt stressed barley (Hordeum vulgare L.) is determined by allele specific increased expression of HKT1;5

High affinity potassium transporters (HKT) are recognized as important genes for crop salt tolerance improvement. In this study, we investigated HvHKT1;5 as a candidate gene for a previously discovered quantitative trait locus that controls shoot Na⁺ and Na⁺/K⁺ ratio in salt-stressed barley lines on a hydroponic system. Two major haplotype groups could be distinguished for this gene in a barley collection of 95 genotypes based on the presence of three intronic insertions; a designated haplotype group A (HGA, same as reference sequence) and haplotype group B (HGB, with insertions). HGB was associated with a much stronger root expression of HKT1;5 compared to HGA, and consequently higher K⁺ and lower Na⁺ and Cl^- concentrations and a lower Na⁺/K⁺ ratio in the shoots three weeks after exposure to 200 mM NaCl. Our experimental results suggest that allelic variation in the promoter region of the HGB gene is linked to the three insertions may be responsible for the observed increase in expression of HvHKT1;5 alleles after one week of salt stress induction. This study shows that in barley - similar to wheat and rice - HKT1;5 is an important contributor to natural variation in shoot Na⁺ exclusion.

An intron-derived motif strongly increases gene expression from transcribed sequences through a splicing independent mechanism in Arabidopsis thaliana

Certain introns significantly increase mRNA accumulation by a poorly understood mechanism. These introns have no effect when located upstream, or more than ~1 Kb downstream, of the start of transcription. We tested the ability of a formerly non-stimulating intron containing 11 copies of the sequence TTNGATYTG, which is over-represented in promoter-proximal introns in Arabidopsis thaliana, to affect expression from various positions. The activity profile of this intron at different locations was similar to that of a natural intron from the UBQ10 gene, suggesting that the motif increases mRNA accumulation by the same mechanism. A series of introns with different numbers of this motif revealed that the effect on expression is linearly dependent on motif copy number up to at least 20, with each copy adding another 1.5-fold increase in mRNA accumulation. Furthermore, 6 copies of the motif stimulated mRNA accumulation to a similar degree from within an intron or when introduced into the 5'-UTR and coding sequences of an intronless construct, demonstrating that splicing is not required for this sequence to boost expression. The ability of this motif to substantially elevate expression from several hundred nucleotides downstream of the transcription start site reveals a novel type of eukaryotic gene regulation.

Arabidopsis mlo3 mutant plants exhibit spontaneous callose deposition and signs of early leaf senescence

Arabidopsis thaliana mlo3 mutant plants are not affected in pathogen infection phenotypes but-reminiscent of mlo2 mutant plants-exhibit spontaneous callose deposition and signs of early leaf senescence. The family of Mildew resistance Locus O (MLO) proteins is best known for its profound effect on the outcome of powdery mildew infections: when the appropriate MLO protein is absent, the plant is fully resistant to otherwise virulent powdery mildew fungi. However, most members of the MLO protein family remain functionally unexplored. Here, we investigate Arabidopsis thaliana MLO3, the closest relative of AtMLO2, AtMLO6 and AtMLO12, which are the Arabidopsis MLO genes implicated in the powdery mildew interaction. The co-expression network of AtMLO3 suggests association of the gene with plant defense-related processes such as salicylic acid homeostasis. Our extensive analysis shows that mlo3 mutants are unaffected regarding their infection phenotype upon challenge with the powdery mildew fungi Golovinomyces orontii and Erysiphe pisi, the oomycete Hyaloperonospora arabidopsidis, and the bacterial pathogen Pseudomonas syringae (the latter both in terms of basal and systemic acquired resistance), indicating that the protein does not play a major role in the response to any of these pathogens. However, mlo3 genotypes display spontaneous callose deposition as well as signs of early senescence in 6- or 7-week-old rosette leaves in the absence of any pathogen challenge, a phenotype that is reminiscent of mlo2 mutant plants. We hypothesize that de-regulated callose deposition in mlo3 genotypes might be the result of a subtle transient aberration of salicylic acid-jasmonic acid homeostasis during development.

Genome-wide identification, expression profiling, and network analysis of AT-hook gene family in maize

AT-hook motif nuclear localized (AHL) genes have diverse but poorly understood biological functions. We identified and analyzed 37 AHL genes in maize. We also discovered four and one additional AHLs in rice and sorghum, respectively, besides those reported earlier. The maize AHLs were classified into two clades (A and B) and three distinct types (I, II, and III) as also reported in Arabidopsis. Phylogenetic and ortholog analyses showed that, while the evolutionary classification was conserved in plants, expansion of the AHL gene family in maize was accompanied with new biological functions. Gene structure analysis showed that, while all but one Type-I AHLs lacked an intron, origin of Type-II and Type-III AHLs was associated with the gain of introns suggesting evolutionarily distinct temporal and spatial expression patterns and, likely, neofunctionalization. Gene duplication analysis revealed that AHLs in maize expanded via dispersive duplication further supporting their functional diversity. To discern these functions, we analyzed 71 transcriptomes from diverse tissues and developmental stages of maize and classified AHLs into eight groups with distinct temporal/spatial expression profiles. Coexpression analysis implicated 5 AHLs and 33 novel genes in networks specific to endosperm, seed, root, leaf, and reproductive tissues indicating their role in the development of these organs. Major processes coregulated by AHLs include pollen development, drought response, senescence, and wound response. We also identified interactions of AHL proteins in coregulating important processes including stress response. These novel insights into the role of AHLs in plant development provide a platform for functional analyses in maize and related grasses.

Analysis of the Functional Relevance of Epigenetic Chromatin Marks in the First Intron Associated with Specific Gene Expression Patterns

We previously showed that the first intron of genes exhibits several interesting characteristics not seen in other introns: 1) it is the longest intron on average in almost all eukaryotes, 2) it presents the highest number of conserved sites, and 3) it exhibits the highest density of regulatory chromatin marks. Here, we expand on our previous study by integrating various multiomics data, leading to further evidence supporting the functionality of sites in the first intron. We first show that trait-associated single-nucleotide polymorphisms (TASs) are significantly enriched in the first intron. We also show that within the first intron, the density of epigenetic chromatin signals is higher near TASs than in distant regions. Furthermore, the distribution of several chromatin regulatory marks is investigated in relation to gene expression specificity (i.e., housekeeping vs. tissue-specific expression), essentiality (essential genes vs. nonessential genes), and levels of gene expression; housekeeping genes or essential genes contain greater proportions of active chromatin marks than tissue-specific genes or nonessential genes, and highly expressed genes exhibit a greater density of chromatin regulatory marks than genes with low expression. Moreover, we observe that genes carrying multiple first-intron TASs interact with each other within a large protein-protein interaction network, ultimately connecting to the UBC protein, a well-established protein involved in ubiquitination. We believe that our results shed light on the functionality of first introns as a genomic entity involved in gene expression regulation.

Combination of the endogenous promoter-intron significantly improves salt and drought tolerance conferred by TdSHN1 transcription factor in transgenic tobacco

Recent years have witnessed a renewed interest in introns as a tool to increase gene expression. We previously isolated TdSHN1 gene encoding a transcription factor in durum wheat. Here we show that TdSHN1 intron contains many CT-stretches and the motif CGATT known to be important for IME. When subjected to bioinformatics analysis using IMEter software, TdSHN1 intron obtained a score of 17.04 which indicates that it can moderately enhance gene expression. TdSHN1 gene including its intron was placed under the control of TdSHN1 endogenous salt and drought-inducible promoter or the constitutive 35S promoter and transferred into tobacco. Transgenic lines were obtained and designated gD (with 35S promoter) and PI (with native promoter). A third construct was also used in which intron-less cDNA was driven by the 35S promoter (cD lines). Results showed that, gD lines exhibited lower stomatal density than cD lines. When subjected to drought and salt stresses, gD lines outperformed intron-less cD lines and WT. Indeed, gD lines exhibited longer roots, higher biomass production, retained more chlorophyll, produced less ROS and MDA and had higher antioxidant activity. qRT-PCR analysis revealed that gD lines had higher TdSHN1 expression levels than cD lines. In addition, expression of ROS-scavengering, stress-related and wax biosynthesis tobacco genes was higher in gD lines compared to cD lines and WT. Interestingly, under stress conditions, PI transgenic lines showed higher TdSHN1 expression levels and outperformed gD lines. These results suggest that TdSHN1 intron enhances gene expression when used alone or in combination with TdSHN1 endogenous promoter.

Introns as Gene Regulators: A Brick on the Accelerator

A picture is beginning to emerge from a variety of organisms that for a subset of genes, the most important sequences that regulate expression are situated not in the promoter but rather are located within introns in the first kilobase of transcribed sequences. The actual sequences involved are difficult to identify either by sequence comparisons or by deletion analysis because they are dispersed, additive, and poorly conserved. However, expression-controlling introns can be identified computationally in species with relatively small introns, based on genome-wide differences in oligomer composition between promoter-proximal and distal introns. The genes regulated by introns are often expressed in most tissues and are among the most highly expressed in the genome. The ability of some introns to strongly stimulate mRNA accumulation from several hundred nucleotides downstream of the transcription start site, even when the promoter has been deleted, reveals that our understanding of gene expression remains incomplete. It is unlikely that any diseases are caused by point mutations or small deletions that reduce the expression of an intron-regulated gene unless splicing is also affected. However, introns may be particularly useful in practical applications such as gene therapy because they strongly activate expression but only affect the transcription unit in which they are located.

The Rapid Methylation of T-DNAs Upon Agrobacterium Inoculation in Plant Leaves

Agrobacterium tumefaciens has been foundational in the development of transgenic plants for both agricultural biotechnology and plant molecular research. However, the transformation efficiency and level of transgene expression obtained for any given construct can be highly variable. These inefficiencies often require screening of many lines to find one with consistent and heritable transgene expression. Transcriptional gene silencing is known to affect transgene expression, and is associated with DNA methylation, especially of cytosines in symmetric CG and CHG contexts. While the specificity, heritability and silencing-associated effects of DNA methylation of transgene sequences have been analyzed in many stably transformed plants, the methylation status of transgene sequences in the T-DNA during the transformation process has not been well-studied. Here we used agro-infiltration of the eGFP reporter gene in Nicotiana benthamiana leaves driven by either an AtEF1α-A4 or a CaMV-35S promoter to study early T-DNA methylation patterns of these promoter sequences. The T-DNA was examined by amplicon sequencing following sodium bisulfite treatment using three different sequencing platforms: Sanger sequencing, Ion Torrent PGM, and the Illumina MiSeq. Rapid DNA methylation was detectable in each promoter region just 2-3 days post-infiltration and the levels continued to rapidly accumulate over the first week, then steadily up to 21 days later. Cytosines in an asymmetric context (CHH) were the most heavily and rapidly methylated. This suggests that early T-DNA methylation may be important in determining the epigenetic and transcriptional fate of integrated transgenes. The Illumina MiSeq platform was the most sensitive and robust way of detecting and following the methylation profiles of the T-DNA promoters. The utility of the methods was then used to show a subtle but significant difference in promoter methylation during intron-mediated enhancement. In addition, the method was able to detect an increase in promoter methylation when the eGFP reporter gene was targeted by siRNAs generated by co-infiltration of a hairpin RNAi construct.

Chloroplast Translation Elongation Factor EF-Tu/SVR11 Is Involved in var2-Mediated Leaf Variegation and Leaf Development in Arabidopsis

Chloroplasts are semiautonomous organelles, retaining their own genomes and gene expression apparatuses but controlled by nucleus genome encoded protein factors during evolution. To analyze the genetic regulatory network of FtsH-mediated chloroplast development in Arabidopsis, a set of suppressor mutants of yellow variegated (var2) have been identified. In this research, we reported the identification of another new var2 suppressor locus, SUPPRESSOR OF VARIEGATION11 (SVR11), which encodes a putative chloroplast-localized prokaryotic type translation elongation factor EF-Tu. SVR11 is likely essential to chloroplast development and plant survival. GUS activity reveals that SVR11 is abundant in the juvenile leaf tissue, lateral roots, and root tips. Interestingly, we found that SVR11 and SVR9 together regulate leaf development, including leaf margin development and cotyledon venation patterns. These findings reinforce the notion that chloroplast translation state triggers retrograde signals regulate not only chloroplast development but also leaf development.

A simple and efficient CRISPR/Cas9 platform for induction of single and multiple, heritable mutations in barley (Hordeum vulgare L.)

BACKGROUND

Genome editing of monocot plants can be accomplished by using the components of the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeat/CRISPR associated Cas9) technology specifically optimized for these types of plants. Here, we present the development of RNA-guided Cas9 system for simplex and multiplex genome editing in barley.

RESULTS

We developed a set of customizable RNA-guided Cas9 binary vectors and sgRNA modules for simplex and multiplex editing in barley. To facilitate the design of RNA-guided Cas9 constructs, the pBract derived binary vectors were adapted to Gateway cloning and only one restriction enzyme was required for construction of the sgRNA. We designed a synthetic, codon optimized Cas9 gene containing the N terminal SV40 nuclear localization signal and the UBQ10 Arabidopsis 1st intron. Two different sgRNAs were constructed for simplex editing and one polycistronic tRNA-gRNA construct (PTG) for multiplex editing using an endogenous tRNA processing system. The RNA-guided Cas9 constructs were validated in transgenic barley plants produced by Agrobacterium-mediated transformation. The highest mutation rate was observed in simplex editing of the cytokinin oxidase/dehydrogenase HvCKX1 gene, where mutations at the hvckx1 locus were detected in 88% of the screened T0 plants. We also proved the efficacy of the PTG construct in the multiplex editing of two CKX genes by obtaining 9 plants (21% of all edited plants) with mutations induced in both HvCKX1 and HvCKX3. Analysis of the T1 lines revealed that mutations in the HvCKX1 gene were transmitted to the next generation of plants. Among 220 screened T1 plants we identified 85 heterozygous and 28 homozygous mutants, most of them bearing frameshift mutations in the HvCKX1 gene. We also observed independent segregation of mutations and the Cas9-sgRNA T-DNA insert in several T1 plants. Moreover, the knockout mutations of the Nud gene generated phenotype mutants with naked grains, and the phenotypic changes were identifiable in T0 plants.

CONCLUSIONS

We demonstrated the effectiveness of an optimized RNA-guided Cas9 system that can be used for generating homozygous knockout mutants in the progeny of transgenic barely plants. This is also the first report of successful multiplex editing in barley using a tRNA processing system.

Expression pattern, antiviral role and regulation analysis of interferon-stimulated gene 15 in black seabream, Acanthopagrus schlegelii

Interferon stimulated gene 15 (ISG15) is an IFN inducible ubiquitin-like protein and plays a critical role in immune response against viral infection. In this study, an ISG15 gene (AsISG15) was cloned and characterized from the marine fish black seabream, Acanthopagrus schlegelii. The full-length cDNA of AsISG15 was 1302 bp and encoded 155 amino acids containing two ubiquitin-like motifs and a LRGG conjugation domain. Multiple alignment and phylogenetic tree showed that AsISG15 shared 31-70% amino acid identity with other known ISG15s and had a closer evolutionary relationship with teleost ISG15s. In vitro, AsISG15 expression was inducible by poly I:C, LPS and red spotted nervous necrosis virus (RGNNV) in cultured black seabream brain cells. In vivo, AsISG15 was ubiquitously expressed in all examined tissues with higher expression levels in eye and gill, and the expression was significantly up-regulated in most tissues post RGNNV infection, especially in liver, spleen and kidney. The testing of antiviral activity showed that silencing AsISG15 significantly increased RGNNV replication in RGNNV infected AsS cells, and the LRGG domain was crucial for the anti-RGNNV activity of AsISG15. By promoter-driven luciferase reporter assay, we demonstrated that two IFN-stimulated response elements within the promoter region of AsISG15 and the promoter-proximal intron were essential for AsISG15 expression. Furthermore, our results showed that the gamma-IFN activation sequence located in the intron was required for the intron mediated enhancement for AsISG15 expression. Our results would provide insights for understanding the underlying regulation mechanism of ISG15 in teleost.

Nucleosome Positioning of Intronless Genes in the Human Genome

Nucleosomes, the basic units of chromatin, are involved in transcription regulation and DNA replication. Intronless genes, which constitute 3 percent of the human genome, differ from intron-containing genes in evolution and function. Our analysis reveals that nucleosome positioning shows a distinct pattern in intronless and intron-containing genes. The nucleosome occupancy upstream of transcription start sites of intronless genes is lower than that of intron-containing genes. In contrast, high occupancy and well positioned nucleosomes are observed along the gene body of intronless genes, which is perfectly consistent with the barrier nucleosome model. Intronless genes have a significantly lower expression level than intron-containing genes and most of them are not expressed in CD4+ T cell lines and GM12878 cell lines, which results from their tissue specificity. However, the highly expressed genes are at the same expression level between the two types of genes. The highly expressed intronless genes require a higher density of RNA Pol II in an elongating state to compensate for the lack of introns. Additionally, 5' and 3' nucleosome depleted regions of highly expressed intronless genes are deeper than those of highly expressed intron-containing genes.

The RNA-Binding Protein NONO Coordinates Hepatic Adaptation to Feeding

The mechanisms by which feeding and fasting drive rhythmic gene expression for physiological adaptation to daily rhythm in nutrient availability are not well understood. Here we show that, upon feeding, the RNA-binding protein NONO accumulates within speckle-like structures in liver cell nuclei. Combining RNA-immunoprecipitation and sequencing (RIP-seq), we find that an increased number of RNAs are bound by NONO after feeding. We further show that NONO binds and regulates the rhythmicity of genes involved in nutrient metabolism post-transcriptionally. Finally, we show that disrupted rhythmicity of NONO target genes has profound metabolic impact. Indeed, NONO-deficient mice exhibit impaired glucose tolerance and lower hepatic glycogen and lipids. Accordingly, these mice shift from glucose storage to fat oxidation, and therefore remain lean throughout adulthood. In conclusion, our study demonstrates that NONO post-transcriptionally coordinates circadian mRNA expression of metabolic genes with the feeding/fasting cycle, thereby playing a critical role in energy homeostasis.

Exploring the role of Inositol 1,3,4-trisphosphate 5/6 kinase-2 (GmITPK2) as a dehydration and salinity stress regulator in Glycine max (L.) Merr. through heterologous expression in E. coli

Phytic acid (PA) is implicative in a spectrum of biochemical and physiological processes involved in plant stress response. Inositol 1,3,4, Tris phosphate 5/6 kinase (ITPK), a polyphosphate kinase that converts Inositol 1,3,4 trisphosphate to Inositol 1,3,4,5/6 tetra phosphate, averting the inositol phosphate pool towards PA biosynthesis, is a key regulator that exists in four different isoforms in soybean. In the present study, in-silico analysis of the promoter region of ITPKs was done and among the four isoforms, promoter region of GmITPK2 showed the presence of two MYB binding elements for drought inducibility and one for ABA response. Expression profiling through qRT-PCR under drought and salinity stress showed higher expression of GmITPK2 isoform compared to the other members of the family. The study revealed GmITPK2 as an early dehydration responsive gene which is also induced by dehydration and exogenous treatment with ABA. To evaluate the osmo-protective role of GmITPK2, attempts were made to assess the bacterial growth on Luria Broth media containing 200 mM NaCl, 16% PEG and 100 μM ABA, individually. The transformed E. coli BL21 (DE3) cells harbouring the GmITPK2 gene depicted better growth on the media compared to the bacterial cells containing the vector alone. Similarly, the growth of the transformed cells in the liquid media containing 200 mM NaCl, 16% PEG and 100 μM ABA showed higher absorbance at 600 nm compared to control, at different time intervals. The GmITPK2 recombinant E. coli cells showing tolerance to drought and salinity thus demonstrated the functional redundancy of the gene across taxa. The purity and specificity of the recombinant protein was assessed and confirmed through PAGE showing a band of ∼35 kDa on western blotting using Anti- Penta His- HRP conjugate antibody. To the best of our knowledge, the present study is the first report exemplifying the role of GmITPK2 isoform in drought and salinity tolerance in soybean.

Patterns of variation in cis-regulatory regions: examining evidence of purifying selection

BACKGROUND

With only 2 % of the human genome consisting of protein coding genes, functionality across the rest of the genome has been the subject of much debate. This has gained further impetus in recent years due to a rapidly growing catalogue of genomic elements, based primarily on biochemical signatures (e.g. the ENCODE project). While the assessment of functionality is a complex task, the presence of selection acting on a genomic region is a strong indicator of importance. In this study, we apply population genetic methods to investigate signals overlaying several classes of regulatory elements.

RESULTS

We disentangle signals of purifying selection acting directly on regulatory elements from the confounding factors of demography and purifying selection linked to e.g. nearby protein coding regions. We confirm the importance of regulatory regions proximal to coding sequence, while also finding differential levels of selection at distal regions. We note differences in purifying selection among transcription factor families. Signals of constraint at some genomic classes were also strongly dependent on their physical location relative to coding sequence. In addition, levels of selection efficacy across genomic classes differed between African and non-African populations.

CONCLUSIONS

In order to assign a valid signal of selection to a particular class of genomic sequence, we show that it is crucial to isolate the signal by accounting for the effects of demography and linked-purifying selection. Our study highlights the intricate interplay of factors affecting signals of selection on functional elements.

Glutamate synthases from conifers: gene structure and phylogenetic studies

BACKGROUND

Plants synthesize glutamate from ammonium by the combined activity of the enzymes glutamine synthetase (GS) and glutamate synthase (GOGAT) through the glutamate synthase cycle. In plants, there are two forms of glutamate synthases that differ in their electron donors, NADH-GOGAT (EC 1.4.1.14) and Fd-GOGAT (EC 1.4.7.1), which have differential roles either in primary ammonia assimilation or in the reassimilation of ammonium from different catabolic processes. Glutamate synthases are complex iron-sulfur flavoproteins containing functional domains involved in the control and coordination of their catalytic activities in annual plants. In conifers, partial cDNA sequences for GOGATs have been isolated and used for gene expression studies. However, knowledge of the gene structure and of phylogenetic relationships with other plant enzymes is quite scant.

RESULTS

Technological advances in conifer megagenomes sequencing have made it possible to obtain full-length cDNA sequences encoding Fd- and NADH-GOGAT from maritime pine, as well as BAC clones containing sequences for NADH-GOGAT and Fd-GOGAT genes. In the current study, we studied the genomic organization of pine GOGAT genes, the size of their exons/introns, copy numbers in the pine genome and relationships with other plant genes. Phylogenetic analysis was performed, and the degree of preservation and dissimilarity of key domains for the catalytic activities of these enzymes in different taxa were determined.

CONCLUSIONS

Fd- and NADH-GOGAT are encoded by single-copy genes in the maritime pine genome. The Fd-GOGAT gene is extremely large spanning more than 330 kb and the presence of very long introns highlights the important contribution of LTR retrotransposons to the gene size in conifers. In contrast, the structure of the NADH-GOGAT gene is similar to the orthologous genes in angiosperms. Our phylogenetic analysis indicates that these two genes had different origins during plant evolution. The results provide new insights into the structure and molecular evolution of these essential genes.

Arabidopsis thaliana MLO genes are expressed in discrete domains during reproductive development

MLOs in Plant Reproduction. The MILDEW RESISTANCE LOCUS-O (MLO) protein family, comprised of 15 members, plays roles in diverse cell-cell communication processes such as powdery mildew susceptibility, root thigmomorphogenesis, and pollen tube reception. The NORTIA (NTA, AtMLO7) gene is expressed in the synergid cells of the female gametophyte where it functions in intercellular communication with the pollen tube. Discrepancies between previously published promoter::GUS and promoter::gene-GUS constructs expression patterns led us to explore the regulation of NTA expression. Here we found via NTA_pro::gNTA-GUS truncations that sequences within the NTA gene negatively regulate its expression in the stomata and carpel walls. This led to the hypothesis that other MLO family members may also have additional regulatory sequences within the gene. MLO_pro::gMLO-GUS constructs were examined for each family member focusing specifically on flowers in order to determine whether other MLOs could play a role in reproductive cell-cell communication. Notably, several MLOs were expressed in the pollen, in the stigma, in the pollinated style, and in the synergids and central cell. These findings indicate that other MLOs in addition to NTA could play a role in reproduction. Previous studies on the MLO family showed that phylogenetically related MLOs had redundant functions in powdery mildew infection and root thigmomorphogenesis; however, MLO expression in reproductive tissues did not strictly follow phylogenetic relationships, indicating that MLOs from different evolutionary origins may have been recruited for function in sexual reproduction.

How introns enhance gene expression

In many eukaryotes, including mammals, plants, yeast, and insects, introns can increase gene expression without functioning as a binding site for transcription factors. This phenomenon was termed 'intron-mediated enhancement'. Introns can increase transcript levels by affecting the rate of transcription, nuclear export, and transcript stability. Moreover, introns can also increase the efficiency of mRNA translation. This review discusses the current knowledge about these processes. The role of splicing in IME and the significance of intron position relative to the sites of transcription and translation initiation are elaborated. Particular emphasis is placed on the question why different introns, present at the same location of the same genes and spliced at a similar high efficiency, can have very different impacts on expression - from almost no effect to considerable stimulation. This situation can be at least partly accounted for by the identification of splicing-unrelated intronic elements with a special ability to enhance mRNA accumulation or translational efficiency. The many factors that could lead to the large variation observed between the impact of introns in different genes and experimental systems are highlighted. It is suggested that there is no sole, definite answer to the question "how do introns enhance gene expression". Rather, each intron-gene combination might undergo its own unique mixture of processes that lead to the perceptible outcome.

Predictive Models of Spatial Transcriptional Response to High Salinity

Plants are exposed to a variety of environmental conditions, and their ability to respond to environmental variation depends on the proper regulation of gene expression in an organ-, tissue-, and cell type-specific manner. Although our knowledge of how stress responses are regulated is accumulating, a genome-wide model of how plant transcription factors (TFs) and cis-regulatory elements control spatially specific stress response has yet to emerge. Using Arabidopsis (Arabidopsis thaliana) as a model, we identified a set of 1,894 putative cis-regulatory elements (pCREs) that are associated with high-salinity (salt) up-regulated genes in the root or the shoot. We used these pCREs to develop computational models that can better predict salt up-regulated genes in the root and shoot compared with models based on known TF binding motifs. In addition, we incorporated TF binding sites identified via large-scale in vitro assays, chromatin accessibility, evolutionary conservation, and pCRE combinatorial relationships in machine learning models and found that only consideration of pCRE combinations led to better performance in salt up-regulation prediction in the root and shoot. Our results suggest that the plant organ transcriptional response to high salinity is regulated by a core set of pCREs and provide a genome-wide view of the cis-regulatory code of plant spatial transcriptional responses to environmental stress.

Intron DNA Sequences Can Be More Important Than the Proximal Promoter in Determining the Site of Transcript Initiation

To more precisely define the positions from which certain intronic regulatory sequences increase mRNA accumulation, the effect of a UBIQUITIN intron on gene expression was tested from six different positions surrounding the transcription start site (TSS) of a reporter gene fusion in Arabidopsis thaliana The intron increased expression from all transcribed positions but had no effect when upstream of the 5'-most TSS. While this implies that the intron must be transcribed to increase expression, the TSS changed when the intron was located in the 5'-untranslated region (UTR), suggesting that the intron affects transcription initiation. Remarkably, deleting 303 nucleotides of the promoter including all known TSSs and all but 18 nucleotides of the 5'-UTR had virtually no effect on the level of gene expression as long as an intron containing stimulatory sequences was included. Instead, transcription was initiated in normally untranscribed sequences the same distance upstream of the intron as when the promoter was intact. These results suggest that certain intronic DNA sequences play unexpectedly large roles in directing transcription initiation and constitute a previously unrecognized type of downstream regulatory element for genes transcribed by RNA polymerase II.

Synthetic introns help identify sequences in the 5' UTR intron of the Glycine max polyubiquitin (Gmubi) promoter that give increased promoter activity

MAIN CONCLUSION

Specific sequences within the leader intron of a soybean polyubiquitin gene stimulated gene expression when placed either within a synthetic intron or upstream of a core promoter. The intron in the 5' untranslated region of the soybean polyubiquitin promoter, Gmubi, seems to contribute to the high activity of this promoter. To identify the stimulatory sequences within the intron, ten different sequential intronic sequences of 40 nt were isolated, cloned as tetrameric repeats and placed upstream of a minimal cauliflower mosaic virus 35S (35S) core promoter, which was used to control expression of the green fluorescent protein. Intron fragment tetramers were also cloned within a modified, native intron, creating a Synthetic INtron Cassette (SINC), which was then placed downstream of Gmubi and 35S core promoters. Intron fragment tetramers and SINC constructs were evaluated using transient expression in lima bean cotyledons and stable expression in soybean hairy roots. Intron fragments, used as tetramers upstream of the 35S core promoter, yielded up to 80 times higher expression than the core promoter in transient expression analyses and ten times higher expression in stably transformed hairy roots. Tetrameric intronic fragments, cloned downstream of the Gmubi and 35S core promoters and within the synthetic intron, also yielded increased transient and stable GFP expression that was up to 4 times higher than Gmubi alone and up to 40 times higher than the 35S core promoter alone. These intron fragments contain sequences that seem to act as promoter regulatory elements and may contribute to the increased expression observed with this native strong promoter. Intron regulatory elements and synthetic introns may provide additional tools for increasing transgene expression in plants.

Plant Actin-Depolymerizing Factors Possess Opposing Biochemical Properties Arising from Key Amino Acid Changes throughout Evolution

Functional divergence in paralogs is an important genetic source of evolutionary innovation. Actin-depolymerizing factors (ADFs) are among the most important actin binding proteins and are involved in generating and remodeling actin cytoskeletal architecture via their conserved F-actin severing or depolymerizing activity. In plants, ADFs coevolved with actin, but their biochemical properties are diverse. Unfortunately, the biochemical function of most plant ADFs and the potential mechanisms of their functional divergence remain unclear. Here, in vitro biochemical analyses demonstrated that all 11 ADF genes in Arabidopsis thaliana exhibit opposing biochemical properties. Subclass III ADFs evolved F-actin bundling (B-type) function from conserved F-actin depolymerizing (D-type) function, and subclass I ADFs have enhanced D-type function. By tracking historical mutation sites on ancestral proteins, several fundamental amino acid residues affecting the biochemical functions of these proteins were identified in Arabidopsis and various plants, suggesting that the biochemical divergence of ADFs has been conserved during the evolution of angiosperm plants. Importantly, N-terminal extensions on subclass III ADFs that arose from intron-sliding events are indispensable for the alteration of D-type to B-type function. We conclude that the evolution of these N-terminal extensions and several conserved mutations produced the diverse biochemical functions of plant ADFs from a putative ancestor.