Plant-specific environmental and developmental signals regulate the mismatch repair protein MSH6 in Arabidopsis thaliana

The DNA mismatch repair (MMR) is a postreplicative system that guarantees genomic stability by correcting mispaired and unpaired nucleotides. In eukaryotic nuclei, MMR is initiated by the binding of heterodimeric MutS homologue (MSH) complexes to the DNA error or lesion. Among these proteins, MSH2-MSH6 is the most abundant heterodimer. Even though the MMR mechanism and proteins are highly conserved throughout evolution, physiological differences between species can lead to different regulatory features. Here, we investigated how light, sugar, and/or hormones modulate Arabidopsis thaliana MSH6 expression pattern. We first characterized the promoter region of MSH6. Phylogenetic shadowing revealed three highly conserved regions. These regions were analyzed by the generation of deletion constructs of the MSH6 full-length promoter fused to the β-glucuronidase (GUS) gene. Combined, our in silico and genetic analyses revealed that a 121-bp promoter fragment was necessary for MSH6 expression and contained potential cis-acting elements involved in light- and hormone-responsive gene expression. Accordingly, light exposure or sugar treatment of four-day old A. thaliana seedlings triggered an upregulation of MSH6 in shoot and root apical meristems. Appropriately, MSH6 was also induced by the stem cell inducer WUSCHEL. Further, the stimulatory effect of light was dependent on the presence of phyA. In addition, treatment of seedlings with auxin or cytokinin also caused an upregulation of MSH6 under darkness. Consistent with auxin signals, MSH6 expression was suppressed in the GATA23 RNAi line compared with the wild type. Our results provide evidence that endogenous factors and environmental signals controlling plant growth and development regulate the MSH6 protein in A. thaliana.

Prediction and analysis of cis-regulatory elements in Dorsal and Ventral patterning genes of Tribolium castaneum and its comparison with Drosophila melanogaster

Insect embryonic development and morphology are characterized by their anterior-posterior and dorsal-ventral (DV) patterning. In Drosophila embryos, DV patterning is mediated by a dorsal protein gradient which activates twist and snail proteins, the important regulators of DV patterning. To activate or repress gene expression, some regulatory proteins bind in clusters to their target gene at sites known as cis-regulatory elements or enhancers. To understand how variations in gene expression in different lineages might lead to different phenotypes, it is necessary to understand enhancers and their evolution. Drosophila melanogaster has been widely studied to understand the interactions between transcription factors and the transcription factor binding sites. Tribolium castaneum is an upcoming model animal which is catching the interest of biologists and the research on the enhancer mechanisms in the insect's axes patterning is still in infancy. Therefore, the current study was designed to compare the enhancers of DV patterning in the two insect species. The sequences of ten proteins involved in DV patterning of D. melanogaster were obtained from Flybase. The protein sequences of T. castaneum orthologous to those obtained from D. melanogaster were acquired from NCBI BLAST, and these were then converted to DNA sequences which were modified by adding 20 kb sequences both upstream and downstream to the gene. These modified sequences were used for further analysis. Bioinformatics tools (Cluster-Buster and MCAST) were used to search for clusters of binding sites (enhancers) in the modified DV genes. The results obtained showed that the transcription factors in Drosophila melanogaster and Tribolium castaneum are nearly identical; however, the number of binding sites varies between the two species, indicating transcription factor binding site evolution, as predicted by two different computational tools. It was observed that dorsal, twist, snail, zelda, and Supressor of Hairless are the transcription factors responsible for the regulation of DV patterning in the two insect species.

Modulation in gene expression and enzyme activity suggested the roles of monodehydroascorbate reductase in development and stress response in bread wheat

Monodehydroascorbate reductase (MDHAR) is a crucial enzymatic antioxidant of the ascorbate-glutathione pathway involved in reactive oxygen species scavenging. Herein, we identified 15 TaMDHAR genes in bread wheat. Phylogenetic analysis revealed their clustering into three groups, which are also related to the subcellular localization in the peroxisome matrix, peroxisome membrane, and chloroplast. Each TaMDHAR protein consisted of two conserved domains; Pyr_redox and Pyr_redox_2 of the pyridine nucleotide disulfide oxidoreductase family. The occurrence of diverse groups of cis-regulatory elements in the promoter region and their interaction with numerous transcription factors suggest assorted functions of TaMDHARs in growth and development and in light, phytohormones, and stress responses. Expression analysis in various tissues further revealed their importance in vegetative and reproductive development. In addition, the differential gene expression and enhanced enzyme activity during drought, heat, and salt treatments exposed their role in abiotic stress response. Interaction of MDHARs with various antioxidant enzymes and biochemicals related to the ascorbate-glutathione cycle exposed their synchronized functioning. Interaction with auxin indicated the probability of cross-talk between antioxidants and auxin signaling. The miR168a, miR169, miR172 and others interaction with various TaMDHARs further directed their association with developmental processes and stress responses. The current study provides extensive information about the importance of TaMDHARs, moreover, the precise role of each gene needs to be established in future studies.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

Methylation-directed regulatory networks determine enhancing and silencing of mutation disease driver genes and explain inter-patient expression variation

BACKGROUND

Common diseases manifest differentially between patients, but the genetic origin of this variation remains unclear. To explore possible involvement of gene transcriptional-variation, we produce a DNA methylation-oriented, driver-gene-wide dataset of regulatory elements in human glioblastomas and study their effect on inter-patient gene expression variation.

RESULTS

In 175 of 177 analyzed gene regulatory domains, transcriptional enhancers and silencers are intermixed. Under experimental conditions, DNA methylation induces enhancers to alter their enhancing effects or convert into silencers, while silencers are affected inversely. High-resolution mapping of the association between DNA methylation and gene expression in intact genomes reveals methylation-related regulatory units (average size = 915.1 base-pairs). Upon increased methylation of these units, their target-genes either increased or decreased in expression. Gene-enhancing and silencing units constitute cis-regulatory networks of genes. Mathematical modeling of the networks highlights indicative methylation sites, which signified the effect of key regulatory units, and add up to make the overall transcriptional effect of the network. Methylation variation in these sites effectively describe inter-patient expression variation and, compared with DNA sequence-alterations, appears as a major contributor of gene-expression variation among glioblastoma patients.

CONCLUSIONS

We describe complex cis-regulatory networks, which determine gene expression by summing the effects of positive and negative transcriptional inputs. In these networks, DNA methylation induces both enhancing and silencing effects, depending on the context. The revealed mechanism sheds light on the regulatory role of DNA methylation, explains inter-individual gene-expression variation, and opens the way for monitoring the driving forces behind deferential courses of cancer and other diseases.

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

Identification of oil palm cis-regulatory elements based on DNA free energy and single nucleotide polymorphism density

Transcription control through cis-regulatory elements (CREs) is one of important regulators of gene expression. This study aimed to identify the location of CREs in oil palm (Elaeis guineensis Jacq.) using the combination of DNA free energy and single nucleotide polymorphism (SNP) density approaches. Promoter region sequences were extracted oil palm genome spanning from 1500 nucleotides (nt) upstream to 1000 nt downstream of every annotated transcription start sites (TSS). Free energy profiles of each promoter region were calculated using PromPredict software. Raw reads from the deep sequencing of 59 oil palm origins were used to calculate SNP density of each promoter region. The result showed that the average free energy (AFE) on the upstream region of TSS is about 1.5 kcal/mol higher compared to the downstream region. Using DNA free energy method, 16,281 regions of CREs were predicted. Most of predicted CREs was located between 1 and 500 nt upstream of TSS. Anti-correlation pattern between free energy and SNP density was observed on the predicted regions of CREs. This anti-correlated pattern was also observed on an experimentally determined promoter of the oil palm metallothionein gene, EgMSP1. Considering the increasing use of promoter information on plant biotechnology, an easy and accurate promoter prediction using the combination of free energy and SNP density method could be recommended.

Frontiers and techniques in plant gene regulation

Understanding plant gene regulation has been a priority for generations of plant scientists. However, due to its complex nature, the regulatory code governing plant gene expression has yet to be deciphered comprehensively. Recently developed methods-often relying on next-generation sequencing technology and state-of-the-art computational approaches-have started to further our understanding of the gene regulatory logic used by plants. In this review, we discuss these methods and the insights into the regulatory code of plants that they can yield.

Genome-wide association study unravels genomic regions associated with chlorophyll fluorescence parameters in wheat (Triticum aestivum L.) under different sowing conditions

KEY MESSAGE

Genome-wide association study identified 205 significant marker-trait associations for chlorophyll fluorescence parameters in wheat. Candidate gene mining, in silico expression, and promoter analyses revealed the potential candidate genes associated with the studied parameters. The present study investigated the effect of varied sowing conditions (viz., early, timely, and late) on different chlorophyll fluorescence parameters in diverse wheat germplasm set comprising of 198 lines over two cropping seasons (2020-2021 and 2021-2022). Further, a genome-wide association study was conducted to identify potential genomic regions associated with these parameters. The results revealed significant impacts of sowing conditions on all fluorescence parameters, with the maximum and minimum effects on FI (26.64%) and FV/FM (2.12%), respectively. Among the 205 marker-trait associations (MTAs) identified, 11 high-confidence MTAs were chosen, exhibiting substantial effects on multiple fluorescence parameters, and each explaining more than 10% of the phenotypic variation. Through gene mining of genomic regions encompassing high-confidence MTAs, we identified a total of 626 unique gene models. In silico expression analysis revealed 42 genes with an expression value exceeding 2 TPM. Among them, 10 genes were identified as potential candidate genes with functional relevance to enhanced photosynthetic efficiency. These genes mainly encoded for the following important proteins/products-ankyrin repeat protein, 2Fe-2S ferredoxin-type iron-sulfur-binding domain, NADH-ubiquinone reductase complex-1 MLRQ subunit, oxidoreductase FAD/NAD(P)-binding, photosystem-I PsaF, and protein kinases. Promoter analysis revealed the presence of light-responsive (viz., GT1-motif, TCCC-motif, I-box, GT1-motif, TCT-motif, and SP-1) and stress-responsive (viz., ABRE, AuxRR-core, GARE-motif, and ARE) cis-regulatory elements, which may be involved in the regulation of identified putative candidate genes. Findings from this study could directly help wheat breeders in selecting lines with favorable alleles for chlorophyll fluorescence, while the identified markers will facilitate marker-assisted selection of potential genomic regions for improved photosynthesis.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

NOMe-HiC: joint profiling of genetic variant, DNA methylation, chromatin accessibility, and 3D genome in the same DNA molecule

Cis-regulatory elements are coordinated to regulate the expression of their targeted genes. However, the joint measurement of cis-regulatory elements' activities and their interactions in spatial proximity is limited by the current sequencing approaches. We describe a method, NOMe-HiC, which simultaneously captures single-nucleotide polymorphisms, DNA methylation, chromatin accessibility (GpC methyltransferase footprints), and chromosome conformation changes from the same DNA molecule, together with the transcriptome, in a single assay. NOMe-HiC shows high concordance with state-of-the-art mono-omic assays across different molecular measurements and reveals coordinated chromatin accessibility at distal genomic segments in spatial proximity and novel types of long-range allele-specific chromatin accessibility.

Smar2C2: A Simple and Efficient Protocol for the Identification of Transcription Start Sites

Promoters and the noncoding sequences that drive their function are fundamental aspects of genes that are critical to their regulation. The transcription preinitiation complex binds and assembles on promoters where it facilitates transcription. The transcription start site (TSS) is located downstream of the promoter sequence and is defined as the location in the genome where polymerase begins transcribing DNA into RNA. Knowing the location of TSSs is useful for annotation of genes, identification of non-coding sequences important to gene regulation, detection of alternative TSSs, and understanding of 5' UTR content. Several existing techniques make it possible to accurately identify TSSs, but are often difficult to perform experimentally, require large amounts of input RNA, or are unable to identify a large number of TSSs from a single sample. Many of these protocols take advantage of template switching reverse transcriptases (TSRTs), which reliably place an adaptor at the 5' end of a first strand synthesis of cDNA. Here, we introduce a protocol that exploits TSRT activity combined with rolling circle amplification to identify TSSs with several unique advantages over existing methods. Sequence adaptors are placed on the 5' and 3' end of the full-length cDNA copy of a transcript. A splint compatible with those adaptors is then used to circularize the full-length cDNA. Linear DNA containing concatemers of the cDNA are generated using rolling circle amplification, and a sequencing library is formed by fragmenting the concatemers. This protocol is straightforward to execute, requiring limited bench time with relatively stable reagents. Using extremely low amounts of RNA input, this protocol produces large numbers of accurate, deduplicated TSSs genome wide. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Splint generation Basic Protocol 2: RNA extraction Basic Protocol 3: cDNA synthesis Basic Protocol 4: cDNA circularization and amplification Basic Protocol 5: Library generation.

Synthetic sub-genomic transcript promoter from Horseradish Latent Virus (HRLV)

We characterized an efficient chimeric sub-genomic transcript promoter from Horseradish Latent Virus, FHS4, active in both dicot and monocot plants, and it could be a potential tool for plant biotechnology. Plant pararetroviruses are a rich source of novel plant promoters widely used for biotechnological applications. Here, we comprehensively characterized a unique sub-genomic transcript (Sgt) promoter of Horseradish Latent Virus (HRLV) and identified a fragment (HS4; - 340 to + 10; 351 bp) that showed the highest expression of reporter genes in both transient and transgenic assays as evidenced by biochemical, histochemical GUS reporter assay and transcript analysis of uidA gene by qRT-PCR. Phylogenetic analysis showed that the HSgt promoter was closely related to the sub-genomic promoter of the Cauliflower Mosaic Virus (CaMV19S). We found that the as-1 element and W-box played an important role in the transcriptional activity of the HS4 promoter. Furthermore, the HS4 promoter was also induced by salicylic acid. Alongside, we enhanced the activity of the HS4 promoter by coupling the enhancer region from Figwort Mosaic Virus (FMV) promoter to the upstream region of it. This hybrid promoter FHS4 was around 1.1 times stronger than the most commonly used promoter, 35S (Cauliflower Mosaic Virus full-length transcript promoter), and was efficient in driving reporter genes in both dicot and monocot plants. Subsequently, transgenic tobacco plants expressing an anti-microbial peptide BrLTP2.1 (Brassica rapa lipid transport protein 2.1), under the control of the FHS4 promoter, were developed. The in vitro anti-fungal assay revealed that the plant-derived BrLTP2.1 protein driven by an FHS4 promoter manifested increased resistance against an important plant fungal pathogen, Alternaria alternata. Finally, we concluded that the FHS4 promoter can be used as an alternative to the 35S promoter and has a high potential to become an efficient tool in plant biotechnology.

Glandular trichome specificity of menthol biosynthesis pathway gene promoters from Mentha × piperita

Several cis-elements including Myb-binding motifs together confer glandular trichome specificity as revealed from heterologous expression and analysis of menthol biosynthesis pathway gene promoters. Glandular Trichomes (GTs) are result of division of epidermal cells that produce diverse metabolites. Species of mint family are important for their essential oil containing many high-value terpenoids, biosynthesized and stored in these GTs. Hence, GTs constitute attractive targets for metabolic engineering and GT-specific promoters are important. In this investigation, the upstream regions of the Mentha × piperita menthol biosynthetic pathway genes (-)-limonene synthase, (-)-P450 limonene-3- hydroxylase, (-)-trans-isopiperitenol dehydrogenase, (-)-Isopiperitenone reductase, ( +)-Pulegone reductase, (-)-Menthone reductase/ (-)-Menthol dehydrogenase and a branched pathway gene ( +)-menthofuran synthase were isolated and characterized. These fragments, fused to β-glucuronidase (GUS) reporter gene of pBI101 binary vector, are able to drive high level gene expression in transgenic tobacco trichomes with strong signals in GTs, except for (-)-Isopiperitenone reductase. The GT-enriched tissue from transformed plants were analysed for GUS enzyme activity and RNA expression which correlates the GUS staining. To characterize the cis-elements responsible for GT-specific expression, a series of 5' deletion constructs for MpPLS and MpPMFS were cloned and analysed in stable transgenic tobacco lines. The specificity of trichome expression was located to -  797 to-  598 bp sequence for (-)-limonene synthase and-  629 to -   530 bp for ( +)-menthofuran synthase promoters containing specific Myb-binding motifs in addition to other unique motifs described for developmental regulation without any defined pattern. All other pathway promoters also recruits specific but different Myb factors as indicated by this analysis.

Mysteries of gene regulation: Promoters are not the sole triggers of gene expression

•

TF binding peaks were widely distributed in nonpromoters, especially downstream regions of transcription termination sites.

•

Exons of non-coding regions were the prominent regions of TF binding.

•

TAD boundaries were colocalized with activating histone marks and TF binding.

•

Genes with distinct functions demonstrated substantially different behaviors in cis-regulation and epigenetic signatures.

Cis-regulatory elements of promoters are essential for gene regulation by transcription factors (TFs). However, the regulatory roles of nonpromoter regions, TFs, and epigenetic marks remain poorly understood in plants. In this study, we characterized the cis-regulatory regions of 53 TFs and 19 histone marks in 328 chromatin immunoprecipitation (ChIP-seq) datasets from Arabidopsis. The genome-wide maps indicated that both promoters and regions around the transcription termination sites of protein-coding genes recruit the most TFs. The maps also revealed a diverse of histone combinations. The analysis suggested that exons play critical roles in the regulation of non-coding genes. Additionally, comparative analysis between heat-stress-responsive and nonresponsive genes indicated that the genes with distinct functions also exhibited substantial differences in cis-regulatory regions, histone regulation, and topologically associating domain (TAD) boundary organization. By integrating multiple high-throughput sequencing datasets, this study generated regulatory models for protein-coding genes, non-coding genes, and TAD boundaries to explain the complexity of transcriptional regulation.

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

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

Dynamic enhancer transcription associates with reprogramming of immune genes during pattern triggered immunity in Arabidopsis

Background

Enhancers are cis-regulatory elements present in eukaryote genomes, which constitute indispensable determinants of gene regulation by governing the spatiotemporal and quantitative expression dynamics of target genes, and are involved in multiple life processes, for instance during development and disease states. The importance of enhancer activity has additionally been highlighted for immune responses in animals and plants; however, the dynamics of enhancer activities and molecular functions in plant innate immunity are largely unknown. Here, we investigated the involvement of distal enhancers in early innate immunity in Arabidopsis thaliana.

Results

A group of putative distal enhancers producing low-abundance transcripts either unidirectionally or bidirectionally are identified. We show that enhancer transcripts are dynamically modulated in plant immunity triggered by microbe-associated molecular patterns and are strongly correlated with open chromatin, low levels of methylated DNA, and increases in RNA polymerase II targeting and acetylated histone marks. Dynamic enhancer transcription is correlated with target early immune gene expression patterns. Cis motifs that are bound by immune-related transcription factors, such as WRKYs and SARD1, are highly enriched within upregulated enhancers. Moreover, a subset of core pattern-induced enhancers are upregulated by multiple patterns from diverse pathogens. The expression dynamics of putative immunity-related enhancers and the importance of WRKY binding motifs for enhancer function were also validated.

Conclusions

Our study demonstrates the general occurrence of enhancer transcription in plants and provides novel information on the distal regulatory landscape during early plant innate immunity, providing new insights into immune gene regulation and ultimately improving the mechanistic understanding of the plant immune system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01362-8.

Cross-talk between enhancers, structural elements and activating transcription factors maintains the 3D architecture and expression of the CFTR gene

Robust protocols to examine 3D chromatin structure have greatly advanced knowledge of gene regulatory mechanisms. Here we focus on the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which provides a paradigm for validating models of gene regulation built upon genome-wide analysis. We examine the mechanisms by which multiple cis-regulatory elements (CREs) at the CFTR gene coordinate its expression in intestinal epithelial cells. Using CRISPR/Cas9 to remove CREs, individually and in tandem, followed by assays of gene expression and higher-order chromatin structure (4C-seq), we reveal the cross-talk and dependency of two cell-specific intronic enhancers. The results suggest a mechanism whereby the locus responds when CREs are lost, which may involve activating transcription factors such as FOXA2. Also, by removing the 5' topologically-associating domain (TAD) boundary, we illustrate its impact on CFTR gene expression and architecture. These data suggest a multi-layered regulatory hierarchy that is highly sensitive to perturbations.

Epigenomics analysis of miRNA cis-regulatory elements in pig muscle and fat tissues

Although large-scale and accurate identification of cis-regulatory elements on pig protein-coding and long non-coding genes has been reported, similar study on pig miRNAs is still lacking. Here, we systematically characterized the cis-regulatory elements of pig miRNAs in muscle and fat by adopting miRNAomes, ChIP-seq, ATAC-seq, RNA-seq and Hi-C data. In total, the cis-regulatory elements of 257 (85.95%) expressed miRNAs including 226 known and 31 novel miRNAs were identified. Especially, the miRNAs associated with super-enhancers, active promoters, and "A" compartment were significantly higher than those associated by typical enhancers, prompters without H3K27ac, and "B" compartment, respectively. The tissue specific transcription factors were the primary determination of core miRNA expression pattern in muscle and fat. Moreover, the miRNA promoters are more evolutionarily conserved than miRNA enhancers, like other type genes. Our study adds additional important information to existing pig epigenetic data and provides essential resource for future in-depth investigation of pig epigenetics.

In Silico and Transcription Analysis of Trehalose-6-phosphate Phosphatase Gene Family of Wheat: Trehalose Synthesis Genes Contribute to Salinity, Drought Stress and Leaf Senescence

Trehalose-6-phosphate phosphatase (TPP) genes take part in trehalose metabolism and also in stress tolerance, which has been well documented in many species but poorly understood in wheat. The present research has identified a family of 31 TPP genes in Triticum aestivum L. through homology searches and classified them into five clades by phylogenetic tree analysis, providing evidence of an evolutionary status with Hordeum vulgare, Brachypodium distachyon and Oryza sativa. The exon-intron distribution revealed a discrete evolutionary history and projected possible gene duplication occurrences. Furthermore, different computational approaches were used to analyze the physical and chemical properties, conserved domains and motifs, subcellular and chromosomal localization, and three-dimensional (3-D) protein structures. Cis-regulatory elements (CREs) analysis predicted that TaTPP promoters consist of CREs related to plant growth and development, hormones, and stress. Transcriptional analysis revealed that the transcription levels of TaTPPs were variable in different developmental stages and organs. In addition, qRT-PCR analysis showed that different TaTPPs were induced under salt and drought stresses and during leaf senescence. Therefore, the findings of the present study give fundamental genomic information and possible biological functions of the TaTPP gene family in wheat and will provide the path for a better understanding of TaTPPs involvement in wheat developmental processes, stress tolerance, and leaf senescence.

Inducible CRISPRa screen identifies putative enhancers

Enhancers are critical cis-regulatory elements that regulate spatiotemporal gene expression and control cell fates. However, the identification of enhancers in native cellular contexts still remains a challenge. Here, we develop an inducible CRISPR activation (CRISPRa) system by transgenic expression of doxycycline (Dox)-inducible dCas9-VPR in mouse embryonic stem cells (iVPR ESC). With this line, a simple introduction of specific guide RNAs targeting promoters or enhancers allows us to realize the effect of CRISPRa in an inducible, reversible, and Dox concentration-dependent manner. Taking advantage of this system, we induce tiled CRISPRa across genomic regions (105 kilobases) surrounding T (Brachyury), one of the key mesodermal development regulator genes. Moreover, we identify several CRISPRa-responsive elements with chromatin features of putative enhancers, including a region the homologous sequence in which humans harbors a body height risk variant. Genetic deletion of this region in ESC does affect subsequent T gene activation and osteogenic differentiation. Therefore, our inducible CRISPRa ESC line provides a convenient platform for high-throughput screens of putative enhancers.

Promoter analysis of the SPATULA (FvSPT) and SPIRAL (FvSPR) genes in the woodland diploid strawberry (Fragaria vesca L.)

The aim of this study was to identify transcription factor (TF) binding sites and cis-regulatory elements (CREs) on the promoters of FvSPR1-like2 (SPIRAL) and FvSPT (SPATULA) genes in the woodland diploid strawberry (Fragaria vesca L.). We identified: (1) MYB59, WRKY25 and WRKY8 TFs which play a role in ethylene signaling; (2) ARF family of TFs which play a role in ARF-mediated auxin signaling on the promoter of FvSPR1-like2 gene; (3) ARR family of TFs which play a role in cytokinin signaling; (4) ERF family of TFs which play a role in ethylene signaling on the promoter of FvSPT. This bioinformatic analysis of TFs and CREs may provide a better understanding of the function of genes involved in, and the mechanism underlying, non-climateric ripening during strawberry fruit maturation.

High promoter sequence variation in subgroup 6 members of R2R3-MYB genes is involved in different floral anthocyanin color patterns in Lilium spp

The spatially and temporally distinct expression of R2R3-MYB positive regulators is among the major mechanisms that create various anthocyanin color patterns in many flowers. However, we do not know how these positive regulators have gained different expression profiles. In the Asiatic hybrid lily 'Lollypop' (derived from the crosses of species belonging to Sinomartagon/Daurolirion section), MYB12 and MYB19S regulate the pigmentation at whole tepals and raised tepal spots, respectively. In the Oriental hybrid lily 'Sorbonne' (derived from the crosses of species belonging to the Archelirion section), MYB12 regulates both whole tepal and raised spot pigmentation. The genes have similar amino acid sequences with similar protein functions but exhibit different expression profiles in lily flowers. As promoters are among the most significant factors affecting gene expression profiles, their promoter sequences were determined in this study. The three genes had very different promoter sequences, and putative cis-regulatory elements were not conserved in numbers or order. To further confirm the promoter functions, tobacco plants were transformed with native promoter-driven MYB12 or MYB19S genes of 'Lollypop.' Expression levels of MYB12 were higher in corolla tubes than in lobes, while those of MYB19S were higher in corolla lobes than in tubes. Thus, the diverse promoter functions were likely to be the leading causes of their different expression profiles and generation of unique color patterns. Finally, the history of R2R3-MYB gene establishment during lily evolution was estimated using sequence data.

Functional divergence of Brassica napus BnaABI1 paralogs in the structurally conserved PP2CA gene subfamily of Brassicaceae

Group A PP2C (PP2CA) genes form a gene subfamily whose members play an important role in regulating many biological processes by dephosphorylation of target proteins. In this study we examined the effects of evolutionary changes responsible for functional divergence of BnaABI1 paralogs in Brassica napus against the background of the conserved PP2CA gene subfamily in Brassicaceae. We performed comprehensive phylogenetic analyses of 192 PP2CA genes in 15 species in combination with protein structure homology modeling. Fundamentally, the number of PP2CA genes remained relatively constant in these taxa, except in the Brassica genus and Camelina sativa. The expansion of this gene subfamily in these species has resulted from whole genome duplication. We demonstrated a high degree of structural conservation of the PP2CA genes, with a few minor variations between the different PP2CA groups. Furthermore, the pattern of conserved sequence motifs in the PP2CA proteins and their secondary and 3D structures revealed strong conservation of the key ion-binding sites. Syntenic analysis of triplicated regions including ABI1 paralogs revealed significant structural rearrangements of the Brassica genomes. The functional and syntenic data clearly show that triplication of BnaABI1 in B. napus has had an impact on its functions, as well as the positions of adjacent genes in the corresponding chromosomal regions. The expression profiling of BnaABI1 genes showed functional divergence, i.e. subfunctionalization, potentially leading to neofunctionalization. These differences in expression are likely due to changes in the promoters of the BnaABI1 paralogs. Our results highlight the complexity of PP2CA gene subfamily evolution in Brassicaceae.

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-β signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

Epigenomic and transcriptomic analysis of chronic inflammatory diseases

Chronic inflammatory diseases (CIDs) have complex pathologies that result from aberrant and persistent immune responses. However, the precise triggers and mechanisms remain elusive. An important aspect of CID research focuses on epigenetics modifications, which regulate gene expression and provide a dynamic transcriptional response to inflammation. In recent years, mounting evidence has demonstrated an association between epigenomic and transcriptomic dysregulation and the phenotypes of CIDs. In particular, epigenetic changes at cis-regulatory elements have provided new insights for immune cell-specific alterations that contribute to disease etiology. Furthermore, the advancements in single-cell genomics provide novel solutions to cell type heterogeneity, which has long posed challenges for CID diagnosis and treatment. In this review, we discuss the current state of epigenomics research of CID and the insights derived from single-cell transcriptomic and epigenomic studies.

Profiling Chromatin Accessibility at Single-cell Resolution

How distinct transcriptional programs are enacted to generate cellular heterogeneity and plasticity, and enable complex fate decisions are important open questions. One key regulator is the cell’s epigenome state that drives distinct transcriptional programs by regulating chromatin accessibility. Genome-wide chromatin accessibility measurements can impart insights into regulatory sequences (in)accessible to DNA-binding proteins at a single-cell resolution. This review outlines molecular methods and bioinformatic tools for capturing cell-to-cell chromatin variation using single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) in a scalable fashion. It also covers joint profiling of chromatin with transcriptome/proteome measurements, computational strategies to integrate multi-omic measurements, and predictive bioinformatic tools to infer chromatin accessibility from single-cell transcriptomic datasets. Methodological refinements that increase power for cell discovery through robust chromatin coverage and integrate measurements from multiple modalities will further expand our understanding of gene regulation during homeostasis and disease.

PeakMatcher facilitates updated Aedes aegypti embryonic cis-regulatory element map

BACKGROUND

The Aedes aegypti mosquito is a threat to human health across the globe. The A. aegypti genome was recently re-sequenced and re-assembled. Due to a combination of long-read PacBio and Hi-C sequencing, the AaegL5 assembly is chromosome complete and significantly improves the assembly in key areas such as the M/m sex-determining locus. Release of the updated genome assembly has precipitated the need to reprocess historical functional genomic data sets, including cis-regulatory element (CRE) maps that had previously been generated for A. aegypti.

RESULTS

We re-processed and re-analyzed the A. aegypti whole embryo FAIRE seq data to create an updated embryonic CRE map for the AaegL5 genome. We validated that the new CRE map recapitulates key features of the original AaegL3 CRE map. Further, we built on the improved assembly in the M/m locus to analyze overlaps of open chromatin regions with genes. To support the validation, we created a new method (PeakMatcher) for matching peaks from the same experimental data set across genome assemblies.

CONCLUSION

Use of PeakMatcher software, which is available publicly under an open-source license, facilitated the release of an updated and validated CRE map, which is available through the NIH GEO. These findings demonstrate that PeakMatcher software will be a useful resource for validation and transferring of previous annotations to updated genome assemblies.

De Novo Assembly of the Northern Cardinal (Cardinalis cardinalis) Genome Reveals Candidate Regulatory Regions for Sexually Dichromatic Red Plumage Coloration

Northern cardinals (Cardinalis cardinalis) are common, mid-sized passerines widely distributed in North America. As an iconic species with strong sexual dichromatism, it has been the focus of extensive ecological and evolutionary research, yet genomic studies investigating the evolution of genotype-phenotype association of plumage coloration and dichromatism are lacking. Here we present a new, highly-contiguous assembly for C. cardinalis We generated a 1.1 Gb assembly comprised of 4,762 scaffolds, with a scaffold N50 of 3.6 Mb, a contig N50 of 114.4 kb and a longest scaffold of 19.7 Mb. We identified 93.5% complete and single-copy orthologs from an Aves dataset using BUSCO, demonstrating high completeness of the genome assembly. We annotated the genomic region comprising the CYP2J19 gene, which plays a pivotal role in the red coloration in birds. Comparative analyses demonstrated non-exonic regions unique to the CYP2J19 gene in passerines and a long insertion upstream of the gene in C. cardinalis Transcription factor binding motifs discovered in the unique insertion region in C. cardinalis suggest potential androgen-regulated mechanisms underlying sexual dichromatism. Pairwise Sequential Markovian Coalescent (PSMC) analysis of the genome reveals fluctuations in historic effective population size between 100,000-250,000 in the last 2 millions years, with declines concordant with the beginning of the Pleistocene epoch and Last Glacial Period. This draft genome of C. cardinalis provides an important resource for future studies of ecological, evolutionary, and functional genomics in cardinals and other birds.

The Tug1 lncRNA locus is essential for male fertility

Background

Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function.

Results

We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function.

Conclusions

Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Transcriptome wide identification and characterization of regulatory genes involved in EAA metabolism and validation through expression analysis in different developmental stages of finger millet spikes

Finger millet is a rich source of seed storage proteins (SSPs). Various regulatory genes play an important role to maintain the quality and accumulation of SSPs in crop seeds. In the present study, nine regulatory genes of EAAs metabolic pathway, i.e., aspartate kinase, homoserine dehydrogenase, threonine synthase, threonine dehydratase, dihydrodipicolinate synthase, cystathionine γ synthase, anthranilate synthase, acetolactate synthase and lysine 2-oxoglutarato reductase/saccharopine dehydrogenase (LOR/SD) were identified from the transcriptomic data of developing spikes of two finger millet genotypes, i.e., GP-45 and GP-1. Results of sequence alignment search and motif/domain analysis showed high similarity of nucleotide sequences of identified regulatory genes with their respective homologs in rice. Results of promoter analysis revealed the presence of various cis-regulatory elements, like nitrogen responsive cis-elements (O2-site and GCN4), light responsive cis-elements, and stress responsive cis-elements. The presence of nine regulatory genes identified from the transcriptomic data of GP-45 and GP-1 was further confirmed by real time expression analysis in high and low protein containing genotypes, i.e., GE-3885 and GE-1437. Results of real time expression analysis showed significantly higher expression (p ≤ 0.01) of regulatory genes in GE-3885 rather than GE-1437 under control and treatment condition. Crude protein content of GE-3885 was found to be significantly higher (p ≤ 0.01) in comparison to GE-1437 under control condition, while under treatment condition GE-1437 was found to be more responsive to KNO₃ treatment rather than GE-3885.

Genome-wide identification and characterization of the fibrillin gene family in Triticum aestivum

Background

The fibrillin (FBN) gene family is highly conserved and widely distributed in the photosynthetic organs of plants. Members of this gene family are involved in the growth and development of plants and their response to biotic and abiotic stresses. Wheat (Triticum aestivum), an important food crop, has a complex genetic background and little progress has occurred in the understanding of its molecular mechanisms.

Methods

In this study, we identified 26 FBN genes in the whole genome of T. aestivum through bioinformatic tools and biotechnological means. These genes were divided into 11 subgroups and were distributed on 11 chromosomes of T. aestivum. Interestingly, most of the TaFBN genes were located on the chromosomes 2A, 2B and 2D. The gene structure of each subgroup of gene family members and the position and number of motifs were highly similar.

Results

The evolutionary analysis results indicated that the affinities of FBNs in monocots were closer together. The tissue-specific analysis revealed that TaFBN genes were expressed in different tissues and developmental stages. In addition, some TaFBNs were involved in one or more biotic and abiotic stress responses. These results provide a basis for further study of the biological function of FBNs.

Genome-wide identification and expression analysis of WRKY transcription factors in pearl millet (Pennisetum glaucum) under dehydration and salinity stress

Background

Plants have developed various sophisticated mechanisms to cope up with climate extremes and different stress conditions, especially by involving specific transcription factors (TFs). The members of the WRKY TF family are well known for their role in plant development, phytohormone signaling and developing resistance against biotic or abiotic stresses. In this study, we performed a genome-wide screening to identify and analyze the WRKY TFs in pearl millet (Pennisetum glaucum; PgWRKY), which is one of the most widely grown cereal crops in the semi-arid regions.

Results

A total number of 97 putative PgWRKY proteins were identified and classified into three major Groups (I-III) based on the presence of WRKY DNA binding domain and zinc-finger motif structures. Members of Group II have been further subdivided into five subgroups (IIa-IIe) based on the phylogenetic analysis. In-silico analysis of PgWRKYs revealed the presence of various cis-regulatory elements in their promoter region like ABRE, DRE, ERE, EIRE, Dof, AUXRR, G-box, etc., suggesting their probable involvement in growth, development and stress responses of pearl millet. Chromosomal mapping evidenced uneven distribution of identified 97 PgWRKY genes across all the seven chromosomes of pearl millet. Synteny analysis of PgWRKYs established their orthologous and paralogous relationship among the WRKY gene family of Arabidopsis thaliana, Oryza sativa and Setaria italica. Gene ontology (GO) annotation functionally categorized these PgWRKYs under cellular components, molecular functions and biological processes. Further, the differential expression pattern of PgWRKYs was noticed in different tissues (leaf, stem, root) and under both drought and salt stress conditions. The expression pattern of PgWRKY33, PgWRKY62 and PgWRKY65 indicates their probable involvement in both dehydration and salinity stress responses in pearl millet.

Conclusion

Functional characterization of identified PgWRKYs can be useful in delineating their role behind the natural stress tolerance of pearl millet against harsh environmental conditions. Further, these PgWRKYs can be employed in genome editing for millet crop improvement.

Enhancer activation by FGF signalling during otic induction

Vertebrate ear progenitors are induced by fibroblast growth factor signalling, however the molecular mechanisms leading to the coordinate activation of downstream targets are yet to be discovered. The ear, like other sensory placodes, arises from the pre-placodal region at the border of the neural plate. Using a multiplex NanoString approach, we determined the response of these progenitors to FGF signalling by examining the changes of more than 200 transcripts that define the otic and other placodes, neural crest and neural plate territories. This analysis identifies new direct and indirect FGF targets during otic induction. Investigating changes in histone marks by ChIP-seq reveals that FGF exposure of pre-placodal cells leads to rapid deposition of active chromatin marks H3K27ac near FGF-response genes, while H3K27ac is depleted in the vicinity of non-otic genes. Genomic regions that gain H3K27ac act as cis-regulatory elements controlling otic gene expression in time and space and define a unique transcription factor signature likely to control their activity. Finally, we show that in response to FGF signalling the transcription factor dimer AP1 recruits the histone acetyl transferase p300 to selected otic enhancers. Thus, during ear induction FGF signalling modifies the chromatin landscape to promote enhancer activation and chromatin accessibility.

Identification and Characterization of Cis-Regulatory Elements for Photoreceptor-Type-Specific Transcription in ZebraFish

Tissue-specific or cell-type-specific transcription of protein-coding genes is controlled by both trans-regulatory elements (TREs) and cis-regulatory elements (CREs). However, it is challenging to identify TREs and CREs, which are unknown for most genes. Here, we describe a protocol for identifying two types of transcription-activating CREs-core promoters and enhancers-of zebrafish photoreceptor type-specific genes. This protocol is composed of three phases: bioinformatic prediction, experimental validation, and characterization of the CREs. To better illustrate the principles and logic of this protocol, we exemplify it with the discovery of the core promoter and enhancer of the mpp5b apical polarity gene (also known as ponli), whose red, green, and blue (RGB) cone-specific transcription requires its enhancer, a member of the rainbow enhancer family. While exemplified with an RGB-cone-specific gene, this protocol is general and can be used to identify the core promoters and enhancers of other protein-coding genes.

Acute salt stress differentially modulates nitrate reductase expression in contrasting salt responsive rice cultivars

Salt stress response includes alteration in the activity of various important enzymes in plants. Nitrate reductase (NR) is one of the known enzyme affected by salt stress. In this study, contrasting salt responsive cultivars (CVS) (IR64-sensitive and CSR 36-tolerant) were considered to study the regulation of NR genes under salt stress conditions. Using Arabidopsis genes Nia1 and Nia2, three different NR genes were identified in rice and their expression study was conducted. Under stress condition, salt-sensitive CVS (IR64) showed a decrease in NR activity under in vitro and in vivo conditions, whereas tolerant CVS showed an increase in NR activity. Different trends for NR activity in contrasting genotype are explained by the variable number of GATA element in the upstream region of the NR gene. This variation of NR activity in contrasting CVS further co-relates with the transcript level of NR genes. The transcript level of three different NR genes also evidenced the effect of CREs in gene regulation. Promoter (1-kb upstream region) of different NR genes contained different abiotic stress-responsive CREs, which explain the differential behavior of these genes towards the abiotic stress. Overall, this study concludes the role of CREs in the regulation of NR gene and indicates the importance of transcriptional control of NR activity under stress condition. This is the first type of report that highlights the role of the regulatory mechanism of NR genes under salt stress condition.

Functional characterization of Gh_A08G1120 (GH3.5) gene reveal their significant role in enhancing drought and salt stress tolerance in cotton

BACKGROUND

Auxins play an important role in plant growth and development; the auxins responsive gene; auxin/indole-3-acetic acid (Aux/IAA), small auxin-up RNAs (SAUR) and Gretchen Hagen3 (GH3) control their mechanisms. The GH3 genes function in homeostasis by the catalytic activities in auxin conjugation and bounding free indole-3-acetic acid (IAA) to amino acids.

RESULTS

In our study, we identified the GH3 genes in three cotton species; Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii, analyzed their chromosomal distribution, phylogenetic relationships, cis-regulatory element function and performed virus induced gene silencing of the novel Gh_A08G1120 (GH3.5) gene. The phylogenetic tree showed four clusters of genes with clade 1, 3 and 4 having mainly members of the GH3 of the cotton species while clade 2 was mainly members belonging to Arabidopsis. There were no paralogous genes, and few orthologous genes were observed between Gossypium and other species. All the GO terms were detected, but only 14 genes were found to have described GO terms in upland cotton, more biological functions were detected, as compared to the other functions. The GH3.17 subfamily harbored the highest number of the cis-regulatory elements, most having promoters towards dehydration-responsiveness. The RNA expression analysis revealed that 10 and 8 genes in drought and salinity stress conditions respectively were upregulated in G. hirsutum. All the genes that were upregulated in plants under salt stress conditions were also upregulated in drought stress; moreover, Gh_A08G1120 (GH3.5) exhibited a significant upregulation across the two stress factors. Functional characterization of Gh_A08G1120 (GH3.5) through virus-induced gene silencing (VIGS) revealed that the VIGS plants ability to tolerate drought and salt stresses was significantly reduced compared to the wild types. The chlorophyll content, relative leaf water content (RLWC), and superoxide dismutase (SOD) concentration level were reduced significantly while malondialdehyde concentration and ion leakage as a measure of cell membrane stability (CMS) increased in VIGS plants under drought and salt stress conditions.

CONCLUSION

This study revealed the significance of the GH3 genes in enabling the plant's adaptation to drought and salt stress conditions as evidenced by the VIGS results and RT-qPCR analysis.

Genome-wide enhancer annotations differ significantly in genomic distribution, evolution, and function

Background

Non-coding gene regulatory enhancers are essential to transcription in mammalian cells. As a result, a large variety of experimental and computational strategies have been developed to identify cis-regulatory enhancer sequences. Given the differences in the biological signals assayed, some variation in the enhancers identified by different methods is expected; however, the concordance of enhancers identified by different methods has not been comprehensively evaluated. This is critically needed, since in practice, most studies consider enhancers identified by only a single method. Here, we compare enhancer sets from eleven representative strategies in four biological contexts.

Results

All sets we evaluated overlap significantly more than expected by chance; however, there is significant dissimilarity in their genomic, evolutionary, and functional characteristics, both at the element and base-pair level, within each context. The disagreement is sufficient to influence interpretation of candidate SNPs from GWAS studies, and to lead to disparate conclusions about enhancer and disease mechanisms. Most regions identified as enhancers are supported by only one method, and we find limited evidence that regions identified by multiple methods are better candidates than those identified by a single method. As a result, we cannot recommend the use of any single enhancer identification strategy in all settings.

Conclusions

Our results highlight the inherent complexity of enhancer biology and identify an important challenge to mapping the genetic architecture of complex disease. Greater appreciation of how the diverse enhancer identification strategies in use today relate to the dynamic activity of gene regulatory regions is needed to enable robust and reproducible results.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5779-x) contains supplementary material, which is available to authorized users.

Allele-specific expression and alternative splicing in horse×donkey and cattle×yak hybrids

Divergence of gene expression and alternative splicing is a crucial driving force in the evolution of species; to date, however the molecular mechanism remains unclear. Hybrids of closely related species provide a suitable model to analyze allele-specific expression (ASE) and allele-specific alternative splicing (ASS). Analysis of ASE and ASS can uncover the differences in cis-regulatory elements between closely related species, while eliminating interference of trans-regulatory elements. Here, we provide a detailed characterization of ASE and ASS from 19 and 10 transcriptome datasets across five tissues from reciprocal-cross hybrids of horse×donkey (mule/hinny) and cattle×yak (dzo), respectively. Results showed that 4.8%-8.7% and 10.8%-16.7% of genes exhibited ASE and ASS, respectively. Notably, lncRNAs and pseudogenes were more likely to show ASE than protein-coding genes. In addition, genes showing ASE and ASS in mule/hinny were found to be involved in the regulation of muscle strength, whereas those of dzo were involved in high-altitude adaptation. In conclusion, our study demonstrated that exploration of genes showing ASE and ASS in hybrids of closely related species is feasible for species evolution research.

Sucrose transporters of resistant grapevine are involved in stress resistance

The whole promoter regions of SUTs in Vitis were firstly isolated. SUTs are involved in the adaptation to biotic and abiotic stresses. The vulnerability of Vitis vinifera to abiotic and biotic stresses limits its yields. In contrast, Vitis amurensis displays resistance to environmental stresses, such as microbial pathogens, low temperatures, and drought. Sucrose transporters (SUTs) are important regulators for plant growth and stress tolerance; however, the role that SUTs play in stress resistance in V. amurensis is not known. Using V. amurensis Ruper. 'Zuoshan-1' and V. vinifera 'Chardonnay', we found that SUC27 was highly expressed in several vegetative organs of Zuoshan-1, SUC12 was weakly expressed or absent in most organs in both the species, and the distribution of SUC11 in source and sink organs was highest in Zuoshan-1. A search for cis-regulatory elements in the promoter sequences of SUTs revealed that they were regulated by light, environmental stresses, physiological correlation, and hormones. The SUTs in Zuoshan-1 mostly show a higher and rapid response than in Chardonnay under the induction by Plasmopara viticola infection, cold, water deficit, and dark conditions. The induction of SUTs was associated with the upregulation of key genes involved in sucrose metabolism and the biosynthesis of plant hormones. These results indicate that stress resistance in Zuoshan-1 is governed by the differential distribution and induction of SUTs by various stimuli, and the subsequent promotion of sucrose metabolism and hormone synthesis.

Quantitative analysis of the ThrbCRM1-centered gene regulatory network

Enhancer activity is determined by both the activity and occupancy of transcription factors as well as the specific sequences they bind. Experimental investigation of this dynamic requires the ability to manipulate components of the system, ideally in as close to an in vivo context as possible. Here we use electroporation of plasmid reporters to define critical parameters of a specific cis-regulatory element, ThrbCRM1, during retinal development. ThrbCRM1 is associated with cone photoreceptor genesis and activated in a subset of developing retinal cells that co-express the Otx2 and Onecut1 (OC1) transcription factors. Variation of reporter plasmid concentration was used to generate dose response curves and revealed an effect of binding site availability on the number and strength of cells with reporter activity. Critical sequence elements of the ThrbCRM1 element were defined using both mutagenesis and misexpression of the Otx2 and OC1 transcription factors in the developing retina. Additionally, these experiments suggest that the ThrbCRM1 element is co-regulated by Otx2 and OC1 even under conditions of sub-optimal binding of OC1.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

In vivo epigenome editing and transcriptional modulation using CRISPR technology

The rapid advancement of CRISPR technology has enabled targeted epigenome editing and transcriptional modulation in the native chromatin context. However, only a few studies have reported the successful editing of the epigenome in adult animals in contrast to the rapidly growing number of in vivo genome editing over the past few years. In this review, we discuss the challenges facing in vivo epigenome editing and new strategies to overcome the huddles. The biggest challenge has been the difficulty in packaging dCas9 fusion proteins required for manipulation of epigenome into the adeno-associated virus (AAV) delivery vehicle. We review the strategies to address the AAV packaging issue, including small dCas9 orthologues, truncated dCas9 mutants, a split-dCas9 system, and potent truncated effector domains. We discuss the dCas9 conjugation strategies to recruit endogenous chromatin modifiers and remodelers to specific genomic loci, and recently developed methods to recruit multiple copies of the dCas9 fusion protein, or to simultaneous express multiple gRNAs for robust epigenome editing or synergistic transcriptional modulation. The use of Cre-inducible dCas9-expressing mice or a genetic cross between dCas9- and sgRNA-expressing flies has also helped overcome the transgene delivery issue. We provide perspective on how a combination use of these strategies can facilitate in vivo epigenome editing and transcriptional modulation.

SMuRF: a novel tool to identify regulatory elements enriched for somatic point mutations

BACKGROUND

Single Nucleotide Variants (SNVs), including somatic point mutations and Single Nucleotide Polymorphisms (SNPs), in noncoding cis-regulatory elements (CREs) can affect gene regulation and lead to disease development. Several approaches have been developed to identify highly mutated regions, but these do not take into account the specific genomic context, and thus likelihood of mutation, of CREs.

RESULTS

Here, we present SMuRF (Significantly Mutated Region Finder), a user-friendly command-line tool to identify these significantly mutated regions from user-defined genomic intervals and SNVs. We demonstrate this using publicly available datasets in which SMuRF identifies 72 significantly mutated CREs in liver cancer, including known mutated gene promoters as well as previously unreported regions.

CONCLUSIONS

SMuRF is a helpful tool to allow the simple identification of significantly mutated regulatory elements. It is open-source and freely available on GitHub ( https://github.com/LupienLab/SMURF ).

Identification and characterization of early photoreceptor cis-regulatory elements and their relation to Onecut1

Background

Cone and rod photoreceptors are two of the primary cell types affected in human retinal disease. Potential strategies to combat these diseases are the use of gene therapy to rescue compromised photoreceptors or to generate new functional photoreceptors to replace those lost in the diseased retina. Cis-regulatory elements specific to cones, rods, or both types of photoreceptors are critical components of successful implementation of these two strategies. The purpose of this study was to identify and characterize the cell type specificity and activity of cis-regulatory elements active in developing photoreceptors.

Methods

Cis-regulatory elements were introduced into the developing chicken and mouse retina by electroporation. Characterization of reporter activity in relation with cell type markers was determined using confocal microscopy. In addition, two high-throughput flow cytometry assay were developed to assess whether these elements were downstream of Onecut1 in the photoreceptor specification network.

Results

The majority of cis-regulatory elements were active in both cone and rod photoreceptors and were largely uninfluenced by a Onecut1 dominant-negative construct. Elements associated with the Thrb, Nr2e3, and Rhodopsin genes showed highly enriched activity in cones or rods, and were affected by interference in Onecut1 signaling. Rhodopsin promoter activity was the most highly influenced by Onecut1 activity and its induction could be modulated by the Maf family transcription factor L-Maf. Nr2e3 elements were observed to have activity in cone photoreceptors and Nr2e3 protein was expressed in developing cone photoreceptors, suggesting a role for this predominant rod gene in cone photoreceptor development.

Conclusions

The analysis presented here provides an experimental framework to determine the specificity and strength of photoreceptor elements within specific genetic networks during development. The Onecut1 transcription factor is one such factor that influences the gene regulatory networks specific to cones and rods, but not those that are common to both.

Electronic supplementary material

The online version of this article (10.1186/s13064-018-0121-x) contains supplementary material, which is available to authorized users.

Biological functions, regulatory mechanisms, and disease relevance of RNA localization pathways

The asymmetric subcellular distribution of RNA molecules from their sites of transcription to specific compartments of the cell is an important aspect of post-transcriptional gene regulation. This involves the interplay of intrinsic cis-regulatory elements within the RNA molecules with trans-acting RNA-binding proteins and associated factors. Together, these interactions dictate the intracellular localization route of RNAs, whose downstream impacts have wide-ranging implications in cellular physiology. In this review, we examine the mechanisms underlying RNA localization and discuss their biological significance. We also review the growing body of evidence pointing to aberrant RNA localization pathways in the development and progression of diseases.

Isolation, characterization, and evaluation of three Citrus sinensis-derived constitutive gene promoters

Regulatory sequences from the citrus constitutive genes cyclophilin (CsCYP), glyceraldehyde-3-phosphate dehydrogenase C2 (CsGAPC2), and elongation factor 1-alpha (CsEF1) were isolated, fused to the uidA gene, and qualitatively and quantitatively evaluated in transgenic sweet orange plants. The 5' upstream region of a gene (the promoter) is the most important component for the initiation and regulation of gene transcription of both native genes and transgenes in plants. The isolation and characterization of gene regulatory sequences are essential to the development of intragenic or cisgenic genetic manipulation strategies, which imply the use of genetic material from the same species or from closely related species. We describe herein the isolation and evaluation of the promoter sequence from three constitutively expressed citrus genes: cyclophilin (CsCYP), glyceraldehyde-3-phosphate dehydrogenase C2 (CsGAPC2), and elongation factor 1-alpha (CsEF1). The functionality of the promoters was confirmed by a histochemical GUS assay in leaves, stems, and roots of stably transformed citrus plants expressing the promoter-uidA construct. Lower uidA mRNA levels were detected when the transgene was under the control of citrus promoters as compared to the expression under the control of the CaMV35S promoter. The association of the uidA gene with the citrus-derived promoters resulted in mRNA levels of up to 60-41.8% of the value obtained with the construct containing CaMV35S driving the uidA gene. Moreover, a lower inter-individual variability in transgene expression was observed amongst the different transgenic lines, where gene constructs containing citrus-derived promoters were used. In silico analysis of the citrus-derived promoter sequences revealed that their activity may be controlled by several putative cis-regulatory elements. These citrus promoters will expand the availability of regulatory sequences for driving gene expression in citrus gene-modification programs.

Identification and molecular characterization of Dof transcription factor gene family preferentially expressed in developing spikes of Eleusine coracana L

We report 48 putative DNA binding with one finger (Dof) TF genes from genome and transcriptome data of finger millet (Eleusine coracana L.; FM), involved in plant developmental process. To characterize seed-specific Dof genes, transcript profiles of 32 EcDof identified from transcriptome data of developing spikes of FM genotypes were further analyzed in different tissues (root, stem, and leaf) and developmental stages of spikes (S₁, S₂, S_3, and S₄) in two FM genotypes [GE1437 (low protein genotype; LPG) and GE3885 (high protein genotype; HPG)]. More than 50% of identified EcDof genes showed expression during seed development processes. Among these, seven genes (EcDof 3, EcDof 5, EcDof 15, EcDof 18, EcDof 22, EcDof 23, and EcDof 31) expressed maximally at specific stages of seed development. Fourteen EcDof genes showed that differential transcript accumulation in vegetative tissue as well as in developing spikes suggests involvement during seed filling and also throughout the plant development. In addition, three EcDof genes (EcDof 9, EcDof 25, and EcDof 28) expressed preferentially at root and stem tissue. The 3D structural prediction of EcDof proteins showed variability in structural attributes. Molecular docking results showed strong binding affinity for seed-specific EcDof-EcO2 with α-prolamine promoters. The identified and characterized EcDof genes will help to dissect the roles of FM seed-specific Dof genes.

Systematic and synthetic approaches to rewire regulatory networks

Microbial gene regulatory networks are composed of cis- and trans-components that in concert act to control essential and adaptive cellular functions. Regulatory components and interactions evolve to adopt new configurations through mutations and network rewiring events, resulting in novel phenotypes that may benefit the cell. Advances in high-throughput DNA synthesis and sequencing have enabled the development of new tools and approaches to better characterize and perturb various elements of regulatory networks. Here, we highlight key recent approaches to systematically dissect the sequence space of cis-regulatory elements and trans-regulators as well as their inter-connections. These efforts yield fundamental insights into the architecture, robustness, and dynamics of gene regulation and provide models and design principles for building synthetic regulatory networks for a variety of practical applications.