Gains and Losses of Transcription Factor Binding Sites in Saccharomyces cerevisiae and Saccharomyces paradoxus

Evolutionary studies of the bHLH transcription factors belonging to MBW complex: their role in seed development

BACKGROUND AND AIMS

The MBW complex consist of proteins belonging to three major families (MYB, bHLH and WDR) involved in various processes throughout plant development: epidermal cell development, mucilage secretory cells and flavonoid biosynthesis. Recently, it has been reported that TT8, encoding a bHLH transcription factor, is involved in the biosynthesis of flavonoids in the seed coat and it also plays a role in bypassing the postzygotic barrier resulting from an unbalance in genetic loads of the parental lines. Here, we focus on the functional evolution, in seed development, of the bHLH proteins that are part of the MBW complex, complemented with a literature review.

METHODS

Phylogenetic analyses performed across seed plants and expression analyses in the reproductive tissues of four selected angiosperms (Arabidopsis thaliana, Brassica napus, Capsella rubella and Solanum lycopersicum) allow us to hypothesize on the evolution of its functions.

KEY RESULTS

TT8 expression in the innermost layer of the seed coat is conserved in the selected angiosperms. However, except for Arabidopsis, TT8 is also expressed in ovules, carpels and fruits. The homologues belonging to the sister clade of TT8, EGL3/GL3, involved in trichome development, are expressed in the outermost layer of the seed coat, suggesting potential roles in mucilage.

CONCLUSIONS

The ancestral function of these genes appears to be flavonoid biosynthesis, and the conservation of TT8 expression patterns in the innermost layer of the seed coat in angiosperms suggests that their function in postzygotic barriers might also be conserved. Moreover, the literature review and the results of the present study suggest a sophisticated association, linking the mechanisms of action of these genes to the cross-communication activity between the different tissues of the seed. Thus, it provides avenues to study the mechanisms of action of TT8 in the postzygotic triploid block, which is crucial because it impacts seed development in unbalanced crosses.

Extensive non-coding sequence divergence between the major human pathogen Aspergillus fumigatus and its relatives

Invasive aspergillosis is a deadly fungal disease; more than 400,000 patients are infected worldwide each year and the mortality rate can be as high as 50-95%. Of the ~450 species in the genus Aspergillus only a few are known to be clinically relevant, with the major pathogen Aspergillus fumigatus being responsible for ~50% of all invasive mold infections. Genomic comparisons between A. fumigatus and other Aspergillus species have historically focused on protein-coding regions. However, most A. fumigatus genes, including those that modulate its virulence, are also present in other pathogenic and non-pathogenic closely related species. Our hypothesis is that differential gene regulation - mediated through the non-coding regions upstream of genes' first codon - contributes to A. fumigatus pathogenicity. To begin testing this, we compared non-coding regions upstream of the first codon of single-copy orthologous genes from the two A. fumigatus reference strains Af293 and A1163 and eight closely related Aspergillus section Fumigati species. We found that these non-coding regions showed extensive sequence variation and lack of homology across species. By examining the evolutionary rates of both protein-coding and non-coding regions in a subset of orthologous genes with highly conserved non-coding regions across the phylogeny, we identified 418 genes, including 25 genes known to modulate A. fumigatus virulence, whose non-coding regions exhibit a different rate of evolution in A. fumigatus. Examination of sequence alignments of these non-coding regions revealed numerous instances of insertions, deletions, and other types of mutations of at least a few nucleotides in A. fumigatus compared to its close relatives. These results show that closely related Aspergillus species that vary greatly in their pathogenicity exhibit extensive non-coding sequence variation and identify numerous changes in non-coding regions of A. fumigatus genes known to contribute to virulence.

Promoter conservation in HDACs points to functional implications

Background

Histone deacetylases (HDACs) are the proteins responsible for removing the acetyl group from lysine residues of core histones in chromosomes, a crucial component of gene regulation. Eleven known HDACs exist in humans and most other vertebrates. While the basic function of HDACs has been well characterized and new discoveries are still being made, the transcriptional regulation of their corresponding genes is still poorly understood.

Results

Here, we conducted a computational analysis of the eleven HDAC promoter sequences in 25 vertebrate species to determine whether transcription factor binding sites (TFBSs) are conserved in HDAC evolution, and if so, whether they provide useful information about HDAC expression and function. Furthermore, we used tissue-specific information of transcription factors to investigate the potential expression patterns of HDACs in different human tissues based on their transcription factor binding sites. We found that the TFBS profiles of most of the HDACs were well conserved in closely related species for all HDAC promoters except HDAC7 and HDAC10. HDAC5 had particularly strong conservation across over half of the species studied, with nearly identical profiles in the primate species. Our comparisons of TFBSs with the tissue specific gene expression profiles of their corresponding TFs showed that most HDACs had the ability to be ubiquitously expressed. A few HDAC promoters exhibited the potential for preferential expression in certain tissues, most notably HDAC11 in gall bladder, while HDAC9 seemed to have less propensity for expression in the nervous system.

Conclusions

In general, we found evolutionary conservation in HDAC promoters that seems to be more prominent for the ubiquitously expressed HDACs. In turn, when conservation did not follow usual phylogeny, human TFBS patterns indicated possible functional relevance. While we found that HDACs appear to uniformly expressed, we confirm that the functional differences in HDACs may be less a matter of location of activity than a question of which proteins and which acetyl groups they may be acting on.

Electronic supplementary material

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

Unrealistic phylogenetic trees may improve phylogenetic footprinting

Motivation

The computational investigation of DNA binding motifs from binding sites is one of the classic tasks in bioinformatics and a prerequisite for understanding gene regulation as a whole. Due to the development of sequencing technologies and the increasing number of available genomes, approaches based on phylogenetic footprinting become increasingly attractive. Phylogenetic footprinting requires phylogenetic trees with attached substitution probabilities for quantifying the evolution of binding sites, but these trees and substitution probabilities are typically not known and cannot be estimated easily.

Results

Here, we investigate the influence of phylogenetic trees with different substitution probabilities on the classification performance of phylogenetic footprinting using synthetic and real data. For synthetic data we find that the classification performance is highest when the substitution probability used for phylogenetic footprinting is similar to that used for data generation. For real data, however, we typically find that the classification performance of phylogenetic footprinting surprisingly increases with increasing substitution probabilities and is often highest for unrealistically high substitution probabilities close to one. This finding suggests that choosing realistic model assumptions might not always yield optimal predictions in general and that choosing unrealistically high substitution probabilities close to one might actually improve the classification performance of phylogenetic footprinting.

Availability and Implementation

The proposed PF is implemented in JAVA and can be downloaded from https://github.com/mgledi/PhyFoo

Supplementary information

Supplementary data are available at Bioinformatics online.