HSDFinder: A BLAST-Based Strategy for Identifying Highly Similar Duplicated Genes in Eukaryotic Genomes

Gene duplication is an important evolutionary mechanism capable of providing new genetic material for adaptive and nonadaptive evolution. However, bioinformatics tools for identifying duplicate genes are often limited to the detection of paralogs in multiple species or to specific types of gene duplicates, such as retrocopies. Here, we present a user-friendly, BLAST-based web tool, called HSDFinder, which can identify, annotate, categorize, and visualize highly similar duplicate genes (HSDs) in eukaryotic nuclear genomes. HSDFinder includes an online heatmap plotting option, allowing users to compare HSDs among different species and visualize the results in different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional categories. The external software requirements are BLAST, InterProScan, and KEGG. The utility of HSDFinder was tested on various model eukaryotic species, including Chlamydomonas reinhardtii, Arabidopsis thaliana, Oryza sativa, and Zea mays as well as the psychrophilic green alga Chlamydomonas sp. UWO241, and was proven to be a practical and accurate tool for gene duplication analyses. The web tool is free to use at http://hsdfinder.com. Documentation and tutorials can be found via the GitHub: https://github.com/zx0223winner/HSDFinder.

The Genetic Profiles of TIF1 and TIF2 Duplicate Genes

<p>As one of the key steps in protein synthesis, translation initiation is subjected to multi-level regulation which is achieved via diverse mechanisms. The cell adjusts protein synthesis accordingly to its status and environment. The degree of contribution of the processes involved in the regulation of translation initiation is still poorly understood. The first part of this study focuses on identifying mechanisms of regulation in a translationally deficient yeast system, impaired by the loss of one or the other of the TIF1/2 duplicate genes, which together code for the eukaryotic initiation factor 4A (eIF4A). A major finding of this research is related to the functional competences associated with the two duplicate members of the gene pair. Although the genetic profile associated with TIF1 highlights a connection with transcriptional process, the majority of transcription-translation inter-talk is allocated with TIF2, along with a dense network of genetic interactions surrounding the SAGA complex. TIF2 is also the only paralog devoted to interactions with a substantial group of functionally related genes involved in early meiotic gene expression. Protein degradation in the global control of protein synthesis represents a fundamental process and accounts for diverse points of control, in particular through ubiquitination/deubiquitination. This research concludes that functional turnover of proteins and the translation/transcription inter-talk emerges as the most significant contributors to the sophistically regulated translational regulation, The second part of this study aims to determine the extent of similarity and divergence between the TIF1 and TIF2 paralogs. Growth of their individual deletion strains was challenged under different chemical and environmental conditions with the intent to explore the relative contributions of each duplicate in response to an extend range of perturbations. The pair of duplicates appeared convincingly robust in coping with these adversities under disparate cellular contexts, thus suggesting a highly conserved and backed-up genetic network. One of the primary treatments made use of lithium, a condition which was hoped to help, along with furthering our understanding of the TIF1 and TIF2 networks, in formulating an explanation on how augmented translation initiation overcomes lithium toxicity. Although this approach did not return results that could be used to address this unresolved topic, evaluation of genetic inhibition and suppression was highly instructive regarding the mechanisms of response triggered upon lithium/galactose stress. Regulation and synchronization of basic cellular processes were affected: emphasis brought on aspects of cell communication highlighted mechanisms articulated by kinase enzymes and the importance of repression of cell cycle progression in control of protein synthesis. Data from the screen also indicated the stress that combined lithium/galactose treatment places on central metabolic pathways, for instance those between the Leloir, gluconeogenesis, and trehalose pathways.</p>

The Genetic Profiles of TIF1 and TIF2 Duplicate Genes

<p>As one of the key steps in protein synthesis, translation initiation is subjected to multi-level regulation which is achieved via diverse mechanisms. The cell adjusts protein synthesis accordingly to its status and environment. The degree of contribution of the processes involved in the regulation of translation initiation is still poorly understood. The first part of this study focuses on identifying mechanisms of regulation in a translationally deficient yeast system, impaired by the loss of one or the other of the TIF1/2 duplicate genes, which together code for the eukaryotic initiation factor 4A (eIF4A). A major finding of this research is related to the functional competences associated with the two duplicate members of the gene pair. Although the genetic profile associated with TIF1 highlights a connection with transcriptional process, the majority of transcription-translation inter-talk is allocated with TIF2, along with a dense network of genetic interactions surrounding the SAGA complex. TIF2 is also the only paralog devoted to interactions with a substantial group of functionally related genes involved in early meiotic gene expression. Protein degradation in the global control of protein synthesis represents a fundamental process and accounts for diverse points of control, in particular through ubiquitination/deubiquitination. This research concludes that functional turnover of proteins and the translation/transcription inter-talk emerges as the most significant contributors to the sophistically regulated translational regulation, The second part of this study aims to determine the extent of similarity and divergence between the TIF1 and TIF2 paralogs. Growth of their individual deletion strains was challenged under different chemical and environmental conditions with the intent to explore the relative contributions of each duplicate in response to an extend range of perturbations. The pair of duplicates appeared convincingly robust in coping with these adversities under disparate cellular contexts, thus suggesting a highly conserved and backed-up genetic network. One of the primary treatments made use of lithium, a condition which was hoped to help, along with furthering our understanding of the TIF1 and TIF2 networks, in formulating an explanation on how augmented translation initiation overcomes lithium toxicity. Although this approach did not return results that could be used to address this unresolved topic, evaluation of genetic inhibition and suppression was highly instructive regarding the mechanisms of response triggered upon lithium/galactose stress. Regulation and synchronization of basic cellular processes were affected: emphasis brought on aspects of cell communication highlighted mechanisms articulated by kinase enzymes and the importance of repression of cell cycle progression in control of protein synthesis. Data from the screen also indicated the stress that combined lithium/galactose treatment places on central metabolic pathways, for instance those between the Leloir, gluconeogenesis, and trehalose pathways.</p>

Paralogs in the PKA Regulon Traveled Different Evolutionary Routes to Divergent Expression in Budding Yeast

Functional divergence of duplicate genes, or paralogs, is an important driver of novelty in evolution. In the model yeast Saccharomyces cerevisiae, there are 547 paralog gene pairs that survive from an interspecies Whole Genome Hybridization (WGH) that occurred ~100MYA. In this work, we report that ~1/6th (110) of these WGH paralogs pairs (or ohnologs) are differentially expressed with a striking pattern upon Protein Kinase A (PKA) inhibition. One member of each pair in this group has low basal expression that increases upon PKA inhibition, while the other has moderate and unchanging expression. For these genes, expression of orthologs upon PKA inhibition in the non-WGH species Kluyveromyces lactis and for PKA-related stresses in other budding yeasts shows unchanging expression, suggesting that lack of responsiveness to PKA was likely the typical ancestral phenotype prior to duplication. Promoter sequence analysis across related budding yeast species further revealed that the subsequent emergence of PKA-dependence took different evolutionary routes. In some examples, regulation by PKA and differential expression appears to have arisen following the WGH, while in others, regulation by PKA appears to have arisen in one of the two parental lineages prior to the WGH. More broadly, our results illustrate the unique opportunities presented by a WGH event for generating functional divergence by bringing together two parental lineages with separately evolved regulation into one species. We propose that functional divergence of two ohnologs can be facilitated through such regulatory divergence.

Adaptive evolution driving the young duplications in six Rosaceae species

Background In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported. Results In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses. Conclusions The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

In situ dissecting the evolution of gene duplication with different histone modification patterns based on high-throughput data analysis in Arabidopsis thaliana

Background Genetic regulation is known to contribute to the divergent expression of duplicate genes; however, little is known about how epigenetic modifications regulate the expression of duplicate genes in plants. Methods The histone modification (HM) profile patterns of different modes of gene duplication, including the whole genome duplication, proximal duplication, tandem duplication and transposed duplication were characterized based on ChIP-chip or ChIP-seq datasets. In this study, 10 distinct HM marks including H2Bub, H3K4me1, H3K4me2, H3K4me3, H3K9ac, H3K9me2, H3K27me1, H3K27me3, H3K36me3 and H3K14ac were analyzed. Moreover, the features of gene duplication with different HM patterns were characterized based on 88 RNA-seq datasets of Arabidopsis thaliana. Results This study showed that duplicate genes in Arabidopsis have a more similar HM pattern than single-copy genes in both their promoters and protein-coding regions. The evolution of HM marks is found to be coupled with coding sequence divergence and expression divergence after gene duplication. We found that functionally selective constraints may impose on epigenetic evolution after gene duplication. Furthermore, duplicate genes with distinct functions have more divergence in histone modification compared with the ones with the same function, while higher expression divergence is found with mutations of chromatin modifiers. This study shows the role of epigenetic marks in regulating gene expression and functional divergence after gene duplication in plants based on sequencing data.

Excision Dominates Pseudogenization During Fractionation After Whole Genome Duplication and in Gene Loss After Speciation in Plants

We take advantage of synteny blocks, the analytical construct enabled at the evolutionary moment of speciation or polyploidization, to follow the independent loss of duplicate genes in two sister species or the loss through fractionation of syntenic paralogs in a doubled genome. By examining how much sequence remains after a contiguous series of genes is deleted, we find that this residue remains at a constant low level independent of how many genes are lost—there are few if any relics of the missing sequence. Pseudogenes are rare or extremely transient in this context. The potential exceptions lie exclusively with a few examples of speciation, where the synteny blocks in some larger genomes tolerate degenerate sequence during genomic divergence of two species, but not after whole genome doubling in the same species where fractionation pressure eliminates virtually all non-coding sequence.

Learning Retention Mechanisms and Evolutionary Parameters of Duplicate Genes from Their Expression Data

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication.