The MEME Suite

Comparative and Expression Analysis of Ubiquitin Conjugating Domain-Containing Genes in Two Pyrus Species

Ripening affects the nutritional contents and quality of fleshy fruits, and it plays an important role during the process of fruit development. Studies have demonstrated that ubiquitin-conjugating (UBC or E2) genes can regulate fruit ripening, but the characterization of UBCs in pear is not well documented. The recently published genome-wide sequences of Pyrus bretschneideri and Pyrus communis have allowed a comprehensive analysis of this important gene family in pear. Using bioinformatics approaches, we identified 83 (PbrUBCs) and 84 (PcpUBCs) genes from P. bretschneideri and P. communis, respectively, which were divided into 13 subfamilies. In total, 198 PbrUBC paralogous, 215 PcpUBC paralogous, and 129 orthologous gene pairs were detected. Some paralogous gene pairs were found to be distributed on the same chromosome, suggesting that these paralogs may be caused by tandem duplications. The expression patterns of most UBC genes were divergent between Pyrus bretschneideri and Pyrus communis during pear fruit development. Remarkably, the transcriptome data showed that UBC genes might play a more important role in fruit ripening for further study. This is the first report on the systematic analysis of two Pyrus UBC gene families, and these data will help further study the role of UBC genes in fruit development and ripening, as well as contribute to the functional verification of UBC genes in pear.

Discovering Patterns From Sequences Using Pattern-Directed Aligned Pattern Clustering

Functional region identification is of fundamental importance for protein sequences analysis. Such knowledge provides better scientific understanding and could assist drug discovery. Up-to-date, domain annotation is one approach, but it needs to leverage existing databases. For de novo discovery, motif discovery locates and aligns locally homologous sub-sequences to obtain a position-weight matrix (PWM), which is a fixed-length representation model, whereas protein functional region size varies. It thus requires computational expensive exhaustive search to obtain a PWM with width of optimal range. This paper presents a new method known as pattern-directed aligned pattern clustering (PD-APCn) to discover and align patterns in conserved protein functional regions. It adopts aligned pattern cluster (APC) with patterns of variable length and strong support to direct the incremental APC expansion. It allows substitution and frame-shift mutations until a robust termination condition is reached. The concept of breakpoint gap is introduced to identify spots of mutations, such as substitution and frame shifts. Experiments on synthetic data sets with different sizes and noise levels showed that PD-APCn outperforms MEME with much higher recall and Fmeasure and computational speed 665 times faster that MEME. When applying to Cytochrome C and Ubiquitin families, it found all key binding sites within the APCs.

A network of epigenomic and transcriptional cooperation encompassing an epigenomic master regulator in cancer

Coordinated experiments focused on transcriptional responses and chromatin states are well-equipped to capture different epigenomic and transcriptomic levels governing the circuitry of a regulatory network. We propose a workflow for the genome-wide identification of epigenomic and transcriptional cooperation to elucidate transcriptional networks in cancer. Gene promoter annotation in combination with network analysis and sequence-resolution of enriched transcriptional motifs in epigenomic data reveals transcription factor families that act synergistically with epigenomic master regulators. By investigating complementary omics levels, a close teamwork of the transcriptional and epigenomic machinery was discovered. The discovered network is tightly connected and surrounds the histone lysine demethylase KDM3A, basic helix-loop-helix factors MYC, HIF1A, and SREBF1, as well as differentiation factors AP1, MYOD1, SP1, MEIS1, ZEB1, and ELK1. In such a cooperative network, one component opens the chromatin, another one recognizes gene-specific DNA motifs, others scaffold between histones, cofactors, and the transcriptional complex. In cancer, due to the ability to team up with transcription factors, epigenetic factors concert mitogenic and metabolic gene networks, claiming the role of a cancer master regulators or epioncogenes. Significantly, specific histone modification patterns are commonly associated with open or closed chromatin states, and are linked to distinct biological outcomes by transcriptional activation or repression. Disruption of patterns of histone modifications is associated with the loss of proliferative control and cancer. There is tremendous therapeutic potential in understanding and targeting histone modification pathways. Thus, investigating cooperation of chromatin remodelers and the transcriptional machinery is not only important for elucidating fundamental mechanisms of chromatin regulation, but also necessary for the design of targeted therapeutics.

A Large-Scale Mutational Analysis of Two-Component Signaling Systems of Lonsdalea quercina Revealed that KdpD-KdpE Regulates Bacterial Virulence Against Host Poplar Trees

Poplar, which is a dominant species in plant communities distributed in the northern hemisphere, is commonly used as a model plant in forestry studies. Poplar production can be inhibited by infections caused by bacteria, including Lonsdalea quercina subsp. populi, which is a gram-negative bacterium responsible for bark canker disease. However, the molecular basis of the pathogenesis remains uncharacterized. In this study, we annotated the two-component signal transduction systems (TCSs) encoded by the L. quercina subsp. populi N-5-1 genome and identified 18 putative histidine kinases and 24 response regulators. A large-scale mutational analysis revealed that 19 TCS genes regulated bacterial virulence against poplar trees. Additionally, the deletion of kdpE or overexpression of kdpD resulted in almost complete loss of bacterial virulence. We observed that kdpE and kdpD formed a bi-cistronic operon. KdpD exhibited autokinase activity and could bind to KdpE (K_d = 5.73 ± 0.64 μM). Furthermore, KdpE is an OmpR family response regulator. A chromatin immunoprecipitation sequencing analysis revealed that KdpE binds to an imperfect palindromic sequence within the promoters of 44 genes, including stress response genes Lqp0434, Lqp3037, and Lqp3270. A comprehensive analysis of TCS functions may help to characterize the regulation of poplar bark canker disease.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

A Zeb2-miR-200c loop controls midbrain dopaminergic neuron neurogenesis and migration

Zeb2 is a homeodomain transcription factor that plays pleiotropic functions during embryogenesis, but its role for midbrain dopaminergic (mDA) neuron development is unknown. Here we report that Zeb2 is highly expressed in progenitor cells in the ventricular zone of the midbrain floor plate and downregulated in postmitotic neuroblasts. Functional experiments show that Zeb2 expression in the embryonic ventral midbrain is dynamically regulated by a negative feedback loop that involves miR-200c. We also find that Zeb2 overexpression reduces the levels of CXCR4, NR4A2, and PITX3 in the developing ventral midbrain in vivo, resulting in migration and mDA differentiation defects. This phenotype was recapitulated by miR-200c knockdown, suggesting that the Zeb2-miR-200c loop prevents the premature differentiation of mDA progenitors into postmitotic cells and their migration. Together, our study establishes Zeb2 and miR-200c as critical regulators that maintain the balance between mDA progenitor proliferation and neurogenesis.

Comprehensive analysis of mitogen-activated protein kinase cascades in chrysanthemum

Background

Mitogen-activated protein kinase (MAPK) cascades, an important type of pathway in eukaryotic signaling networks, play a key role in plant defense responses, growth and development.

Methods

Phylogenetic analysis and conserved motif analysis of the MKK and MPK families in Arabidopsis thaliana, Helianthus annuus and Chrysanthemum morifolium classified MKK genes and MPK genes. qRT-PCR was used for the expression patterns of CmMPK and CmMKK genes, and yeast two-hybrid assay was applied to clear the interaction between CmMPKs and CmMKKs.

Results

We characterized six MKK genes and 11 MPK genes in chrysanthemum based on transcriptomic sequences and classified these genes into four groups. qRT-PCR analysis demonstrated that CmMKKs and CmMPKs exhibited various expression patterns in different organs of chrysanthemum and in response to abiotic stresses and phytohormone treatments. Furthermore, a yeast two-hybrid assay was applied to analyze the interaction between CmMKKs and CmMPKs and reveal the MAPK cascades in chrysanthemum.

Discussion

Our data led us to propose that CmMKK4-CmMPK13 and CmMKK2-CmMPK4 may be involved in regulating salt resistance and in the relationship between CmMKK9 and CmMPK6 and temperature stress.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male–female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Sex-specific meiotic crossover (CO) landscapes have been identified in multiple species. Here, the authors show that male and female meioses in maize have similar CO landscapes, and differences between COs in the two sexes only exists in their location relative to transcription start sites and some chromatin marks.

Novel Insights from Comparative In Silico Analysis of Green Microalgal Cellulases

The assumption that cellulose degradation and assimilation can only be carried out by heterotrophic organisms was shattered in 2012 when it was discovered that the unicellular green alga, Chlamydomonas reinhardtii (Cr), can utilize cellulose for growth under CO₂-limiting conditions. Publications of genomes/transcriptomes of the colonial microalgae, Gonium pectorale (Gp) and Volvox carteri (Vc), between 2010⁻2016 prompted us to look for cellulase genes in these algae and to compare them to cellulases from bacteria, fungi, lower/higher plants, and invertebrate metazoans. Interestingly, algal catalytic domains (CDs), belonging to the family GH9, clustered separately and showed the highest (33⁻42%) and lowest (17⁻36%) sequence identity with respect to cellulases from invertebrate metazoans and bacteria, respectively, whereas the identity with cellulases from plants was only 27⁻33%. Based on comparative multiple alignments and homology models, the domain arrangement and active-site architecture of algal cellulases are described in detail. It was found that all algal cellulases are modular, consisting of putative novel cysteine-rich carbohydrate-binding modules (CBMs) and proline/serine-(PS) rich linkers. Two genes were found to encode a protein with a putative Ig-like domain and a cellulase with an unknown domain, respectively. A feature observed in one cellulase homolog from Gp and shared by a spinach cellulase is the existence of two CDs separated by linkers and with a C-terminal CBM. Dockerin and Fn-3-like domains, typically found in bacterial cellulases, are absent in algal enzymes. The targeted gene expression analysis shows that two Gp cellulases consisting, respectively, of a single and two CDs were upregulated upon filter paper addition to the medium.

Asymmetric distribution of biomolecules of maternal origin in the Xenopus laevis egg and their impact on the developmental plan

Asymmetric cell division is a ubiquitous feature during the development of higher organisms. Asymmetry is achieved by differential localization or activities of biological molecules such as proteins, and coding and non-coding RNAs. Here, we present subcellular transcriptomic and proteomic analyses along the animal-vegetal axis of Xenopus laevis eggs. More than 98% of the maternal mRNAs could be categorized into four localization profile groups: animal, vegetal, extremely vegetal, and a newly described group of mRNAs that we call extremely animal, which are mRNAs enriched in the animal cortex region. 3′UTRs of localized mRNAs were analyzed for localization motifs. Several putative motifs were discovered for vegetal and extremely vegetal mRNAs, while no distinct conserved motifs for the extremely animal mRNAs were identified, suggesting different localization mechanisms. Asymmetric profiles were also found for proteins, with correlation to those of corresponding mRNAs. Based on unexpected observation of the profiles of the homoeologous genes exd2 we propose a possible mechanism of genetic evolution.

Putative bovine topological association domains and CTCF binding motifs can reduce the search space for causative regulatory variants of complex traits

BACKGROUND

Topological association domains (TADs) are chromosomal domains characterised by frequent internal DNA-DNA interactions. The transcription factor CTCF binds to conserved DNA sequence patterns called CTCF binding motifs to either prohibit or facilitate chromosomal interactions. TADs and CTCF binding motifs control gene expression, but they are not yet well defined in the bovine genome. In this paper, we sought to improve the annotation of bovine TADs and CTCF binding motifs, and assess whether the new annotation can reduce the search space for cis-regulatory variants.

RESULTS

We used genomic synteny to map TADs and CTCF binding motifs from humans, mice, dogs and macaques to the bovine genome. We found that our mapped TADs exhibited the same hallmark properties of those sourced from experimental data, such as housekeeping genes, transfer RNA genes, CTCF binding motifs, short interspersed elements, H3K4me3 and H3K27ac. We showed that runs of genes with the same pattern of allele-specific expression (ASE) (either favouring paternal or maternal allele) were often located in the same TAD or between the same conserved CTCF binding motifs. Analyses of variance showed that when averaged across all bovine tissues tested, TADs explained 14% of ASE variation (standard deviation, SD: 0.056), while CTCF explained 27% (SD: 0.078). Furthermore, we showed that the quantitative trait loci (QTLs) associated with gene expression variation (eQTLs) or ASE variation (aseQTLs), which were identified from mRNA transcripts from 141 lactating cows' white blood and milk cells, were highly enriched at putative bovine CTCF binding motifs. The linearly-furthermost, and most-significant aseQTL and eQTL for each genic target were located within the same TAD as the gene more often than expected (Chi-Squared test P-value < 0.001).

CONCLUSIONS

Our results suggest that genomic synteny can be used to functionally annotate conserved transcriptional components, and provides a tool to reduce the search space for causative regulatory variants in the bovine genome.

Evolutionary and genomic analysis of the caleosin/peroxygenase (CLO/PXG) gene/protein families in the Viridiplantae

Bioinformatics analyses of caleosin/peroxygenases (CLO/PXG) demonstrated that these genes are present in the vast majority of Viridiplantae taxa for which sequence data are available. Functionally active CLO/PXG proteins with roles in abiotic stress tolerance and lipid droplet storage are present in some Trebouxiophycean and Chlorophycean green algae but are absent from the small number of sequenced Prasinophyceaen genomes. CLO/PXG-like genes are expressed during dehydration stress in Charophyte algae, a sister clade of the land plants (Embryophyta). CLO/PXG-like sequences are also present in all of the >300 sequenced Embryophyte genomes, where some species contain as many as 10–12 genes that have arisen via selective gene duplication. Angiosperm genomes harbour at least one copy each of two distinct CLO/PX isoforms, termed H (high) and L (low), where H-forms contain an additional C-terminal motif of about 30–50 residues that is absent from L-forms. In contrast, species in other Viridiplantae taxa, including green algae, non-vascular plants, ferns and gymnosperms, contain only one (or occasionally both) of these isoforms per genome. Transcriptome and biochemical data show that CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression. CLO/PXG proteins can associate with cytosolic lipid droplets and/or bilayer membranes. Many of the analysed isoforms also have peroxygenase activity and are involved in oxylipin metabolism. The distribution of CLO/PXG-like genes is consistent with an origin >1 billion years ago in at least two of the earliest diverging groups of the Viridiplantae, namely the Chlorophyta and the Streptophyta, after the Viridiplantae had already diverged from other Archaeplastidal groups such as the Rhodophyta and Glaucophyta. While algal CLO/PXGs have roles in lipid packaging and stress responses, the Embryophyte proteins have a much wider spectrum of roles and may have been instrumental in the colonisation of terrestrial habitats and the subsequent diversification as the major land flora.

The Sucrose Synthase Gene Family in Chinese Pear (Pyrus bretschneideri Rehd.): Structure, Expression, and Evolution

Sucrose synthase (SS) is a key enzyme involved in sucrose metabolism that is critical in plant growth and development, and particularly quality of the fruit. Sucrose synthase gene families have been identified and characterized in plants various plants such as tobacco, grape, rice, and Arabidopsis. However, there is still lack of detailed information about sucrose synthase gene in pear. In the present study, we performed a systematic analysis of the pear (Pyrus bretschneideri Rehd.) genome and reported 30 sucrose synthase genes. Subsequently, gene structure, phylogenetic relationship, chromosomal localization, gene duplications, promoter regions, collinearity, RNA-Seq data and qRT-PCR were conducted on these sucrose synthase genes. The transcript analysis revealed that 10 PbSSs genes (30%) were especially expressed in pear fruit development. Additionally, qRT-PCR analysis verified the RNA-seq data and shown that PbSS30, PbSS24, and PbSS15 have a potential role in the pear fruit development stages. This study provides important insights into the evolution of sucrose synthase gene family in pear and will provide assistance for further investigation of sucrose synthase genes functions in the process of fruit development, fruit quality and resistance to environmental stresses.

Viromic Analysis of Wastewater Input to a River Catchment Reveals a Diverse Assemblage of RNA Viruses

Detection of viruses in the environment is heavily dependent on PCR-based approaches that require reference sequences for primer design. While this strategy can accurately detect known viruses, it will not find novel genotypes or emerging and invasive viral species. In this study, we investigated the use of viromics, i.e., high-throughput sequencing of the biosphere's viral fraction, to detect human-/animal-pathogenic RNA viruses in the Conwy river catchment area in Wales, United Kingdom. Using a combination of filtering and nuclease treatment, we extracted the viral fraction from wastewater and estuarine river water and sediment, followed by high-throughput RNA sequencing (RNA-Seq) analysis on the Illumina HiSeq platform, for the discovery of RNA virus genomes. We found a higher richness of RNA viruses in wastewater samples than in river water and sediment, and we assembled a complete norovirus genotype GI.2 genome from wastewater effluent, which was not contemporaneously detected by conventional reverse transcription-quantitative PCR (qRT-PCR). The simultaneous presence of diverse rotavirus signatures in wastewater indicated the potential for zoonotic infections in the area and suggested runoff from pig farms as a possible origin of these viruses. Our results show that viromics can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided appropriate and rigorous controls are included. IMPORTANCE Enteric viruses cause gastrointestinal illness and are commonly transmitted through the fecal-oral route. When wastewater is released into river systems, these viruses can contaminate the environment. Our results show that we can use viromics to find the range of potentially pathogenic viruses that are present in the environment and identify prevalent genotypes. The ultimate goal is to trace the fate of these pathogenic viruses from origin to the point where they are a threat to human health, informing reference-based detection methods and water quality management.

Viromic Analysis of Wastewater Input to a River Catchment Reveals a Diverse Assemblage of RNA Viruses

Detection of viruses in the environment is heavily dependent on PCR-based approaches that require reference sequences for primer design. While this strategy can accurately detect known viruses, it will not find novel genotypes or emerging and invasive viral species. In this study, we investigated the use of viromics, i.e., high-throughput sequencing of the biosphere's viral fraction, to detect human-/animal-pathogenic RNA viruses in the Conwy river catchment area in Wales, United Kingdom. Using a combination of filtering and nuclease treatment, we extracted the viral fraction from wastewater and estuarine river water and sediment, followed by high-throughput RNA sequencing (RNA-Seq) analysis on the Illumina HiSeq platform, for the discovery of RNA virus genomes. We found a higher richness of RNA viruses in wastewater samples than in river water and sediment, and we assembled a complete norovirus genotype GI.2 genome from wastewater effluent, which was not contemporaneously detected by conventional reverse transcription-quantitative PCR (qRT-PCR). The simultaneous presence of diverse rotavirus signatures in wastewater indicated the potential for zoonotic infections in the area and suggested runoff from pig farms as a possible origin of these viruses. Our results show that viromics can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided appropriate and rigorous controls are included. IMPORTANCE Enteric viruses cause gastrointestinal illness and are commonly transmitted through the fecal-oral route. When wastewater is released into river systems, these viruses can contaminate the environment. Our results show that we can use viromics to find the range of potentially pathogenic viruses that are present in the environment and identify prevalent genotypes. The ultimate goal is to trace the fate of these pathogenic viruses from origin to the point where they are a threat to human health, informing reference-based detection methods and water quality management.

VRTN is Required for the Development of Thoracic Vertebrae in Mammals

Vertnin (VRTN) variants are associated with thoracic vertebral number (TVN) in pigs. However, the biological function of VRTN remains poorly understood. Here we first conducted a range of experiments to demonstrate that VRTN is a responsible gene for TVN and two causative variants in the regulatory region of VRTN additively regulate TVN. Then, we show that VRTN is a novel DNA-binding transcription factor as it localizes exclusively in the nucleus, binds to DNA on a genome-wide scale and regulates the transcription of a set of genes that harbor VRTN binding motifs. Next, we illustrate that VRTN is essential for the development of thoracic vertebrae. Vrtn-null embryos display somitogenesis defect with the failure of axial rotation and fewer somites at the thoracic somite stage. Half of Vrtn heterozygous mice show abnormal spinal development with fewer thoracic vertebrae and ribs than their wild-type littermates. Lastly, we reveal that VRTN could modulate somite segmentation via the Notch signaling pathway. The findings advance our understanding of the mechanisms underlying the development of thoracic vertebrate in mammals, and VRTN causative variants provide a robust tool to improve pork production by selecting the alleles increasing the number of thoracic vertebrae and ribs.

Host-mediated RNAi of a Notch-like receptor gene in Meloidogyne incognita induces nematode resistance

GLP-1 (abnormal germline proliferation) is a Notch-like receptor protein that plays an essential role in pharyngeal development. In this study, an orthologue of Caenorhabditis elegans glp-1 was identified in Meloidogyne incognita. A computational analysis revealed that the orthologue contained almost all the domains present in the C. elegans gene: specifically, the LIN-12/Notch repeat, the ankyrin repeat, a transmembrane domain and different ligand-binding motifs were present in orthologue, but the epidermal growth factor-like motif was not observed. An expression analysis showed differential expression of glp-1 throughout the life cycle of M. incognita, with relatively higher expression in the egg stage. To evaluate the silencing efficacy of Mi-glp-1, transgenic Arabidopsis plants carrying double-stranded RNA constructs of glp-1 were generated, and infection of these plants with M. incognita resulted in a 47-50% reduction in the numbers of galls, females and egg masses. Females obtained from the transgenic RNAi lines exhibited 40-60% reductions in the transcript levels of the targeted glp-1 gene compared with females isolated from the control plants. Second-generation juveniles (J2s), which were descendants of the infected females from the transgenic lines, showed aberrant phenotypes. These J2s exhibited a significant decrease in the overall distance from the stylet to the metacorpus region, and this effect was accompanied by disruption around the metacorporeal bulb of the pharynx. The present study suggests a role for this gene in organ (pharynx) development during embryogenesis in M. incognita and its potential use as a target in the management of nematode infestations in plants.

70ProPred: a predictor for discovering sigma70 promoters based on combining multiple features

BACKGROUND

Promoter is an important sequence regulation element, which is in charge of gene transcription initiation. In prokaryotes, σ70 promoters regulate the transcription of most genes. The promoter recognition has been a crucial part of gene structure recognition. It's also the core issue of constructing gene transcriptional regulation network. With the successfully completion of genome sequencing from an increasing number of microbe species, the accurate identification of σ70 promoter regions in DNA sequence is not easy.

RESULTS

In order to improve the prediction accuracy of sigma70 promoters in prokaryote, a promoter recognition model 70ProPred was established. In this work, two sequence-based features, including position-specific trinucleotide propensity based on single-stranded characteristic (PSTNPss) and electron-ion potential values for trinucleotides (PseEIIP), were assessed to build the best prediction model. It was found that 79 features of PSTNPSS combined with 64 features of PseEIIP obtained the best performance for sigma70 promoter identification, with a promising accuracy and the Matthews correlation coefficient (MCC) at 95.56% and 0.90, respectively.

CONCLUSION

The jackknife tests showed that 70ProPred outperforms the existing sigma70 promoter prediction approaches in terms of accuracy and stability. Additionally, this approach can also be extended to predict promoters of other species. In order to facilitate experimental biologists, an online web server for the proposed method was established, which is freely available at http://server.malab.cn/70ProPred/ .

Molecular evolution and diversification of the SMXL gene family

Strigolactones (SLs) are a relatively recent addition to the list of plant hormones that control different aspects of plant development. SL signalling is perceived by an α/β hydrolase, DWARF 14 (D14). A close homolog of D14, KARRIKIN INSENSTIVE2 (KAI2), is involved in perception of an uncharacterized molecule called karrikin (KAR). Recent studies in Arabidopsis identified the SUPPRESSOR OF MAX2 1 (SMAX1) and SMAX1-LIKE 7 (SMXL7) to be potential SCF-MAX2 complex-mediated proteasome targets of KAI2 and D14, respectively. Genetic studies on SMXL7 and SMAX1 demonstrated distinct developmental roles for each, but very little is known about these repressors in terms of their sequence features. In this study, we performed an extensive comparative analysis of SMXLs and determined their phylogenetic and evolutionary history in the plant lineage. Our results show that SMXL family members can be sub-divided into four distinct phylogenetic clades/classes, with an ancient SMAX1. Further, we identified the clade-specific motifs that have evolved and that might act as determinants of SL-KAR signalling specificity. These specificities resulted from functional diversities among the clades. Our results suggest that a gradual co-evolution of SMXL members with their upstream receptors D14/KAI2 provided an increased specificity to both the SL perception and response in land plants.

Genome Wide Identification, Evolutionary, and Expression Analysis of VQ Genes from Two Pyrus Species

The VQ motif-containing gene, a member of the plant-specific genes, is involved in the plant developmental process and various stress responses. The VQ motif-containing gene family has been studied in several plants, such as rice (Oryza sativa), maize (Zea mays), and Arabidopsis (Arabidopsis thaliana). However, no systematic study has been performed in Pyrus species, which have important economic value. In our study, we identified 41 and 28 VQ motif-containing genes in Pyrus bretschneideri and Pyrus communis, respectively. Phylogenetic trees were calculated using A. thaliana and O. sativa VQ motif-containing genes as a template, allowing us to categorize these genes into nine subfamilies. Thirty-two and eight paralogous of VQ motif-containing genes were found in P. bretschneideri and P. communis, respectively, showing that the VQ motif-containing genes had a more remarkable expansion in P. bretschneideri than in P. communis. A total of 31 orthologous pairs were identified from the P. bretschneideri and P. communis VQ motif-containing genes. Additionally, among the paralogs, we found that these duplication gene pairs probably derived from segmental duplication/whole-genome duplication (WGD) events in the genomes of P. bretschneideri and P. communis, respectively. The gene expression profiles in both P. bretschneideri and P. communis fruits suggested functional redundancy for some orthologous gene pairs derived from a common ancestry, and sub-functionalization or neo-functionalization for some of them. Our study provided the first systematic evolutionary analysis of the VQ motif-containing genes in Pyrus, and highlighted the diversification and duplication of VQ motif-containing genes in both P. bretschneideri and P. communis.

A B-ARR-mediated cytokinin transcriptional network directs hormone cross-regulation and shoot development

Cytokinin fulfills its diverse roles in planta through a series of transcriptional responses. We identify the in vivo DNA binding site profiles for three genetically redundant type-B ARABIDOPSIS RESPONSE REGULATORS (B-ARRs): ARR1, ARR10, and ARR12. The expression and genome-wide DNA binding locations of the three B-ARRs extensively overlap. Constructing a primary cytokinin response transcriptional network reveals a recurring theme of widespread cross-regulation between the components of the cytokinin pathway and other plant hormone pathways. The B-ARRs are found to have similar DNA binding motifs, though sequences flanking the core motif were degenerate. Cytokinin treatments amalgamate the three different B-ARRs motifs to identical DNA binding signatures (AGATHY, H(a/t/c), Y(t/c)) which suggests cytokinin may regulate binding activity of B-ARR family members. Furthermore, we find that WUSCHEL, a key gene required for apical meristem maintenance, is a cytokinin-dependent B-ARR target gene, demonstrating the importance of the cytokinin transcription factor network in shoot development.

Systematic analysis and comparison of the PHD-Finger gene family in Chinese pear (Pyrus bretschneideri) and its role in fruit development

PHD-finger proteins, which belongs to the type of zinc finger family, and that play an important role in the regulation of both transcription and the chromatin state in eukaryotes. Currently, PHD-finger proteins have been well studied in animals, while few studies have been carried out on their function in plants. In the present study, 129 non-redundant PHD-finger genes were identified from 5 Rosaceae species (pear, apple, strawberry, mei, and peach); among them, 31 genes were identified in pear. Subsequently, we carried out a bioinformatics analysis of the PHD-finger genes. Thirty-one PbPHD genes were divided into 7 subfamilies based on the phylogenetic analysis, which are consistent with the intron-exon and conserved motif analyses. In addition, we identified five segmental duplication events, implying that the segmental duplications might be a crucial role in the expansion of the PHD-finger gene family in pear. The microsynteny analysis of five Rosaceae species showed that there were independent duplication events in addition to the genome-wide duplication of the pear genome. Subsequently, ten expressed PHD-finger genes of pear fruit were identified using qRT-PCR, and one of these genes, PbPHD10, was identified as an important candidate gene for the regulation of lignin synthesis. Our research provides useful information for the further analysis of the function of PHD-finger gene family in pear.

Gcn4 Binding in Coding Regions Can Activate Internal and Canonical 5' Promoters in Yeast

Gcn4 is a yeast transcriptional activator induced by amino acid starvation. ChIP-seq analysis revealed 546 genomic sites occupied by Gcn4 in starved cells, representing ∼30% of Gcn4-binding motifs. Surprisingly, only ∼40% of the bound sites are in promoters, of which only ∼60% activate transcription, indicating extensive negative control over Gcn4 function. Most of the remaining ∼300 Gcn4-bound sites are within coding sequences (CDSs), with ∼75 representing the only bound sites near Gcn4-induced genes. Many such unconventional sites map between divergent antisense and sub-genic sense transcripts induced within CDSs adjacent to induced TBP peaks, consistent with Gcn4 activation of cryptic bidirectional internal promoters. Mutational analysis confirms that Gcn4 sites within CDSs can activate sub-genic and full-length transcripts from the same or adjacent genes, showing that functional Gcn4 binding is not confined to promoters. Our results show that internal promoters can be regulated by an activator that functions at conventional 5'-positioned promoters.

Interactions between genetic variation and cellular environment in skeletal muscle gene expression

From whole organisms to individual cells, responses to environmental conditions are influenced by genetic makeup, where the effect of genetic variation on a trait depends on the environmental context. RNA-sequencing quantifies gene expression as a molecular trait, and is capable of capturing both genetic and environmental effects. In this study, we explore opportunities of using allele-specific expression (ASE) to discover cis-acting genotype-environment interactions (GxE)—genetic effects on gene expression that depend on an environmental condition. Treating 17 common, clinical traits as approximations of the cellular environment of 267 skeletal muscle biopsies, we identify 10 candidate environmental response expression quantitative trait loci (reQTLs) across 6 traits (12 unique gene-environment trait pairs; 10% FDR per trait) including sex, systolic blood pressure, and low-density lipoprotein cholesterol. Although using ASE is in principle a promising approach to detect GxE effects, replication of such signals can be challenging as validation requires harmonization of environmental traits across cohorts and a sufficient sampling of heterozygotes for a transcribed SNP. Comprehensive discovery and replication will require large human transcriptome datasets, or the integration of multiple transcribed SNPs, coupled with standardized clinical phenotyping.

Massive GGAAs in genomic repetitive sequences serve as a nuclear reservoir of NF-κB

Nuclear factor κB (NF-κB) is a DNA-binding transcription factor. Characterizing its genomic binding sites is crucial for understanding its gene regulatory function and mechanism in cells. This study characterized the binding sites of NF-κB RelA/p65 in the tumor neurosis factor-α (TNFα) stimulated HeLa cells by a precise chromatin immunoprecipitation-sequencing (ChIP-seq). The results revealed that NF-κB binds nontraditional motifs (nt-motifs) containing conserved GGAA quadruplet. Moreover, nt-motifs mainly distribute in the peaks nearby centromeres that contain a larger number of repetitive elements such as satellite, simple repeats and short interspersed nuclear elements (SINEs). This intracellular binding pattern was then confirmed by the in vitro detection, indicating that NF-κB dimers can bind the nontraditional κB (nt-κB) sites with low affinity. However, this binding hardly activates transcription. This study thus deduced that NF-κB binding nt-motifs may realize functions other than gene regulation as NF-κB binding traditional motifs (t-motifs). To testify the deduction, many ChIP-seq data of other cell lines were then analyzed. The results indicate that NF-κB binding nt-motifs is also widely present in other cells. The ChIP-seq data analysis also revealed that nt-motifs more widely distribute in the peaks with low-fold enrichment. Importantly, it was also found that NF-κB binding nt-motifs is mainly present in the resting cells, whereas NF-κB binding t-motifs is mainly present in the stimulated cells. Astonishingly, no known function was enriched by the gene annotation of nt-motif peaks. Based on these results, this study proposed that the nt-κB sites that extensively distribute in larger numbers of repeat elements function as a nuclear reservoir of NF-κB. The nuclear NF-κB proteins stored at nt-κB sites in the resting cells may be recruited to the t-κB sites for regulating its target genes upon stimulation.

Global dynamics of stage-specific transcription factor binding during thymocyte development

In vertebrates, multiple transcription factors (TFs) bind to gene regulatory elements (promoters, enhancers, and silencers) to execute developmental expression changes. ChIP experiments are often used to identify where TFs bind to regulatory elements in the genome, but the requirement of TF-specific antibodies hampers analyses of tens of TFs at multiple loci. Here we tested whether TF binding predictions using ATAC-seq can be used to infer the identity of TFs that bind to functionally validated enhancers of the Cd4, Cd8, and Gata3 genes in thymocytes. We performed ATAC-seq at four distinct stages of development in mouse thymus, probing the chromatin accessibility landscape in double negative (DN), double positive (DP), CD4 single positive (SP4) and CD8 SP (SP8) thymocytes. Integration of chromatin accessibility with TF motifs genome-wide allowed us to infer stage-specific occupied TF binding sites within known and potentially novel regulatory elements. Our results provide genome-wide stage-specific T cell open chromatin profiles, and allow the identification of candidate TFs that drive thymocyte differentiation at each developmental stage.

Sex- and Tissue-Specific Expression Profiles of Odorant Binding Protein and Chemosensory Protein Genes in Bradysia odoriphaga (Diptera: Sciaridae)

Bradysia odoriphaga is an agricultural pest insect affecting the production of Chinese chive and other liliaceous vegetables in China, and it is significantly attracted by sex pheromones and the volatiles derived from host plants. Despite verification of this chemosensory behavior, however, it is still unknown how B. odoriphaga recognizes these volatile compounds on the molecular level. Many of odorant binding proteins (OBPs) and chemosensory proteins (CSPs) play crucial roles in olfactory perception. Here, we identified 49 OBP and 5 CSP genes from the antennae and body transcriptomes of female and male adults of B. odoriphaga, respectively. Sequence alignment and phylogenetic analysis among Dipteran OBPs and CSPs were analyzed. The sex- and tissue-specific expression profiles of 54 putative chemosensory genes among different tissues were investigated by quantitative real-time PCR (qRT-PCR). qRT-PCR analysis results suggested that 22 OBP and 3 CSP genes were enriched in the antennae, indicating they might be essential for detection of general odorants and pheromones. Among these antennae-enriched genes, nine OBPs (BodoOBP2/4/6/8/12/13/20/28/33) were enriched in the male antennae and may play crucial roles in the detection of sex pheromones. Moreover, some OBP and CSP genes were enriched in non-antennae tissues, such as in the legs (BodoOBP3/9/19/21/34/35/38/39/45 and BodoCSP1), wings (BodoOBP17/30/32/37/44), abdomens and thoraxes (BodoOBP29/36), and heads (BodoOBP14/23/31 and BodoCSP2), suggesting that these genes might be involved in olfactory, gustatory, or other physiological processes. Our findings provide a starting point to facilitate functional research of these chemosensory genes in B. odoriphaga at the molecular level.

Structural and functional characterisation of two novel durum wheat annexin genes in response to abiotic stress

Abiotic stress results in massive loss of crop productivity throughout the world. Understanding the plant gene regulatory mechanisms involved in stress responses is very important. Annexins are a conserved multigene family of Ca-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signalling during plant growth and development. Annexins function to counteract oxidative stress, maintain cell redox homeostasis and enhance drought tolerance. A full-length cDNA of two genes (TdAnn6 and TdAnn12) encoding annexin proteins were isolated and characterised from Tunisian durum wheat varieties (Triticum turgidum L. subsp. durum cv. Mahmoudi). Analyses of the deduced proteins encoded by annexin cDNAs (TdAnn6 and TdAnn12) indicate the presence of the characteristic four repeats of 70-75 amino acids and the motifs proposed to be involved in Ca2+ binding. Gene expression patterns obtained by real-time PCR revealed differential temporal and spatial regulation of the two annexin genes in durum wheat under different abiotic stress conditions such as salt (NaCl 150mM), osmotic (10% polyethylene glycol 8000), ionic (LiCl 10mM), oxidative (H2O2), ABA (100µM), salicylic acid (10mM), cold (4°C) and heat (37°C) stress. The two annexin genes were not regulated by heavy metal stress (CdCl2 150µM). Moreover, heterologous expression of TdAnn6 and TdAnn12 in yeast improves its tolerance to abiotic stresses, suggesting annexin's involvement in theses stress tolerance mechanisms. Taken together, our results show that the two newly isolated wheat annexin might play an active role in modulating plant cell responses to abiotic stress responses.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

Bioinformatic characterization of the Anoctamin Superfamily of Ca2+-activated ion channels and lipid scramblases

Our laboratory has developed bioinformatic strategies for identifying distant phylogenetic relationships and characterizing families and superfamilies of transport proteins. Results using these tools suggest that the Anoctamin Superfamily of cation and anion channels, as well as lipid scramblases, includes three functionally characterized families: the Anoctamin (ANO), Transmembrane Channel (TMC) and Ca²⁺-permeable Stress-gated Cation Channel (CSC) families; as well as four families of functionally uncharacterized proteins, which we refer to as the Anoctamin-like (ANO-L), Transmembrane Channel-like (TMC-L), and CSC-like (CSC-L1 and CSC-L2) families. We have constructed protein clusters and trees showing the relative relationships among the seven families. Topological analyses suggest that the members of these families have essentially the same topologies. Comparative examination of these homologous families provides insight into possible mechanisms of action, indicates the currently recognized organismal distributions of these proteins, and suggests drug design potential for the disease-related channel proteins.

In Silico Selection Approach to Develop DNA Aptamers for a Stem-like Cell Subpopulation of Non-small Lung Cancer Adenocarcinoma Cell Line A549

Background

Detection of circulating lung cancer cells with cancer-stem like characteristics would represent an improved tool for disease prognosis. However, current antibodies based methods have some disadvantages and therefore cell SELEX (Systematic Evolution of Ligands by Exponential Enrichment) was used to develop DNA aptamers, recognizing cell surface markers of non-small lung carcinoma (NSLC) cells.

Materials and methods

The human adenocarcinoma cell line A549 was used for selection in seven cell SELEX cycles. We used human blood leukocytes for negative selection, and lung stem cell protein marker CD90 antibody binding A549 cells for positive selection.

Results

The obtained oligonucleotide sequences after the seventh SELEX cycle were subjected to in silico selection analysis based on three independent types of bioinformatics approaches, selecting two closely related aptamer candidates in terms of consensus sequences, structural motifs, binding affinity (Kd) and stability (ΔG). We selected and identified the aptamer A155_18 with very good binding characteristics to A459 cells, selected for CD90 antibody binding. The calculated phylogenetic tree showed that aptamers A155_18 and the known A549 cell aptamer S6 have a close structural relationship. MEME sequence analysis showed that they share two unique motifs, not present in other sequences.

Conclusions

The novel aptamer A155_18 has strong binding affinity for A549 lung carcinoma cell line subpopulation that is expressing stem cell marker CD90, indicating a possible stemness, characteristic for the A459 line, or a subpopulation present within this cell line. This aptamer can be applied as diagnostic tool, identifying NSLC circulating cells.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Bivalent Chromatin Domains in Glioblastoma Reveal a Subtype-Specific Signature of Glioma Stem Cells

Glioblastoma multiforme (GBM) can be clustered by gene expression into four main subtypes associated with prognosis and survival, but enhancers and other gene-regulatory elements have not yet been identified in primary tumors. Here, we profiled six histone modifications and CTCF binding as well as gene expression in primary gliomas and identified chromatin states that define distinct regulatory elements across the tumor genome. Enhancers in mesenchymal and classical tumor subtypes drove gene expression associated with cell migration and invasion, whereas enhancers in proneural tumors controlled genes associated with a less aggressive phenotype in GBM. We identified bivalent domains marked by activating and repressive chromatin modifications. Interestingly, the gene interaction network from common (subtype-independent) bivalent domains was highly enriched for homeobox genes and transcription factors and dominated by SHH and Wnt signaling pathways. This subtype-independent signature of early neural development may be indicative of poised dedifferentiation capacity in glioblastoma and could provide potential targets for therapy.Significance: Enhancers and bivalent domains in glioblastoma are regulated in a subtype-specific manner that resembles gene regulation in glioma stem cells. Cancer Res; 78(10); 2463-74. ©2018 AACR.

The histone demethylase KDM3A regulates the transcriptional program of the androgen receptor in prostate cancer cells

The lysine demethylase 3A (KDM3A, JMJD1A or JHDM2A) controls transcriptional networks in a variety of biological processes such as spermatogenesis, metabolism, stem cell activity, and tumor progression. We matched transcriptomic and ChIP-Seq profiles to decipher a genome-wide regulatory network of epigenetic control by KDM3A in prostate cancer cells. ChIP-Seq experiments monitoring histone 3 lysine 9 (H3K9) methylation marks show global histone demethylation effects of KDM3A. Combined assessment of histone demethylation events and gene expression changes presented major transcriptional activation suggesting that distinct oncogenic regulators may synergize with the epigenetic patterns by KDM3A. Pathway enrichment analysis of cells with KDM3A knockdown prioritized androgen signaling indicating that KDM3A plays a key role in regulating androgen receptor activity. Matched ChIP-Seq and knockdown experiments of KDM3A in combination with ChIP-Seq of the androgen receptor resulted in a gain of H3K9 methylation marks around androgen receptor binding sites of selected transcriptional targets in androgen signaling including positive regulation of KRT19, NKX3-1, KLK3, NDRG1, MAF, CREB3L4, MYC, INPP4B, PTK2B, MAPK1, MAP2K1, IGF1, E2F1, HSP90AA1, HIF1A, and ACSL3. The cancer systems biology analysis of KDM3A-dependent genes identifies an epigenetic and transcriptional network in androgen response, hypoxia, glycolysis, and lipid metabolism. Genome-wide ChIP-Seq data highlights specific gene targets and the ability of epigenetic master regulators to control oncogenic pathways and cancer progression.

Genome-Wide Analysis Characterization and Evolution of SBP Genes in Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume

The SQUAMOSA promoter binding protein (SBP)-box proteins are plant-specific transcriptional factors in plants. SBP TFs are known to play important functions in a diverse development process and also related in the process of evolutionary novelties. SBP gene family has been characterized in several plant species, but little is known about molecular evolution, functional divergence and comprehensive study of SBP gene family in Rosacea. We carried out genome-wide investigations and identified 14, 32, 17, and 17 SBP genes from four Rosacea species (Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume, respectively). According to phylogenetic analysis arranged the SBP protein sequences in seven groups. Localization of SBP genes presented an uneven distribution on corresponding chromosomes of Rosacea species. Our analyses designated that the SBP genes duplication events (segmental and tandem) and divergence. In addition, due to highly conserved structure pattern of SBP genes, recommended that highly conserved region of microsyneteny in the Rosacea species. Type I and II functional divergence was detected among various amino acids in SBP proteins, while there was no positive selection according to substitutional model analysis using PMAL software. These results recommended that the purifying selection might be leading force during the evolution process and dominate conservation of SBP genes in Rosacea species according to environmental selection pressure analysis. Our results will provide basic understanding and foundation for future research insights on the evolution of the SBP genes in Rosacea.

MarR-Family Transcription Factor HpaR Controls Expression of the vgrR-vgrS Operon of Xanthomonas campestris pv. campestris

MarR (multiple antibiotic resistance regulator)-family transcription factors (TFs), which regulate the expression of virulence factors and other physiological pathways in pathogenic bacteria, are regarded as ideal molecular targets for the development of novel antimicrobial strategies. In the plant bacterial pathogen Xanthomonas campestris pv. campestris, HpaR, a typical MarR-family TF, is associated with bacterial virulence, but its mechanism of virulence regulation remains unclear. Here, we dissected the HpaR regulon using high-throughput RNA sequencing and chromatin immunoprecipitation sequencing. HpaR directly or indirectly controls the expression of approximately 448 genes; it acts both as a transcriptional activator and a repressor to control the expression of downstream genes by directly binding to their promoter regions. The consensus HpaR-binding DNA motifs contain imperfect palindromic sequences similar to [G/T]CAACAATT[C/T]TTG. In-depth analysis revealed that HpaR positively modulates transcription level of the vgrR-vgrS operon that encodes an important two-component signal transduction system to sense iron depletion and regulate bacterial virulence. Epistasis analysis demonstrated that vgrR-vgrS is a core downstream component of HpaR regulation, as overexpression of vgrR restored the phenotypic deficiencies caused by a hpaR mutation. This dissection of the HpaR regulon should facilitate future studies focused on the activating mechanism of HpaR during bacterial infection.

Whole cell-SELEX of aptamers with a tyrosine-like side chain against live bacteria

In an effort to expand the binding and recognition capabilities of aptamers, a nucleoside triphosphate modified with a phenol that mimics the side chain of tyrosine was used in the selection of DNA aptamers against live bacteria. Of multiple modified aptamers that were isolated against Escherichia coli DH5α cells, one aptamer displays high selectivity and affinity for the target cells and is greatly enriched for phenol-modified dU nucleotides (dU^y, 47.5%). When the same sequences are synthesized with TTP, no binding is observed. Taken together, these findings highlight the value of using modified nucleotide triphosphates in aptamer selections and portends success in SELEX against an array of whole cells as targets.

Allele-Specific Chromatin Recruitment and Therapeutic Vulnerabilities of ESR1 Activating Mutations

Estrogen receptor α (ER) ligand-binding domain (LBD) mutations are found in a substantial number of endocrine treatment-resistant metastatic ER-positive (ER⁺) breast cancers. We investigated the chromatin recruitment, transcriptional network, and genetic vulnerabilities in breast cancer models harboring the clinically relevant ER mutations. These mutants exhibit both ligand-independent functions that mimic estradiol-bound wild-type ER as well as allele-specific neomorphic properties that promote a pro-metastatic phenotype. Analysis of the genome-wide ER binding sites identified mutant ER unique recruitment mediating the allele-specific transcriptional program. Genetic screens identified genes that are essential for the ligand-independent growth driven by the mutants. These studies provide insights into the mechanism of endocrine therapy resistance engendered by ER mutations and potential therapeutic targets.

Identification of AUXIN RESPONSE FACTOR gene family from Prunus sibirica and its expression analysis during mesocarp and kernel development

BACKGROUND

Auxin response factors (ARFs) in auxin signaling pathway are an important component that can regulate the transcription of auxin-responsive genes involved in almost all aspects of plant growth and development. To our knowledge, the comprehensive and systematic characterization of ARF genes has never been reported in Prunus sibirica, a novel woody biodiesel feedstock in China.

RESULTS

In this study, we identified 14 PsARF genes with a perfect open reading frame (ORF) in P. sibirica by using its previous transcriptomic data. Conserved motif analysis showed that all identified PsARF proteins had typical DNA-binding and ARF domain, but 5 members (PsARF3, 8 10, 16 and 17) lacked the dimerization domain. Phylogenetic analysis of the ARF proteins generated from various plant species indicated that ARFs could be categorized into 4 major groups (Class I, II, III and IV), in which all identified ARFs from P. sibirica showed a closest relationship with those from P. mume. Comparison of the expression profiles of 14 PsARF genes in different developmental stages of Siberian apricot mesocarp (SAM) and kernel (SAK) reflected distinct temporal or spatial expression patterns for PsARF genes. Additionally, based on the expressed data from fruit and seed development of multiple plant species, we identified 1514 ARF-correlated genes using weighted gene co-expression network analysis (WGCNA). And the major portion of ARF-correlated gene was characterized to be involved in protein, nucleic acid and carbohydrate metabolic, transport and regulatory processes.

CONCLUSIONS

In summary, we systematically and comprehensively analyzed the structure, expression pattern and co-expression network of ARF gene family in P. sibirica. All our findings provide theoretical foundation for the PsARF gene family and will pave the way for elucidating the precise role of PsARF genes in SAM and SAK development.

Genome-Wide Methylation Patterns in Androgen-Independent Prostate Cancer Cells: A Comprehensive Analysis Combining MeDIP-Bisulfite, RNA, and microRNA Sequencing Data

This study aimed to investigate the mechanisms underlying the development of the androgen-independent phenotype in prostate cancer. Methylation patterns were detected in androgen-independent and androgen-dependent lymph node carcinoma of the prostate (LNCaP) prostate carcinoma cells based on methylated DNA immunoprecipitation-bisulfite sequencing data and differentially methylated regions (DMRs) were identified. Differentially expressed genes (DEGs) and micro RNAs (miRNAs) with DMRs (named MDEGs and MDEmiRNAs) were identified by combining transcriptome and methylation data, and transcription factor (TF)-DEGs with DMRs in promoter (PMDEGs) and MDEmiRNA-MDEGs networks were constructed. Furthermore, a time-course analysis of gene transcription during androgen deprivation was performed based on microarray data and DMRs, MDEGs, and DEmiRNAs were validated. In total, 18,447 DMRs, 3369 MDEGs, 850 PMDEGs, and 1 MDEmiRNA (miR-429) were identified. A TF-target network (94 PMDEGs and 5 TFs) and a miRNA–target network (172 MDEGs and miR-429) were constructed. Based on the time-course analysis of genes in the networks, NEDD4L and PBX3 were targeted by SOX5, while GNAQ, ANLN, and KIF11 were targeted by miR-429. The expression levels of these genes and miR-429 were confirmed by quantitative real-time polymerase chain reaction. Additionally, 109 DMRs were confirmed using additional public datasets. The regulatory pathways SOX5-NEDD4L/PBX3, miR429-GNAQ/ANLN—RHOA, and miR429-ANLN—KIF11 may participate in the progression of the androgen-independent phenotype in prostate cancer.

Methods to Study Solo/Orphan Quorum-Sensing Receptors

LuxR solos/orphans are very widespread among Proteobacteria; however they are surprisingly understudied given that they are likely to play a major role in cell-cell communication in bacteria. Here we describe three simple methodologies/approaches that can be used in order to begin to study this subgroup of quorum sensing-related LuxR receptors.

Alternative polyadenylation is involved in auxin-based plant growth and development

Auxin is widely involved in plant growth and development. However, the molecular mechanism on how auxin carries out this work is unclear. In particular, the effect of auxin on pre-mRNA post-transcriptional regulation is mostly unknown. By using a poly(A) tag (PAT) sequencing approach, mRNA alternative polyadenylation (APA) profiles after auxin treatment were revealed. We showed that hundreds of poly(A) site clusters (PACs) are affected by auxin at the transcriptome level, where auxin reduces PAC distribution in 5'-untranslated region (UTR), but increases in the 3'UTR. APA site usage frequencies of 42 genes were switched by auxin, suggesting that auxin affects the choice of poly(A) sites. Furthermore, poly(A) signal selection was altered after auxin treatment. For example, a mutant of poly(A) signal binding protein CPSF30 showed altered sensitivity to auxin treatment, indicating interactions between auxin and the poly(A) signal recognition machinery. We also found that auxin activity on lateral root development is likely mediated by altered expression of ARF7, ARF19 and IAA14 through poly(A) site switches. Our results shed light on the molecular mechanisms of auxin responses relative to its interactions with mRNA polyadenylation.

Integrated omics analysis of root-preferred genes across diverse rice varieties including Japonica and indica cultivars

Plant root systems play essential roles in developmental processes, such as the absorption of water and inorganic nutrients, and structural support. Gene expression is affected by growth conditions and the genetic background of plants. To identify highly conserved root-preferred genes in rice across diverse growth conditions and varieties, we used two independent meta-anatomical expression profiles based on a large collection of Affymetrix and Agilent 44K microarray data sets available for public use. We then identified 684 loci with root-preferred expression, which were validated with in silico analysis using both meta-expression profiles. The expression patterns of four candidate genes were confirmed in vivo by monitoring expression of β-glucuronidase under control of the candidate-gene promoters, providing new tools to manipulate agronomic traits associated with roots. We also utilized real-time PCR to examine the root-preferential expression of 14 genes across four rice varieties, including japonica and indica cultivars. Using a database of rice genes with known functions, we identified the reported functions of 39 out of the 684 candidate genes. Sixteen genes are directly involved in root development, while the remaining are involved in processes indirectly related to root development (i.e., soil-stress tolerance or growth retardation). This indicates the importance of our candidate genes for studies on root development and function. Gene ontology enrichment analysis in the 'biological processes' category revealed that root-preferred genes in rice are closely associated with nutrient transport-related genes, indicating that the primary role of roots is the uptake of nutrients from soil. In addition, predicted protein-protein interaction analysis suggested a molecular network for root development composed of 215 interactions associated with 44 root-preferred or root development-related genes. Taken together, our data provide an important foundation for future research on root development in rice.

Considerations on Experimental Design and Data Analysis of Chromatin Immunoprecipitation Experiments

Arguably one of the most valuable techniques to study chromatin organization, ChIP is the method of choice to map the contacts established between proteins and genomic DNA. Ever since its inception, more than 30 years ago, ChIP has been constantly evolving, improving, and expanding its capabilities and reach. Despite its widespread use by many laboratories across a wide variety of disciplines, ChIP assays can be sometimes challenging to design, and are often sensitive to variations in practical implementation.In this chapter, we provide a general overview of the ChIP method and its most common variations, with a special focus on ChIP-seq. We try to address some of the most important aspects that need to be taken into account in order to design and perform experiments that generate the most reproducible, high-quality data. Some of the main topics covered include the use of properly characterized antibodies, alternatives to chromatin preparation, the need for proper controls, and some recommendations about ChIP-seq data analysis.

NAC Family Transcription Factors in Tobacco and Their Potential Role in Regulating Leaf Senescence

The NAC family is one of the largest families of plant-specific transcription factors (TFs) and NAC proteins play important regulatory roles in a variety of developmental and stress response processes in plants. Members of the NAC family TFs have been shown to be important regulators of leaf senescence in a number of plant species. Here we report the identification of the NAC family in tobacco (Nicotiana tabacum) and characterization of the potential role of some of the tobacco NAC TFs in regulating leaf senescence. A total of 154 NAC genes (NtNACs) were identified and clustered together with the Arabidopsis NAC family into fifteen groups (a-o). Transcriptome data analysis followed by qRT-PCR validation showed that the majority of the senescence-up-regulated NtNACs fall into subgroups NAC-b and f. A number of known senescence regulators from Arabidopsis also belong to these two subgroups. Among these senescence-up-regulated NtNACs, NtNAC080, a close homolog of AtNAP, is a positive regulator of leaf senescence. Overexpression of NtNAC080 caused early senescence in Arabidopsis leaves and NtNAC080 mutation induced by Cas9/gRNA in tobacco led to delayed leaf senescence.

TsNAC1 Is a Key Transcription Factor in Abiotic Stress Resistance and Growth

NAC proteins constitute one of the largest families of plant-specific transcription factors, and a number of these proteins participate in the regulation of plant development and responses to abiotic stress. T. HALOPHILA STRESS RELATED NAC1 (TsNAC1), cloned from the halophyte Thellungiella halophila, is a NAC transcription factor gene, and its overexpression can improve abiotic stress resistance, especially in salt stress tolerance, in both T. halophila and Arabidopsis (Arabidopsis thaliana) and retard the growth of these plants. In this study, the transcriptional activation activity of TsNAC1 and RD26 from Arabidopsis was compared with the target genes' promoter regions of TsNAC1 from T. halophila, and the results showed that the transcriptional activation activity of TsNAC1 was higher in tobacco (Nicotiana tabacum) and yeast. The target sequence of the promoter from the target genes also was identified, and TsNAC1 was shown to target the positive regulators of ion transportation, such as T. HALOPHILA H⁺-PPASE1, and the transcription factors MYB HYPOCOTYL ELONGATION-RELATED and HOMEOBOX12 In addition, TsNAC1 negatively regulates the expansion of cells, inhibits LIGHT-DEPENDENT SHORT HYPOCOTYLS1 and UDP-XYLOSYLTRANSFERASE2, and directly controls the expression of MULTICOPY SUPPRESSOR OF IRA14 Based on these results, we propose that TsNAC1 functions as an important upstream regulator of plant abiotic stress responses and vegetative growth.

Genome-wide identification and cadmium induced expression profiling of sulfate transporter (SULTR) genes in sorghum (Sorghum bicolor L.)

Sulfur is an essential element for all living organisms. Plants can convert inorganic sulfur into organic sulfur compounds by complex enzymatic steps. In this study, we conducted a genome-wide analysis of sulfate transporter genes (SULTRs) in the sorghum (Sorghum bicolor) genome and examined expression profiles of SbSULTR genes under 200 µM cadmium (Cd) exposure. As a result of sorghum genome analysis, 11 SULTR genes were identified, including SbSULTR1;1, SbSULTR1;2, SbSULTR1;3, SbSULTR2;1, SbSULTR2;2, SbSULTR3;1, SbSULTR3;2, SbSULTR3;3, SbSULTR3;4, SbSULTR3;5, and SbSULTR4. Given names are based on phylogeny and chromosomal locations. Except SbSULTR4, all SbSULTR proteins contained Sulfate_transp (PF00916), STAS (PF01740) domains and 12 trans-membrane domains. Phylogenetic analysis revealed that four major groups were identified such as SULTR1, 2, 3, and 4 groups and SULTR4 group was separated to other SULTR groups. In promotor sequences of SbSULTR genes, many diverse cis-acting elements were found mainly related with physiological processes such as light, stress and hormone responsiveness. The expression profiles of SbSULTR genes showed that SULTR1;2, 1;3, 3;3, and 3;5 genes up-regulated in root, while expression level of SULTR4 decreased under 200 µM Cd exposure. The predicted 3D structures of SULTR proteins showed some conformational changes, suggesting functional diversities of SbSULTRs. Finally, results of this study may contribute towards understanding SbSULTR genes and their regulations and roles in Cd stress in sorghum.

Kiwifruit R2R3-MYB transcription factors and contribution of the novel AcMYB75 to red kiwifruit anthocyanin biosynthesis

Red kiwifruit (Actinidia chinensis) is a popular fresh fruit with a high market value due to its unique color, caused by anthocyanin accumulation. The R2R3-MYB transcription factors (TFs) have important roles in plant development and anthocyanin metabolism. In this first comprehensive study of R2R3-MYBs in kiwifruit, a total of 93 R2R3-MYB genes, including five novel previously unannotated AcMYBs, were identified. Their phylogenic relationship, exon-intron structures, and conserved motifs were analyzed. Based on transcriptome data, 60 AcMYBs were expressed (FPKM > 1) across seven developmental stages of kiwifruit, revealing five expression patterns. One of the 5 newly identified R2R3 TFs, AcMYB75, showed an anthocyanin accumulation-linked expression pattern during fruit development. AcMYB75 localized to the nucleus and has an active transactivation domain, verifying it as a transcription factor. AcMYB75 protein specifically bound the promoter of the anthocyanin biosynthesis gene ANS in yeast one-hybrid system and in vivo. In 35 S:AcMYB75 Arabidopsis plants, anthocyanin significantly accumulated in leaves, and the expression of anthocyanin biosynthetic genes was greatly up-regulated. Together, these results suggest that AcMYB75 is involved in anthocyanin biosynthesis in kiwifruit. These findings will increase our understanding of AcMYBs involved in anthocyanin biosynthesis, and also benefit further functional characterization of R2R3-MYB genes in kiwifruit.