CoCoView - A Codon Conservation Viewer via Sequence Logos

Sequence logos are a simple way to display a set of aligned sequences, and they are useful to identify conserved patterns. Since their introduction, several tools have been developed for generating these representations at the single residue level (amino acids or nucleotides). We have developed a tool to build sequence logos of protein-coding sequences at the codon level, allowing more accurate analysis of coding-sequences as they represent synonymous and non-synonymous changes instead of showing only changes that imply on amino acid substitutions. We built CoCoView on top of the Logomaker Python API. It creates codon sequence logos from a multiple sequence alignment of protein-coding sequences. Some properties of the data and the generated logos can be controlled by the end-users, such as data redundancy, plot type and alphabet color.

• Split aligned sequences into codon positions;

• For each position compute codon frequency and information content;

• Use the computed information to plot the graphic.

Bioinformatic analyses to uncover genes involved in trehalose metabolism in the polyploid sugarcane

Trehalose-6-phosphate (T6P) is an intermediate of trehalose biosynthesis that plays an essential role in plant metabolism and development. Here, we comprehensively analyzed sequences from enzymes of trehalose metabolism in sugarcane, one of the main crops used for bioenergy production. We identified protein domains, phylogeny, and in silico expression levels for all classes of enzymes. However, post-translational modifications and residues involved in catalysis and substrate binding were analyzed only in trehalose-6-phosphate synthase (TPS) sequences. We retrieved 71 putative full-length TPS, 93 trehalose-6-phosphate phosphatase (TPP), and 3 trehalase (TRE) of sugarcane, showing all their conserved domains, respectively. Putative TPS (Classes I and II) and TPP sugarcane sequences were categorized into well-known groups reported in the literature. We measured the expression levels of the sequences from one sugarcane leaf transcriptomic dataset. Furthermore, TPS Class I has specific N-glycosylation sites inserted in conserved motifs and carries catalytic and binding residues in its TPS domain. Some of these residues are mutated in TPS Class II members, which implies loss of enzyme activity. Our approach retrieved many homo(eo)logous sequences for genes involved in trehalose metabolism, paving the way to discover the role of T6P signaling in sugarcane.

The experience of teaching introductory programming skills to bioscientists in Brazil

Computational biology has gained traction as an independent scientific discipline over the last years in South America. However, there is still a growing need for bioscientists, from different backgrounds, with different levels, to acquire programming skills, which could reduce the time from data to insights and bridge communication between life scientists and computer scientists. Python is a programming language extensively used in bioinformatics and data science, which is particularly suitable for beginners. Here, we describe the conception, organization, and implementation of the Brazilian Python Workshop for Biological Data. This workshop has been organized by graduate and undergraduate students and supported, mostly in administrative matters, by experienced faculty members since 2017. The workshop was conceived for teaching bioscientists, mainly students in Brazil, on how to program in a biological context. The goal of this article was to share our experience with the 2020 edition of the workshop in its virtual format due to the Coronavirus Disease 2019 (COVID-19) pandemic and to compare and contrast this year's experience with the previous in-person editions. We described a hands-on and live coding workshop model for teaching introductory Python programming. We also highlighted the adaptations made from in-person to online format in 2020, the participants' assessment of learning progression, and general workshop management. Lastly, we provided a summary and reflections from our personal experiences from the workshops of the last 4 years. Our takeaways included the benefits of the learning from learners' feedback (LLF) that allowed us to improve the workshop in real time, in the short, and likely in the long term. We concluded that the Brazilian Python Workshop for Biological Data is a highly effective workshop model for teaching a programming language that allows bioscientists to go beyond an initial exploration of programming skills for data analysis in the medium to long term.

A Reductionist Approach Using Primary and Metastatic Cell-Derived Extracellular Vesicles Reveals Hub Proteins Associated with Oral Cancer Prognosis

Oral squamous cell carcinoma (OSCC) has high mortality rates that are largely associated with lymph node metastasis. However, the molecular mechanisms that drive OSCC metastasis are unknown. Extracellular vesicles (EVs) are membrane-bound particles that play a role in intercellular communication and impact cancer development and progression. Thus, profiling EVs would be of great significance to decipher their role in OSCC metastasis. For that purpose, we used a reductionist approach to map the proteomic, miRNA, metabolomic, and lipidomic profiles of EVs derived from human primary tumor (SCC-9) cells and matched lymph node metastatic (LN1) cells. Distinct omics profiles were associated with the metastatic phenotype, including 670 proteins, 217 miRNAs, 26 metabolites, and 63 lipids differentially abundant between LN1 cell– and SCC-9 cell–derived EVs. A multi-omics integration identified 11 ‘hub proteins’ significantly decreased at the metastatic site compared with primary tumor–derived EVs. We confirmed the validity of these findings with analysis of data from multiple public databases and found that low abundance of seven ‘hub proteins’ in EVs from metastatic lymph nodes (ALDH7A1, CAD, CANT1, GOT1, MTHFD1, PYGB, and SARS) is correlated with reduced survival and tumor aggressiveness in patients with cancer. In summary, this multi-omics approach identified proteins transported by EVs that are associated with metastasis and which may potentially serve as prognostic markers in OSCC.

•

Proteomic, miRNA, metabolomic, and lipidomic profiles were mapped in oral cancer EVs.

•

The molecular profile of EVs was associated with the lymph node metastatic phenotype.

•

A multi-omics integrative analysis revealed 11 highly connected ‘hub proteins.’

•

‘Hub proteins’ from EVs are candidates as prognostic markers in oral cancer.

Overexpression of a Sugarcane BAHD Acyltransferase Alters Hydroxycinnamate Content in Maize Cell Wall

The purification of hydroxycinnamic acids [p-coumaric acid (pCA) and ferulic acid (FA)] from grass cell walls requires high-cost processes. Feedstocks with increased levels of one hydroxycinnamate in preference to the other are therefore highly desirable. We identified and conducted expression analysis for nine BAHD acyltransferase ScAts genes from sugarcane. The high conservation of AT10 proteins, together with their similar gene expression patterns, supported a similar role in distinct grasses. Overexpression of ScAT10 in maize resulted in up to 75% increase in total pCA content. Mild hydrolysis and derivatization followed by reductive cleavage (DFRC) analysis showed that pCA increase was restricted to the hemicellulosic portion of the cell wall. Furthermore, total FA content was reduced up to 88%, resulting in a 10-fold increase in the pCA/FA ratio. Thus, we functionally characterized a sugarcane gene involved in pCA content on hemicelluloses and generated a C4 plant that is promising for valorizing pCA production in biorefineries.

Shedding Light on the Dynamic Role of the "Target of Rapamycin" Kinase in the Fast-Growing C4 Species Setaria viridis, a Suitable Model for Biomass Crops

The Target of Rapamycin (TOR) kinase pathway integrates energy and nutrient availability into metabolism promoting growth in eukaryotes. The overall higher efficiency on nutrient use translated into faster growth rates in C₄ grass plants led to the investigation of differential transcriptional and metabolic responses to short-term chemical TOR complex (TORC) suppression in the model Setaria viridis. In addition to previously described responses to TORC inhibition (i.e., general growth arrest, translational repression, and primary metabolism reprogramming) in Arabidopsis thaliana (C₃), the magnitude of changes was smaller in S. viridis, particularly regarding nutrient use efficiency and C allocation and partitioning that promote biosynthetic growth. Besides photosynthetic differences, S. viridis and A. thaliana present several specificities that classify them into distinct lineages, which also contribute to the observed alterations mediated by TOR. Indeed, cell wall metabolism seems to be distinctly regulated according to each cell wall type, as synthesis of non-pectic polysaccharides were affected in S. viridis, whilst assembly and structure in A. thaliana. Our results indicate that the metabolic network needed to achieve faster growth seems to be less stringently controlled by TORC in S. viridis.

Analysis of NAC Domain Transcription Factor Genes of Tectona grandis L.f. Involved in Secondary Cell Wall Deposition

NAC proteins are one of the largest families of plant-specific transcription factors (TFs). They regulate diverse complex biological processes, including secondary xylem differentiation and wood formation. Recent genomic and transcriptomic studies of Tectona grandis L.f. (teak), one of the most valuable hardwood trees in the world, have allowed identification and analysis of developmental genes. In the present work, T. grandis NAC genes were identified and analyzed regarding to their evolution and expression profile during wood formation. We analyzed the recently published T. grandis genome, and identified 130 NAC proteins that are coded by 107 gene loci. These proteins were classified into 23 clades of the NAC family, together with Populus, Eucalyptus, and Arabidopsis. Data on transcript expression revealed specific temporal and spatial expression patterns for the majority of teak NAC genes. RT-PCR indicated expression of VND genes (Tg11g04450-VND2 and Tg15g08390-VND4) related to secondary cell wall formation in xylem vessels of 16-year-old juvenile trees. Our findings open a way to further understanding of NAC transcription factor genes in T. grandis wood biosynthesis, while they are potentially useful for future studies aiming to improve biomass and wood quality using biotechnological approaches.

Gene regulatory networks on transfer entropy (GRNTE): a novel approach to reconstruct gene regulatory interactions applied to a case study for the plant pathogen Phytophthora infestans

Background

The increasing amounts of genomics data have helped in the understanding of the molecular dynamics of complex systems such as plant and animal diseases. However, transcriptional regulation, although playing a central role in the decision-making process of cellular systems, is still poorly understood. In this study, we linked expression data with mathematical models to infer gene regulatory networks (GRN). We present a simple yet effective method to estimate transcription factors’ GRNs from transcriptional data.

Method

We defined interactions between pairs of genes (edges in the GRN) as the partial mutual information between these genes that takes into account time and possible lags in time from one gene in relation to another. We call this method Gene Regulatory Networks on Transfer Entropy (GRNTE) and it corresponds to Granger causality for Gaussian variables in an autoregressive model. To evaluate the reconstruction accuracy of our method, we generated several sub-networks from the GRN of the eukaryotic yeast model, Saccharomyces cerevisae. Then, we applied this method using experimental data of the plant pathogen Phytophthora infestans. We evaluated the transcriptional expression levels of 48 transcription factors of P. infestans during its interaction with one moderately resistant and one susceptible cultivar of yellow potato (Solanum tuberosum group Phureja), using RT-qPCR. With these data, we reconstructed the regulatory network of P. infestans during its interaction with these hosts.

Results

We first evaluated the performance of our method, based on the transfer entropy (GRNTE), on eukaryotic datasets from the GRNs of the yeast S. cerevisae. Results suggest that GRNTE is comparable with the state-of-the-art methods when the parameters for edge detection are properly tuned. In the case of P. infestans, most of the genes considered in this study, showed a significant change in expression from the onset of the interaction (0 h post inoculum - hpi) to the later time-points post inoculation. Hierarchical clustering of the expression data discriminated two distinct periods during the infection: from 12 to 36 hpi and from 48 to 72 hpi for both the moderately resistant and susceptible cultivars. These distinct periods could be associated with two phases of the life cycle of the pathogen when infecting the host plant: the biotrophic and necrotrophic phases.

Conclusions

Here we presented an algorithmic solution to the problem of network reconstruction in time series data. This analytical perspective makes use of the dynamic nature of time series data as it relates to intrinsically dynamic processes such as transcription regulation, were multiple elements of the cell (e.g., transcription factors) act simultaneously and change over time. We applied the algorithm to study the regulatory network of P. infestans during its interaction with two hosts which differ in their level of resistance to the pathogen. Although the gene expression analysis did not show differences between the two hosts, the results of the GRN analyses evidenced rewiring of the genes’ interactions according to the resistance level of the host. This suggests that different regulatory processes are activated in response to different environmental cues. Applications of our methodology showed that it could reliably predict where to place edges in the transcriptional networks and sub-networks. The experimental approach used here can help provide insights on the biological role of these interactions on complex processes such as pathogenicity. The code used is available at https://github.com/jccastrog/GRNTE under GNU general public license 3.0.

Electronic supplementary material

The online version of this article (10.1186/s12976-019-0103-7) contains supplementary material, which is available to authorized users.

Genome sequence of the H2-producing Clostridium beijerinckii strain Br21 isolated from a sugarcane vinasse treatment plant

We report on the nearly complete genome sequence of Clostridium beijerinckii strain Br21, formerly isolated from a sugarcarne vinasse wastewater treatment plant. The resulting genome is ca. 5.9 Mbp in length and resembles the size of previously published C. beijerinckii genomes. We annotated the genome sequence and predicted a total of 5323 genes. Strain Br21 has a genetic toolkit that allows it to exploit diverse sugars that are often found after lignocellulosic biomass pretreatment to yield products of commercial interest. Besides the whole set of genes encoding for enzymes underlying hydrogen production, the genome of the new strain includes genes that enable carbon sources conversion into butanol, ethanol, acetic acid, butyric acid, and the chemical block 1,3-propanediol, which is used to obtain polymers. Moreover, the genome of strain Br21 has a higher number of ORFs with predicted beta-glucosidase activity as compared to other C. beijerinckii strains described in the KEGG database. These characteristics make C. beijerinckii strain Br21 a remarkable candidate for direct use in biotechnological processes and attest that it is a potential biocatalyst supplier.

Gene Co-expression Network Reveals Potential New Genes Related to Sugarcane Bagasse Degradation in Trichoderma reesei RUT-30

The biomass-degrading fungus Trichoderma reesei has been considered a model for cellulose degradation, and it is the primary source of the industrial enzymatic cocktails used in second-generation (2G) ethanol production. However, although various studies and advances have been conducted to understand the cellulolytic system and the transcriptional regulation of T. reesei, the whole set of genes related to lignocellulose degradation has not been completely elucidated. In this study, we inferred a weighted gene co-expression network analysis based on the transcriptome dataset of the T. reesei RUT-C30 strain aiming to identify new target genes involved in sugarcane bagasse breakdown. In total, ~70% of all the differentially expressed genes were found in 28 highly connected gene modules. Several cellulases, sugar transporters, and hypothetical proteins coding genes upregulated in bagasse were grouped into the same modules. Among them, a single module contained the most representative core of cellulolytic enzymes (cellobiohydrolase, endoglucanase, β-glucosidase, and lytic polysaccharide monooxygenase). In addition, functional analysis using Gene Ontology (GO) revealed various classes of hydrolytic activity, cellulase activity, carbohydrate binding and cation:sugar symporter activity enriched in these modules. Several modules also showed GO enrichment for transcription factor activity, indicating the presence of transcriptional regulators along with the genes involved in cellulose breakdown and sugar transport as well as other genes encoding proteins with unknown functions. Highly connected genes (hubs) were also identified within each module, such as predicted transcription factors and genes encoding hypothetical proteins. In addition, various hubs contained at least one DNA binding site for the master activator Xyr1 according to our in silico analysis. The prediction of Xyr1 binding sites and the co-expression with genes encoding carbohydrate active enzymes and sugar transporters suggest a putative role of these hubs in bagasse cell wall deconstruction. Our results demonstrate a vast range of new promising targets that merit additional studies to improve the cellulolytic potential of T. reesei strains and to decrease the production costs of 2G ethanol.

Transcriptome and secretome analysis of Aspergillus fumigatus in the presence of sugarcane bagasse

BACKGROUND

Sugarcane bagasse has been proposed as a lignocellulosic residue for second-generation ethanol (2G) produced by breaking down biomass into fermentable sugars. The enzymatic cocktails for biomass degradation are mostly produced by fungi, but low cost and high efficiency can consolidate 2G technologies. A. fumigatus plays an important role in plant biomass degradation capabilities and recycling. To gain more insight into the divergence in gene expression during steam-exploded bagasse (SEB) breakdown, this study profiled the transcriptome of A. fumigatus by RNA sequencing to compare transcriptional profiles of A. fumigatus grown on media containing SEB or fructose as the sole carbon source. Secretome analysis was also performed using SDS-PAGE and LC-MS/MS.

RESULTS

The maximum activities of cellulases (0.032 U mL-1), endo-1,4-β--xylanase (10.82 U mL-1) and endo-1,3-β glucanases (0.77 U mL-1) showed that functional CAZymes (carbohydrate-active enzymes) were secreted in the SEB culture conditions. Correlations between transcriptome and secretome data identified several CAZymes in A. fumigatus. Particular attention was given to CAZymes related to lignocellulose degradation and sugar transporters. Genes encoding glycoside hydrolase classes commonly expressed during the breakdown of cellulose, such as GH-5, 6, 7, 43, 45, and hemicellulose, such as GH-2, 10, 11, 30, 43, were found to be highly expressed in SEB conditions. Lytic polysaccharide monooxygenases (LPMO) classified as auxiliary activity families AA9 (GH61), CE (1, 4, 8, 15, 16), PL (1, 3, 4, 20) and GT (1, 2, 4, 8, 20, 35, 48) were also differentially expressed in this condition. Similarly, the most important enzymes related to biomass degradation, including endoxylanases, xyloglucanases, β-xylosidases, LPMOs, α-arabinofuranosidases, cellobiohydrolases, endoglucanases and β-glucosidases, were also identified in the secretome.

CONCLUSIONS

This is the first report of a transcriptome and secretome experiment of Aspergillus fumigatus in the degradation of pretreated sugarcane bagasse. The results suggest that this strain employs important strategies for this complex degradation process. It was possible to identify a set of genes and proteins that might be applied in several biotechnology fields. This knowledge can be exploited for the improvement of 2G ethanol production by the rational design of enzymatic cocktails.

ploidyNGS: visually exploring ploidy with Next Generation Sequencing data

Summary

ploidyNGS is a model-free, open source tool to visualize and explore ploidy levels in a newly sequenced genome, exploiting short read data. We tested ploidyNGS using both simulated and real NGS data of the model yeast Saccharomyces cerevisiae. ploidyNGS allows the identification of the ploidy level of a newly sequenced genome in a visual way.

Availability and Implementation

ploidyNGS is available under the GNU General Public License (GPL) at https://github.com/diriano/ploidyNGS. ploidyNGS is implemented in Python and R.

Contact

diriano@gmail.com.

Structural and functional characterization of a highly secreted α-l-arabinofuranosidase (GH62) from Aspergillus nidulans grown on sugarcane bagasse

Carbohydrate-Active Enzymes are key enzymes for biomass-to-bioproducts conversion. α-l-Arabinofuranosidases that belong to the Glycoside Hydrolase family 62 (GH62) have important applications in biofuel production from plant biomass by hydrolyzing arabinoxylans, found in both the primary and secondary cell walls of plants. In this work, we identified a GH62 α-l-arabinofuranosidase (AnAbf62A_wt) that was highly secreted when Aspergillus nidulans was cultivated on sugarcane bagasse. The gene AN7908 was cloned and transformed in A. nidulans for homologous production of AnAbf62A_wt, and we confirmed that the enzyme is N-glycosylated at asparagine 83 by mass spectrometry analysis. The enzyme was also expressed in Escherichia coli and the studies of circular dichroism showed that the melting temperature and structural profile of AnAbf62A_wt and the non-glycosylated enzyme from E. coli (AnAbf62A_deglyc) were highly similar. In addition, the designed glycomutant AnAbf62A_N83Q presented similar patterns of secretion and activity to the AnAbf62A_wt, indicating that the N-glycan does not influence the properties of this enzyme. The crystallographic structure of AnAbf62A_deglyc was obtained and the 1.7Å resolution model showed a five-bladed β-propeller fold, which is conserved in family GH62. Mutants AnAbf62A_Y312F and AnAbf62A_Y312S showed that Y312 was an important substrate-binding residue. Molecular dynamics simulations indicated that the loop containing Y312 could access different conformations separated by moderately low energy barriers. One of these conformations, comprising a local minimum, is responsible for placing Y312 in the vicinity of the arabinose glycosidic bond, and thus, may be important for catalytic efficiency.

Comparative transcriptome analysis reveals different strategies for degradation of steam-exploded sugarcane bagasse by Aspergillus niger and Trichoderma reesei

Background

Second generation (2G) ethanol is produced by breaking down lignocellulosic biomass into fermentable sugars. In Brazil, sugarcane bagasse has been proposed as the lignocellulosic residue for this biofuel production. The enzymatic cocktails for the degradation of biomass-derived polysaccharides are mostly produced by fungi, such as Aspergillus niger and Trichoderma reesei. However, it is not yet fully understood how these microorganisms degrade plant biomass. In order to identify transcriptomic changes during steam-exploded bagasse (SEB) breakdown, we conducted a RNA-seq comparative transcriptome profiling of both fungi growing on SEB as carbon source.

Results

Particular attention was focused on CAZymes, sugar transporters, transcription factors (TFs) and other proteins related to lignocellulose degradation. Although genes coding for the main enzymes involved in biomass deconstruction were expressed by both fungal strains since the beginning of the growth in SEB, significant differences were found in their expression profiles. The expression of these enzymes is mainly regulated at the transcription level, and A. niger and T. reesei also showed differences in TFs content and in their expression. Several sugar transporters that were induced in both fungal strains could be new players on biomass degradation besides their role in sugar uptake. Interestingly, our findings revealed that in both strains several genes that code for proteins of unknown function and pro-oxidant, antioxidant, and detoxification enzymes were induced during growth in SEB as carbon source, but their specific roles on lignocellulose degradation remain to be elucidated.

Conclusions

This is the first report of a time-course experiment monitoring the degradation of pretreated bagasse by two important fungi using the RNA-seq technology. It was possible to identify a set of genes that might be applied in several biotechnology fields. The data suggest that these two microorganisms employ different strategies for biomass breakdown. This knowledge can be exploited for the rational design of enzymatic cocktails and 2G ethanol production improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3857-5) contains supplementary material, which is available to authorized users.

Draft genome sequencing of the sugarcane hybrid SP80-3280

Sugarcane commercial cultivar SP80-3280 has been used as a model for genomic analyses in Brazil. Here we present a draft genome sequence employing Illumina TruSeq Synthetic Long reads. The dataset is available from NCBI BioProject with accession PRJNA272769.

Draft genome sequencing of the sugarcane hybrid SP80-3280

Sugarcane commercial cultivar SP80-3280 has been used as a model for genomic analyses in Brazil. Here we present a draft genome sequence employing Illumina TruSeq Synthetic Long reads. The dataset is available from NCBI BioProject with accession PRJNA272769.

Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus Aspergillus

BACKGROUND

The fungal genus Aspergillus is of critical importance to humankind. Species include those with industrial applications, important pathogens of humans, animals and crops, a source of potent carcinogenic contaminants of food, and an important genetic model. The genome sequences of eight aspergilli have already been explored to investigate aspects of fungal biology, raising questions about evolution and specialization within this genus.

RESULTS

We have generated genome sequences for ten novel, highly diverse Aspergillus species and compared these in detail to sister and more distant genera. Comparative studies of key aspects of fungal biology, including primary and secondary metabolism, stress response, biomass degradation, and signal transduction, revealed both conservation and diversity among the species. Observed genomic differences were validated with experimental studies. This revealed several highlights, such as the potential for sex in asexual species, organic acid production genes being a key feature of black aspergilli, alternative approaches for degrading plant biomass, and indications for the genetic basis of stress response. A genome-wide phylogenetic analysis demonstrated in detail the relationship of the newly genome sequenced species with other aspergilli.

CONCLUSIONS

Many aspects of biological differences between fungal species cannot be explained by current knowledge obtained from genome sequences. The comparative genomics and experimental study, presented here, allows for the first time a genus-wide view of the biological diversity of the aspergilli and in many, but not all, cases linked genome differences to phenotype. Insights gained could be exploited for biotechnological and medical applications of fungi.

Draft genome sequencing of the sugarcane hybrid SP80-3280

Sugarcane commercial cultivar SP80-3280 has been used as a model for genomic analyses in Brazil. Here we present a draft genome sequence employing Illumina TruSeq Synthetic Long reads. The dataset is available from NCBI BioProject with accession PRJNA272769.

RNAseq reveals hydrophobins that are involved in the adaptation of Aspergillus nidulans to lignocellulose

BACKGROUND

Sugarcane is one of the world's most profitable crops. Waste steam-exploded sugarcane bagasse (SEB) is a cheap, abundant, and renewable lignocellulosic feedstock for the next-generation biofuels. In nature, fungi seldom exist as planktonic cells, similar to those found in the nutrient-rich environment created within an industrial fermenter. Instead, fungi predominantly form biofilms that allow them to thrive in hostile environments.

RESULTS

In turn, we adopted an RNA-sequencing approach to interrogate how the model fungus, Aspergillus nidulans, adapts to SEB, revealing the induction of carbon starvation responses and the lignocellulolytic machinery, in addition to morphological adaptations. Genetic analyses showed the importance of hydrophobins for growth on SEB. The major hydrophobin, RodA, was retained within the fungal biofilm on SEB fibres. The StuA transcription factor that regulates fungal morphology was up-regulated during growth on SEB and controlled hydrophobin gene induction. The absence of the RodA or DewC hydrophobins reduced biofilm formation. The loss of a RodA or a functional StuA reduced the retention of the hydrolytic enzymes within the vicinity of the fungus. Hence, hydrophobins promote biofilm formation on SEB, and may enhance lignocellulose utilisation via promoting a compact substrate-enzyme-fungus structure.

CONCLUSION

This novel study highlights the importance of hydrophobins to the formation of biofilms and the efficient deconstruction of lignocellulose.

Physiological and transcriptional analyses of developmental stages along sugarcane leaf

BACKGROUND

Sugarcane is one of the major crops worldwide. It is cultivated in over 100 countries on 22 million ha. The complex genetic architecture and the lack of a complete genomic sequence in sugarcane hamper the adoption of molecular approaches to study its physiology and to develop new varieties. Investments on the development of new sugarcane varieties have been made to maximize sucrose yield, a trait dependent on photosynthetic capacity. However, detailed studies on sugarcane leaves are scarce. In this work, we report the first molecular and physiological characterization of events taking place along a leaf developmental gradient in sugarcane.

RESULTS

Photosynthetic response to CO2 indicated divergence in photosynthetic capacity based on PEPcase activity, corroborated by activity quantification (both in vivo and in vitro) and distinct levels of carbon discrimination on different segments along leaf length. Additionally, leaf segments had contrasting amount of chlorophyll, nitrogen and sugars. RNA-Seq data indicated a plethora of biochemical pathways differentially expressed along the leaf. Some transcription factors families were enriched on each segment and their putative functions corroborate with the distinct developmental stages. Several genes with higher expression in the middle segment, the one with the highest photosynthetic rates, were identified and their role in sugarcane productivity is discussed. Interestingly, sugarcane leaf segments had a different transcriptional behavior compared to previously published data from maize.

CONCLUSION

This is the first report of leaf developmental analysis in sugarcane. Our data on sugarcane is another source of information for further studies aiming to understand and/or improve C4 photosynthesis. The segments used in this work were distinct in their physiological status allowing deeper molecular analysis. Although limited in some aspects, the comparison to maize indicates that all data acquired on one C4 species cannot always be easily extrapolated to other species. However, our data indicates that some transcriptional factors were segment-specific and the sugarcane leaf undergoes through the process of suberizarion, photosynthesis establishment and senescence.

Insights into immune responses in oral cancer through proteomic analysis of saliva and salivary extracellular vesicles

The development and progression of oral cavity squamous cell carcinoma (OSCC) involves complex cellular mechanisms that contribute to the low five-year survival rate of approximately 20% among diagnosed patients. However, the biological processes essential to tumor progression are not completely understood. Therefore, detecting alterations in the salivary proteome may assist in elucidating the cellular mechanisms modulated in OSCC and improve the clinical prognosis of the disease. The proteome of whole saliva and salivary extracellular vesicles (EVs) from patients with OSCC and healthy individuals were analyzed by LC-MS/MS and label-free protein quantification. Proteome data analysis was performed using statistical, machine learning and feature selection methods with additional functional annotation. Biological processes related to immune responses, peptidase inhibitor activity, iron coordination and protease binding were overrepresented in the group of differentially expressed proteins. Proteins related to the inflammatory system, transport of metals and cellular growth and proliferation were identified in the proteome of salivary EVs. The proteomics data were robust and could classify OSCC with 90% accuracy. The saliva proteome analysis revealed that immune processes are related to the presence of OSCC and indicate that proteomics data can contribute to determining OSCC prognosis.

Identification of transcription factor genes and their correlation with the high diversity of stramenopiles

The biological diversity among Stramenopiles is striking; they range from large multicellular seaweeds to tiny unicellular species, they embrace many ecologically important autothrophic (e.g., diatoms, brown algae), and heterotrophic (e.g., oomycetes) groups. Transcription factors (TFs) and other transcription regulators (TRs) regulate spatial and temporal gene expression. A plethora of transcriptional regulatory proteins have been identified and classified into families on the basis of sequence similarity. The purpose of this work is to identify the TF and TR complement in diverse species belonging to Stramenopiles in order to understand how these regulators may contribute to their observed diversity. We identified and classified 63 TF and TR families in 11 species of Stramenopiles. In some species we found gene families with high relative importance. Taking into account the 63 TF and TR families identified, 28 TF and TR families were established to be positively correlated with specific traits like number of predicted proteins, number of flagella and number of cell types during the life cycle. Additionally, we found gains and losses in TF and TR families specific to some species and clades, as well as, two families with high abundance specific to the autotrophic species and three families with high abundance specific to the heterotropic species. For the first time, there is a systematic search of TF and TR families in Stramenopiles. The attempts to uncover relationships between these families and the complexity of this group may be of great impact, considering that there are several important pathogens of plants and animals, as well as, important species involved in carbon cycling. Specific TF and TR families identified in this work appear to be correlated with particular traits in the Stramenopiles group and may be correlated with the high complexity and diversity in Stramenopiles.

Pseudozyma brasiliensis sp. nov., a xylanolytic, ustilaginomycetous yeast species isolated from an insect pest of sugarcane roots

A novel ustilaginomycetous yeast isolated from the intestinal tract of an insect pest of sugarcane roots in Ribeirão Preto, São Paulo State, Brazil, represents a novel species of the genus Pseudozyma based on molecular analyses of the D1/D2 rDNA large subunit and the internal transcribed spacer (ITS1+ITS2) regions. The name Pseudozyma brasiliensis sp. nov. is proposed for this species, with GHG001(T) ( = CBS 13268(T) = UFMG-CM-Y307(T)) as the type strain. P. brasiliensis sp. nov. is a sister species of Pseudozyma vetiver, originally isolated from leaves of vetiver grass and sugarcane in Thailand. P. brasiliensis sp. nov. is able to grow well with xylan as the sole carbon source and produces high levels of an endo-1,4-xylanase that has a higher specific activity in comparison with other eukaryotic xylanases. This enzyme has a variety of industrial applications, indicating the great biotechnological potential of P. brasiliensis.

Annotation of a hybrid partial genome of the coffee rust (Hemileia vastatrix) contributes to the gene repertoire catalog of the Pucciniales

Coffee leaf rust caused by the fungus Hemileia vastatrix is the most damaging disease to coffee worldwide. The pathogen has recently appeared in multiple outbreaks in coffee producing countries resulting in significant yield losses and increases in costs related to its control. New races/isolates are constantly emerging as evidenced by the presence of the fungus in plants that were previously resistant. Genomic studies are opening new avenues for the study of the evolution of pathogens, the detailed description of plant-pathogen interactions and the development of molecular techniques for the identification of individual isolates. For this purpose we sequenced 8 different H. vastatrix isolates using NGS technologies and gathered partial genome assemblies due to the large repetitive content in the coffee rust hybrid genome; 74.4% of the assembled contigs harbor repetitive sequences. A hybrid assembly of 333 Mb was built based on the 8 isolates; this assembly was used for subsequent analyses. Analysis of the conserved gene space showed that the hybrid H. vastatrix genome, though highly fragmented, had a satisfactory level of completion with 91.94% of core protein-coding orthologous genes present. RNA-Seq from urediniospores was used to guide the de novo annotation of the H. vastatrix gene complement. In total, 14,445 genes organized in 3921 families were uncovered; a considerable proportion of the predicted proteins (73.8%) were homologous to other Pucciniales species genomes. Several gene families related to the fungal lifestyle were identified, particularly 483 predicted secreted proteins that represent candidate effector genes and will provide interesting hints to decipher virulence in the coffee rust fungus. The genome sequence of Hva will serve as a template to understand the molecular mechanisms used by this fungus to attack the coffee plant, to study the diversity of this species and for the development of molecular markers to distinguish races/isolates.

Draft Genome Sequence of Pseudozyma brasiliensis sp. nov. Strain GHG001, a High Producer of Endo-1,4-Xylanase Isolated from an Insect Pest of Sugarcane

Here, we present the nuclear and mitochondrial genome sequences of Pseudozyma brasiliensis sp. nov. strain GHG001. P. brasiliensis sp. nov. is the closest relative of Pseudozyma vetiver. P. brasiliensis sp. nov. is capable of growing on xylose or xylan as a sole carbon source and has great biotechnological potential.

Genome-wide identification of regulatory elements and reconstruction of gene regulatory networks of the green alga Chlamydomonas reinhardtii under carbon deprivation

The unicellular green alga Chlamydomonas reinhardtii is a long-established model organism for studies on photosynthesis and carbon metabolism-related physiology. Under conditions of air-level carbon dioxide concentration [CO2], a carbon concentrating mechanism (CCM) is induced to facilitate cellular carbon uptake. CCM increases the availability of carbon dioxide at the site of cellular carbon fixation. To improve our understanding of the transcriptional control of the CCM, we employed FAIRE-seq (formaldehyde-assisted Isolation of Regulatory Elements, followed by deep sequencing) to determine nucleosome-depleted chromatin regions of algal cells subjected to carbon deprivation. Our FAIRE data recapitulated the positions of known regulatory elements in the promoter of the periplasmic carbonic anhydrase (Cah1) gene, which is upregulated during CCM induction, and revealed new candidate regulatory elements at a genome-wide scale. In addition, time series expression patterns of 130 transcription factor (TF) and transcription regulator (TR) genes were obtained for cells cultured under photoautotrophic condition and subjected to a shift from high to low [CO2]. Groups of co-expressed genes were identified and a putative directed gene-regulatory network underlying the CCM was reconstructed from the gene expression data using the recently developed IOTA (inner composition alignment) method. Among the candidate regulatory genes, two members of the MYB-related TF family, Lcr1 (Low-CO 2 response regulator 1) and Lcr2 (Low-CO2 response regulator 2), may play an important role in down-regulating the expression of a particular set of TF and TR genes in response to low [CO2]. The results obtained provide new insights into the transcriptional control of the CCM and revealed more than 60 new candidate regulatory genes. Deep sequencing of nucleosome-depleted genomic regions indicated the presence of new, previously unknown regulatory elements in the C. reinhardtii genome. Our work can serve as a basis for future functional studies of transcriptional regulator genes and genomic regulatory elements in Chlamydomonas.

Molecular Evolution of Slow and Quick Anion Channels (SLACs and QUACs/ALMTs)

Electrophysiological analyses conducted about 25 years ago detected two types of anion channels in the plasma membrane of guard cells. One type of channel responds slowly to changes in membrane voltage while the other responds quickly. Consequently, they were named SLAC, for SLow Anion Channel, and QUAC, for QUick Anion Channel. Recently, genes SLAC1 and QUAC1/ALMT12, underlying the two different anion current components, could be identified in the model plant Arabidopsis thaliana. Expression of the gene products in Xenopus oocytes confirmed the quick and slow current kinetics. In this study we provide an overview on our current knowledge on slow and quick anion channels in plants and analyze the molecular evolution of ALMT/QUAC-like and SLAC-like channels. We discovered fingerprints that allow screening databases for these channel types and were able to identify 192 (177 non-redundant) SLAC-like and 422 (402 non-redundant) ALMT/QUAC-like proteins in the fully sequenced genomes of 32 plant species. Phylogenetic analyses provided new insights into the molecular evolution of these channel types. We also combined sequence alignment and clustering with predictions of protein features, leading to the identification of known conserved phosphorylation sites in SLAC1-like channels along with potential sites that have not been yet experimentally confirmed. Using a similar strategy to analyze the hydropathicity of ALMT/QUAC-like channels, we propose a modified topology with additional transmembrane regions that integrates structure and function of these membrane proteins. Our results suggest that cross-referencing phylogenetic analyses with position-specific protein properties and functional data could be a very powerful tool for genome research approaches in general.

Phylogenetic analysis of k(+) transporters in bryophytes, lycophytes, and flowering plants indicates a specialization of vascular plants

As heritage from early evolution, potassium (K(+)) is absolutely necessary for all living cells. It plays significant roles as stabilizer in metabolism and is important for enzyme activation, stabilization of protein synthesis, and neutralization of negative charges on cellular molecules as proteins and nucleic acids. Land plants even enlarged this spectrum of K(+) utilization after having gone ashore, despite the fact that K(+) is far less available in their new oligotrophic habitats than in sea water. Inevitably, plant cells had to improve and to develop unique transport systems for K(+) accumulation and distribution. In the past two decades a manifold of K(+) transporters from flowering plants has been identified at the molecular level. The recently published genome of the fern ally Selaginella moellendorffii now helps in providing a better understanding on the molecular changes involved in the colonization of land and the development of the vasculature and the seeds. In this article we present an inventory of K(+) transporters of this lycophyte and pigeonhole them together with their relatives from the moss Physcomitrella patens, the monocotyledon Oryza sativa, and two dicotyledonous species, the herbaceous plant Arabidopsis thaliana, and the tree Populus trichocarpa. Interestingly, the transition of green plants from an aqueous to a dry environment coincides with a dramatic reduction in the diversity of voltage-gated potassium channels followed by a diversification on the basis of one surviving K(+) channel class. The first appearance of K(+) release (K(out)) channels in S. moellendorffii that were shown in Arabidopsis to be involved in xylem loading and guard cell closure coincides with the specialization of vascular plants and may indicate an important adaptive step.

Combined transcription factor profiling, microarray analysis and metabolite profiling reveals the transcriptional control of metabolic shifts occurring during tomato fruit development

Maturation of fleshy fruits such as tomato (Solanum lycopersicum) is subject to tight genetic control. Here we describe the development of a quantitative real-time PCR platform that allows accurate quantification of the expression level of approximately 1000 tomato transcription factors. In addition to utilizing this novel approach, we performed cDNA microarray analysis and metabolite profiling of primary and secondary metabolites using GC-MS and LC-MS, respectively. We applied these platforms to pericarp material harvested throughout fruit development, studying both wild-type Solanum lycopersicum cv. Ailsa Craig and the hp1 mutant. This mutant is functionally deficient in the tomato homologue of the negative regulator of the light signal transduction gene DDB1 from Arabidopsis, and is furthermore characterized by dramatically increased pigment and phenolic contents. We choose this particular mutant as it had previously been shown to have dramatic alterations in the content of several important fruit metabolites but relatively little impact on other ripening phenotypes. The combined dataset was mined in order to identify metabolites that were under the control of these transcription factors, and, where possible, the respective transcriptional regulation underlying this control. The results are discussed in terms of both programmed fruit ripening and development and the transcriptional and metabolic shifts that occur in parallel during these processes.

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

Proteome-wide survey of phosphorylation patterns affected by nuclear DNA polymorphisms in Arabidopsis thaliana

Background

Protein phosphorylation is an important post-translational modification influencing many aspects of dynamic cellular behavior. Site-specific phosphorylation of amino acid residues serine, threonine, and tyrosine can have profound effects on protein structure, activity, stability, and interaction with other biomolecules. Phosphorylation sites can be affected in diverse ways in members of any species, one such way is through single nucleotide polymorphisms (SNPs). The availability of large numbers of experimentally identified phosphorylation sites, and of natural variation datasets in Arabidopsis thaliana prompted us to analyze the effect of non-synonymous SNPs (nsSNPs) onto phosphorylation sites.

Results

From the analyses of 7,178 experimentally identified phosphorylation sites we found that: (i) Proteins with multiple phosphorylation sites occur more often than expected by chance. (ii) Phosphorylation hotspots show a preference to be located outside conserved domains. (iii) nsSNPs affected experimental phosphorylation sites as much as the corresponding non-phosphorylated amino acid residues. (iv) Losses of experimental phosphorylation sites by nsSNPs were identified in 86 A. thaliana proteins, among them receptor proteins were overrepresented.

These results were confirmed by similar analyses of predicted phosphorylation sites in A. thaliana. In addition, predicted threonine phosphorylation sites showed a significant enrichment of nsSNPs towards asparagines and a significant depletion of the synonymous substitution. Proteins in which predicted phosphorylation sites were affected by nsSNPs (loss and gain), were determined to be mainly receptor proteins, stress response proteins and proteins involved in nucleotide and protein binding. Proteins involved in metabolism, catalytic activity and biosynthesis were less affected.

Conclusions

We analyzed more than 7,100 experimentally identified phosphorylation sites in almost 4,300 protein-coding loci in silico, thus constituting the largest phosphoproteomics dataset for A. thaliana available to date. Our findings suggest a relatively high variability in the presence or absence of phosphorylation sites between different natural accessions in receptor and other proteins involved in signal transduction. Elucidating the effect of phosphorylation sites affected by nsSNPs on adaptive responses represents an exciting research goal for the future.

PlnTFDB: updated content and new features of the plant transcription factor database

The Plant Transcription Factor Database (PlnTFDB; http://plntfdb.bio.uni-potsdam.de/v3.0/) is an integrative database that provides putatively complete sets of transcription factors (TFs) and other transcriptional regulators (TRs) in plant species (sensu lato) whose genomes have been completely sequenced and annotated. The complete sets of 84 families of TFs and TRs from 19 species ranging from unicellular red and green algae to angiosperms are included in PlnTFDB, representing >1.6 billion years of evolution of gene regulatory networks. For each gene family, a basic description is provided that is complemented by literature references, and multiple sequence alignments of protein domains. TF or TR gene entries include information of expressed sequence tags, 3D protein structures of homologous proteins, domain architecture and cross-links to other computational resources online. Moreover, the different species in PlnTFDB are linked to each other by means of orthologous genes facilitating cross-species comparisons.

QuantPrime--a flexible tool for reliable high-throughput primer design for quantitative PCR

Background

Medium- to large-scale expression profiling using quantitative polymerase chain reaction (qPCR) assays are becoming increasingly important in genomics research. A major bottleneck in experiment preparation is the design of specific primer pairs, where researchers have to make several informed choices, often outside their area of expertise. Using currently available primer design tools, several interactive decisions have to be made, resulting in lengthy design processes with varying qualities of the assays.

Results

Here we present QuantPrime, an intuitive and user-friendly, fully automated tool for primer pair design in small- to large-scale qPCR analyses. QuantPrime can be used online through the internet or on a local computer after download; it offers design and specificity checking with highly customizable parameters and is ready to use with many publicly available transcriptomes of important higher eukaryotic model organisms and plant crops (currently 295 species in total), while benefiting from exon-intron border and alternative splice variant information in available genome annotations. Experimental results with the model plant Arabidopsis thaliana, the crop Hordeum vulgare and the model green alga Chlamydomonas reinhardtii show success rates of designed primer pairs exceeding 96%.

Conclusion

QuantPrime constitutes a flexible, fully automated web application for reliable primer design for use in larger qPCR experiments, as proven by experimental data. The flexible framework is also open for simple use in other quantification applications, such as hydrolyzation probe design for qPCR and oligonucleotide probe design for quantitative in situ hybridization. Future suggestions made by users can be easily implemented, thus allowing QuantPrime to be developed into a broad-range platform for the design of RNA expression assays.

GabiPD: the GABI primary database--a plant integrative 'omics' database

The GABI Primary Database, GabiPD (http://www.gabipd.org/), was established in the frame of the German initiative for Genome Analysis of the Plant Biological System (GABI). The goal of GabiPD is to collect, integrate, analyze and visualize primary information from GABI projects. GabiPD constitutes a repository and analysis platform for a wide array of heterogeneous data from high-throughput experiments in several plant species. Data from different ‘omics’ fronts are incorporated (i.e. genomics, transcriptomics, proteomics and metabolomics), originating from 14 different model or crop species. We have developed the concept of GreenCards for text-based retrieval of all data types in GabiPD (e.g. clones, genes, mutant lines). All data types point to a central Gene GreenCard, where gene information is integrated from genome projects or NCBI UniGene sets. The centralized Gene GreenCard allows visualizing ESTs aligned to annotated transcripts as well as displaying identified protein domains and gene structure. Moreover, GabiPD makes available interactive genetic maps from potato and barley, and protein 2DE gels from Arabidopsis thaliana and Brassica napus. Gene expression and metabolic-profiling data can be visualized through MapManWeb. By the integration of complex data in a framework of existing knowledge, GabiPD provides new insights and allows for new interpretations of the data.

Transcription factors regulating leaf senescence in Arabidopsis thaliana

Senescence is a highly regulated process, eventually leading to cell and tissue disintegration: a physiological process associated with nutrient (e.g. nitrogen) redistribution from leaves to reproductive organs. Senescence is not observed in young leaves, indicating that repressors efficiently act to suppress cell degradation during early leaf development and/or that senescence activators are switched on when a leaf ages. Thus, massive regulatory network re-wiring likely constitutes an important component of the pre-senescence process. Transcription factors (TFs) have been shown to be central elements of such regulatory networks. Here, we used quantitative real-time polymerase chain reaction (qRT-PCR) analysis to study the expression of 1880 TF genes during pre-senescence and early-senescence stages of leaf development, using Arabidopsis thaliana as a model. We show that the expression of 185 TF genes changes when leaves develop from half to fully expanded leaves and finally enter partial senescence. Our analysis identified 41 TF genes that were gradually up-regulated as leaves progressed through these developmental stages. We also identified 144 TF genes that were down-regulated during senescence. A considerable number of the senescence-regulated TF genes were found to respond to abiotic stress, and salt stress appeared to be the major factor controlling their expression. Our data indicate a peculiar fine-tuning of developmental shifts during late-leaf development that is controlled by TFs.

The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes

Background

Transcription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms.

Methodology/Principal Findings

We identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy.

Conclusions/Significance

Our data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments.

The Chlamydomonas genome reveals the evolution of key animal and plant functions

Chlamydomonas reinhardtii is a unicellular green alga whose lineage diverged from land plants over 1 billion years ago. It is a model system for studying chloroplast-based photosynthesis, as well as the structure, assembly, and function of eukaryotic flagella (cilia), which were inherited from the common ancestor of plants and animals, but lost in land plants. We sequenced the approximately 120-megabase nuclear genome of Chlamydomonas and performed comparative phylogenomic analyses, identifying genes encoding uncharacterized proteins that are likely associated with the function and biogenesis of chloroplasts or eukaryotic flagella. Analyses of the Chlamydomonas genome advance our understanding of the ancestral eukaryotic cell, reveal previously unknown genes associated with photosynthetic and flagellar functions, and establish links between ciliopathy and the composition and function of flagella.

The Chlamydomonas genome reveals the evolution of key animal and plant functions

Chlamydomonas reinhardtii is a unicellular green alga whose lineage diverged from land plants over 1 billion years ago. It is a model system for studying chloroplast-based photosynthesis, as well as the structure, assembly, and function of eukaryotic flagella (cilia), which were inherited from the common ancestor of plants and animals, but lost in land plants. We sequenced the approximately 120-megabase nuclear genome of Chlamydomonas and performed comparative phylogenomic analyses, identifying genes encoding uncharacterized proteins that are likely associated with the function and biogenesis of chloroplasts or eukaryotic flagella. Analyses of the Chlamydomonas genome advance our understanding of the ancestral eukaryotic cell, reveal previously unknown genes associated with photosynthetic and flagellar functions, and establish links between ciliopathy and the composition and function of flagella.

Genome-wide analysis of ABA-responsive elements ABRE and CE3 reveals divergent patterns in Arabidopsis and rice

Background

In plants, complex regulatory mechanisms are at the core of physiological and developmental processes. The phytohormone abscisic acid (ABA) is involved in the regulation of various such processes, including stomatal closure, seed and bud dormancy, and physiological responses to cold, drought and salinity stress. The underlying tissue or plant-wide control circuits often include combinatorial gene regulatory mechanisms and networks that we are only beginning to unravel with the help of new molecular tools. The increasing availability of genomic sequences and gene expression data enables us to dissect ABA regulatory mechanisms at the individual gene expression level. In this paper we used an in-silico-based approach directed towards genome-wide prediction and identification of specific features of ABA-responsive elements. In particular we analysed the genome-wide occurrence and positional arrangements of two well-described ABA-responsive cis-regulatory elements (CREs), ABRE and CE3, in thale cress (Arabidopsis thaliana) and rice (Oryza sativa).

Results

Our results show that Arabidopsis and rice use the ABA-responsive elements ABRE and CE3 distinctively. Earlier reports for various monocots have identified CE3 as a coupling element (CE) associated with ABRE. Surprisingly, we found that while ABRE is equally abundant in both species, CE3 is practically absent in Arabidopsis. ABRE-ABRE pairs are common in both genomes, suggesting that these can form functional ABA-responsive complexes (ABRCs) in Arabidopsis and rice. Furthermore, we detected distinct combinations, orientation patterns and DNA strand preferences of ABRE and CE3 motifs in rice gene promoters.

Conclusion

Our computational analyses revealed distinct recruitment patterns of ABA-responsive CREs in upstream sequences of Arabidopsis and rice. The apparent absence of CE3s in Arabidopsis suggests that another CE pairs with ABRE to establish a functional ABRC capable of interacting with transcription factors. Further studies will be needed to test whether the observed differences are extrapolatable to monocots and dicots in general, and to understand how they contribute to the fine-tuning of the hormonal response. The outcome of our investigation can now be used to direct future experimentation designed to further dissect the ABA-dependent regulatory networks.

PlnTFDB: an integrative plant transcription factor database

BACKGROUND

Transcription factors (TFs) are key regulatory proteins that enhance or repress the transcriptional rate of their target genes by binding to specific promoter regions (i.e. cis-acting elements) upon activation or de-activation of upstream signaling cascades. TFs thus constitute master control elements of dynamic transcriptional networks. TFs have fundamental roles in almost all biological processes (development, growth and response to environmental factors) and it is assumed that they play immensely important functions in the evolution of species. In plants, TFs have been employed to manipulate various types of metabolic, developmental and stress response pathways. Cross-species comparison and identification of regulatory modules and hence TFs is thought to become increasingly important for the rational design of new plant biomass. Up to now, however, no computational repository is available that provides access to the largely complete sets of transcription factors of sequenced plant genomes.

DESCRIPTION

PlnTFDB is an integrative plant transcription factor database that provides a web interface to access large (close to complete) sets of transcription factors of several plant species, currently encompassing Arabidopsis thaliana (thale cress), Populus trichocarpa (poplar), Oryza sativa (rice), Chlamydomonas reinhardtii and Ostreococcus tauri. It also provides an access point to its daughter databases of a species-centered representation of transcription factors (OstreoTFDB, ChlamyTFDB, ArabTFDB, PoplarTFDB and RiceTFDB). Information including protein sequences, coding regions, genomic sequences, expressed sequence tags (ESTs), domain architecture and scientific literature is provided for each family.

CONCLUSION

We have created lists of putatively complete sets of transcription factors and other transcriptional regulators for five plant genomes. They are publicly available through http://plntfdb.bio.uni-potsdam.de. Further data will be included in the future when the sequences of other plant genomes become available.

Orphan transcripts in Arabidopsis thaliana: identification of several hundred previously unrecognized genes

Expressed sequence tags (ESTs) represent a huge resource for the discovery of previously unknown genetic information and functional genome assignment. In this study we screened a collection of 178 292 ESTs from Arabidopsis thaliana by testing them against previously annotated genes of the Arabidopsis genome. We identified several hundreds of new transcripts that match the Arabidopsis genome at so far unassigned loci. The transcriptional activity of these loci was independently confirmed by comparison with the Salk Whole Genome Array Data. To a large extent, the newly identified transcriptionally active genomic regions do not encode 'classic' proteins, but instead generate non-coding RNAs and/or small peptide-coding RNAs of presently unknown biological function. More than 560 transcripts identified in this study are not represented by the Affymetrix GeneChip arrays currently widely used for expression profiling in A. thaliana. Our data strongly support the hypothesis that numerous previously unknown genes exist in the Arabidopsis genome.