Chungoo Park

Differential requirements for mRNA folding partially explain why highly expressed proteins evolve slowly

Significance The expression level of a gene is a leading determinant of its rate of protein sequence evolution, but the underlying mechanisms are unclear. We show that as the mRNA concentration increases, natural selection for mRNA folding intensifies, resulting in larger fractions of mutations deleterious to mRNA folding and lower rates of protein evolution. Counterintuitively, selection for mRNA folding also impacts the nonsynonymous-to-synonymous nucleotide substitution rate ratio, requiring a revision of the current interpretation of this ratio as a measure of protein-level selection. These findings demonstrate a prominent role of selection at the mRNA level in molecular evolution. The cause of the tremendous among-protein variation in the rate of sequence evolution is a central subject of molecular evolution. Expression level has been identified as a leading determinant of this variation among genes encoded in the same genome, but the underlying mechanisms are not fully understood. We here propose and demonstrate that a requirement for stronger folding of more abundant mRNAs results in slower evolution of more highly expressed genes and proteins. Specifically, we show that: (i) the higher the expression level of a gene, the greater the selective pressure for its mRNA to fold; (ii) random mutations are more likely to decrease mRNA folding when occurring in highly expressed genes than in lowly expressed genes; and (iii) amino acid substitution rate is negatively correlated with mRNA folding strength, with or without the control of expression level. Furthermore, synonymous (dS) and nonsynonymous (dN) nucleotide substitution rates are both negatively correlated with mRNA folding strength. However, counterintuitively, dS and dN are differentially constrained by selection for mRNA folding, resulting in a significant correlation between mRNA folding strength and dN/dS, even when gene expression level is controlled. The direction and magnitude of this correlation is determined primarily by the G+C frequency at third codon positions. Together, these findings explain why highly expressed genes evolve slowly, demonstrate a major role of natural selection at the mRNA level in constraining protein evolution, and reveal a previously unrecognized and unexpected form of nonprotein-level selection that impacts dN/dS.

High Expression Hampers Horizontal Gene Transfer

Horizontal gene transfer (HGT), the movement of genetic material from one species to another, is a common phenomenon in prokaryotic evolution. Although the rate of HGT is known to vary among genes, our understanding of the cause of this variation, currently summarized by two rules, is far from complete. The first rule states that informational genes, which are involved in DNA replication, transcription, and translation, have lower transferabilities than operational genes. The second rule asserts that protein interactivity negatively impacts gene transferability. Here, we hypothesize that high expression hampers HGT, because the fitness cost of an HGT to the recipient, arising from the 1) energy expenditure in transcription and translation, 2) cytotoxic protein misfolding, 3) reduction in cellular translational efficiency, 4) detrimental protein misinteraction, and 5) disturbance of the optimal protein concentration or cell physiology, increases with the expression level of the transferred gene. To test this hypothesis, we examined laboratory and natural HGTs to Escherichia coli. We observed lower transferabilities of more highly expressed genes, even after controlling the confounding factors from the two established rules and the genic GC content. Furthermore, expression level predicts gene transferability better than all other factors examined. We also confirmed the significant negative impact of gene expression on the rate of HGTs to 127 of 133 genomes of eubacteria and archaebacteria. Together, these findings establish the gene expression level as a major determinant of horizontal gene transferability. They also suggest that most successful HGTs are initially slightly deleterious, fixed because of their negligibly low costs rather than high benefits to the recipient.

Genomic evidence for elevated mutation rates in highly expressed genes

Reporter gene assays have demonstrated both transcription‐associated mutagenesis (TAM) and transcription‐coupled repair, but the net impact of transcription on mutation rate remains unclear, especially at the genomic scale. Using comparative genomics of related species as well as mutation accumulation lines, we show in yeast that the rate of point mutation in a gene increases with the expression level of the gene. Transcription induces mutagenesis on both DNA strands, indicating simultaneous actions of several TAM mechanisms. A significant positive correlation is also detected between the human germline mutation rate and expression level. These results indicate that transcription is overall mutagenic.

Strong Purifying Selection at Genes Escaping X Chromosome Inactivation

To achieve dosage balance of X-linked genes between mammalian males and females, one female X chromosome becomes inactivated. However, approximately 15% of genes on this inactivated chromosome escape X chromosome inactivation (XCI). Here, using a chromosome-wide analysis of primate X-linked orthologs, we test a hypothesis that such genes evolve under a unique selective pressure. We find that escape genes are subject to stronger purifying selection than inactivated genes and that positive selection does not significantly affect the evolution of these genes. The strength of selection does not differ between escape genes with similar versus different expression levels in males versus females. Intriguingly, escape genes possessing Y homologs evolve under the strongest purifying selection. We also found evidence of stronger conservation in gene expression levels in escape than inactivated genes. We hypothesize that divergence in function and expression between X and Y gametologs is driving such strong purifying selection for escape genes.

Coding region structural heterogeneity and turnover of transcription start sites contribute to divergence in expression between duplicate genes.

BackgroundGene expression divergence is one manifestation of functional differences between duplicate genes. Although rapid accumulation of expression divergence between duplicate gene copies has been observed, the driving mechanisms behind this phenomenon have not been explored in detail.ResultsWe examine which factors influence expression divergence between human duplicate genes, utilizing the latest genome-wide data sets. We conclude that the turnover of transcription start sites between duplicate genes occurs rapidly after gene duplication and that gene pairs with shared transcription start sites have significantly higher expression similarity than those without shared transcription start sites. Moreover, we find that most (55%) duplicate gene pairs do not retain the same coding sequence structure between the two duplicate copies and this also contributes to divergence in their expression. Furthermore, the proportion of aligned sequences in cis-regulatory regions between the two copies is positively correlated with expression similarity. Surprisingly, we find no effect of copy-specific transposable element insertions on the divergence of duplicate gene expression.ConclusionsOur results suggest that turnover of transcription start sites, structural heterogeneity of coding sequences, and divergence of cis-regulatory regions between copies play a pivotal role in determining the expression divergence of duplicate genes.

Human papillomavirus genotyping using HPV DNA chip analysis in Korean women.

This study was designed to investigate the genotypes of human papillomavirus (HPV) in Korean women who had abnormal cervical cytology and to evaluate the clinical accuracy of HPV DNA chip analysis for the diagnosis of cervical neoplasia. Liquid-based cytology preparations, HPV DNA chip analysis, and cervical biopsy were performed in 2358 women. High-risk HPV was identified in 23.5% of 1650 histologically confirmed normal samples (including cervicitis and squamous metaplasia) and in 81.8% of 708 samples with cervical intraepithelial neoplasia (CIN) and carcinoma (P< 0.01). The major prevalent high-risk HPV genotypes in 381 samples of CIN II/III were HPV-16, -58, -33, and -31, in order of prevalence rate (average overall, 78.0%), and HPV-16, -18, -58, and -33 (average overall, 81.2%) in 133 samples of squamous cell carcinoma (SCC). The infection rate of HPV-16 was significantly higher than that of other high-risk HPV genotypes in all normal, CIN, and SCC cases (P< 0.01) and increased with more advanced squamous cervical lesions (P< 0.01). The detection accuracy of high-risk HPV using HPV DNA chip analysis for CIN II or worse was as follows: sensitivity 84% (81–87%), specificity 72% (70–74%), positive predictive value 47% (44–50%), and negative predictive value 94% (92–95%). These results suggest that HPV DNA chip analysis may be a reliable diagnostic tool for the detection of cervical neoplasia and that there are geographic differences in the distribution of high-risk HPV genotypes.

Genome-wide evolutionary conservation of N-glycosylation sites

Although posttranslational protein modifications are generally thought to perform important cellular functions, recent studies showed that a large fraction of phosphorylation sites are not evolutionarily conserved. Whether the same is true for other protein modifications, such as N-glycosylation is an open question. N-glycosylation is a form of cotranslational and posttranslational modification that occurs by enzymatic addition of a polysaccharide, or glycan, to an asparagine (N) residue of a protein. Examining a large set of experimentally determined mouse N-glycosylation sites, we find that the evolutionary rate of glycosylated asparagines is significantly lower than that of nonglycosylated asparagines of the same proteins. We further confirm that the conservation of glycosylated asparagines is accompanied by the conservation of the canonical motif sequence for glycosylation, suggesting that the above substitution rate difference is related to glycosylation. Interestingly, when solvent accessibility is considered, the substitution rate disparity between glycosylated and nonglycosylated asparagines is highly significant at solvent accessible sites but not at solvent inaccessible sites. Thus, although the solvent inaccessible glycosylation sites were experimentally identified, they are unlikely to be genuine or physiologically important. For solvent accessible asparagines, our analysis reveals a widespread and strong functional constraint on glycosylation, unlike what has been observed for phosphorylation sites in most studies, including our own analysis. Because the majority of N-glycosylation occurs at solvent accessible sites, our results show an overall functional importance for N-glycosylation.

Draft genome of the sea cucumber Apostichopus japonicus and genetic polymorphism among color variants

Abstract The Japanese sea cucumber (Apostichopus japonicus Selenka 1867) is an economically important species as a source of seafood and ingredient in traditional medicine. It is mainly found off the coasts of northeast Asia. Recently, substantial exploitation and widespread biotic diseases in A. japonicus have generated increasing conservation concern. However, the genomic knowledge base and resources available for researchers to use in managing this natural resource and to establish genetically based breeding systems for sea cucumber aquaculture are still in a nascent stage. A total of 312 Gb of raw sequences were generated using the Illumina HiSeq 2000 platform and assembled to a final size of 0.66 Gb, which is about 80.5% of the estimated genome size (0.82 Gb). We observed nucleotide-level heterozygosity within the assembled genome to be 0.986%. The resulting draft genome assembly comprising 132 607 scaffolds with an N50 value of 10.5 kb contains a total of 21 771 predicted protein-coding genes. We identified 6.6–14.5 million heterozygous single nucleotide polymorphisms in the assembled genome of the three natural color variants (green, red, and black), resulting in an estimated nucleotide diversity of 0.00146. We report the first draft genome of A. japonicus and provide a general overview of the genetic variation in the three major color variants of A. japonicus. These data will help provide a comprehensive view of the genetic, physiological, and evolutionary relationships among color variants in A. japonicus, and will be invaluable resources for sea cucumber genomic research.

Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome of Rhigonema thysanophora (Nematoda: Chromadorea)

We determined the complete mitochondrial genome sequence of Rhigonema thysanophora, the first representative of Rhigonematomorpha, and used this sequence along with 57 other nematode species for phylogenetic analyses. The R. thysanophora mtDNA is 15 015 bp and identical to all other chromadorean nematode mtDNAs published to date in that it contains 36 genes (lacking atp8) encoded in the same direction. Phylogenetic analyses of nucleotide and amino acid sequence data for the 12 protein‐coding genes recovered Rhigonematomorpha as the sister group to the heterakoid species, Ascaridia columbae (Ascaridomorpha). The organization of R. thysanophora mtDNA resembles the most common pattern for the Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha clade in gene order, but with some substantial gene rearrangements. This similarity in gene order is in agreement with the sequence‐based analyses that indicate a close relationship between Rhigonematomorpha and Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha. These results are consistent with certain analyses of nuclear SSU rDNA for R. thysanophora and some earlier classification systems that asserted phylogenetic affinity between Rhigonematomorpha and Ascaridomorpha, but inconsistent with morphology‐based phylogenetic hypotheses that suggested a close (taxonomic) relationship between rhigonematomorphs and oxyuridomorphs (pinworms). These observations must be tempered by noting that few rhigonematomorph species have been sequenced and included in phylogenetic analyses, and preliminary studies based on SSU rDNA suggest the group is not monophyletic. Additional mitochondrial genome sequences of rhigonematids are needed to characterize their phylogenetic relationships within Chromadorea, and to increase understanding of mitochondrial genome evolution.

A computational approach to candidate gene prioritization for X-linked mental retardation using annotation-based binary filtering and motif-based linear discriminatory analysis

AbstractBackgroundSeveral computational candidate gene selection and prioritization methods have recently been developed. These in silico selection and prioritization techniques are usually based on two central approaches - the examination of similarities to known disease genes and/or the evaluation of functional annotation of genes. Each of these approaches has its own caveats. Here we employ a previously described method of candidate gene prioritization based mainly on gene annotation, in accompaniment with a technique based on the evaluation of pertinent sequence motifs or signatures, in an attempt to refine the gene prioritization approach. We apply this approach to X-linked mental retardation (XLMR), a group of heterogeneous disorders for which some of the underlying genetics is known.ResultsThe gene annotation-based binary filtering method yielded a ranked list of putative XLMR candidate genes with good plausibility of being associated with the development of mental retardation. In parallel, a motif finding approach based on linear discriminatory analysis (LDA) was employed to identify short sequence patterns that may discriminate XLMR from non-XLMR genes. High rates (>80%) of correct classification was achieved, suggesting that the identification of these motifs effectively captures genomic signals associated with XLMR vs. non-XLMR genes. The computational tools developed for the motif-based LDA is integrated into the freely available genomic analysis portal Galaxy (http://main.g2.bx.psu.edu/). Nine genes (APLN, ZC4H2, MAGED4, MAGED4B, RAP2C, FAM156A, FAM156B, TBL1X, and UXT) were highlighted as highly-ranked XLMR methods.ConclusionsThe combination of gene annotation information and sequence motif-orientated computational candidate gene prediction methods highlight an added benefit in generating a list of plausible candidate genes, as has been demonstrated for XLMR. Reviewers: This article was reviewed by Dr Barbara Bardoni (nominated by Prof Juergen Brosius); Prof Neil Smalheiser and Dr Dustin Holloway (nominated by Prof Charles DeLisi).