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

Retracing the evolution of Pneumocystis species, with a focus on the human pathogen Pneumocystis jirovecii

SUMMARYEvery human being is presumed to be infected by the fungus Pneumocystis jirovecii at least once in his or her lifetime. This fungus belongs to a large group of species that appear to exclusively infect mammals, with P. jirovecii being the only one known to cause disease in humans. The mystery of P. jirovecii origin and speciation is just beginning to unravel. Here, we provide a review of the major steps of P. jirovecii evolution. The Pneumocystis genus likely originated from soil or plant-associated organisms during the period of Cretaceous ~165 million years ago and successfully shifted to mammals. The transition coincided with a substantial loss of genes, many of which are related to the synthesis of nutrients that can be scavenged from hosts or cell wall components that could be targeted by the mammalian immune system. Following the transition, the Pneumocystis genus cospeciated with mammals. Each species specialized at infecting its own host. Host specialization is presumably built at least partially upon surface glycoproteins, whose protogene was acquired prior to the genus formation. P. jirovecii appeared at ~65 million years ago, overlapping with the emergence of the first primates. P. jirovecii and its sister species P. macacae, which infects macaques nowadays, may have had overlapping host ranges in the distant past. Clues from molecular clocks suggest that P. jirovecii did not cospeciate with humans. Molecular evidence suggests that Pneumocystis speciation involved chromosomal rearrangements and the mounting of genetic barriers that inhibit gene flow among species.

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.

Novel islands of GGC and GCC repeats coincide with human evolution

BACKGROUND

Because of high mutation rate, overrepresentation in genic regions, and link with various neurological, neurodegenerative, and movement disorders, GGC and GCC short tandem repeats (STRs) are prone to natural selection. Among a number of lacking data, the 3-repeats of these STRs remain widely unexplored.

RESULTS

In a genome-wide search in human, here we mapped GGC and GCC STRs of ≥3-repeats, and found novel islands of up to 45 of those STRs, populating spans of 1 to 2 kb of genomic DNA. RGPD4 and NOC4L harbored the densest (GGC)3 (probability 3.09061E-71) and (GCC)3 (probability 1.72376E-61) islands, respectively, and were human-specific. We also found prime instances of directional incremented density of STRs at specific loci in human versus other species, including the FOXK2 and SKI GGC islands. The genes containing those islands significantly diverged in expression in human versus other species, and the proteins encoded by those genes interact closely in a physical interaction network, consequence of which may be human-specific characteristics such as higher order brain functions.

CONCLUSION

We report novel islands of GGC and GCC STRs of evolutionary relevance to human. The density, and in some instances, periodicity of these islands support them as a novel genomic entity, which need to be further explored in evolutionary, mechanistic, and functional platforms.

Characterization of PYL gene family and identification of HaPYL genes response to drought and salt stress in sunflower

In the context of global climate change, drought and soil salinity are some of the most devastating abiotic stresses affecting agriculture today. PYL proteins are essential components of abscisic acid (ABA) signaling and play critical roles in responding to abiotic stressors, including drought and salt stress. Although PYL genes have been studied in many species, their roles in responding to abiotic stress are still unclear in the sunflower. In this study, 19 HaPYL genes, distributed on 15 of 17 chromosomes, were identified in the sunflower. Fragment duplication is the main cause of the expansion of PYL genes in the sunflower genome. Based on phylogenetic analysis, HaPYL genes were divided into three subfamilies. Members in the same subfamily share similar protein motifs and gene exon-intron structures, except for the second subfamily. Tissue expression patterns suggested that HaPYLs serve different functions when responding to developmental and environmental signals in the sunflower. Exogenous ABA treatment showed that most HaPYLs respond to an increase in the ABA level. Among these HaPYLs, HaPYL2a, HaPYL4d, HaPYL4g, HaPYL8a, HaPYL8b, HaPYL8c, HaPYL9b, and HaPYL9c were up-regulated with PEG6000 treatment and NaCl treatment. This indicates that they may play a role in resisting drought and salt stress in the sunflower by mediating ABA signaling. Our findings provide some clues to further explore the functions of PYL genes in the sunflower, especially with regards to drought and salt stress resistance.

Dyads of GGC and GCC form hotspot colonies that coincide with the evolution of human and other great apes

BACKGROUND

GGC and GCC short tandem repeats (STRs) are of various evolutionary, biological, and pathological implications. However, the fundamental two-repeats (dyads) of these STRs are widely unexplored.

RESULTS

On a genome-wide scale, we mapped (GGC)2 and (GCC)2 dyads in human, and found monumental colonies (distance between each dyad < 500 bp) of extraordinary density, and in some instances periodicity. The largest (GCC)2 and (GGC)2 colonies were intergenic, homogeneous, and human-specific, consisting of 219 (GCC)2 on chromosome 2 (probability < 1.545E-219) and 70 (GGC)2 on chromosome 9 (probability = 1.809E-148). We also found that several colonies were shared in other great apes, and directionally increased in density and complexity in human, such as a colony of 99 (GCC)2 on chromosome 20, that specifically expanded in great apes, and reached maximum complexity in human (probability 1.545E-220). Numerous other colonies of evolutionary relevance in human were detected in other largely overlooked regions of the genome, such as chromosome Y and pseudogenes. Several of the genes containing or nearest to those colonies were divergently expressed in human.

CONCLUSION

In conclusion, (GCC)2 and (GGC)2 form unprecedented genomic colonies that coincide with the evolution of human and other great apes. The extent of the genomic rearrangements leading to those colonies support overlooked recombination hotspots, shared across great apes. The identified colonies deserve to be studied in mechanistic, evolutionary, and functional platforms.

The presence of an intron relieves gene repression caused by promoter-proximal four-bp specific sequences in yeast

Introns can enhance gene expression in eukaryotic cells in a process called intron-mediated enhancement (IME). The levels of enhancement are affected not only by the intron sequence but also by coding sequences (CDSs). However, the parts of CDSs responsible for mediating IME have not yet been identified. In this study, we identified an IME-mediating sequence by analyzing three pairs of IME-sensitive and -insensitive CDSs in Saccharomyces cerevisiae. Expression of the CDSs yCLuc, yRoGLU1, and KmBGA1 was enhanced by the presence of an intron (i.e., they were IME sensitive), but the expression of each corresponding codon-changed CDS, which encoded the identical amino acid sequence, was not enhanced (i.e., they were IME insensitive). Interestingly, the IME-insensitive CDSs showed higher expression levels that were like intron-enhanced expression of IME-sensitive CDSs, suggesting that expression of IME-sensitive CDSs was repressed. A four-nucleotide sequence (TCTT) located in the promoter-proximal position of either the untranslated or coding region was found to be responsible for repression in IME-sensitive CDSs, and repression caused by the TCTT sequence was relieved by the presence of an intron. Further, it was found that the expression of intron-containing yeast-native genes, UBC4 and MPT5, was repressed by TCTT in the CDS but relieved by the introns. These results indicate that TCTT sequences in promoter-proximal positions repress gene expression and that introns play a role in relieving gene repression caused by sequences such as TCTT.

Intron-mediated enhancement of DIACYLGLYCEROL ACYLTRANSFERASE1 expression in energycane promotes a step change for lipid accumulation in vegetative tissues

BACKGROUND

Metabolic engineering for hyperaccumulation of lipids in vegetative tissues is a novel strategy for enhancing energy density and biofuel production from biomass crops. Energycane is a prime feedstock for this approach due to its high biomass production and resilience under marginal conditions. DIACYLGLYCEROL ACYLTRANSFERASE (DGAT) catalyzes the last and only committed step in the biosynthesis of triacylglycerol (TAG) and can be a rate-limiting enzyme for the production of TAG.

RESULTS

In this study, we explored the effect of intron-mediated enhancement (IME) on the expression of DGAT1 and resulting accumulation of TAG and total fatty acid (TFA) in leaf and stem tissues of energycane. To maximize lipid accumulation these evaluations were carried out by co-expressing the lipogenic transcription factor WRINKLED1 (WRI1) and the TAG protect factor oleosin (OLE1). Including an intron in the codon-optimized TmDGAT1 elevated the accumulation of its transcript in leaves by seven times on average based on 5 transgenic lines for each construct. Plants with WRI1 (W), DGAT1 with intron (Di), and OLE1 (O) expression (WDiO) accumulated TAG up to a 3.85% of leaf dry weight (DW), a 192-fold increase compared to non-modified energycane (WT) and a 3.8-fold increase compared to the highest accumulation under the intron-less gene combination (WDO). This corresponded to TFA accumulation of up to 8.4% of leaf dry weight, a 2.8-fold or 6.1-fold increase compared to WDO or WT, respectively. Co-expression of WDiO resulted in stem accumulations of TAG up to 1.14% of DW or TFA up to 2.08% of DW that exceeded WT by 57-fold or 12-fold and WDO more than twofold, respectively. Constitutive expression of these lipogenic "push pull and protect" factors correlated with biomass reduction.

CONCLUSIONS

Intron-mediated enhancement (IME) of the expression of DGAT resulted in a step change in lipid accumulation of energycane and confirmed that under our experimental conditions it is rate limiting for lipid accumulation. IME should be applied to other lipogenic factors and metabolic engineering strategies. The findings from this study may be valuable in developing a high biomass feedstock for commercial production of lipids and advanced biofuels.

Genome-wide identification and characterization of ABA receptor pyrabactin resistance 1-like protein (PYL) family in oat

Abscisic acid (ABA) is a phytohormone that plays an important role in plant growth and development. Meanwhile, ABA also plays a key role in the plant response to abiotic stressors such as drought and high salinity. The pyrabactin resistance 1-like (PYR/PYL) protein family of ABA receptors is involved in the initial step of ABA signal transduction. However, no systematic studies of the PYL family in "Avena sativa, a genus Avena in the grass family Poaceae," have been conducted to date. Thus, in this study, we performed a genome-wide screening to identify PYL genes in oat and characterized their responses to drought stress. A total of 12 AsPYL genes distributed on nine chromosomes were identified. The phylogenetic analysis divided these AsPYLs into three subfamilies, based on structural and functional similarities. Gene and motif structure analysis of AsPYLs revealed that members of each subfamily share similar gene and motif structure. Segmental duplication appears to be the driving force for the expansion of PYLs, Furthermore, stress-responsive AsPYLs were detected through RNA-seq analysis. The qRT-PCR analysis of 10 AsPYL genes under drought, salt, and ABA stress revealed that AsPYL genes play an important role in stress response. These data provide a reference for further studies on the oat PYL gene family and its function.

The development of DNA vaccines against SARS-CoV-2

BACKGROUND

The COVID-19 pandemic exerted significant impacts on public health and global economy. Research efforts to develop vaccines at warp speed against SARS-CoV-2 led to novel mRNA, viral vectored, and inactivated vaccines being administered. The current COVID-19 vaccines incorporate the full S protein of the SARS-CoV-2 Wuhan strain but rapidly emerging variants of concern (VOCs) have led to significant reductions in protective efficacies. There is an urgent need to develop next-generation vaccines which could effectively prevent COVID-19.

METHODS

PubMed and Google Scholar were systematically reviewed for peer-reviewed papers up to January 2023.

RESULTS

A promising solution to the problem of emerging variants is a DNA vaccine platform since it can be easily modified. Besides expressing whole protein antigens, DNA vaccines can also be constructed to include specific nucleotide genes encoding highly conserved and immunogenic epitopes from the S protein as well as from other structural/non-structural proteins to develop effective vaccines against VOCs. DNA vaccines are associated with low transfection efficiencies which could be enhanced by chemical, genetic, and molecular adjuvants as well as delivery systems.

CONCLUSIONS

The DNA vaccine platform offers a promising solution to the design of effective vaccines. The challenge of limited immunogenicity in humans might be solved through the use of genetic modifications such as the addition of nuclear localization signal (NLS) peptide gene, strong promoters, MARs, introns, TLR agonists, CD40L, and the development of appropriate delivery systems utilizing nanoparticles to increase uptake by APCs in enhancing the induction of potent immune responses.

Cell Fate Programming by Transcription Factors and Epigenetic Machinery in Stomatal Development

The development of multi-cellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic-helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling - the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knock-down of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

Plant Promoters and Terminators for High-Precision Bioengineering

High-precision bioengineering and synthetic biology require fine-tuning gene expression at both transcriptional and posttranscriptional levels. Gene transcription is tightly regulated by promoters and terminators. Promoters determine the timing, tissues and cells, and levels of the expression of genes. Terminators mediate transcription termination of genes and affect mRNA levels posttranscriptionally, e.g., the 3'-end processing, stability, translation efficiency, and nuclear to cytoplasmic export of mRNAs. The promoter and terminator combination affects gene expression. In the present article, we review the function and features of plant core promoters, proximal and distal promoters, and terminators, and their effects on and benchmarking strategies for regulating gene expression.

Internode elongation in energy cane shows remarkable clues on lignocellulosic biomass biosynthesis in Saccharum hybrids

Energy cane is a dedicated crop to high biomass production and selected during Saccharum breeding programs to fit specific industrial needs for 2G bioethanol production. Internode elongation is one of the most important characteristics in Saccharum hybrids due to its relationship with crop yield. In this study, we selected the third internode elongation of the energy cane. To characterize this process, we divided the internode into five sections and performed a detailed transcriptome analysis (RNA-Seq) and cell wall characterization. The histological analyses revealed a remarkable gradient that spans from cell division and protoxylem lignification to the internode maturation and complete vascular bundle lignification. RNA-Seq analysis revealed more than 11,000 differentially expressed genes between the sections internal. Gene ontology analyzes showed enriched categories in each section, as well as the most expressed genes in each section, presented different biological processes. We found that the internode elongation and division zones have a large number of unique genes. Evaluated the specific profile of genes related to primary and secondary cell wall formation, cellulose synthesis, hemicellulose, lignin, and growth-related genes. For each section these genes presented different profiles along the internode in elongation in energy cane. The results of this study provide an overview of the regulation of gene expression of an internode elongation in energy cane. Gene expression analysis revealed promising candidates for transcriptional regulation of energy cane lignification and evidence key genes for the regulation of internode development, which can serve as a basis for understanding the molecular regulatory mechanisms that support the growth and development of plants in the Saccahrum complex.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

The Expanding Role of Alternative Splicing in Vascular Smooth Muscle Cell Plasticity

Vascular smooth muscle cells (VSMCs) display extraordinary phenotypic plasticity. This allows them to differentiate or dedifferentiate, depending on environmental cues. The ability to ‘switch’ between a quiescent contractile phenotype to a highly proliferative synthetic state renders VSMCs as primary mediators of vascular repair and remodelling. When their plasticity is pathological, it can lead to cardiovascular diseases such as atherosclerosis and restenosis. Coinciding with significant technological and conceptual innovations in RNA biology, there has been a growing focus on the role of alternative splicing in VSMC gene expression regulation. Herein, we review how alternative splicing and its regulatory factors are involved in generating protein diversity and altering gene expression levels in VSMC plasticity. Moreover, we explore how recent advancements in the development of splicing-modulating therapies may be applied to VSMC-related pathologies.

CRY1 Gene Polymorphism and Racing Performance of Homing Pigeons

Cryptochromes (CRY) are the family of proteins proposed as the putative magnetoreceptor molecules. In birds, among others in pigeons, CRY1 is widely expressed in a retina. Homing pigeons are known for their navigational abilities, and pigeon racing is a popular sport. So, the purpose of this study was to analyze the variability of the nucleotide sequence of the homing pigeon CRY1 gene, spanning the region coding the two amino acids W320 and W374 of Trp-triad, and estimate the relationship between genotypes and the racing performance. Investigations were carried out on 129 pigeons. Analysis of sequencing results indicated the AG to TT change within the seventh intron of CRY1 gene. Genotypes were determined by the forced PCR-RFLP method. The influence of detected polymorphism on the results of racing pigeons in 100-400 km flights was shown. The AG/TT individuals achieved significantly higher (p ≤ 0.05) mean values of ace points (AP) than the AG/AG ones. Regarding the detected nucleotide change localization, the polymorphism may be involved in CRY1 gene expression modulation. The AG to TT change in CRY1 gene may be considered as a potential genetic marker of racing performance in homing pigeons.

Increased expression of ANAC017 primes for accelerated senescence

Recent studies in Arabidopsis (Arabidopsis thaliana) have reported conflicting roles for NAC DOMAIN CONTAINING PROTEIN 17 (ANAC017), a transcription factor regulating mitochondria-to-nuclear signaling, and its closest paralog NAC DOMAIN CONTAINING PROTEIN 16 (ANAC016), in leaf senescence. By synchronizing senescence in individually darkened leaves of knockout and overexpressing mutants from these contrasting studies, we demonstrate that elevated ANAC017 expression consistently causes accelerated senescence and cell death. A time-resolved transcriptome analysis revealed that senescence-associated pathways such as autophagy are not constitutively activated in ANAC017 overexpression lines, but require a senescence-stimulus to trigger accelerated induction. ANAC017 transcript and ANAC017-target genes are constitutively upregulated in ANAC017 overexpression lines, but surprisingly show a transient "super-induction" 1 d after senescence induction. This induction of ANAC017 and its target genes is observed during the later stages of age-related and dark-induced senescence, indicating the ANAC017 pathway is also activated in natural senescence. In contrast, knockout mutants of ANAC017 showed lowered senescence-induced induction of ANAC017 target genes during the late stages of dark-induced senescence. Finally, promoter binding analyses show that the ANAC016 promoter sequence is directly bound by ANAC017, so ANAC016 likely acts downstream of ANAC017 and is directly transcriptionally controlled by ANAC017 in a feed-forward loop during late senescence.

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.

Efficient expression and function of a receptor-like kinase in wheat powdery mildew defence require an intron-located MYB binding site

The LRK10-like receptor kinases (LRK10L-RLKs) are ubiquitously present in higher plants, but knowledge of their expression and function is still limited. Here, we report expression and functional analysis of TtdLRK10L-1, a typical LRK10L-RLK in durum wheat (Triticum turgidum L. ssp. durum). The introns of TtdLRK10L-1 contained multiple kinds of predicted cis-elements. To investigate the potential effect of these cis-elements on TtdLRK10L-1 expression and function, two types of transgenic wheat lines were prepared, which expressed a GFP-tagged TtdLRK10L-1 protein (TtdLRK10L-1:GFP) from the cDNA or genomic DNA (gDNA) sequence of TtdLRK10L-1 under the native promoter. TtdLRK10L-1:GFP expression was up-regulated by the powdery mildew pathogen Blumeria graminis f. sp. tritici (Bgt) in both types of transgenic plants, with the scale of the elevation being much stronger in the gDNA lines. Both types of transgenic plants exhibited enhanced resistance to Bgt infection relative to wild type control. Notably, the Bgt defence activated in the gDNA lines was significantly stronger than that in the cDNA lines. Further analysis revealed that a putative MYB transcription factor binding site (MYB-BS, CAGTTA) located in TtdLRK10L-1 intron I was critical for the efficient expression and function of TtdLRK10L-1 in Bgt defence. This MYB-BS could also increase the activity of a superpromoter widely used in ectopic gene expression studies in plants. Together, our results deepen the understanding of the expression and functional characteristics of LRK10L-RLKs. TtdLRK10L-1 is likely useful for further dissecting the molecular processes underlying wheat defence against Bgt and for developing Bgt resistant wheat crops.

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.

Genome-wide identification and characterization of ABA receptor PYL gene family in rice

BACKGROUND

Abscisic acid (ABA), a key phytohormone that controls plant growth and stress responses, is sensed by the pyrabactin resistance 1(PYR1)/PYR1-like (PYL)/regulatory components of the ABA receptor (RCAR) family of proteins. Comprehensive information on evolution and function of PYL gene family in rice (Oryza sativa) needs further investigation. This study made detailed analysis on evolutionary relationship between PYL family members, collinearity, synteny, gene structure, protein motifs, cis-regulatory elements (CREs), SNP variations, miRNAs targeting PYLs and expression profiles in different tissues and stress responses.

RESULTS

Based on sequence homology with Arabidopsis PYL proteins, we identified a total of 13 PYLs in rice (BOP clade) and maize (PACCMAD clade), while other members of BOP (wheat - each diploid genome, barley and Brachypodium) and PACCMAD (sorghum and foxtail millet) have 8-9 PYLs. The phylogenetic analysis divided PYLs into three subfamilies that are structurally and functionally conserved across species. Gene structure and motif analysis of OsPYLs revealed that members of each subfamily have similar gene and motif structure. Segmental duplication appears be the driving force for the expansion of PYLs, and the majority of the PYLs underwent evolution under purifying selection in rice. 32 unique potential miRNAs that might target PYLs were identified in rice. Thus, the predicted regulation of PYLs through miRNAs in rice is more elaborate as compared with B. napus. Further, the miRNAs identified to in this study were also regulated by stresses, which adds additional layer of regulation of PYLs. The frequency of SAPs identified was higher in indica cultivars and were predominantly located in START domain that participate in ABA binding. The promoters of most of the OsPYLs have cis-regulatory elements involved in imparting abiotic stress responsive expression. In silico and q-RT-PCR expression analyses of PYL genes revealed multifaceted role of ABARs in shaping plant development as well as abiotic stress responses.

CONCLUSION

The predicted miRNA mediated regulation of OsPYLs and stress regulated expression of all OsPYLs, at least, under one stress, lays foundation for further validation and fine tuning ABA receptors for stress tolerance without yield penalty in rice.

Comparative sequence analysis across Brassicaceae, regulatory diversity in KCS5 and KCS6 homologs from Arabidopsis thaliana and Brassica juncea, and intronic fragment as a negative transcriptional regulator

Intra- and epicuticular-waxes primarily comprising of very long chain aliphatic lipid (VLCFA), terpenoids and secondary metabolites such as sterol and flavonoids played a major role in successful colonization of terrestrial ecosystem by aquatic plants and are thus considered as a key evolutionary innovation. The key rate limiting step of Fatty Acid (FA) biosynthesis of condensation/elongation are catalyzed by the enzyme, β-ketoacyl coenzyme A synthase (KCS), part of FAE (Fatty Acid Elongase) complex. KCS6 has been shown to be responsible for elongation using C22 fatty acid as substrate and is considered essential for synthesis of VLCFA for cuticular waxes. Earlier studies have established KCS5 as a close paralog of KCS6 in Arabidopsis thaliana, albeit with non-redundant function. We subsequently established segmental duplication responsible for origin of KCS6-KCS5 paralogy which is exclusive to Brassicaceae. In the present study, we aim to understand impact of duplication on regulatory diversification and evolution, through sequence and functional analysis of cis-regulatory element of KCS5 and KCS6. High level of sequence variation leading to conservation of only the proximal end of the promoter corresponding to the core promoter was observed among Brassicaceae members; such high diversity was also revealed when sliding window analysis revealed only two to three phylogenetic footprints. Profiling of transcription factor binding sites (TFBS) across Brassicaceae shows presence of light, hormone and stress responsive motifs; a few motifs involved in tissue specific expression (Skn-1; endosperm) were also detected. Functional characterization using transcriptional fusion constructs revealed regulatory diversification when promoter activity of homologs from A. thaliana and Brassica juncea were compared. When subjected to 5-Azacytidine, altered promoter activity was observed, implying role of DNA methylation in transcriptional regulation. Finally, investigation of the role of an 87 bp fragment from first intron that is retained in a splice variant, revealed it to be a transcriptional repressor. This is a first report on comparative sequence and functional analysis of transcriptional regulation of KCS5 and KCS6; further studies are required before manipulation of cuticular waxes as a strategy for mitigating stress.

The Role of Polymorphisms in Vitamin D-Related Genes in Response to Vitamin D Supplementation

Background. Vitamin D deficiency represents a major healthcare problem. Vitamin D status is influenced by genetic and environmental determinants. Several observational studies have evaluated the association of single-nucleotide polymorphisms (SNPs) in vitamin D-related genes and vitamin D levels. Nevertheless, little is known about the role of these SNPs in the response to vitamin D supplementation. We conducted an interventional study to define the association between SNPs in vitamin D-related genes and the response to vitamin D supplementation in 100 self-reported healthy women of Arab ancestry for the majority. Methods. A total of 100 healthy female subjects received a weekly oral dose of 50,000 IU vitamin D for 12 weeks. Serum vitamin D concentration and metabolic profiles were measured at baseline and 12 weeks post-vitamin D supplementation. The genotypes of 37 SNPs selected from previously reported vitamin D-related genes have been assessed by Fluidigm genotyping assay. Results. Rs731236 (VDR gene) and rs7116978 (CYP2R1 gene) showed a significant association with vitamin D status. The rs731236 GG genotype and the rs7116978 CC genotype were associated with a “vitamin D sufficiency” state. Rs731236 GG and rs7116978 CC genotypes showed a higher response to vitamin D supplementation. Transcription factor binding site prediction analysis showed altered binding sites for transcription factors according to the different rs7116978 alleles. Interestingly, the 37 SNPs previously established to play a role in vitamin D-related pathways explained very little of the response to vitamin D supplementation in our cohort, suggesting the existence of alternative loci whose number and effect size need to be investigated in future studies. Conclusion. In this paper, we present novel data on vitamin D-related SNPs and response to vitamin D supplementation demonstrating the feasibility of applying functional genomic approaches in interventional studies to assess individual-level responses to vitamin D supplementation.

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.

Genome-Wide Identification and Characterization of UTR-Introns of Citrus sinensis

Introns exist not only in coding sequences (CDSs) but also in untranslated regions (UTRs) of a gene. Recent studies in animals and model plants such as Arabidopsis have revealed that the UTR-introns (UIs) are widely presented in most genomes and involved in regulation of gene expression or RNA stability. In the present study, we identified introns at both 5'UTRs (5UIs) and 3'UTRs (3UIs) of sweet orange genes, investigated their size and nucleotide distribution characteristics, and explored the distribution of cis-elements in the UI sequences. Functional category of genes with predicted UIs were further analyzed using GO, KEGG, and PageMan enrichment. In addition, the organ-dependent splicing and abundance of selected UI-containing genes in root, leaf, and stem were experimentally determined. Totally, we identified 825 UI- and 570 3UI-containing transcripts, corresponding to 617 and 469 genes, respectively. Among them, 74 genes contain both 5UI and 3UI. Nucleotide distribution analysis showed that 5UI distribution is biased at both ends of 5'UTR whiles 3UI distribution is biased close to the start site of 3'UTR. Cis- elements analysis revealed that 5UI and 3UI sequences were rich of promoter-enhancing related elements, indicating that they might function in regulating the expression through them. Function enrichment analysis revealed that genes containing 5UI are significantly enriched in the RNA transport pathway. While, genes containing 3UI are significantly enriched in splicesome. Notably, many pentatricopeptide repeat-containing protein genes and the disease resistance genes were identified to be 3UI-containing. RT-PCR result confirmed the existence of UIs in the eight selected gene transcripts whereas alternative splicing events were found in some of them. Meanwhile, qRT-PCR result showed that UIs were differentially expressed among organs, and significant correlation was found between some genes and their UIs, for example: The expression of VPS28 and its 3UI was significantly negative correlated. This is the first report about the UIs in sweet orange from genome-wide level, which could provide evidence for further understanding of the role of UIs in gene expression regulation.