Identification of promoter regions and regulatory sites

Targeted bisulfite sequencing of Scots pine adaptation-related genes

During environmental changes, epigenetic processes can enable adaptive responses faster than natural selection. In plants, very little is known about the role of DNA methylation during long-term adaptation. Scots pine is a widely distributed coniferous species which must adapt to different environmental conditions throughout its long lifespan. Thus, epigenetic modifications may contribute towards this direction. We provide bisulfite next-generation sequencing data from the putative promoters and exons of eight adaptation-related genes (A3IP2, CCA1, COL1, COL2, FTL2, MFT1, PHYO, and ZTL) in three Scots pine populations located in northern and southern parts of Finland. DNA methylation levels were studied in the two seed tissues: the maternal megagametophyte which contributes to embryo viability, and the biparental embryo which represents the next generation. In most genes, differentially methylated cytosines (DMCs) were in line with our previously demonstrated gene expression differences found in the same Scots pine populations. In addition, we found a strong correlation of total methylation levels between the embryo and megagametophyte tissues of a given individual tree, which indicates that DNA methylation can be inherited from the maternal parent. In conclusion, our results imply that DNA methylation differences may contribute to the adaptation of Scots pine populations in different climatic conditions.

R2R3-MYB transcription factor CsMYB60 controls mature fruit skin color by regulating flavonoid accumulation in cucumber

Skin color is an important trait that determines the cosmetic appearance and quality of fruits. In cucumber, the skin color ranges from white to brown in mature fruits. However, the genetic basis for this important trait remains unclear. We conducted a genome-wide association study of natural cucumber populations, along with map-based cloning techniques, on an F2 population resulting from a cross between Pepino (with yellow-brown fruit skin) and Zaoer-N (with creamy fruit skin). We identified CsMYB60 as a candidate gene responsible for skin coloration in mature cucumber fruits. In cucumber accessions with white to pale yellow skin color, a premature stop mutation (C to T) was found in the second exon region of CsMYB60, whereas light yellow cucumber accessions exhibited splicing premature termination caused by an intronic mutator-like element insertion in CsMYB60. Transgenic CsMYB60c cucumber plants displayed a yellow-brown skin color by promoting accumulation of flavonoids, especially hyperoside, a yellow-colored flavonol. CsMYB60c encodes a nuclear protein that primarily acts as a transcriptional activator through its C-terminal activation motif. RNA sequencing and DNA affinity purification sequencing assays revealed that CsMYB60c promotes skin coloration by directly binding to the YYTACCTAMYT motif in the promoter regions of flavonoid biosynthetic genes, including CsF3'H, which encodes flavonoid 3'-hydroxylase. The findings of our study not only offer insight into the function of CsMYB60 as dominantly controlling fruit coloration, but also highlight that intronic DNA mutations can have a similar phenotypic impact as exonic mutations, which may be valuable in future cucumber breeding programs.

Multi-Allelic Mitochondrial DNA Deletions in an Adult Dog with Chronic Weakness, Exercise Intolerance and Lactic Acidemia

(1) Background: An adult dog was presented to a board-certified veterinary neurologist for evaluation of chronic weakness, exercise intolerance and lactic acidemia. (2) Methods: A mitochondrial myopathy was diagnosed based on the histological and histochemical phenotype of numerous COX-negative muscle fibers. Whole-genome sequencing established the presence of multiple extended deletions in the mitochondrial DNA (mtDNA), with the highest prevalence between the 1-11 kb positions of the approximately 16 kb mitochondrial chromosome. Such findings are typically suggestive of an underlying nuclear genome variant affecting mitochondrial replication, repair, or metabolism. (3) Results: Numerous variants in the nuclear genome unique to the case were identified in the whole-genome sequence data, and one, the insertion of a DYNLT1 retrogene, whose parent gene is a regulator of the mitochondrial voltage-dependent anion channel (VDAC), was considered a plausible causal variant. (4) Conclusions: Here, we add mitochondrial deletion disorders to the spectrum of myopathies affecting adult dogs.

APETALA2-like Floral Homeotic Protein Up-Regulating FaesAP1_2 Gene Involved in Floral Development in Long-Homostyle Common Buckwheat

In the rosid species Arabidopsis thaliana, the AP2-type AP2 transcription factor (TF) is required for specifying the sepals and petals identities and confers a major A-function to antagonize the C-function in the outer floral whorls. In the asterid species Petunia, the AP2-type ROB TFs are required for perianth and pistil development, as well as repressing the B-function together with TOE-type TF BEN. In Long-homostyle (LH) Fagopyrum esculentum, VIGS-silencing showed that FaesAP2 is mainly involved in controlling filament and style length, but FaesTOE is mainly involved in regulating filament length and pollen grain development. Both FaesAP2 (AP2-type) and FaesTOE (TOE-type) are redundantly involved in style and/or filament length determination instead of perianth development. However, neither FaesAP2 nor FaesTOE could directly repress the B and/or C class genes in common buckwheat. Moreover, the FaesAP1_2 silenced flower showed tepal numbers, and filament length decreased obviously. Interestingly, yeast one-hybrid (Y1H) and dual-luciferase reporter (DR) further suggested that FaesTOE directly up-regulates FaesAP1_2 to be involved in filament length determination in LH common buckwheat. Moreover, the knockdown of FaesTOE expression could result in expression down-regulation of the directly target FaesAP1_2 in the FaesTOE-silenced LH plants. Our findings uncover a stamen development pathway in common buckwheat and offer deeper insight into the functional evolution of AP2 orthologs in the early-diverging core eudicots.

DNA Methylation Analysis Reveals Potential Mechanism in Takifugu rubripes Against Cryptocaryon irritans Infection

Takifugu rubripes (T. rubripes) is a valuable commercial fish, and Cryptocaryon irritans (C. irritans) has a significant impact on its aquaculture productivity. DNA methylation is one of the earliest discovered ways of gene epigenetic modification and also an important form of modification, as well as an essential type of alteration that regulates gene expression, including immune response. To further explore the anti-infection mechanism of T. rubripes in inhibiting this disease, we determined genome-wide DNA methylation profiles in the gill of T. rubripes using whole-genome bisulfite sequencing (WGBS) and combined with RNA sequence (RNA-seq). A total of 4659 differentially methylated genes (DMGs) in the gene body and 1546 DMGs in the promoter between the infection and control group were identified. And we identified 2501 differentially expressed genes (DEGs), including 1100 upregulated and 1401 downregulated genes. After enrichment analysis, we identified DMGs and DEGs of immune-related pathways including MAPK, Wnt, ErbB, and VEGF signaling pathways, as well as node genes prkcb, myca, tp53, and map2k2a. Based on the RNA-Seq results, we plotted a network graph to demonstrate the relationship between immune pathways and functional related genes, in addition to gene methylation and expression levels. At the same time, we predicted the CpG island and transcription factor of four immune-related key genes prkcb and mapped the gene structure. These unique discoveries could be helpful in the understanding of C. irritans pathogenesis, and the candidate genes screened may serve as optimum methylation-based biomarkers that can be utilized for the correct diagnosis and therapy T. rubripes in the development of the ability to resist C. irritans infection.

Genotype-specific expression of uncle fester suggests a role in allorecognition education in a basal chordate

Histocompatibility is the ability to discriminate between self and non-self tissues, and has been described in species throughout the metazoa. Despite its universal presence, histocompatibility genes utilized by different phyla are unique- those found in sponges, cnidarians, ascidians and vertebrates are not orthologous. Thus, the origins of these sophisticated recognition systems, and any potential functional commonalities between them are not understood. A well-studied histocompatibility system exists in the botryllid ascidians, members of the chordate subphylum, Tunicata, and provides an opportunity to do so. Histocompatibility in the botryllids occurs at the tips of an extracorporeal vasculature that come into contact when two individuals grow into proximity. If compatible, the vessels will fuse, forming a parabiosis between the two individuals. If incompatible, the two vessels will reject- an inflammatory reaction that results in melanin scar formation at the point of contact, blocking anastomosis. Compatibility is determined by a single, highly polymorphic locus called the fuhc with the following rules: individuals that share one or both fuhc alleles will fuse, while those who share neither will reject. The fuhc locus encodes multiple proteins with roles in allorecognition, including one called uncle fester, which is necessary and sufficient to initiate the rejection response. Here we report the existence of genotype-specific expression levels of uncle fester, differing by up to 8-fold at the mRNA-level, and that these expression levels are constant and maintained for the lifetime of an individual. We also found that these differences had functional consequences: the expression level of uncle fester correlated with the speed and severity of the rejection response. These findings support previous conclusions that uncle fester levels modulate the rejection response, and may be responsible for controlling the variation observed in the timing and intensity of the reaction. The maintenance of genotype specific expression of uncle fester is also evidence of an education process reminiscent of that which occurs in mammalian Natural Killer (NK) cells. In turn, this suggests that while histocompatibility receptors and ligands evolve via convergent evolution, they may utilize conserved intracellular machinery to interpret binding events at the cell surface.

Promoter Prediction in Agrobacterium tumefaciens Strain C58 by Using Artificial Intelligence Strategies

Promoters are the genomic regions upstream of genes that RNA polymerase binds in order to initiate gene transcription. Understanding the regulation of gene expression depends on being able to identify promoters, because they are the most important component of gene expression. Agrobacterium tumefaciens (A. tumefaciens) strain C58 was the subject of this study with the goal of creating a machine learning-based model to predict promoters. In this study, nucleotide density (ND), k-mer, and one-hot were used to encode the promoter sequence. Support vector machine (SVM) on fivefold cross-validation with incremental feature selection (IFS) was used to optimize the generated features. These improved characteristics were then used to distinguish promoter sequences by feeding them into the random forest (RF) classifier. Tenfold cross-validation (CV) analysis revealed that the projected model has the ability to produce an accuracy of 84.22%.

Epigenetic reprogramming around IFN1 and IFNy2 promoters in rainbow trout cells inoculated with infectious pancreatic necrosis virus (IPNV)

Infectious pancreatic necrosis virus (IPNV) has proven to effectively evade the host antiviral responses. This study clarifies whether the modulation of the antiviral immune response exerted by IPNV involves epigenetic mechanisms. An in-silico characterization of the rainbow trout IFN1 and IFNγ2 promoters was performed, identifying the islands or sequences rich in CpG dinucleotides and the putative transcription factor binding sites (TBS) for both gene promoters. RTS11 cells (rainbow trout monocyte/macrophage) were infected with IPNV, and the course of viral infection was followed up to 48 h post infection (hpi). Infected cells showed increased IFN1 and IFNγ2 transcriptional expression at 6 and 24 hpi, respectively. IPNV infection caused increases and decreases in global IFNγ2 promoter methylation at 6 and 24 hpi, respectively. The CpG dinucleotides at positions -392 and + 38 of this promoter were the most sensitive to methylation changes. The IFN1 promoter remained fully unmethylated during the course of the infection, similar to the control. The changes in the methylation pattern observed for the IFNγ2 promoter were coincident with the changes in DNA methyltransferase (DNMT) expression levels, increasing at 6 hpi and decreasing below basal level at 24 hpi. Similarly, the H4 histones associated with the IFN1 and IFNγ2 promoters were hyperacetylated at 6 hpi, subsequently decreasing their acetylation below basal levels at 24 hpi, in both promoters. Coincidentally with the above, overexpression of histone acetyltransferase (HAT) was observed at 6 hpi and of histone deacetylase (HDAC) at 24 hpi, with return to baseline of HAT. These results suggest that IPNV would epigenetically modulate the expression of IFN1 by changing acetylation levels of the histones H4 associated with its promoter. Also, the modulation of the expression of IFNy2 would be by switching methylation/demethylation levels of its promoter, in addition to changes in acetylation levels of histones H4 associated with this promoter. This study is the first to demonstrate the effect of epigenetic reprogramming after IPNV infection in salmonid cells, demonstrating that promoter methylation/demethylation level and changes in the histone code associated with promoters may play a role in the modulation of the immune response induced by the virus.

Classification of Promoter Sequences from Human Genome

We have developed a new method for promoter sequence classification based on a genetic algorithm and the MAHDS sequence alignment method. We have created four classes of human promoters, combining 17,310 sequences out of the 29,598 present in the EPD database. We searched the human genome for potential promoter sequences (PPSs) using dynamic programming and position weight matrices representing each of the promoter sequence classes. A total of 3,065,317 potential promoter sequences were found. Only 1,241,206 of them were located in unannotated parts of the human genome. Every other PPS found intersected with either true promoters, transposable elements, or interspersed repeats. We found a strong intersection between PPSs and Alu elements as well as transcript start sites. The number of false positive PPSs is estimated to be 3 × 10-8 per nucleotide, which is several orders of magnitude lower than for any other promoter prediction method. The developed method can be used to search for PPSs in various eukaryotic genomes.

Codon adaptation by synonymous mutations impacts the functional properties of the estrogen receptor-alpha protein in breast cancer cells

Oestrogen receptor-alpha (ERα) positivity is intimately associated with the development of hormone-dependent breast cancers. A major challenge in the treatment of these cancers is to understand and overcome the mechanisms of endocrine resistance. Recently, two distinct translation programmes using specific transfer RNA (tRNA) repertoires and codon usage frequencies were evidenced during cell proliferation and differentiation. Considering the phenotype switch of cancer cells to more proliferating and less-differentiated states, we can speculate that the changes in the tRNA pool and codon usage that likely occur make the ERα coding sequence no longer adapted, impacting translational rate, co-translational folding and the resulting functional properties of the protein. To verify this hypothesis, we generated an ERα synonymous coding sequence whose codon usage was optimized to the frequencies observed in genes expressed specifically in proliferating cells and then investigated the functional properties of the encoded receptor. We demonstrate that such a codon adaptation restores ERα activities to levels observed in differentiated cells, including: (a) an enhanced contribution exerted by transactivation function 1 (AF1) in ERα transcriptional activity; (b) enhanced interactions with nuclear receptor corepressor 1 and 2 [NCoR1 and NCoR2 (also known as SMRT) respectively], promoting repressive capability; and (c) reduced interactions with SRC proto-oncogene, non-receptor tyrosine kinase (Src) and phosphoinositide 3-kinase (PI3K) p85 kinases, inhibiting MAPK and AKT signalling pathway.

The aryl hydrocarbon receptor cell intrinsically promotes resident memory CD8+ T cell differentiation and function

The Aryl hydrocarbon receptor (Ahr) regulates the differentiation and function of CD4+ T cells; however, its cell-intrinsic role in CD8+ T cells remains elusive. Herein we show that Ahr acts as a promoter of resident memory CD8+ T cell (TRM) differentiation and function. Genetic ablation of Ahr in mouse CD8+ T cells leads to increased CD127-KLRG1+ short-lived effector cells and CD44+CD62L+ T central memory cells but reduced granzyme-B-producing CD69+CD103+ TRM cells. Genome-wide analyses reveal that Ahr suppresses the circulating while promoting the resident memory core gene program. A tumor resident polyfunctional CD8+ T cell population, revealed by single-cell RNA-seq, is diminished upon Ahr deletion, compromising anti-tumor immunity. Human intestinal intraepithelial CD8+ T cells also highly express AHR that regulates in vitro TRM differentiation and granzyme B production. Collectively, these data suggest that Ahr is an important cell-intrinsic factor for CD8+ T cell immunity.

Motif and conserved module analysis in DNA (promoters, enhancers) and RNA (lncRNA, mRNA) using AlModules

Nucleic acid motifs consist of conserved and variable nucleotide regions. For functional action, several motifs are combined to modules. The tool AIModules allows identification of such motifs including combinations of them and conservation in several nucleic acid stretches. AIModules recognizes conserved motifs and combinations of motifs (modules) allowing a number of interesting biological applications such as analysis of promoter and transcription factor binding sites (TFBS), identification of conserved modules shared between several gene families, e.g. promoter regions, but also analysis of shared and conserved other DNA motifs such as enhancers and silencers, in mRNA (motifs or regulatory elements e.g. for polyadenylation) and lncRNAs. The tool AIModules presented here is an integrated solution for motif analysis, offered as a Web service as well as downloadable software. Several nucleotide sequences are queried for TFBSs using predefined matrices from the JASPAR DB or by using one's own matrices for diverse types of DNA or RNA motif discovery. Furthermore, AIModules can find TFBSs common to two or more sequences. Demanding high or low conservation, AIModules outperforms other solutions in speed and finds more modules (specific combinations of TFBS) than alternative available software. The application also searches RNA motifs such as polyadenylation site or RNA-protein binding motifs as well as DNA motifs such as enhancers as well as user-specified motif combinations ( https://bioinfo-wuerz.de/aimodules/ ; alternative entry pages: https://aimodules.heinzelab.de or https://www.biozentrum.uni-wuerzburg.de/bioinfo/computing/aimodules ). The application is free and open source whether used online, on-site, or locally.

Moderate DNA methylation changes associated with nitrogen remobilization and leaf senescence in Arabidopsis

The lifespan of plants is restricted by environmental and genetic components. Following the transition to reproductive growth, leaf senescence ends cellular life in monocarpic plants to remobilize nutrients to storage organs. In Arabidopsis, we initially observed altered leaf to seed ratios, faster senescence progression, altered leaf nitrogen recovery after transient nitrogen removal, and ultimately enhanced nitrogen remobilization from the leaves in two methylation mutants (ros1 and the triple dmr1/2 cmt3 knockout). Analysis of the DNA methylome in wild type Col-0 leaves identified an initial moderate decline of cytosine methylation with progressing leaf senescence, predominantly in the CG context. Late senescence was associated with moderate de novo methylation of cytosines, primarily in the CHH context. Relatively few differentially methylated regions, including one in the ROS1 promoter linked to down-regulation of ROS1, were present, but these were unrelated to known senescence-associated genes. Differential methylation patterns were identified in transcription factor binding sites, such as the W-boxes that are targeted by WRKYs. Methylation in artificial binding sites impaired transcription factor binding in vitro. However, it remains unclear how moderate methylome changes during leaf senescence are linked with up-regulated genes during senescence.

De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Database of Potential Promoter Sequences in the Capsicum annuum Genome

In this study, we used a mathematical method for the multiple alignment of highly divergent sequences (MAHDS) to create a database of potential promoter sequences (PPSs) in the Capsicum annuum genome. To search for PPSs, 20 statistically significant classes of sequences located in the range from −499 to +100 nucleotides near the annotated genes were calculated. For each class, a position–weight matrix (PWM) was computed and then used to identify PPSs in the C. annuum genome. In total, 825,136 PPSs were detected, with a false positive rate of 0.13%. The PPSs obtained with the MAHDS method were tested using TSSFinder, which detects transcription start sites. The databank of the found PPSs provides their coordinates in chromosomes, the alignment of each PPS with the PWM, and the level of statistical significance as a normal distribution argument, and can be used in genetic engineering and biotechnology.

The Epigenetic Regulation in Plant Specialized Metabolism: DNA Methylation Limits Paclitaxel in vitro Biotechnological Production

Environmental conditions are key factors in the modulation of the epigenetic mechanisms regulating gene expression in plants. Specifically, the maintenance of cell cultures in optimal in vitro conditions alters methylation patterns and, consequently, their genetic transcription and metabolism. Paclitaxel production in Taxus x media cell cultures is reduced during its maintenance in in vitro conditions, compromising the biotechnological production of this valuable anticancer agent. To understand how DNA methylation influences taxane production, the promoters of three genes (GGPPS, TXS, and DBTNBT) involved in taxane biosynthesis have been studied, comparing the methylation patterns between a new line and one of ~14 years old. Our work revealed that while the central promoter of the GGPPS gene is protected from cytosine methylation accumulation, TXS and DBTNBT promoters accumulate methylation at different levels. The DBTNBT promoter of the old line is the most affected, showing a 200 bp regulatory region where all the cytosines were methylated. This evidence the existence of specific epigenetic regulatory mechanisms affecting the last steps of the pathway, such as the DBTNBT promoter. Interestingly, the GGPPS promoter, a regulatory sequence of a non-specific taxane biosynthetic gene, was not affected by this mechanism. In addition, the relationship between the detected methylation points and the predicted transcription factor binding sites (TFBS) showed that the action of TFs would be compromised in the old line, giving a further explanation for the production reduction in in vitro cell cultures. This knowledge could help in designing novel strategies to enhance the biotechnological production of taxanes over time.

A New Homotetramer Hemoglobin in the Pulmonary Surfactant of Plateau Zokors (Myospalax Baileyi)

The plateau zokor (Myospalax baileyi) is a native species to the Qinghai-Tibetan Plateau, inhabiting hypoxia and hypercapnia sealed subterranean burrows that pose several unique physiological challenges. In this study, we observed a novel heme-containing protein in the pulmonary surfactant (PS) of plateau zokor, identified the encoding gene of the protein, predicted its origination and structure, verified its expression in alveolar epithelial cells, and determined the protein’s affinity to oxygen and its effect on the oxygen-dissolving capability in the PS of plateau zokors. The protein is an unusual homotetramer hemoglobin consisting of four γ-like subunits, and the subunit is encoded by a paralog gene of γ, that is γ-like. The divergence time of γ-like from γ is estimated by the molecular clock to be about 2.45 Mya. The generation of γ-like in plateau zokors might well relate to long-time stress of the high land hypoxia. Unlike γ, the γ-like has a hypoxia response element (HRE) and a lung tissue-specific enhancer in its upstream region, and it is expressed specifically in lung tissues and up-regulated by hypoxia. The protein is named as γ₄-like which is expressed specifically in Alveolar epithelial type II (ATII) cells and secreted into the alveolar cavities through the osmiophilic multilamellar body (LBs). The γ₄-like has a higher affinity to oxygen, and that increases significantly oxygen-dissolving capability in the PS of plateau zokors by its oxygenation function, which might be beneficial for the plateau zokors to obtain oxygen from the severe hypoxia environments by facilitating oxygen diffusion from alveoli to blood.

IMGT® Biocuration and Analysis of the Rhesus Monkey IG Loci

The adaptive immune system, along with the innate immune system, are the two main biological processes that protect an organism from pathogens. The adaptive immune system is characterized by the specificity and extreme diversity of its antigen receptors. These antigen receptors are the immunoglobulins (IG) or antibodies of the B cells and the T cell receptors (TR) of the T cells. The IG are proteins that have a dual role in immunity: they recognize antigens and trigger elimination mechanisms, to rid the body of foreign cells. The synthesis of the immunoglobulin heavy and light chains requires gene rearrangements at the DNA level in the IGH, IGK, and IGL loci. The rhesus monkey (Macaca mulatta) is one of the most widely used nonhuman primate species in biomedical research. In this manuscript, we provide a thorough analysis of the three IG loci of the Mmul_10 assembly of rhesus monkey, integrating IMGT previously existing data. Detailed characterization of IG genes includes their localization and position in the loci, the determination of the allele functionality, and the description of the regulatory elements of their promoters as well as the sequences of the conventional recombination signals (RS). This complete annotation of the genomic IG loci of Mmul_10 assembly and the highly detailed IG gene characterization could be used as a model, in additional rhesus monkey assemblies, for the analysis of the IG allelic polymorphism and structural variation, which have been described in rhesus monkeys.

ZmERF21 directly regulates hormone signaling and stress-responsive gene expression to influence drought tolerance in maize seedlings

Drought stress adversely impacts crop development and yield. Maize frequently encounters drought stress during its life cycle. Improvement of drought tolerance is a priority of maize breeding programs. Here, we identified a novel transcription factor encoding gene, APETALA2 (AP2)/Ethylene response factor (ERF), which is tightly associated with drought tolerance in maize seedlings. ZmERF21 is mainly expressed in the root and leaf and it can be highly induced by polyethylene glycol treatment. Genetic analysis showed that the zmerf21 mutant plants displayed a reduced drought tolerance phenotype, accompanied by phenotypical and physiological changes that are commonly observed in drought conditions. Overexpression of ZmERF21 in maize significantly increased the chlorophyll content and activities of antioxidant enzymes under drought conditions. RNA-Seq and DNA affinity purification sequencing analysis further revealed that ZmERF21 may directly regulate the expression of genes related to hormone (ethylene, abscisic acid) and Ca signaling as well as other stress-response genes through binding to the promoters of potential target genes. Our results thereby provided molecular evidence of ZmERF21 is involved in the drought stress response of maize.

Characterization of Differentially Expressed Genes under Salt Stress in Olive

Climate change, currently taking place worldwide and also in the Mediterranean area, is leading to a reduction in water availability and to groundwater salinization. Olive represents one of the most efficient tree crops to face these scenarios, thanks to its natural ability to tolerate moderate salinity and drought. In the present work, four olive cultivars (Koroneiki, Picual, Royal de Cazorla and Fadak86) were exposed to high salt stress conditions (200 mM of NaCl) in greenhouse, in order to evaluate their tolerance level and to identify key genes involved in salt stress response. Molecular and physiological parameters, as well as plant growth and leaves’ ions Na⁺ and K⁺ content were measured. Results of the physiological measurements showed Royal de Cazorla as the most tolerant cultivar, and Fadak86 and Picual as the most susceptible ones. Ten candidate genes were analyzed and their complete genomic, CDS and protein sequences were identified. The expression analysis of their transcripts through reverse transcriptase quantitative PCR (RT-qPCR) demonstrated that only OeNHX7, OeP5CS, OeRD19A and OePetD were upregulated in tolerant cultivars, thus suggesting their key role in the activation of a salt tolerance mechanism.

The 5'-upstream region of WRKY18 transcription factor from banana is a stress-inducible promoter with strong expression in guard cells

Increasing crop productivity in an ever-changing environmental scenario is a major challenge for maintaining the food supply worldwide. Generation of crops having broad-spectrum pathogen resistance with the ability to cope with water scarcity is the only solution to feed the expanding world population. Stomatal closure has implications on pathogen colonization and drought tolerance. Recent studies have provided novel insights into networks involved in stomatal closure which is being used in biotechnological applications for improving crop endurance. Despite that genetic engineering of stomata requires guard cell preferred or specific regulatory regions to avoid undesirable side effects. In the present study, we describe the 5'-upstream regulatory region of the WRKY18 transcription factor of banana and functionally analyzed its stress meditated activation and strong guard cell preferred activity. Expression of MusaWRKY18 is augmented in leaves of banana cultivars Karibale Monthan, Rasthali and Grand Nain under multiple stress conditions suggesting its role in stress responses of banana plants. Transgenic tobacco lines harboring P_MusaWRKY18 -β-D-glucuronidase (GUS) were regenerated and GUS staining demonstrated substantial GUS expression in guard cells which corroborates with multiple Dof1 binding cis-elements in P_MusaWRKY18 . Fluorescent β-galactosidase assay demonstrated the stress-mediated strong induction profiles of P_MusaWRKY18 at different time points in transgenic tobacco lines exposed to drought, high-salinity, cold, and applications of abscisic acid, salicylic acid, methyl jasmonate, and ethephon. This study sheds novel insights into guard cell preferred expression of WRKY genes under stress and confirm the utility of P_MusaWRKY18 for exploring guard cell functions and guard cell engineering.

Exploring the diversity of promoter and 5'UTR sequences in ancestral, historic and modern wheat

A data set of promoter and 5'UTR sequences of homoeo-alleles of 459 wheat genes that contribute to agriculturally important traits in 95 ancestral and commercial wheat cultivars is presented here. The high-stringency myBaits technology used made individual capture of homoeo-allele promoters possible, which is reported here for the first time. Promoters of most genes are remarkably conserved across the 83 hexaploid cultivars used with <7 haplotypes per promoter and 21% being identical to the reference Chinese Spring. InDels and many high-confidence SNPs are located within predicted plant transcription factor binding sites, potentially changing gene expression. Most haplotypes found in the Watkins landraces and a few haplotypes found in Triticum monococcum, germplasms hitherto not thought to have been used in modern wheat breeding, are already found in many commercial hexaploid wheats. The full data set which is useful for genomic and gene function studies and wheat breeding is available at https://rrescloud.rothamsted.ac.uk/index.php/s/DMCFDu5iAGTl50u/authenticate.

Glutamate Dehydrogenase (GdhA) of Streptococcus pneumoniae Is Required for High Temperature Adaptation

During its progression from the nasopharynx to other sterile and nonsterile niches of its human host, Streptococcus pneumoniae must cope with changes in temperature. We hypothesized that the temperature adaptation is an important facet of pneumococcal survival in the host. Here, we evaluated the effect of temperature on pneumococcus and studied the role of glutamate dehydrogenase (GdhA) in thermal adaptation associated with virulence and survival. Microarray analysis revealed a significant transcriptional response to changes in temperature, affecting the expression of 252 genes in total at 34°C and 40°C relative to at 37°C. One of the differentially regulated genes was gdhA, which is upregulated at 40°C and downregulated at 34°C relative to 37°C. Deletion of gdhA attenuated the growth, cell size, biofilm formation, pH survival, and biosynthesis of proteins associated with virulence in a temperature-dependent manner. Moreover, deletion of gdhA stimulated formate production irrespective of temperature fluctuation. Finally, ΔgdhA grown at 40°C was less virulent than other temperatures or the wild type at the same temperature in a Galleria mellonella infection model, suggesting that GdhA is required for pneumococcal virulence at elevated temperature.

TSSFinder-fast and accurate ab initio prediction of the core promoter in eukaryotic genomes

Promoter annotation is an important task in the analysis of a genome. One of the main challenges for this task is locating the border between the promoter region and the transcribing region of the gene, the transcription start site (TSS). The TSS is the reference point to delimit the DNA sequence responsible for the assembly of the transcribing complex. As the same gene can have more than one TSS, so to delimit the promoter region, it is important to locate the closest TSS to the site of the beginning of the translation. This paper presents TSSFinder, a new software for the prediction of the TSS signal of eukaryotic genes that is significantly more accurate than other available software. We currently are the only application to offer pre-trained models for six different eukaryotic organisms: Arabidopsis thaliana, Drosophila melanogaster, Gallus gallus, Homo sapiens, Oryza sativa and Saccharomyces cerevisiae. Additionally, our software can be easily customized for specific organisms using only 125 DNA sequences with a validated TSS signal and corresponding genomic locations as a training set. TSSFinder is a valuable new tool for the annotation of genomes. TSSFinder source code and docker container can be downloaded from http://tssfinder.github.io. Alternatively, TSSFinder is also available as a web service at http://sucest-fun.org/wsapp/tssfinder/.

E2F1-Induced lncRNA BAIAP2-AS1 Overexpression Contributes to the Malignant Progression of Hepatocellular Carcinoma via miR-361-3p/SOX4 Axis

Currently, plenty of researches have revealed that long noncoding RNAs (lncRNAs) can act as crucial roles during the progression of various tumors, including hepatocellular carcinoma (HCC). Here, we measured the expression of lncRNA BAIAP2 antisense RNA 1(BAIAP2-AS1) as well as its contribution to the developments of HCC. In this study, the expressions of BAIAP2-AS1 and SOX4 were distinctly upregulated in HCC cells and tissues, and high BAIAP2-AS1 may be a novel biomarker for HCC. E2F1 activated BAIAP2-AS1 expression. The silence of BAIAP2-AS1 inhibited the proliferation and metastasis of HepG2 and PLC5 cells. Assays for relationship verification showed that BAIAP2-AS1 regulated the expression of SOX4 and miR-361-3p. Rescue experiments further confirmed the positive interaction between miR-361-3p and BAIAP2-AS1 as well as between miR-361-3p and SOX4. Overall, BAIAP2-AS1 modulated the miR-361-3p/SOX4 axis to promote the development of HCC. Thus, our study offers a potential therapeutic target for treating HCC.

Convolutional neural networks (CNNs): concepts and applications in pharmacogenomics

Convolutional neural networks (CNNs) have been used to extract information from various datasets of different dimensions. This approach has led to accurate interpretations in several subfields of biological research, like pharmacogenomics, addressing issues previously faced by other computational methods. With the rising attention for personalized and precision medicine, scientists and clinicians have now turned to artificial intelligence systems to provide them with solutions for therapeutics development. CNNs have already provided valuable insights into biological data transformation. Due to the rise of interest in precision and personalized medicine, in this review, we have provided a brief overview of the possibilities of implementing CNNs as an effective tool for analyzing one-dimensional biological data, such as nucleotide and protein sequences, as well as small molecular data, e.g., simplified molecular-input line-entry specification, InChI, binary fingerprints, etc., to categorize the models based on their objective and also highlight various challenges. The review is organized into specific research domains that participate in pharmacogenomics for a more comprehensive understanding. Furthermore, the future intentions of deep learning are outlined.

Disparate Dynamics of Gene Body and cis-Regulatory Element Evolution Illustrated for the Senescence-Associated Cysteine Protease Gene SAG12 of Plants

Gene regulation networks precisely orchestrate the expression of genes that are closely associated with defined physiological and developmental processes such as leaf senescence in plants. The Arabidopsis thaliana senescence-associated gene 12 (AtSAG12) encodes a cysteine protease that is (i) involved in the degradation of chloroplast proteins and (ii) almost exclusively expressed during senescence. Transcription factors, such as WRKY53 and WRKY45, bind to W-boxes in the promoter region of AtSAG12 and play key roles in its activation. Other transcription factors, such as bZIPs, might have accessory functions in their gene regulation, as several A-boxes have been identified and appear to be highly overrepresented in the promoter region compared to the whole genome distribution but are not localized within the regulatory regions driving senescence-associated expression. To address whether these two regulatory elements exhibiting these different properties are conserved in other closely related species, we constructed phylogenetic trees of the coding sequences of orthologs of AtSAG12 and screened their respective 2000 bp promoter regions for the presence of conserved cis-regulatory elements, such as bZIP and WRKY binding sites. Interestingly, the functional relevant upstream located W-boxes were absent in plant species as closely related as Arabidopsis lyrata, whereas an A-box cluster appeared to be conserved in the Arabidopsis species but disappeared in Brassica napus. Several orthologs were present in other species, possibly because of local or whole genome duplication events, but with distinct cis-regulatory sites in different locations. However, at least one gene copy in each family analyzed carried one W-box and one A-box in its promoter. These gene differences in SAG12 orthologs are discussed in the framework of cis- and trans-regulatory factors, of promoter and gene evolution, of genetic variation, and of the enhancement of the adaptability of plants to changing environmental conditions.

A Small Regulatory RNA Generated from the malK 5' Untranslated Region Targets Gluconeogenesis in Vibrio Species

Vsr217 is a small RNA from Vibrio tasmaniensis LGP32, a pathogen associated with mortality events affecting juvenile oysters. The vsr217 gene is located within the 5' untranslated region (UTR) of malK, encoding the ATPase component of the maltose importer, and is conserved within the genus Vibrio. In the presence of maltose, vsr217 is regulated by MalT, the positive regulator of the maltose regulon. vsr217 is required in cis for the full expression of malK. In addition, Vsr217 acts in trans to downregulate the expression of fbp encoding fructose-1,6-bisphosphatase, an enzyme involved in gluconeogenesis. Thus, in the presence of maltose, the induction of Vsr217 is expected to promote glycolysis by negatively regulating the expression of a key enzyme of gluconeogenesis. IMPORTANCE Juvenile pacific oysters have been subject in recent years to summer mortality episodes with deep economic consequences. The pathogen Vibrio tasmaniensis has been associated with such mortality events. For bacterial pathogens, survival within the host requires profound metabolic adaptations according to available resources. All kinds of regulatory elements, including noncoding RNAs, orchestrate this response. Oysters are rich in glycogen, a precursor of maltose, and we previously reported that V. tasmaniensis maltose-regulated genes are strongly induced during oyster infection. Here, we report the dual mechanism by which a small regulatory RNA, generated from the 5' untranslated region of a gene belonging to the maltose regulon, acts both in cis and trans. In cis, it stimulates growth on maltose, and in trans, it downregulates the expression of a gene associated with gluconeogenesis, thus coordinating maltose utilization with central carbon metabolism.

Mathematical Algorithm for Identification of Eukaryotic Promoter Sequences

Identification of promoter sequences in the eukaryotic genome, by computer methods, is an important task of bioinformatics. However, this problem has not been solved since the best algorithms have a false positive probability of 10−3–10−4 per nucleotide. As a result of full genome analysis, there may be more false positives than annotated gene promoters. The probability of a false positive should be reduced to 10−6–10−8 to reduce the number of false positives and increase the reliability of the prediction. The method for multi alignment of the promoter sequences was developed. Then, mathematical methods were developed for calculation of the statistically important classes of the promoter sequences. Five promoter classes, from the rice genome, were created. We developed promoter classes to search for potential promoter sequences in the rice genome with a false positive number less than 10−8 per nucleotide. Five classes of promoter sequences contain 1740, 222, 199, 167 and 130 promoters, respectively. A total of 145,277 potential promoter sequences (PPSs) were identified. Of these, 18,563 are promoters of known genes, 87,233 PPSs intersect with transposable elements, and 37,390 PPSs were found in previously unannotated sequences. The number of false positives for a randomly mixed rice genome is less than 10−8 per nucleotide. The method developed for detecting PPSs was compared with some previously used approaches. The developed mathematical method can be used to search for genes, transposable elements, and transcript start sites in eukaryotic genomes.

A 598-bp InDel Variation in the Promoter Region of Bna.SOC1.A05 Is Predominantly Present in Winter Type Rapeseeds

During rapeseed domestication and breeding, genetic diversity allowed to adapt it to different eco-geographical regions and to shape its useful traits. Structural variations (SVs), including presence/absence variations (PAVs), are thought to play a major role in the genetic diversity and phenotypic plasticity of rapeseed. In this study, we detected a 598-bp PAV within the promoter region of an Arabidopsis ortholog of a major flowering time gene and a downstream target of FLC, SOC1, which is one of the first genes that are upregulated in rapeseed during vernalization. Further analysis showed that the insertion is present predominantly in winter types while absent in spring types. The 589-bp sequence is present only in the A sub-genome indicating that it originated from Brassica rapa. Since the genomic region around Bna.SOC1.A05 showed a strong reduction in nucleotide diversity, the insertion might represent a larger selected sweep for rapeseed adaptation. Cis-element analysis showed that the insertion contains an ACGTG box, which is the strongest binding motif for the HY5 transcription factor in Arabidopsis. In addition, expression analyses showed that mRNA levels of Bna.SOC1.A05 were lower in accessions carrying the insertion compared to the ones that had no insertion.

Use of Mathematical Methods for the Biosafety Assessment of Agricultural Crops

In Russia and around the world, there are important questions regarding the potential threats to national and biological safety created by genetic technologies and the need to improve or introduce new, justified, and adequate measures for their control, regulation, and prevention. The article shows that a significant volume of the global market is occupied by five major transgenic crops, and producers are ready to switch to crops with an edited genome that has been approved in the United States, Argentina, and other countries. We propose a qualitatively new approach to the risk assessment of edited plants, "Safe Design," and we have also developed an extremely important, fundamentally new approach to the development of methods that combine next-generation sequencing (NGS) and Bioinformatics for the assessment of the crop import biosafety. The proposed mathematical approach provides a detailed analysis of the possible insertions of DNA fragments into the genome of edited crops and a clarification of their biological significance. The developed method can be used in the rapid screening of plants for the presence of potentially dangerous genes, viral sequences, and nonspecific promoter sequences.

Multiple Alignment of Promoter Sequences from the Arabidopsis thaliana L. Genome

In this study, we developed a new mathematical method for performing multiple alignment of highly divergent sequences (MAHDS), i.e., sequences that have on average more than 2.5 substitutions per position (x). We generated sets of artificial DNA sequences with x ranging from 0 to 4.4 and applied MAHDS as well as currently used multiple sequence alignment algorithms, including ClustalW, MAFFT, T-Coffee, Kalign, and Muscle to these sets. The results indicated that most of the existing methods could produce statistically significant alignments only for the sets with x < 2.5, whereas MAHDS could operate on sequences with x = 4.4. We also used MAHDS to analyze a set of promoter sequences from the Arabidopsis thaliana genome and discovered many conserved regions upstream of the transcription initiation site (from -499 to +1 bp); a part of the downstream region (from +1 to +70 bp) also significantly contributed to the obtained alignments. The possibilities of applying the newly developed method for the identification of promoter sequences in any genome are discussed. A server for multiple alignment of nucleotide sequences has been created.

First insights into the molecular basis association between promoter polymorphisms of the IL1B gene and Helicobacter pylori infection in the Sudanese population: computational approach

BACKGROUND

Helicobacter pylori (H. pylori) infects nearly half of the world's population with a variation in incidence among different geographic regions. Genetic variants in the promoter regions of the IL1B gene can affect cytokine expression and creates a condition of hypoacidity which favors the survival and colonization of H. pylori. Therefore, the aim of this study was to characterize the polymorphic sites in the 5'- region [-687_ + 297] of IL1B in H. pylori infection using in silico tools.

RESULTS

A total of five nucleotide variations were detected in the 5'-regulatory region [-687_ + 297] of IL1B which led to the addition or alteration of transcription factor binding sites (TFBSs) or composite regulatory elements (CEs). Genotyping of IL1B - 31 C > T revealed a significant association between -31 T and susceptibility to H. pylori infection in the studied population (P = 0.0363). Comparative analysis showed conservation rates of IL1B upstream [-368_ + 10] region above 70% in chimpanzee, rhesus monkey, a domesticated dog, cow and rat.

CONCLUSIONS

In H. pylori-infected patients, three detected SNPs (- 338, - 155 and - 31) located in the IL1B promoter were predicted to alter TFBSs and CE, which might affect the gene expression. These in silico predictions provide insight for further experimental in vitro and in vivo studies of the regulation of IL1B expression and its relationship to H. pylori infection. However, the recognition of regulatory motifs by computer algorithms is fundamental for understanding gene expression patterns.

Prediction of Rice Transcription Start Sites Using TransPrise: A Novel Machine Learning Approach

As the interest in genetic resequencing increases, so does the need for effective mathematical, computational, and statistical approaches. One of the difficult problems in genome annotation is determination of precise positions of transcription start sites. In this paper, we present TransPrise-an efficient deep learning tool for predicting positions of eukaryotic transcription start sites. TransPrise offers significant improvement over existing promoter-prediction methods. To illustrate this, we compared predictions of TransPrise with the TSSPlant approach for well-annotated genome of Oryza sativa. Using a computer with a graphics processing unit, the run time of TransPrise is 250 min on a genome of 374 Mb long.We provide the full basis for the comparison and encourage users to freely access a set of our computational tools to facilitate and streamline their own analyses. The ready-to-use Docker image with all the necessary packages, models, and code as well as the source code of the TransPrise algorithm are available at http://compubioverne.group/ . The source code is ready to use and to be customized to predict TSS in any eukaryotic organism.

A Defect in the Twin-Arginine Translocation Pathway Decreases the Tolerance of Xanthomonas campestris pv. campestris to Phenazines

Phenazine-1-carboxylic acid (PCA), a member of phenazines secreted by microorganisms, inhibits the growth of many bacteria and fungi. Xanthomonas campestris pv. campestris is the causal agent of black rot, the most important disease of cruciferous crops worldwide, and is more tolerant to PCA than other Xanthomonas species. Previous studies reported that reactive oxygen species (ROS) scavenging ability is involved in regulating the PCA tolerance of Xanthomonas species. Additionally, the cytochrome c maturation (CCM) system has been found to play a more important role in tolerance to phenazines than the ROS scavenging system. In this study, a highly PCA-sensitive insertion mutant of X. campestris pv. campestris, X-5, was identified and studied. The insertion site of X-5 was found to be in tatB gene (XC_4183), which encodes a subunit of the twin-arginine translocation (TAT) complex. Disruption of the three genes of TAT pathway resulted in decreased biological fitness and reduced tolerance to phenazines in comparison with the wild-type strain 8004. These results imply that the tolerance mechanism of the TAT pathway to phenazines is related to the CCM system, but not due to the ROS scavenging system. Furthermore, respiration-related characteristic tests and peptide analysis suggested that disruption of the TAT complex causes a defect in the cytochrome bc1 complex, which may be involved in the tolerance to phenazines. In summary, this study sheds new light on the critical role of the TAT pathway in influencing the fitness and phenazines tolerance of Xanthomonas species.

Heterologous Gene Expression System Using the Cold-Inducible CnAFP Promoter in Chlamydomonas reinhardtii

To increase the availability of microalgae as producers of valuable compounds, it is necessary to develop novel systems for gene expression regulation. Among the diverse expression systems available in microalgae, none are designed to induce expression by low temperature. In this study, we explored a cold-inducible system using the antifreeze protein (AFP) promoter from a polar diatom, Chaetoceros neogracile. A vector containing the CnAFP promoter (pCnAFP) was generated to regulate nuclear gene expression, and reporter genes (Gaussia luciferase (GLuc) and mVenus fluorescent protein (mVenus)) were successfully expressed in the model microalga, Chlamydomonas reinhardtii. In particular, under the control of pCnAFP, the expression of these genes was increased at low temperature, unlike pAR1, a promoter that is widely used for gene expression in C. reinhardtii. Promoter truncation assays showed that cold inducibility was still present even when pCnAFP was shortened to 600 bp, indicating the presence of a low-temperature response element between -600 and -477 bp. Our results show the availability of new heterologous gene expression systems with cold-inducible promoters and the possibility to find novel low-temperature response factors in microalgae. Through further improvement, this cold-inducible promoter could be used to develop more efficient expression tools.

Serovar-dependent differences in Hfq-regulated phenotypes inActinobacillus pleuropneumoniae

The RNA chaperone Hfq regulates diverse processes in numerous bacteria. In this study, we compared phenotypes (growth rate, adherence, response to different stress conditions and virulence in Galleria mellonella) of wild-type (WT) and isogenic hfq mutants of three serovars (1, 8 and 15) of the porcine pathogen Actinobacillus pleuropneumoniae. Similar growth in rich broth was seen for all strains except Ap1∆hfq, which showed slightly reduced growth throughout the 24 h time course, and the complemented Ap8∆hfqC mutant had a prolonged lag phase. Differences were seen between the three serovar WT strains regarding adherence, stress response and virulence in G. mellonella, and deletion of hfq affected some, but not all of these phenotypes, depending on serovar. Complementation by expression of cloned hfq from an endogenous promoter only restored some WT phenotypes, indicating that complex regulatory networks may be involved, and that levels of Hfq may be as important as presence/absence of the protein regarding its contribution to gene regulation. Our results support that Hfq is a pleiotropic global regulator in A. pleuropneumoniae, but serovar-related differences exist. These results highlight the importance of testing multiple strains/serovars within a given species when determining contributions of global regulators, such as Hfq, to expression of complex phenotypes.

Ginsenosides regulate adventitious root formation in Panax ginseng via a CLE45-WOX11 regulatory module

Adventitious root branching is vital to plant growth and regeneration, but the regulation of this process remains unclear. We therefore investigated how ginsenosides regulate adventitious root branching in Panax ginseng. Cell proliferation and adventitious root branching were decreased in the presence of ginsenoside Rb1 and a high concentration of ginsenoside Re, but increased when treating with a low concentration of Re. Moreover, the exogenous application of a synthetic dodeca-amino acid peptide that has a CLAVATA3/EMBRYO SURROUNDING REGION-related (CLE) motif corresponding to PgCLE45 retarded root growth in both ginseng and Arabidopsis. The root Re levels and the expression of the DDS, CYP716A47, and CYP716A53 genes that encode enzymes involved in ginsenoside synthesis were decreased in the presence of PgCLE45. The expression profiles of PgWOX and PgCLE genes were determined to further investigate the CLE-WOX signaling pathway. The levels of PgWOX11 transcripts showed an inverse pattern to PgCLE45 transcripts. Using yeast one-hybrid assay, EMSA, and ChIP assay, we showed that PgWOX11 bound to the PgCLE45 promoter, which contained the HD motif. Transient expression assay showed that PgWOX11 induced the expression of PgCLE45 in adventitious roots, while PgCLE45 suppressed the expression of PgWOX11. These results suggest that there is a negative feedback regulation between PgCLE45 and PgWOX11. Taken together, these data show that ginsenosides regulate adventitious root branching via a novel PgCLE45-PgWOX11 regulatory loop, providing a potential mechanism for the regulation of adventitious root branching.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max

BACKGROUND

To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region.

RESULTS

RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes.

CONCLUSIONS

We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

Intronic regions of the human coagulation factor VIII gene harboring transcription factor binding sites with a strong bias towards the short-interspersed elements

Increasing data show that intronic derived regulatory elements, such as transcription factor binding sites (TFBs), play key roles in gene regulation, and malfunction. Accordingly, characterizing the sequence context of the intronic regions of the human coagulation factor VIII (hFVIII) gene can be important. In this study, the intronic regions of the hFVIII gene were scrutinized based on in-silico methods. The results disclosed that these regions harbor a rich array of functional elements such as repetitive elements (REs), splicing sites, and transcription factor binding sites (TFBs). Among these elements, TFBs and REs showed a significant distribution and correlation to each other. This survey indicated that 31% of TFBs are localized in the intronic regions of the gene. Moreover, TFBs indicate a strong bias in the regions far from splice sites of introns with mapping to different REs. Accordingly, TFBs showed highly bias toward Short Interspersed Elements (SINEs), which in turn they covering about 12% of the total of REs. However, the distribution pattern of TFBs-REs showed different bias in the intronic regions, spatially into the Introns 13 and 25. The rich array of SINE-TFBs and CR1-TFBs were situated within 5′UTR of the gene that may be an important driving force for regulatory innovation of the hFVIII gene. Taken together, these data may lead to revealing intronic regions with the capacity to renewing gene regulatory networks of the hFVIII gene. On the other hand, these correlations might provide the novel idea for a new hypothesis of molecular evolution of the FVIII gene, and treatment of Hemophilia A which should be considered in future studies.

Genome-wide evolutionary characterization and expression analysis of SIAMESE-RELATED family genes in maize

Background

The SIAMESE (SIM) locus is a cell-cycle kinase inhibitor (CKI) gene that has to date been identified only in plants; it encodes a protein that promotes transformation from mitosis to endoreplication. Members of the SIAMESE-RELATED (SMR) family have similar functions, and some are related to cell-cycle responses and abiotic stresses. However, the functions of SMRs are poorly understood in maize (Zea mays L.).

Results

In the present study, 12 putative SMRs were identified throughout the entire genome of maize, and these were clustered into six groups together with the SMRs from seven other plant species. Members of the ZmSMR family were divided into four groups according to their protein sequences. Various cis-acting elements in the upstream sequences of ZmSMRs responded to abiotic stresses. Expression analyses revealed that all ZmSMRs were upregulated at 5, 20, 25, and 35 days after pollination. In addition, we found that ZmSMR9/11/12 may have regulated the initiation of endoreplication in endosperm central cells. Additionally, ZmSMR2/10 may have been primarily responsible for the endoreplication regulation of outer endosperm or aleurone cells. The relatively high expression levels of almost all ZmSMRs in the ears and tassels also implied that these genes may function in seed development. The effects of treatments with ABA, heat, cold, salt, and drought on maize seedlings and expression of ZmSMR genes suggested that ZmSMRs were strongly associated with response to abiotic stresses.

Conclusion

The present study is the first to conduct a genome-wide analysis of members of the ZmSMR family by investigating their locations in chromosomes, identifying regulatory elements in their promoter regions, and examining motifs in their protein sequences. Expression analysis of different endosperm developmental periods, tissues, abiotic stresses, and hormonal treatments suggests that ZmSMR genes may function in endoreplication and regulate the development of reproductive organs. These results may provide valuable information for future studies of the functions of the SMR family in maize.

Mapping and Quantification of Non-Coding RNA Originating from the rDNA in Human Glioma Cells

Ribosomal DNA is one of the most conserved parts of the genome, especially in its rRNA coding regions, but some puzzling pieces of its noncoding repetitive sequences harbor secrets of cell growth and development machinery. Disruptions in the neat mechanisms of rDNA orchestrating the cell functioning result in malignant conversion. In cancer cells, the organization of rRNA coding genes and their transcription somehow differ from that of normal cells, but little is known about the particular mechanism for this switch. In this study, we demonstrate that the region ~2 kb upstream of the rDNA promoter is transcriptionally active in one type of the most malignant human brain tumors, and we compare its expression rate to that of healthy human tissues and cell cultures. Sense and antisense non-coding RNA transcripts were detected and mapped, but their secondary structure and functions remain to be elucidated. We propose that the transcripts may relate to a new class of so-called promoter-associated RNAs (pRNAs), or have some other regulatory functions. We also hope that the expression of these non-coding RNAs can be used as a marker in glioma diagnostics and prognosis.

Control of the neuroprotective Lipocalin Apolipoprotein D expression by alternative promoter regions and differentially expressed mRNA 5' UTR variants

The Lipocalin Apolipoprotein D (ApoD) is one of the few genes consistently overexpressed in the aging brain, and in most neurodegenerative and psychiatric diseases. Its functions include metabolism regulation, myelin management, neuroprotection, and longevity regulation. Knowledge of endogenous regulatory mechanisms controlling brain disease-triggered ApoD expression is relevant if we want to boost pharmacologically its neuroprotecting potential. In addition to classical transcriptional control, Lipocalins have a remarkable variability in mRNA 5’UTR-dependent translation efficiency. Using bioinformatic analyses, we uncover strong selective pressures preserving ApoD 5’UTR properties, indicating unexpected functional conservation. PCR amplifications demonstrate the production of five 5’UTR variants (A-E) in mouse ApoD, with diverse expression levels across tissues and developmental stages. Importantly, Variant E is specifically expressed in the oxidative stress-challenged brain. Predictive analyses of 5’UTR secondary structures and enrichment in elements restraining translation, point to Variant E as a tight regulator of ApoD expression. We find two genomic regions conserved in human and mouse ApoD: a canonical (α) promoter region and a previously unknown region upstream of Variant E that could function as an alternative mouse promoter (β). Luciferase assays demonstrate that both α and β promoter regions can drive expression in cultured mouse astrocytes, and that Promoter β activity responds proportionally to incremental doses of the oxidative stress generator Paraquat. We postulate that Promoter β works in association with Variant E 5’UTR as a regulatory tandem that organizes ApoD gene expression in the nervous system in response to oxidative stress, the most common factor in aging and neurodegeneration.

Expression patterns of NbrgsCaM family genes in Nicotiana benthamiana and their potential roles in development and stress responses

rgsCaM has been reported as a calmodulin-like (CML) factor induced by viral infection in Nicotiana. There are three CMLs that belong to the rgsCaM family in Arabidopsis thaliana. In this study, we found a total of 5 NbrgsCaM coding sequences in N. benthamiana genome. We analyzed transcription patterns of NbrgsCaMs in transgenic plants expressing a β-glucuronidase (GUS) under the promoter of NbrgsCaMs by histochemistry staining and RT-qPCR. Similar to their Arabidopsis homologs, most NbrgsCaMs have an overlapping but distinct expression pattern in response to developmental and environmental changes. Specifically, the NbrgsCaM4 promoter exhibited robust activity and showed distinct regulatory response to viral infection, developmental stages and other abiotic stimuli. Overall, these findings provide clues for further understanding of the NbrgsCaM family genes in regulating plant growth and development under biotic stress and environmental stimulation.

Thapsigargin promotes colorectal cancer cell migration through upregulation of lncRNA MALAT1

Colorectal cancer (CRC) is the third most common tumor in the world; however, the role and mechanism of endoplasmic reticulum (ER) stress in CRC metastasis remains largely unclear. Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is a long non-coding RNA (lncRNA), which has previously been associated with CRC metastasis. It has been suggested that ER stress pathways regulate lncRNA expression; however, the effect of ER stress on MALAT1 expression in cancer is unknown. The present study aimed to investigate the relationship between ER stress pathways, MALAT1 expression and cell migration in CRC cells. ER stress was induced by thapsigargin (TG); low dose TG induced the migration of HT29 and HCT116 cells, but not SW1116 and SW620 cells. This effect was associated with increased expression levels of MALAT1, as the knockdown of MALAT1 prevented TG-induced cell migration. TG-induced MALAT1 expression was associated with inositol-requiring enzyme 1 (IRE1) expression and activation of the protein kinase R (PKR)-like ER kinase (PERK) signaling pathway. X-box-binding protein 1 (XBP1) and activating transcription factor 4 (ATF4) binding sites were predicted to be located in the MALAT1 gene promoter regions and the expression of MALAT1 was positively associated with XBP1 and ATF4 expression levels in CRC tissue samples. Thus, these findings indicated that ER stress may promote the migration of CRC cells and contribute to the progression of CRC through the activation of the IRE1/XBP1 and PERK/eIF2α/ATF4 signaling pathways. In conclusion, to the best of our knowledge, this study is the first report that lncRNA MALAT1 expression is regulated by the IRE1/XBP1 pathway in CRC.

Arabidopsis MYB4 plays dual roles in flavonoid biosynthesis

Flavonoids are major secondary metabolites derived from the plant phenylpropanoid pathway that play important roles in plant development and also have benefits for human health. So-called MBW ternary complexes involving R2R3-MYB and basic helix-loop-helix (bHLH) transcription factors along with WD-repeat proteins have been reported to regulate expression of the biosynthetic genes in the flavonoid pathway. MYB4 and its closest homolog MYB7 have been suggested to function as repressors of phenylpropanoid metabolism. However, the detailed mechanism by which they act has not been fully elucidated. Here, we show that Arabidopsis thaliana MYB4 and its homologs MYB7 and MYB32 interact with the bHLH transcription factors TT8, GL3 and EGL3 and thereby interfere with the transcriptional activity of the MBW complexes. In addition, MYB4 can also inhibit flavonoid accumulation by repressing expression of the gene encoding Arogenate Dehydratase 6 (ADT6), which catalyzes the final step in the biosynthesis of phenylalanine, the precursor for flavonoid biosynthesis. MYB4 potentially represses not only the conventional ADT6 encoding the plastidial enzyme but also the alternative isoform encoding the cytosolic enzyme. We suggest that MYB4 plays dual roles in modulating the flavonoid biosynthetic pathway in Arabidopsis.