Human-specific nonsense mutations identified by genome sequence comparisons

Association of Polygenetic Risk Scores Related to Immunity and Inflammation with Hyperthyroidism Risk and Interactions between the Polygenetic Scores and Dietary Factors in a Large Cohort

Graves's disease and thyroiditis induce hyperthyroidism, the causes of which remain unclear, although they are involved with genetic and environmental factors. We aimed to evaluate polygenetic variants for hyperthyroidism risk and their interaction with metabolic parameters and nutritional intakes in an urban hospital-based cohort. A genome-wide association study (GWAS) of participants with (cases; n = 842) and without (controls, n = 38,799) hyperthyroidism was used to identify and select genetic variants. In clinical and lifestyle interaction with PRS, 312 participants cured of hyperthyroidism were excluded. Single nucleotide polymorphisms (SNPs) associated with gene-gene interactions were selected by hyperthyroidism generalized multifactor dimensionality reduction. Polygenic risk scores (PRSs) were generated by summing the numbers of selected SNP risk alleles. The best gene-gene interaction model included tumor-necrosis factor (TNF)_rs1800610, mucin 22 (MUC22)_rs1304322089, tribbles pseudokinase 2 (TRIB2)_rs1881145, cytotoxic T-lymphocyte-associated antigen 4 (CTLA4)_rs231775, lipoma-preferred partner (LPP)_rs6780858, and human leukocyte antigen (HLA)-J_ rs767861647. The PRS of the best model was positively associated with hyperthyroidism risk by 1.939-fold (1.317-2.854) after adjusting for covariates. PRSs interacted with age, metabolic syndrome, and dietary inflammatory index (DII), while hyperthyroidism risk interacted with energy, calcium, seaweed, milk, and coffee intake (P < 0.05). The PRS impact on hyperthyroidism risk was observed in younger (<55 years) participants and adults without metabolic syndrome. PRSs were positively associated with hyperthyroidism risk in participants with low dietary intakes of energy (OR = 2.74), calcium (OR = 2.84), seaweed (OR = 3.43), milk (OR = 2.91), coffee (OR = 2.44), and DII (OR = 3.45). In conclusion, adults with high PRS involved in inflammation and immunity had a high hyperthyroidism risk exacerbated under low intakes of energy, calcium, seaweed, milk, or coffee. These results can be applied to personalized nutrition in a clinical setting.

Specificity of the Metallothionein-1 Response by Cadmium-Exposed Normal Human Urothelial Cells

Occupational and environmental exposure to cadmium is associated with the development of urothelial cancer. The metallothionein (MT) family of genes encodes proteins that sequester metal ions and modulate physiological processes, including zinc homeostasis. Little is known about the selectivity of expression of the different MT isoforms. Here, we examined the effect of cadmium exposure on MT gene and isoform expression by normal human urothelial (NHU) cell cultures. Baseline and cadmium-induced MT gene expression was characterized by next-generation sequencing and RT-PCR; protein expression was assessed by Western blotting using isoform-specific antibodies. Expression of the zinc transporter-1 (SLC30A1) gene was also assessed. NHU cells displayed transcription of MT-2A, but neither MT-3 nor MT-4 genes. Most striking was a highly inducer-specific expression of MT-1 genes, with cadmium inducing transcription of MT-1A, MT-1G, MT-1H, and MT-1M. Whereas MT-1G was also induced by zinc and nickel ions and MT-1H by iron, both MT-1A and MT-1M were highly cadmium-specific, which was confirmed for protein using isoform-specific antibodies. Protein but not transcript endured post-exposure, probably reflecting sequestration. SLC30A1 transcription was also affected by cadmium ion exposure, potentially reflecting perturbation of intracellular zinc homeostasis. We conclude that human urothelium displays a highly inductive profile of MT-1 gene expression, with two isoforms identified as highly specific to cadmium, providing candidate transcript and long-lived protein biomarkers of cadmium exposure.

Nε-lysine acetylation in the endoplasmic reticulum - a novel cellular mechanism that regulates proteostasis and autophagy

Protein post-translational modifications (PTMs) take many shapes, have many effects and are necessary for cellular homeostasis. One of these PTMs, Nε-lysine acetylation, was thought to occur only in the mitochondria, cytosol and nucleus, but this paradigm was challenged in the past decade with the discovery of lysine acetylation in the lumen of the endoplasmic reticulum (ER). This process is governed by the ER acetylation machinery: the cytosol:ER-lumen acetyl-CoA transporter AT-1 (also known as SLC33A1), and the ER-resident lysine acetyltransferases ATase1 and ATase2 (also known as NAT8B and NAT8, respectively). This Review summarizes the more recent biochemical, cellular and mouse model studies that underscore the importance of the ER acetylation process in maintaining protein homeostasis and autophagy within the secretory pathway, and its impact on developmental and age-associated diseases.

MOXD2, a Gene Possibly Associated with Olfaction, Is Frequently Inactivated in Birds

Vertebrate MOXD2 encodes a monooxygenase DBH-like 2 protein that could be involved in neurotransmitter metabolism, potentially during olfactory transduction. Loss of MOXD2 in apes and whales has been proposed to be associated with evolution of olfaction in these clades. We analyzed 57 bird genomes to identify MOXD2 sequences and found frequent loss of MOXD2 in 38 birds. Among the 57 birds, 19 species appeared to have an intact MOXD2 that encoded a full-length protein; 32 birds had a gene with open reading frame-disrupting point mutations and/or exon deletions; and the remaining 6 species did not show any MOXD2 sequence, suggesting a whole-gene deletion. Notably, among 10 passerine birds examined, 9 species shared a common genomic deletion that spanned several exons, implying the gene loss occurred in a common ancestor of these birds. However, 2 closely related penguin species, each of which had an inactive MOXD2, did not share any mutation, suggesting an independent loss after their divergence. Distribution of the 38 birds without an intact MOXD2 in the bird phylogenetic tree clearly indicates that MOXD2 loss is widespread and independent in bird lineages. We propose that widespread MOXD2 loss in some bird lineages may be implicated in the evolution of olfactory perception in these birds.

Loss of gene function and evolution of human phenotypes

Humans have acquired many distinct evolutionary traits after the human-chimpanzee divergence. These phenotypes have resulted from genetic changes that occurred in the human genome and were retained by natural selection. Comparative primate genome analyses reveal that loss-of-function mutations are common in the human genome. Some of these gene inactivation events were revealed to be associated with the emergence of advantageous phenotypes and were therefore positively selected and fixed in modern humans (the "less-ismore" hypothesis). Representative cases of human gene inactivation and their functional implications are presented in this review. Functional studies of additional inactive genes will provide insight into the molecular mechanisms underlying acquisition of various human-specific traits.

CEP89 is required for mitochondrial metabolism and neuronal function in man and fly

It is estimated that the human mitochondrial proteome consists of 1000-1500 distinct proteins. The majority of these support the various biochemical pathways that are active in these organelles. Individuals with an oxidative phosphorylation disorder of unknown cause provide a unique opportunity to identify novel genes implicated in mitochondrial biology. We identified a homozygous deletion of CEP89 in a patient with isolated complex IV deficiency, intellectual disability and multisystemic problems. CEP89 is a ubiquitously expressed and highly conserved gene of unknown function. Immunocytochemistry and cellular fractionation experiments showed that CEP89 is present both in the cytosol and in the mitochondrial intermembrane space. Furthermore, we ascertained in vitro that downregulation of CEP89 resulted in a severe decrease in complex IV in-gel activity and altered mobility, suggesting that the complex is aberrantly formed. Two-dimensional BN-SDS gel analysis revealed that CEP89 associates with a high-molecular weight complex. Together, these data confirm a role for CEP89 in mitochondrial metabolism. In addition, we modeled CEP89 loss of function in Drosophila. Ubiquitous knockdown of fly Cep89 decreased complex IV activity and resulted in complete lethality. Furthermore, Cep89 is required for mitochondrial integrity, membrane depolarization and synaptic transmission of photoreceptor neurons, and for (sub)synaptic organization of the larval neuromuscular junction. Finally, we tested neuronal Cep89 knockdown flies in the light-off jump reflex habituation assay, which revealed its role in learning. We conclude that CEP89 proteins play an important role in mitochondrial metabolism, especially complex IV activity, and are required for neuronal and cognitive function across evolution.

Lysine acetylation in the lumen of the ER: a novel and essential function under the control of the UPR

The N(ε)-amino group of lysine residues can be transiently modified by the addition of an acetyl group. Recognized functions of N(ε)-lysine acetylation include regulation of activity, molecular stabilization and conformational assembly of a protein. For more than forty years lysine acetylation was thought to occur only in the cytosol and nucleus. Targets included cytoskeletal-associated proteins as well as transcription factors, histone proteins and proteins involved in DNA recombination and repair. However, in 2007 we reported that a type I membrane protein involved in the pathogenesis of Alzheimer's disease was transiently acetylated on the ε amino group of seven lysine residues while transiting along the secretory pathway. Surprisingly, the acetylation occurred in the lumen of the endoplasmic reticulum (ER) forcing us to reconsider old paradigms. Indeed, if lysine acetylation can occur in the lumen of the ER, then all the essential biochemical elements of the reaction must be available in the lumen of the organelle. Follow-up studies revealed the existence of ER-based acetyl-CoA:lysine acetyltransferases as well as a membrane transporter that translocates acetyl-CoA from the cytosol into the ER lumen. Large-scale proteomics showed that the list of substrates of the ER-based acetylation machinery includes both transiting and resident proteins. Finally, genetic studies revealed that this machinery is tightly linked to human diseases. Here, we describe these exciting findings as well as recent biochemical and cellular advances, and discuss possible impact on both human physiology and pathology.

Inactivation of the MSLNL gene encoding mesothelin-like protein during African great ape evolution

Loss of gene function is implicated in the emergence of novel phenotypes during organism evolution. Here, we report the inactivation of the MSLNL gene encoding mesothelin-like protein in African great ape evolution. Human MSLNL has a nonsense mutation in exon 10 and two polymorphic mutations: a frameshift in exon 3 and a nonsense codon in exon 8. The gorilla gene also shows multiple deleterious mutations, including a premature stop codon, a deletion, and a splice site mutation. Molecular evolutionary analysis indicated relaxed selection pressure on MSLNL in African great ape lineages, which suggested that MSLNL might have become inactivated before the divergence of human, chimpanzee and gorilla. The mouse Mslnl gene is highly expressed in olfactory epithelium and moderately expressed in several other tissues. We propose that the loss of MSLNL may be associated with the evolution of the olfactory system in African great apes including human.

Up-regulation of metallothionein gene expression in parkinsonian astrocytes

The role of glial cells in Parkinson's disease (PD) is unclear. We have previously reported a striking up-regulation of DnaJB6 heat shock protein in PD substantia nigra astrocytes. Whole genome transcriptome analysis also indicated increased expression of metallothionein genes in substantia nigra and cortex of sporadic PD cases. Metallothioneins are metal-binding proteins in the CNS that are released by astrocytes and associated with neuroprotection. Metallothionein expression was investigated in 18 PD cases and 15 non-PD controls using quantitative real-time polymerase chain reaction (qRT-PCR), in situ hybridisation (ISH) and immunocytochemistry (ICC). We observed a strong increase in the expression of metallothioneins MT1E, MT1F, MT1G, MT1H, MT1M, MT1X and MT2A in both PD nigra and frontal cortex. Expression of LRP2 (megalin), the neuronal metallothionein receptor was also significantly increased. qRT-PCR confirmed metallothionein up-regulation. Astrocytes were found to be the main source of metallothioneins 1 and 2 based on ISH results, and this finding was confirmed by ICC. Our findings demonstrate metallothionein expression by reactive astrocytes in PD nigra and support a neuroprotective role for these cells. The traditional view that nigral astrocytes are non-reactive in PD is clearly incorrect. However, it is possible that astrocytes are themselves affected by the disease process which may explain their comparatively modest and previously overlooked response.

Identification of novel phosphorylation modification sites in human proteins that originated after the human-chimpanzee divergence

MOTIVATION

Phosphorylation modifications of specific protein residues are involved in a wide range of biological processes such as modulation of intracellular signal networks. Here, we present the development and application of a bioinformatics procedure for systematic identification of human-specific phosphorylation sites in proteins that may have occurred after the human-chimpanzee divergence.

RESULTS

We collected annotated human phosphorylation sites and compared each site to orthologous mammalian proteins across taxa including chimpanzee, orangutan, rhesus macaque, marmoset, mouse, dog, cow, elephant, opossum and platypus. We identified 37 human-specific gains of annotated phosphorylation sites in 35 proteins: 22 serines, 12 threonines and 3 tyrosines. The novel phosphorylation sites are situated in highly conserved segments of the protein. Proteins with novel phosphorylation sites are involved in crucial biological processes such as cell division (AURKB, CASC5, MKI67 and PDCD4) and chromatin remodeling (HIRA, HIRIP3, HIST1H1T, NAP1L4 and LRWD1). Modified phosphorylatable residues produce novel target sites for protein kinases such as cyclin-dependent kinases and casein kinases, possibly resulting in rewiring and fine-tuning of phosphorylation regulatory networks. The potential human-specific phosphorylation sites identified in this study are useful as candidates for functional analysis to identify novel phenotypes in humans.

CONTACT

hahny@cau.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gains, losses and changes of function after gene duplication: study of the metallothionein family

Metallothioneins (MT) are small proteins involved in heavy metal detoxification and protection against oxidative stress and cancer. The mammalian MT family originated through a series of duplication events which generated four major genes (MT1 to MT4). MT1 and MT2 encode for ubiquitous proteins, while MT3 and MT4 evolved to accomplish specific roles in brain and epithelium, respectively. Herein, phylogenetic, transcriptional and polymorphic analyses are carried out to expose gains, losses and diversification of functions that characterize the evolutionary history of the MT family. The phylogenetic analyses show that all four major genes originated through a single duplication event prior to the radiation of mammals. Further expansion of the MT1 gene has occurred in the primate lineage reaching in humans a total of 13 paralogs, five of which are pseudogenes. In humans, the reading frame of all five MT1 pseudogenes is reconstructed by sequence homology with a functional duplicate revealing that loss of invariant cysteines is the most frequent event accounting for pseudogeneisation. Expression analyses based on EST counts and RT-PCR experiments show that, as for MT1 and MT2, human MT3 is also ubiquitously expressed while MT4 transcripts are present in brain, testes, esophagus and mainly in thymus. Polymorphic variation reveals two deleterious mutations (Cys30Tyr and Arg31Trp) in MT4 with frequencies reaching about 30% in African and Asian populations suggesting the gene is inactive in some individuals and physiological compensation for its loss must arise from a functional equivalent. Altogether our findings provide novel data on the evolution and diversification of MT gene duplicates, a valuable resource for understanding the vast set of biological processes in which these proteins are involved.

Identification of human-specific transcript variants induced by DNA insertions in the human genome

MOTIVATION

Many genes in the human genome produce a wide variety of transcript variants resulting from alternative exon splicing, differential promoter usage, or altered polyadenylation site utilization that may function differently in human cells. Here, we present a bioinformatics method for the systematic identification of human-specific novel transcript variants that might have arisen after the human-chimpanzee divergence.

RESULTS

The procedure involved collecting genomic insertions that are unique to the human genome when compared with orthologous chimpanzee and rhesus macaque genomic regions, and that are expressed in the transcriptome as exons evidenced by mRNAs and/or expressed sequence tags (ESTs). Using this procedure, we identified 112 transcript variants that are specific to humans; 74 were associated with known genes and the remaining transcripts were located in unannotated genomic loci. The original source of inserts was mostly transposable elements including L1, Alu, SVA, and human endogenous retroviruses (HERVs). Interestingly, some non-repetitive genomic segments were also involved in the generation of novel transcript variants. Insert contributions to the transcripts included promoters, terminal exons and insertions in exons, splice donors and acceptors and complete exon cassettes. Comparison of personal genomes revealed that at least seven loci were polymorphic in humans. The exaptation of human-specific genomic inserts as novel transcript variants may have increased human gene versatility or affected gene regulation.

Divergence, demography and gene loss along the human lineage

Genomic DNA sequences are an irreplaceable source for reconstructing the vanished past of living organisms. Based on updated sequence data, this paper summarizes our studies on species divergence time, ancient population size and functional loss of genes in the primate lineage leading to modern humans (Homo sapiens sapiens). The inter- and intraspecific comparisons of DNA sequences suggest that the human lineage experienced a rather severe bottleneck in the Middle Pleistocene, throughout which period the subdivided African population played a predominant role in shaping the genetic architecture of modern humans. Also, published and newly identified human-specific pseudogenes (HSPs) are enumerated in order to infer their significance for human evolution. Of the 121 candidate genes obtained, authentic HSPs turn out to comprise only 25 olfactory receptor genes, four T cell receptor genes and nine other genes. The fixation of HSPs has been too rare over the past 6–7 Myr to account for species differences between humans and chimpanzees.

An exon-based comparative variant analysis pipeline to study the scale and role of frameshift and nonsense mutation in the human-chimpanzee divergence

Chimpanzees and humans are closely related but differ in many deadly human diseases and other characteristics in physiology, anatomy, and pathology. In spite of decades of extensive research, crucial questions about the molecular mechanisms behind the differences are yet to be understood. Here I report ExonVar, a novel computational pipeline for Exon-based human-chimpanzee comparative Variant analysis. The objective is to comparatively analyze mutations specifically those that caused the frameshift and nonsense mutations and to assess their scale and potential impacts on human-chimpanzee divergence. Genomewide analysis of human and chimpanzee exons with ExonVar identified a number of species-specific, exon-disrupting mutations in chimpanzees but much fewer in humans. Many were found on genes involved in important biological processes such as T cell lineage development, the pathogenesis of inflammatory diseases, and antigen induced cell death. A “less-is-more” model was previously established to illustrate the role of the gene inactivation and disruptions during human evolution. Here this analysis suggested a different model where the chimpanzee-specific exon-disrupting mutations may act as additional evolutionary force that drove the human-chimpanzee divergence. Finally, the analysis revealed a number of sequencing errors in the chimpanzee and human genome sequences and further illustrated that they could be corrected without resequencing.

Genome-wide analysis of chimpanzee genes with premature termination codons

Background

Premature termination codons (PTCs) cause mRNA degradation or a truncated protein and thereby contribute to the transcriptome and proteome divergence between species. Here we present the first genome-wide study of PTCs in the chimpanzee. By comparing the human and chimpanzee genome sequences we identify and characterize genes with PTCs, in order to understand the contribution of these mutations to the transcriptome diversity between the species.

Results

We have studied a total of 13,487 human-chimpanzee gene pairs and found that ~8% were affected by PTCs in the chimpanzee. A majority (764/1,109) of PTCs were caused by insertions or deletions and the remaining part was caused by substitutions. The distribution of PTC genes varied between chromosomes, with Y having the highest proportion. Furthermore, the density of PTC genes varied on a megabasepair scale within chromosomes and we found the density to be correlated both with indel divergence and proximity to the telomere. Within genes, PTCs were more common close to the 5' and 3' ends of the amino acid sequence. Gene Ontology classification revealed that olfactory receptor genes were over represented among the PTC genes.

Conclusion

Our results showed that the density of PTC genes fluctuated across the genome depending on the local genomic context. PTCs were preferentially located in the terminal parts of the transcript, which generally have a lower frequency of functional domains, indicating that selection was operating against PTCs at sites central to protein function. The enrichment of GO terms associated with olfaction suggests that PTCs may have influenced the difference in the repertoire of olfactory genes between humans and chimpanzees. In summary, 8% of the chimpanzee genes were affected by PTCs and this type of variation is likely to have an important effect on the transcript and proteomic divergence between humans and chimpanzees.

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.

N-acetyltransferase 8, a positional candidate for blood pressure and renal regulation: resequencing, association and in silico study

Background

Kidneys have an important function in blood pressure (BP) regulation and elevated BP may lead to kidney failure. Chr2p12-p13 region linked to BP traits in multiple studies harbours a potential candidate for BP and renal function, N-acetyltransferase 8 (NAT8) expressed in embryonic and adult kidney and associated with nephrotoxicity response.

Methods/Results

We report the first study exploring NAT8 as a potential candidate gene for blood pressure and kidney function. The resequencing (n = 42, random Estonian samples) identified 15 NAT8 polymorphisms, including 6 novel variants. The diversity of NAT8 5' upstream region (π/bp = 0.00320) exceeded up to 10 times the variation in the NAT8 genic region (π/bp = 0.00037) as well as the average variation (π/bp = 0.00040) for the promoters of 29 reference genes associated with hypertension. We suggest that a potential source for such high variation could be an active gene conversion process from NAT8B duplicate gene to NAT8. Similarly to NAT8, several reference genes with the most variable upstream regions have also duplicate copies. The NAT8 promoter SNPs were targeted with pilot quantitative association studies for blood pressure (n = 137, healthy unrelated individuals) and for the index of kidney function – estimated glomerular filtration rate (eGFR; n = 157 hypertensives with and without nephropathy). Minor alleles of these polymorphisms revealed a significant protective effect against elevated systolic BP as well as kidney failure in hypertension patients (p < 0.05; linear regression model, addictive effect).

Conclusion

The full resequencing and pilot association study of a novel positional candidate gene for blood pressure and renal function, human N-acetyltransferase 8, suggested a contribution of highly variable NAT8 promoter polymorphisms in determination of systolic blood pressure and eGFR. Based on in silico analysis, we raise the hypothesis that the alternative SNP alleles of the NAT8 upstream region may have differential effect on gene expression.

Large-scale identification of human genes implicated in epidermal barrier function

Identification of genes expressed in epidermal granular keratinocytes by ORESTES, including a number that are highly specific for these cells.

Background

During epidermal differentiation, keratinocytes progressing through the suprabasal layers undergo complex and tightly regulated biochemical modifications leading to cornification and desquamation. The last living cells, the granular keratinocytes (GKs), produce almost all of the proteins and lipids required for the protective barrier function before their programmed cell death gives rise to corneocytes. We present here the first analysis of the transcriptome of human GKs, purified from healthy epidermis by an original approach.

Results

Using the ORESTES method, 22,585 expressed sequence tags (ESTs) were produced that matched 3,387 genes. Despite normalization provided by this method (mean 4.6 ORESTES per gene), some highly transcribed genes, including that encoding dermokine, were overrepresented. About 330 expressed genes displayed less than 100 ESTs in UniGene clusters and are most likely to be specific for GKs and potentially involved in barrier function. This hypothesis was tested by comparing the relative expression of 73 genes in the basal and granular layers of epidermis by quantitative RT-PCR. Among these, 33 were identified as new, highly specific markers of GKs, including those encoding a protease, protease inhibitors and proteins involved in lipid metabolism and transport. We identified filaggrin 2 (also called ifapsoriasin), a poorly characterized member of the epidermal differentiation complex, as well as three new lipase genes clustered with paralogous genes on chromosome 10q23.31. A new gene of unknown function, C1orf81, is specifically disrupted in the human genome by a frameshift mutation.

Conclusion

These data increase the present knowledge of genes responsible for the formation of the skin barrier and suggest new candidates for genodermatoses of unknown origin.

Inactivation of MOXD2 and S100A15A by exon deletion during human evolution

We devised a bioinformatics method for systematic identification of putative human-specific exon-deletion mutations that occurred after the divergence of human and chimpanzee and experimentally verified 2 of the predicted mutations in MOXD2 and S100A15A genes. MOXD2 gene encodes a monooxygenase that is highly conserved in mammals and is mostly expressed in the olfactory epithelium in mouse. The presence of a deletion of the last 2 exons and a polymorphic nonsense mutation in exon 6 suggests that MOXD2 gene is inactive in humans. S100A15A is a member of the S100 family of calcium-binding proteins, the mouse ortholog of which is expressed during epidermal maturation. Human S100A15A gene is likely to be inactive because the start codon-bearing exon is deleted in human. We propose that modification or inactivation of MOXD2 and S100A15A genes have contributed to the loss of certain smell sense in humans and to the development of human skin.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.