Lengths of coding and noncoding regions of a gene correlate with gene essentiality and rates of evolution

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

Contradiction in Star-Allele Nomenclature of Pharmacogenes between Common Haplotypes and Rare Variants

The nomenclature of star alleles has been widely used in pharmacogenomics to enhance treatment outcomes, predict drug response variability, and reduce adverse reactions. However, the discovery of numerous rare functional variants through genome sequencing introduces complexities into the star-allele system. This study aimed to assess the nature and impact of the rapid discovery of numerous rare functional variants in the traditional haplotype-based star-allele system. We developed a new method to construct haplogroups, representing a common ancestry structure, by iteratively excluding rare and functional variants of the 25 representative pharmacogenes using the 2504 genomes from the 1000 Genomes Project. In total, 192 haplogroups and 288 star alleles were identified, with an average of 7.68 ± 4.2 cross-ethnic haplogroups per gene. Most of the haplogroups (70.8%, 136/192) were highly aligned with their corresponding classical star alleles (VI = 1.86 ± 0.78), exhibiting higher genetic diversity than the star alleles. Approximately 41.3% (N = 119) of the star alleles in the 2504 genomes did not belong to any of the haplogroups, and most of them (91.3%, 105/116) were determined by a single variant according to the allele-definition table provided by CPIC. These functional single variants had low allele frequency (MAF < 1%), high evolutionary conservation, and variant deleteriousness, which suggests significant negative selection. It is suggested that the traditional haplotype-based naming system for pharmacogenetic star alleles now needs to be adjusted by balancing both traditional haplotyping and newly emerging variant-sequencing approaches to reduce naming complexity.

Genome-wide analysis of glyoxalase-like gene families in grape (Vitis vinifera L.) and their expression profiling in response to downy mildew infection

BACKGROUND

The glyoxalase system usually comprises two enzymes, glyoxalase I (GLYI) and glyoxalase II (GLYII). This system converts cytotoxic methylglyoxal (MG) into non-toxic D-lactate in the presence of reduced glutathione (GSH) in two enzymatic steps. Recently, a novel type of glyoxalase III (GLYIII) activity has observed in Escherichia coli that can detoxify MG into D-lactate directly, in one step, without a cofactor. Investigation of the glyoxalase enzymes of a number of plant species shows the importance of their roles in response both to abiotic and to biotic stresses. Until now, glyoxalase gene families have been identified in the genomes of four plants, Arabidopsis, Oryza sativa, Glycine max and Medicago truncatula but no similar study has been done with the grapevine Vitis vinifera L.

RESULTS

In this study, four GLYI-like, two GLYII-like and three GLYIII-like genes are identified from the genome database of grape. All these genes were analysed in detail, including their chromosomal locations, phylogenetic relationships, exon-intron distributions, protein domain organisations and the presence of conserved binding sites. Using quantitative real-time PCR analysis (qRT-PCR), the expression profiles of these genes were analysed in different tissues of grape, and also when under infection stress from downy mildew (Plasmopara viticola). The study reveals that most VvGLY-like genes had higher expressions in stem, leaf, tendril and ovule but lower expressions in the flower. In addition, most of the VvGLY-like gene members were P. viticola responsive with high expressions 6-12 h and 96-120 h after inoculation. However, VvGLYI-like1 was highly expressed 48 h after inoculation, similar to VvPR1 and VvNPR1 which are involved in the defence response.

CONCLUSIONS

This study identified the GLYI-like, GLYII-like and GLYIII-like full gene families of the grapevine. Based on a phylogenetic analysis and the presence of conserved binding sites, we speculate that these glyoxalase-like genes in grape encode active glyoxalases. Moreover, our study provides a basis for discussing the roles of VvGLYI-like, VvGLYII-like and VvGLYIII-like genes in grape's response to downy mildew infection. Our results shed light on the selection of candidate genes for downy mildew tolerance in grape and lay the foundation for further functional investigations of these glyoxalase genes.

Correlates of evolutionary rates in the murine sperm proteome

Background

Protein-coding genes expressed in sperm evolve at different rates. To gain deeper insight into the factors underlying this heterogeneity we examined the relative importance of a diverse set of previously described rate correlates in determining the evolution of murine sperm proteins.

Results

Using partial rank correlations we detected several major rate indicators: Phyletic gene age, numbers of protein-protein interactions, and survival essentiality emerged as particularly important rate correlates in murine sperm proteins. Tissue specificity, numbers of paralogs, and untranslated region lengths also correlate significantly with sperm genes’ evolutionary rates, albeit to a lesser extent. Multifunctionality, coding sequence or average intron lengths, and mean expression level have insignificant or virtually no independent effects on evolutionary rates in murine sperm genes. Gene ontology enrichment analyses of three equally sized murine sperm protein groups classified based on their evolutionary rates indicate strongest sperm-specific functional specialization in the most quickly evolving gene class.

Conclusions

We propose a model according to which slowly evolving murine sperm proteins tend to be constrained by factors such as survival essentiality, network connectivity, and/or broad expression. In contrast, evolutionary change may arise especially in less constrained sperm proteins, which might, moreover, be prone to specialize to reproduction-related functions. Our results should be taken into account in future studies on rate variations of reproductive genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1157-6) contains supplementary material, which is available to authorized users.

Characterization of OglDREB2A gene from African rice (Oryza glaberrima), comparative analysis and its transcriptional regulation under salinity stress

In this study, AP2 DNA-binding domain-containing transcription factor, OglDREB2A, was cloned from the African rice (Oryza glaberrima) and compared with 3000 rice genotypes. Further, the phylogenetic and various structural analysis was performed using in silico approaches. Further, to understand its allelic variation in rice, SNPs and indels were detected among the 3000 rice genotypes which indicated that while coding region is highly conserved, yet noncoding regions such as UTR and intron contained most of the variation. Phylogenetic analysis of the OglDREB2A sequence in different Oryza as well as in diverse eudicot species revealed that DREB from various Oryza species were diversed much earlier than other genes. Further, structural features and in silico analyses provided insights into different properties of OglDREB2A protein. The neutrality test on the coding region of OglDREB2A from different genotypes of O. glaberrima showed the lack of selection in this gene. Among the different developmental stages, it was upregulated at tillering and flag leaf under salinity treatment indicating its positive role in seedling and reproductive stage tolerance. Real-time PCR analysis also indicated the conserve expression pattern of this gene under salinity stress across the three different Oryza species having different degree of salinity tolerance.

Disease-associated mitochondrial mutations and the evolution of primate mitogenomes

Several human diseases have been associated with mutations in mitochondrial genes comprising a set of confirmed and reported mutations according to the MITOMAP database. An analysis of complete mitogenomes across 139 primate species showed that most confirmed disease-associated mutations occurred in aligned codon positions and gene regions under strong purifying selection resulting in a strong evolutionary conservation. Only two confirmed variants (7.1%), coding for the same amino acids accounting for severe human diseases, were identified without apparent pathogenicity in non-human primates, like the closely related Bornean orangutan. Conversely, reported disease-associated mutations were not especially concentrated in conserved codon positions, and a large fraction of them occurred in highly variable ones. Additionally, 88 (45.8%) of reported mutations showed similar variants in several non-human primates and some of them have been present in extinct species of the genus Homo. Considering that recurrent mutations leading to persistent variants throughout the evolutionary diversification of primates are less likely to be severely damaging to fitness, we suggest that these 88 mutations are less likely to be pathogenic. Conversely, 69 (35.9%) of reported disease-associated mutations occurred in extremely conserved aligned codon positions which makes them more likely to damage the primate mitochondrial physiology.

Essentiality Is a Strong Determinant of Protein Rates of Evolution during Mutation Accumulation Experiments in Escherichia coli

The Neutral Theory of Molecular Evolution is considered the most powerful theory to understand the evolutionary behavior of proteins. One of the main predictions of this theory is that essential proteins should evolve slower than dispensable ones owing to increased selective constraints. Comparison of genomes of different species, however, has revealed only small differences between the rates of evolution of essential and nonessential proteins. In some analyses, these differences vanish once confounding factors are controlled for, whereas in other cases essentiality seems to have an independent, albeit small, effect. It has been argued that comparing relatively distant genomes may entail a number of limitations. For instance, many of the genes that are dispensable in controlled lab conditions may be essential in some of the conditions faced in nature. Moreover, essentiality can change during evolution, and rates of protein evolution are simultaneously shaped by a variety of factors, whose individual effects are difficult to isolate. Here, we conducted two parallel mutation accumulation experiments in Escherichia coli, during 5,500–5,750 generations, and compared the genomes at different points of the experiments. Our approach (a short-term experiment, under highly controlled conditions) enabled us to overcome many of the limitations of previous studies. We observed that essential proteins evolved substantially slower than nonessential ones during our experiments. Strikingly, rates of protein evolution were only moderately affected by expression level and protein length.

Expression Divergence Is Correlated with Sequence Evolution but Not Positive Selection in Conifers

The evolutionary and genomic determinants of sequence evolution in conifers are poorly understood, and previous studies have found only limited evidence for positive selection. Using RNAseq data, we compared gene expression profiles to patterns of divergence and polymorphism in 44 seedlings of lodgepole pine (Pinus contorta) and 39 seedlings of interior spruce (Picea glauca × engelmannii) to elucidate the evolutionary forces that shape their genomes and their plastic responses to abiotic stress. We found that rapidly diverging genes tend to have greater expression divergence, lower expression levels, reduced levels of synonymous site diversity, and longer proteins than slowly diverging genes. Similar patterns were identified for the untranslated regions, but with some exceptions. We found evidence that genes with low expression levels had a larger fraction of nearly neutral sites, suggesting a primary role for negative selection in determining the association between evolutionary rate and expression level. There was limited evidence for differences in the rate of positive selection among genes with divergent versus conserved expression profiles and some evidence supporting relaxed selection in genes diverging in expression between the species. Finally, we identified a small number of genes that showed evidence of site-specific positive selection using divergence data alone. However, estimates of the proportion of sites fixed by positive selection (α) were in the range of other plant species with large effective population sizes suggesting relatively high rates of adaptive divergence among conifers.