Segmental duplications in euchromatic regions of human chromosome 5: a source of evolutionary instability and transcriptional innovation

Comprehensive In Silico Analysis of Retrotransposon Insertions within the Survival Motor Neuron Genes Involved in Spinal Muscular Atrophy

Transposable elements (TEs) are interspersed repetitive and mobile DNA sequences within the genome. Better tools for evaluating TE-derived sequences have provided insights into the contribution of TEs to human development and disease. Spinal muscular atrophy (SMA) is an autosomal recessive motor neuron disease that is caused by deletions or mutations in the Survival Motor Neuron 1 (SMN1) gene but retention of its nearly perfect orthologue SMN2. Both genes are highly enriched in TEs. To establish a link between TEs and SMA, we conducted a comprehensive, in silico analysis of TE insertions within the SMN1/2 loci of SMA, carrier and healthy genomes. We found an Alu insertion in the promoter region and one L1 element in the 3'UTR that may play an important role in alternative promoter as well as in alternative transcriptional termination. Additionally, several intronic Alu repeats may influence alternative splicing via RNA circularization and causes the presence of new alternative exons. These Alu repeats present throughout the genes are also prone to recombination events that could lead to SMN1 exons deletions and, ultimately, SMA. TE characterization of the SMA genomic region could provide for a better understanding of the implications of TEs on human disease and genomic evolution.

Fusion Genes in Prostate Cancer: A Comparison in Men of African and European Descent

Simple Summary

Men of African origin have a 2–3 times greater chance of developing prostate cancer than those of European origin, and of patients that are diagnosed with the disease, men of African descent are 2 times more likely to die compared to white men. Men of African origin are still greatly underrepresented in genetic studies and clinical trials. This, unfortunately, means that new discoveries in cancer treatment are missing key information on the group with a greater chance of mortality. A fusion gene is a hybrid gene formed from two previously independent genes. Fusion genes have been found to be common in all main types of human cancer. The objective of this study was to increase our knowledge of fusion genes in prostate cancer using computational approaches and to compare fusion genes between men of African and European origin. This identified novel gene fusions unique to men of African origin and suggested that this group has a greater number of fusion genes.

Abstract

Prostate cancer is one of the most prevalent cancers worldwide, particularly affecting men living a western lifestyle and of African descent, suggesting risk factors that are genetic, environmental, and socioeconomic in nature. In the USA, African American (AA) men are disproportionately affected, on average suffering from a higher grade of the disease and at a younger age compared to men of European descent (EA). Fusion genes are chimeric products formed by the merging of two separate genes occurring as a result of chromosomal structural changes, for example, inversion or trans/cis-splicing of neighboring genes. They are known drivers of cancer and have been identified in 20% of cancers. Improvements in genomics technologies such as RNA-sequencing coupled with better algorithms for prediction of fusion genes has added to our knowledge of specific gene fusions in cancers. At present AA are underrepresented in genomic studies of prostate cancer. The primary goal of this study was to examine molecular differences in predicted fusion genes in a cohort of AA and EA men in the context of prostate cancer using computational approaches. RNA was purified from prostate tissue specimens obtained at surgery from subjects enrolled in the study. Fusion gene predictions were performed using four different fusion gene detection programs. This identified novel putative gene fusions unique to AA and suggested that the fusion gene burden was higher in AA compared to EA men.

Genomic Variability in the Survival Motor Neuron Genes (SMN1 and SMN2): Implications for Spinal Muscular Atrophy Phenotype and Therapeutics Development

Spinal muscular atrophy (SMA) is a leading genetic cause of infant death worldwide that is characterized by loss of spinal motor neurons leading to muscle weakness and atrophy. SMA results from the loss of survival motor neuron 1 (SMN1) gene but retention of its paralog SMN2. The copy numbers of SMN1 and SMN2 are variable within the human population with SMN2 copy number inversely correlating with SMA severity. Current therapeutic options for SMA focus on increasing SMN2 expression and alternative splicing so as to increase the amount of SMN protein. Recent work has demonstrated that not all SMN2, or SMN1, genes are equivalent and there is a high degree of genomic heterogeneity with respect to the SMN genes. Because SMA is now an actionable disease with SMN2 being the primary target, it is imperative to have a comprehensive understanding of this genomic heterogeneity with respect to hybrid SMN1–SMN2 genes generated by gene conversion events as well as partial deletions of the SMN genes. This review will describe this genetic heterogeneity in SMA and its impact on disease phenotype as well as therapeutic efficacy.

Temperature-sensitive spinal muscular atrophy-causing point mutations lead to SMN instability, locomotor defects and premature lethality in Drosophila

Spinal muscular atrophy (SMA) is the leading genetic cause of death in young children, arising from homozygous deletion or mutation of the survival motor neuron 1 (SMN1) gene. SMN protein expressed from a paralogous gene, SMN2, is the primary genetic modifier of SMA; small changes in overall SMN levels cause dramatic changes in disease severity. Thus, deeper insight into mechanisms that regulate SMN protein stability should lead to better therapeutic outcomes. Here, we show that SMA patient-derived missense mutations in the Drosophila SMN Tudor domain exhibit a pronounced temperature sensitivity that affects organismal viability, larval locomotor function and adult longevity. These disease-related phenotypes are domain specific and result from decreased SMN stability at elevated temperature. This system was utilized to manipulate SMN levels during various stages of Drosophila development. Owing to a large maternal contribution of mRNA and protein, Smn is not expressed zygotically during embryogenesis. Interestingly, we find that only baseline levels of SMN are required during larval stages, whereas high levels of the protein are required during pupation. This previously uncharacterized period of elevated SMN expression, during which the majority of adult tissues are formed and differentiated, could be an important and translationally relevant developmental stage in which to study SMN function. Taken together, these findings illustrate a novel in vivo role for the SMN Tudor domain in maintaining SMN homeostasis and highlight the necessity for high SMN levels at crucial developmental time points that are conserved from Drosophila to humans.

Comprehensive molecular screening strategy of OCLN in band-like calcification with simplified gyration and polymicrogyria

Occludin (OCLN) is an important component of the tight junction complex, providing apical intercellular connections between adjacent cells in endothelial and epithelial tissue. In 2010 O'Driscoll et al reported mutations in OCLN to cause band-like calcification with simplified gyration and polymicrogyria (BLC-PMG). BLC-PMG is a rare autosomal recessive syndrome, characterized by early onset seizures, progressive microcephaly, severe developmental delay and deep cortical gray matter and basal ganglia calcification with symmetrical, predominantly fronto-parietal, polymicrogyria. Here we report 4 additional cases of BLC-PMG with novel OCLN mutations, and provide a summary of the published mutational spectrum. More generally, we describe a comprehensive molecular screening strategy taking into account the technical challenges associated with the genetic architecture of OCLN, which include the presence of a pseudo-gene and copy number variants.

Monkey-based research on human disease: the implications of genetic differences

Assertions that the use of monkeys to investigate human diseases is valid scientifically are frequently based on a reported 90-93% genetic similarity between the species. Critical analyses of the relevance of monkey studies to human biology, however, indicate that this genetic similarity does not result in sufficient physiological similarity for monkeys to constitute good models for research, and that monkey data do not translate well to progress in clinical practice for humans. Salient examples include the failure of new drugs in clinical trials, the highly different infectivity and pathology of SIV/HIV, and poor extrapolation of research on Alzheimer's disease, Parkinson's disease and stroke. The major molecular differences underlying these inter-species phenotypic disparities have been revealed by comparative genomics and molecular biology - there are key differences in all aspects of gene expression and protein function, from chromosome and chromatin structure to post-translational modification. The collective effects of these differences are striking, extensive and widespread, and they show that the superficial similarity between human and monkey genetic sequences is of little benefit for biomedical research. The extrapolation of biomedical data from monkeys to humans is therefore highly unreliable, and the use of monkeys must be considered of questionable value, particularly given the breadth and potential of alternative methods of enquiry that are currently available to scientists.

Telomeric Region of the Spinal Muscular Atrophy Locus Is Susceptible to Structural Variations

BACKGROUND

Most patients with spinal muscular atrophy lack the survival motor neuron 1 gene (SMN1) in the telomeric region of the spinal muscular atrophy locus on chromosome 5q13. On the other hand, the copy number of SMN2, a centromeric homolog of SMN1, is increased in many of these patients. This study aimed to clarify the mechanism underlying these structural variations.

METHODS

We determined the copy numbers of telomeric and centromeric genes in the spinal muscular atrophy locus of 86 patients and 22 control subjects using multiplex ligation-dependent probe amplification analysis. Then, we chose 74 patients lacking SMN1 exons 7 and 8, and compared their dataset with that of 22 control subjects retaining SMN1 exons 7 and 8.

RESULTS

The SMN2 copy number was shown to vary widely and to correlate with the disease severity of the patients. Interestingly, telomeric NAIP and telomeric GTF2H2 showed similar tendencies. We also noted positive correlations among the copy number of SMN2 and the telomeric genes of the spinal muscular atrophy locus. However, the copy numbers of centromeric NAIP and centromeric GTF2H2 were stable among the patients, with both approximating a value of two.

CONCLUSION

Our findings suggested that the telomeric region of the spinal muscular atrophy locus appears to be susceptible to structural variation, whereas the centromeric region is stable. Moreover, according to our results, new SMN2 copies may be generated in the telomeric region of the spinal muscular atrophy locus, supporting the SMN1-to-SMN2 gene conversion theory.

Copy Number Variations in the Survival Motor Neuron Genes: Implications for Spinal Muscular Atrophy and Other Neurodegenerative Diseases

Proximal spinal muscular atrophy (SMA), a leading genetic cause of infant death worldwide, is an early-onset, autosomal recessive neurodegenerative disease characterized by the loss of spinal α-motor neurons. This loss of α-motor neurons is associated with muscle weakness and atrophy. SMA can be classified into five clinical grades based on age of onset and severity of the disease. Regardless of clinical grade, proximal SMA results from the loss or mutation of SMN1 (survival motor neuron 1) on chromosome 5q13. In humans a large tandem chromosomal duplication has lead to a second copy of the SMN gene locus known as SMN2. SMN2 is distinguishable from SMN1 by a single nucleotide difference that disrupts an exonic splice enhancer in exon 7. As a result, most of SMN2 mRNAs lack exon 7 (SMNΔ7) and produce a protein that is both unstable and less than fully functional. Although only 10–20% of the SMN2 gene product is fully functional, increased genomic copies of SMN2 inversely correlates with disease severity among individuals with SMA. Because SMN2 copy number influences disease severity in SMA, there is prognostic value in accurate measurement of SMN2 copy number from patients being evaluated for SMA. This prognostic value is especially important given that SMN2 copy number is now being used as an inclusion criterion for SMA clinical trials. In addition to SMA, copy number variations (CNVs) in the SMN genes can affect the clinical severity of other neurological disorders including amyotrophic lateral sclerosis (ALS) and progressive muscular atrophy (PMA). This review will discuss how SMN1 and SMN2 CNVs are detected and why accurate measurement of SMN1 and SMN2 copy numbers is relevant for SMA and other neurodegenerative diseases.

Tandem RNA chimeras contribute to transcriptome diversity in human population and are associated with intronic genetic variants

Chimeric RNAs originating from two or more different genes are known to exist not only in cancer, but also in normal tissues, where they can play a role in human evolution. However, the exact mechanism of their formation is unknown. Here, we use RNA sequencing data from 462 healthy individuals representing 5 human populations to systematically identify and in depth characterize 81 RNA tandem chimeric transcripts, 13 of which are novel. We observe that 6 out of these 81 chimeras have been regarded as cancer-specific. Moreover, we show that a prevalence of long introns at the fusion breakpoint is associated with the chimeric transcripts formation. We also find that tandem RNA chimeras have lower abundances as compared to their partner genes. Finally, by combining our results with genomic data from the same individuals we uncover intronic genetic variants associated with the chimeric RNA formation. Taken together our findings provide an important insight into the chimeric transcripts formation and open new avenues of research into the role of intronic genetic variants in post-transcriptional processing events.

Solving the puzzle of spinal muscular atrophy: what are the missing pieces?

Spinal muscular atrophy (SMA) is an autosomal recessive, lower motor neuron disease. Clinical heterogeneity is pervasive: three infantile (type I-III) and one adult-onset (type IV) forms are recognized. Type I SMA is the most common genetic cause of death in infancy and accounts for about 50% of all patients with SMA. Most forms of SMA are caused by mutations of the survival motor neuron (SMN1) gene. A second gene that is 99% identical to SMN1 (SMN2) is located in the same region. The only functionally relevant difference between the two genes identified to date is a C → T transition in exon 7 of SMN2, which determines an alternative spliced isoform that predominantly excludes exon 7. Thus, SMN2 genes do not produce sufficient full length SMN protein to prevent the onset of the disease. Since the identification of the causative mutation, biomedical research of SMA has progressed by leaps and bounds: from clues on the function of SMN protein, to the development of different models of the disease, to the identification of potential treatments, some of which are currently in human trials. The aim of this review is to elucidate the current state of knowledge, emphasizing how close we are to the solution of the puzzle that is SMA, and, more importantly, to highlight the missing pieces of this puzzle. Filling in these gaps in our knowledge will likely accelerate the development and delivery of efficient treatments for SMA patients and be a prerequisite towards achieving our final goal, the cure of SMA.

Genetic analysis of cardiac-specific transcription factors reveals novel insights into molecular causes of congenital heart disease

Heart development is complex and requires the sequential and timely interplay of regulatory master proteins, notably several transcription factors. Germline mutations in the human transcription factor genes have been associated with congenital heart disease (CHD), but familial cases studied so far have different mutations. There is no strict genotype-phenotype correlation and mutations in transcription factor genes are rare in unrelated patients. Most cases of CHD come from unaffected family members. The study of archived, but morphologically well-characterized malformed hearts for DNA alterations provides important clues regarding cardiogenic transcription factor genes that would lead to loss-of-function of the protein. Identification of tissue-restricted multiple mutations and multiple haplotypes suggests that somatic mutation and mosaicism are linked to cardiac anomalies. Altogether, somatic mutations and genomic instability in the diseased cardiac tissues of patients with CHD provide a novel mechanism of disease.

Two neighboring microdeletions of 5q13.2 in a child with oculo-auriculo-vertebral spectrum

We describe a patient with multiple congenital anomalies, including hemifacial microsomia, asymmetric macrostomia, dysplastic mandible, multiple preauricular tags, atresia of the external auricular canal, and vertebral anomalies, which coincide with oculo-auriculo-vertebral spectrum. G-banding ( approximately 850 band level) showed a normal 46, XY karyotype. A genome-wide screen for copy number variations (CNVs) using single nucleotide polymorphism (SNP) arrays revealed a 1Mb and a 167 kb deletion both on chromosome 5q13.2, which were absent in the parents and in 27 controls. Sixteen genes were located in the deleted region, including BIR1C and OCLN, which are involved in apoptosis. Haploinsufficiency of these genes may be contributing to the phenotype in this patient. To our knowledge, there are no previous reports of this 5q13.2 deletion in a patient with oculo-auriculo-vertebral spectrum.

Identification of human specific gene duplications relative to other primates by array CGH and quantitative PCR

In order to identify human lineage specific (HLS) copy number differences (CNDs) compared to other primates, we performed pair wise comparisons (human vs. chimpanzee, gorilla and orangutan) by using cDNA array comparative genomic hybridization (CGH). A set of 23 genes with HLS duplications were identified, as well as other lineage differences in gene copy number specific of chimpanzee, gorilla and orangutan. Each species has gained more copies of specific genes rather than losing gene copies. Eleven of the 23 genes have only been observed to have undergone HLS duplication in Fortna et al. (2004) and in the present study. Then, seven of these 11 genes were analyzed by quantitative PCR in chimpanzee, gorilla and orangutan, as well as in other six primate species (Hylobates lar, Cercopithecus aethiops, Papio hamadryas, Macaca mulatta, Lagothrix lagothricha, and Saimiri sciureus). Six genes confirmed array CGH data, and four of them appeared to have bona fide HLS duplications (ABCB10, E2F6, CDH12, and TDG genes). We propose that these gene duplications have a potential to contribute to specific human phenotypes.

The evolution of human segmental duplications and the core duplicon hypothesis

Duplicated sequences are important sources of genetic instability and in the evolution of new gene function within species. Hominids have a preponderance of intrachromosomal duplications organized in an interspersed fashion, as opposed to tandem duplications, which are common in other mammalian genomes such as mouse, dog, and cow. Multiple lines of evidence, including sequence divergence, comparative primate genomes, and fluorescence in situ hybridization (FISH) analyses, point to an excess of segmental duplications in the common ancestor of humans and African great apes. We find that much of the interspersed human duplication architecture within chromosomes is focused around common sequence elements referred to as "core duplicons." These cores correspond to the expansion of gene families, some of which show signatures of positive selection and lack orthologs present in other mammalian species. This genomic architecture predisposes apes and humans not only to extensive genetic diversity, but also to large-scale structural diversity mediated by nonallelic homologous recombination. In humans, many de novo large-scale genomic changes mediated by these duplications are associated with neuropsychiatric and neurodevelopmental disease. We propose that the disadvantage of a high rate of new mutations is offset by the selective advantage of newly minted genes within the cores.

Primate-specific spliced PMCHL RNAs are non-protein coding in human and macaque tissues

Background

Brain-expressed genes that were created in primate lineage represent obvious candidates to investigate molecular mechanisms that contributed to neural reorganization and emergence of new behavioural functions in Homo sapiens. PMCHL1 arose from retroposition of a pro-melanin-concentrating hormone (PMCH) antisense mRNA on the ancestral human chromosome 5p14 when platyrrhines and catarrhines diverged. Mutations before divergence of hylobatidae led to creation of new exons and finally PMCHL1 duplicated in an ancestor of hominids to generate PMCHL2 at the human chromosome 5q13. A complex pattern of spliced and unspliced PMCHL RNAs were found in human brain and testis.

Results

Several novel spliced PMCHL transcripts have been characterized in human testis and fetal brain, identifying an additional exon and novel splice sites. Sequencing of PMCHL genes in several non-human primates allowed to carry out phylogenetic analyses revealing that the initial retroposition event took place within an intron of the brain cadherin (CDH12) gene, soon after platyrrhine/catarrhine divergence, i.e. 30–35 Mya, and was concomitant with the insertion of an AluSg element. Sequence analysis of the spliced PMCHL transcripts identified only short ORFs of less than 300 bp, with low (VMCH-p8 and protein variants) or no evolutionary conservation. Western blot analyses of human and macaque tissues expressing PMCHL RNA failed to reveal any protein corresponding to VMCH-p8 and protein variants encoded by spliced transcripts.

Conclusion

Our present results improve our knowledge of the gene structure and the evolutionary history of the primate-specific chimeric PMCHL genes. These genes produce multiple spliced transcripts, bearing short, non-conserved and apparently non-translated ORFs that may function as mRNA-like non-coding RNAs.

Segmental duplications arise from Pol32-dependent repair of broken forks through two alternative replication-based mechanisms

The propensity of segmental duplications (SDs) to promote genomic instability is of increasing interest since their involvement in numerous human genomic diseases and cancers was revealed. However, the mechanism(s) responsible for their appearance remain mostly speculative. Here, we show that in budding yeast, replication accidents, which are most likely transformed into broken forks, play a causal role in the formation of SDs. The Pol32 subunit of the major replicative polymerase Polδ is required for all SD formation, demonstrating that SDs result from untimely DNA synthesis rather than from unequal crossing-over. Although Pol32 is known to be required for classical (Rad52-dependant) break-induced replication, only half of the SDs can be attributed to this mechanism. The remaining SDs are generated through a Rad52-independent mechanism of template switching between microsatellites or microhomologous sequences. This new mechanism, named microhomology/microsatellite-induced replication (MMIR), differs from all known DNA double-strand break repair pathways, as MMIR-mediated duplications still occur in the combined absence of homologous recombination, microhomology-mediated, and nonhomologous end joining machineries. The interplay between these two replication-based pathways explains important features of higher eukaryotic genomes, such as the strong, but not strict, association between SDs and transposable elements, as well as the frequent formation of oncogenic fusion genes generating protein innovations at SD junctions.

Duplications of long segments of chromosomes are frequently observed in multicellular organisms (∼5% of our genome, for instance). They appear as a fundamental trait of the recent genome evolution in great apes and are often associated with chromosomal instability, capable of increasing genetic polymorphism among individuals, but also having dramatic consequences as a source of diseases and cancer. Despite their importance, the molecular mechanisms of formation of segmental duplications remain unclear. Using a specifically designed experimental system in the baker's yeast Saccharomyces cerevisiae, hundreds of naturally occurring segmental duplications encompassing dozens of genes were selected. With the help of modern molecular methods coupled to detailed genetic analysis, we show that such duplication events are frequent and result from untimely DNA synthesis accidents produced by two distinct molecular mechanisms: the well-known break-induced replication and a novel mechanism of template switching between low-complexity or microhomologous sequences. These two mechanisms, rather than unequal recombination events, contribute in comparable proportions to duplication formation, the latter being prone to create novel gene fusions at chromosomal junctions. The mechanisms identified in yeast could explain the origin of a variety of genetic diseases in human, such as hemophilia A, Pelizaeus-Merzbacher disease, or some neurological disorders.

Construction, characterization, and chromosomal mapping of a fosmid library of the white-cheeked gibbon (Nomascus leucogenys)

Gibbons have experienced extensive karyotype rearrangements during evolution and represent an ideal model for studying the underlying molecular mechanism of evolutionary chromosomal rearrangements. It is anticipated that the cloning and sequence characterization of evolutionary chromosomal breakpoints will provide vital insights into the molecular force that has driven such a radical karyotype reshuffle in gibbons. We constructed and characterized a high-quality fosmid library of the white-cheeked gibbon (Nomascus leucogenys) containing 192,000 non-redundant clones with an average insert size of 38 kb and 2.5-fold genome coverage. By end sequencing of 100 randomly selected fosmid clones, we generated 196 sequence tags for the library. These end-sequenced fosmid clones were then mapped onto the chromosomes of the white-cheeked gibbon by fluorescence in situ hybridization, and no spurious chimeric clone was detected. BLAST search against the human genome showed a good correlation between the number of hit clones and the number of chromosomes, an indication of unbiased chromosomal distribution of the fosmid library. The chromosomal distribution of the mapped clones is also consistent with the BLAST search result against human and white-cheeked gibbon genomes. The fosmid library and the mapped clones will serve as a valuable resource for further studying gibbons’ chromosomal rearrangements and the underlying molecular mechanism as well as for comparative genomic study in the lesser apes.

Segmental duplications and evolutionary plasticity at tumor chromosome break-prone regions

We have previously found that the borders of evolutionarily conserved chromosomal regions often coincide with tumor-associated deletion breakpoints within human 3p12-p22. Moreover, a detailed analysis of a frequently deleted region at 3p21.3 (CER1) showed associations between tumor breaks and gene duplications. We now report on the analysis of 54 chromosome 3 breaks by multipoint FISH (mpFISH) in 10 carcinoma-derived cell lines. The centromeric region was broken in five lines. In lines with highly complex karyotypes, breaks were clustered near known fragile sites, FRA3B, FRA3C, and FRA3D (three lines), and in two other regions: 3p12.3-p13 ( approximately 75 Mb position) and 3q21.3-q22.1 ( approximately 130 Mb position) (six lines). All locations are shown based on NCBI Build 36.1 human genome sequence. The last two regions participated in three of four chromosome 3 inversions during primate evolution. Regions at 75, 127, and 131 Mb positions carry a large ( approximately 250 kb) segmental duplication (tumor break-prone segmental duplication [TBSD]). TBSD homologous sequences were found at 15 sites on different chromosomes. They were located within bands frequently involved in carcinoma-associated breaks. Thirteen of them have been involved in inversions during primate evolution; 10 were reused by breaks during mammalian evolution; 14 showed copy number polymorphism in man. TBSD sites showed an increase in satellite repeats, retrotransposed sequences, and other segmental duplications. We propose that the instability of these sites stems from specific organization of the chromosomal region, associated with location at a boundary between different CG-content isochores and with the presence of TBSDs and "instability elements," including satellite repeats and retroviral sequences.

Evolution of man in the light of molecular genetics: a review. Part I. Our evolutionary history and genomics

The discovery in the mid 1970s of efficient methods of DNA sequencing and their subsequent development into more and more rapid procedures followed by sequencing the genomes of many species, including man in 2001, revolutionised the whole of biology. Remarkably, new light could be cast on the evolutionary relations of different species, and the tempo and mode of evolution within a given species, notably man, could quantitatively be illuminated including ongoing evolution possibly involving also the size of the brains. This review is a short summary of the results of the molecular genetic investigations of human evolution including the time and place of the formation of our species, our evolutionary relation to the closest living species relatives as well as extinct forms of the genus Homo. The nature and amount of genetic polymorphism in man is also considered with special emphasis on the causes of this variation, and the role of natural selection in human evolution. A consensus about the mosaic nature of our genome and the rather dynamic structure of our ancestral population is gradually emerging. The modern gene pool has most likely been contributed to several different ancestral demes either before or after the emergence of the anatomically modern human phenotype in the extent that even the nature of the evolutionary lineage leading to the anatomically modern man as a distinct biological species is disputable. Regulation of the function of genes, as well as the evolution of brains will be dealt with in the second part of this review.

Hominoid lineage specific amplification of low-copy repeats on 22q11.2 (LCR22s) associated with velo-cardio-facial/digeorge syndrome

Segmental duplications or low-copy repeats (LCRs) constitute approximately 5% of the sequenced portion of the human genome and are associated with many human congenital anomaly disorders. The low-copy repeats on chromosome 22q11.2 (LCR22s) mediate chromosomal rearrangements resulting in deletions, duplications and translocations. The evolutionary mechanisms leading to LCR22 formation is unknown. Four genes, USP18, BCR, GGTLA and GGT, map adjacent to the LCR22s and pseudogene copies are located within them. It has been hypothesized that gene duplication occurred during primate evolution, followed by recombination events, forming pseudogene copies. We investigated whether gene duplication could be detected in non-human hominoid species. FISH mapping was performed using probes to the four functional gene loci. There was evidence for a single copy in humans but additional copies in hominoid species. We then compared LCR22 copy number using LCR22 FISH probes. Lineage specific LCR22 variation was detected in the hominoid species supporting the hypothesis. To independently validate initial findings, real time PCR, and screening of gorilla BAC library filters were performed. This was compared to array comparative genome hybridization data available. The most striking finding was a dramatic amplification of LCR22s in the gorilla. The LCR22s localized to the telomeric or subtelomeric bands of gorilla chromosomes. The most parsimonious explanation is that the LCR22s became amplified by inter-chromosomal recombination between telomeric bands. In summary, our results are consistent with a lineage specific coupling between gene and LCR22 duplication events. The LCR22s thus serve as an important model for evolution of genome variation.

A Drosophila melanogaster model of spinal muscular atrophy reveals a function for SMN in striated muscle

Mutations in human survival motor neurons 1 (SMN1) cause spinal muscular atrophy (SMA) and are associated with defects in assembly of small nuclear ribonucleoproteins (snRNPs) in vitro. However, the etiological link between snRNPs and SMA is unclear. We have developed a Drosophila melanogaster system to model SMA in vivo. Larval-lethal Smn-null mutations show no detectable snRNP reduction, making it unlikely that these animals die from global snRNP deprivation. Hypomorphic mutations in Smn reduce dSMN protein levels in the adult thorax, causing flightlessness and acute muscular atrophy. Mutant flight muscle motoneurons display pronounced axon routing and arborization defects. Moreover, Smn mutant myofibers fail to form thin filaments and phenocopy null mutations in Act88F, which is the flight muscle-specific actin isoform. In wild-type muscles, dSMN colocalizes with sarcomeric actin and forms a complex with alpha-actinin, the thin filament crosslinker. The sarcomeric localization of Smn is conserved in mouse myofibrils. These observations suggest a muscle-specific function for SMN and underline the importance of this tissue in modulating SMA severity.

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.

Modeling non-random deletions in cancer

Chromosome deletions do abound in cancer and are detected in certain regions in a non-random manner. Although their relevance remains elusive, it is a general agreement that segmental losses provide the cell with selective growth advantage. Consequently these may contain genes and/or regulatory sequences that control normal growth and inhibit malignancy. We have developed a monochromosomal hybrid based experimental model for the generation and functional analysis of deletions, that is called "elimination test" (Et). Focused on human chromosome 3 - that was known to carry multiple 3p deletions - the Et was expected to restrict a 3p tumor suppressor region to a sufficiently small segment that permits the selection of a critically important candidate gene. Surprisingly, we detected three regions that were lost in all or majority of tumors: CER1 (3p21.3, Mb: 43.32-45.74), CER2 (3p22, Mb: 37.83-39.06) and FER (3p14.3-p21.2, Mb: 50.12-58.03). In contrast a 3q26-qter region (CRR) was regularly retained. CER1 - our main focus - contains multiple genes that may inhibit tumor growth, but 3 genes, RIS1, LF (LTF) and LIMD1 have already the necessary experimental support to be considered bona fide tumor suppressors. Tumor suppressor region borders display instability features including: (1) they break in evolution and in tumors, (2) they evolve horizontally, and (3) they are enriched with pseudogene insertions. The most remarkable features at the breakpoint cluster regions were segmental duplications that drive horizontal evolution and contribute to cancer associated instability.

Complex patterns of copy number variation at sites of segmental duplications: an important category of structural variation in the human genome

The structural diversity of the human genome is much higher than previously assumed although its full extent remains unknown. To investigate the association between segmental duplications that display constitutive copy number differences (CNDs) between humans and the great apes and those which exhibit polymorphic copy number variations (CNVs) between humans, we analysed a BAC array enriched with segmental duplications displaying such CNDs. This study documents for the first time that in addition to human-specific gains common to all humans, these duplication clusters (DCs) also exhibit polymorphic CNVs > 40 kb. Segmental duplication is known to have been a frequent event during human genome evolution. Importantly, among the CNV-associated genes identified here, those involved in transcriptional regulation were found to be significantly overrepresented. Complex patterns of variation were evident at sites of DCs, manifesting as inter-individual differentially sized copy number alterations at the same genomic loci. Thus, CNVs associated with segmental duplications do not simply represent insertion/deletion polymorphisms, but rather constitute a wide variety of rearrangements involving differential amplification and partial gains and losses with high inter-individual variability. Although the number of CNVs was not found to differ between Africans and Caucasians/Asians, the average number of variant patterns per locus was significantly lower in Africans. Thus, complex variation patterns characterizing segmental duplications result from relatively recent genomic rearrangements. The high number of these rearrangements, some of which are potentially recurrent, together with differences in population size and expansion dynamics, may account for the greater diversity of CNV in Caucasians/Asians as compared with Africans.

Primate segmental duplications: crucibles of evolution, diversity and disease

Compared with other mammals, the genomes of humans and other primates show an enrichment of large, interspersed segmental duplications (SDs) with high levels of sequence identity. Recent evidence has begun to shed light on the origin of primate SDs, pointing to a complex interplay of mechanisms and indicating that distinct waves of duplication took place during primate evolution. There is also evidence for a strong association between duplication, genomic instability and large-scale chromosomal rearrangements. Exciting new findings suggest that SDs have not only created novel primate gene families, but might have also influenced current human genic and phenotypic variation on a previously unappreciated scale. A growing number of examples link natural human genetic variation of these regions to susceptibility to common disease.

A preliminary comparative analysis of primate segmental duplications shows elevated substitution rates and a great-ape expansion of intrachromosomal duplications

Compared with other sequenced animal genomes, human segmental duplications appear larger, more interspersed, and disproportionately represented as high-sequence identity alignments. Global sequence divergence estimates of human duplications have suggested an expansion relatively recently during hominoid evolution. Based on primate comparative sequence analysis of 37 unique duplication-transition regions, we establish a molecular clock for their divergence that shows a significant increase in their effective substitution rate when compared with unique genomic sequence. Fluorescent in situ hybridization (FISH) analyses from 1053 random nonhuman primate BACs indicate that great-ape species have been enriched for interspersed segmental duplications compared with representative Old World and New World monkeys. These findings support computational analyses that show a 12-fold excess of recent (>98%) intrachromosomal duplications when compared with duplications between nonhomologous chromosomes. These architectural shifts in genomic structure and elevated substitution rates have important implications for the emergence of new genes, gene-expression differences, and structural variation among humans and great apes.

Breakpoint analysis of the pericentric inversion distinguishing human chromosome 4 from the homologous chromosome in the chimpanzee (Pan troglodytes)

The study of breakpoints that occurred during primate evolution promises to yield valuable insights into the mechanisms underlying chromosome rearrangements in both evolution and pathology. Karyotypic differences between humans and chimpanzees include nine pericentric inversions, which may have potentiated the parapatric speciation of hominids and chimpanzees 5-6 million years ago. Detailed analysis of the respective chromosomal breakpoints is a prerequisite for any assessment of the genetic consequences of these inversions. The breakpoints of the inversion that distinguishes human chromosome 4 (HSA4) from its chimpanzee counterpart were identified by fluorescence in situ hybridization (FISH) and comparative sequence analysis. These breakpoints, at HSA4p14 and 4q21.3, do not disrupt the protein coding region of a gene, although they occur in regions with an abundance of LINE and LTR-elements. At 30 kb proximal to the breakpoint in 4q21.3, we identified an as yet unannotated gene, C4orf12, that lacks an homologous counterpart in rodents and is expressed at a 33-fold higher level in human fibroblasts as compared to chimpanzee. Seven out of 11 genes that mapped to the breakpoint regions have been previously analyzed using oligonucleotide-microarrays. One of these genes, WDFY3, exhibits a three-fold difference in expression between human and chimpanzee. To investigate whether the genomic architecture might have facilitated the inversion, comparative sequence analysis was used to identify an approximately 5-kb inverted repeat in the breakpoint regions. This inverted repeat is inexact and comprises six subrepeats with 78 to 98% complementarity. (TA)-rich repeats were also noted at the breakpoints. These findings imply that genomic architecture, and specifically high-copy repetitive elements, may have made a significant contribution to hominoid karyotype evolution, predisposing specific genomic regions to rearrangements.

A molecular cytogenetic study of chromosome evolution in chimpanzee

We applied multitude multicolor banding (mMCB) in combination with a novel FISH DNA probe set including subcentromeric, subtelomeric and whole chromosome painting probes (subCTM) to characterize a Pan paniscus (PPA) cell line. These powerful techniques allowed us to refine the breakpoints of a pericentric inversion on chimpanzee chromosome 4, and discovered a novel cryptic pericentric inversion in chimpanzee chromosome 11. mMCB provided a starting point for mapping and high resolution analysis of breakpoints on PPA chromosome 4, which are within a long terminal repeat (LTR) and surrounded by segmental duplications, as well as the integration/expansion sites of the interstitial heterochromatin on chimpanzee chromosomes 6 and 14. Moreover, we found evidence at hand for different types of heterochromatin in the chimpanzee genome. Finally, shedding new light on the human/chimpanzee speciation, karyotypes of three members of the genus Pan were studied by mMCB and no cytogenetic differences were found although the phylogenetic distance between these subspecies is suggested to be 2.5 million years.

Evidence for widespread reticulate evolution within human duplicons

Approximately 5% of the human genome consists of segmental duplications that can cause genomic mutations and may play a role in gene innovation. Reticulate evolutionary processes, such as unequal crossing-over and gene conversion, are known to occur within specific duplicon families, but the broader contribution of these processes to the evolution of human duplications remains poorly characterized. Here, we use phylogenetic profiling to analyze multiple alignments of 24 human duplicon families that span >8 Mb of DNA. Our results indicate that none of them are evolving independently, with all alignments showing sharp discontinuities in phylogenetic signal consistent with reticulation. To analyze these results in more detail, we have developed a quartet method that estimates the relative contribution of nucleotide substitution and reticulate processes to sequence evolution. Our data indicate that most of the duplications show a highly significant excess of sites consistent with reticulate evolution, compared with the number expected by nucleotide substitution alone, with 15 of 30 alignments showing a >20-fold excess over that expected. Using permutation tests, we also show that at least 5% of the total sequence shares 100% sequence identity because of reticulation, a figure that includes 74 independent tracts of perfect identity >2 kb in length. Furthermore, analysis of a subset of alignments indicates that the density of reticulation events is as high as 1 every 4 kb. These results indicate that phylogenetic relationships within recently duplicated human DNA can be rapidly disrupted by reticulate evolution. This finding has important implications for efforts to finish the human genome sequence, complicates comparative sequence analysis of duplicon families, and could profoundly influence the tempo of gene-family evolution.

Evolutionary implications of pericentromeric gene expression in humans

Human pericentromeric sequences are enriched for recent sequence duplications. The continual creation and shuffling of these duplications can create novel intron-exon structures and it has been suggested that these regions have a function as gene nurseries. However, these sequences are also rich in satellite repeats which can repress transcription, and analyses of chromosomes 10 and 21 have suggested that they are transcript poor. Here, we investigate the relationship between pericentromeric duplication and transcription by analyzing the in silico transcriptional profiles within the proximal 1.5 Mb of genomic sequence on all human chromosome arms in relation to duplication status. We identify an approximately 5x excess of transcripts specific to cancer and/or testis in pericentromeric duplications compared to surrounding single copy sequence, with the expression of >50% of all transcripts in duplications being restricted to these tissues. We also identify an approximately 5x excess of transcripts in duplications which contain large quantities of interspersed repeats. These results indicate that the transcriptional profiles of duplicated and single copy sequences within pericentromeric DNA are distinct, suggesting that pericentromeric instability is unlikely to represent a common route for gene creation but may have a disproportionate effect upon genes whose function is restricted to the germ line.

Molecular characterisation of the pericentric inversion that distinguishes human chromosome 5 from the homologous chimpanzee chromosome

Human and chimpanzee karyotypes differ by virtue of nine pericentric inversions that serve to distinguish human chromosomes 1, 4, 5, 9, 12, 15, 16, 17, and 18 from their chimpanzee orthologues. In this study, we have analysed the breakpoints of the pericentric inversion characteristic of chimpanzee chromosome 4, the homologue of human chromosome 5. Breakpoint-spanning BAC clones were identified from both the human and chimpanzee genomes by fluorescence in situ hybridisation, and the precise locations of the breakpoints were determined by sequence comparisons. In stark contrast to some other characterised evolutionary rearrangements in primates, this chimpanzee-specific inversion appears not to have been mediated by either gross segmental duplications or low-copy repeats, although micro-duplications were found adjacent to the breakpoints. However, alternating purine-pyrimidine (RY) tracts were detected at the breakpoints, and such sequences are known to adopt non-B DNA conformations that are capable of triggering DNA breakage and genomic rearrangements. Comparison of the breakpoint region of human chromosome 5q15 with the orthologous regions of the chicken, mouse, and rat genomes, revealed similar but non-identical syntenic disruptions in all three species. The clustering of evolutionary breakpoints within this chromosomal region, together with the presence of multiple pathological breakpoints in the vicinity of both 5p15 and 5q15, is consistent with the non-random model of chromosomal evolution and suggests that these regions may well possess intrinsic features that have served to mediate a variety of genomic rearrangements, including the pericentric inversion in chimpanzee chromosome 4.

Serial segmental duplications during primate evolution result in complex human genome architecture

The human genome is particularly rich in low-copy repeats (LCRs) or segmental duplications (5%-10%), and this characteristic likely distinguishes us from lower mammals such as rodents. How and why the complex human genome architecture consisting of multiple LCRs has evolved remains an open question. Using molecular and computational analyses of human and primate genomic regions, we analyzed the structure and evolution of LCRs that resulted in complex architectural features of the human genome in proximal 17p. We found that multiple LCRs of different origins are situated adjacent to one another, whereas each LCR changed at different time points between >25 to 3-7 million years ago (Mya) during primate evolution. Evolutionary studies in primates suggested communication between the LCRs by gene conversion. The DNA transposable element MER1-Charlie3 and retroviral ERVL elements were identified at the breakpoint of the t(4;19) chromosome translocation in Gorilla gorilla, suggesting a potential role for transpositions in evolution of the primate genome. Thus, a series of consecutive segmental duplication events during primate evolution resulted in complex genome architecture in proximal 17p. Some of the more recent events led to the formation of novel genes that in human are expressed primarily in the brain. Our observations support the contention that serial segmental duplication events might have orchestrated primate evolution by the generation of novel fusion/fission genes as well as potentially by genomic inversions associated with decreased recombination rates facilitating gene divergence.

The structure and evolution of centromeric transition regions within the human genome

An understanding of how centromeric transition regions are organized is a critical aspect of chromosome structure and function; however, the sequence context of these regions has been difficult to resolve on the basis of the draft genome sequence. We present a detailed analysis of the structure and assembly of all human pericentromeric regions (5 megabases). Most chromosome arms (35 out of 43) show a gradient of dwindling transcriptional diversity accompanied by an increasing number of interchromosomal duplications in proximity to the centromere. At least 30% of the centromeric transition region structure originates from euchromatic gene-containing segments of DNA that were duplicatively transposed towards pericentromeric regions at a rate of six-seven events per million years during primate evolution. This process has led to the formation of a minimum of 28 new transcripts by exon exaptation and exon shuffling, many of which are primarily expressed in the testis. The distribution of these duplicated segments is nonrandom among pericentromeric regions, suggesting that some regions have served as preferential acceptors of euchromatic DNA.

The domains of apoptosis: a genomics perspective

Apoptosis plays important roles in many facets of normal physiology in animal species, including programmed cell death associated with fetal development or metamorphosis, tissue homeostasis, immune cell education, and some aspects of aging. Defects in the regulation of apoptosis contribute to multiple diseases associated with either inappropriate cell loss or pathological cell accumulation. Host-pathogen interactions have additionally provided evolutionary pressure for apoptosis as a defense mechanism against viruses and microbes, sometimes linking apoptosis mechanisms with inflammatory responses. To a large extent, the apoptosis machinery can be viewed as a network, with different nodes connected by physical interactions of evolutionarily conserved domains. These domains can serve as signatures for identification of proteins involved in the network. In particular, the caspase recruitment domains (CARDs); death effector domains (DEDs); death domains (DDs); BIR (baculovirus IAP repeat) domains of inhibitor of apoptosis proteins (IAPs); Bcl-2 family proteins; caspase protease domains; and endonuclease-associated CIDE (cell death-inducing DFF45-like effector) domains are found in common in proteins involved in apoptosis. In the genomes of mammals, genes encoding proteins that carry one or more of these signature domains are often present in multiple copies, making up diverse gene families that permit tissue-specific and highly regulated control of cell life and death decisions through combinations of stimulus-specific gene expression and complex protein interaction networks. In this Review, we organize the repertoire of apoptosis proteins of humans into domain families, drawing comparisons with homologs in other vertebrate and invertebrate animal species, and discuss some of the functional implications of these findings.

The SPANX gene family of cancer/testis-specific antigens: rapid evolution and amplification in African great apes and hominids

Human sperm protein associated with the nucleus on the X chromosome (SPANX) genes comprise a gene family with five known members (SPANX-A1, -A2, -B, -C, and -D), encoding cancer/testis-specific antigens that are potential targets for cancer immunotherapy. These highly similar paralogous genes cluster on the X chromosome at Xq27. We isolated and sequenced primate genomic clones homologous to human SPANX. Analysis of these clones and search of the human genome sequence revealed an uncharacterized group of genes, SPANX-N, which are present in all primates as well as in mouse and rat. In humans, four SPANX-N genes comprise a series of tandem duplicates at Xq27; a fifth member of this subfamily is located at Xp11. Similarly to SPANX-A/D, human SPANX-N genes are expressed in normal testis and some melanoma cell lines; testis-specific expression of SPANX is also conserved in mouse. Analysis of the taxonomic distribution of the long and short forms of the intron indicates that SPANX-N is the ancestral form, from which the SPANX-A/D subfamily evolved in the common ancestor of the hominoid lineage. Strikingly, the coding sequences of the SPANX genes evolved much faster than the intron and the 5' untranslated region. There is a strong correlation between the rates of evolution of synonymous and nonsynonymous codon positions, both of which are accelerated 2-fold or more compared to the noncoding sequences. Thus, evolution of the SPANX family appears to have involved positive selection that affected not only the protein sequence but also the synonymous sites in the coding sequence.

Shuffling of genes within low-copy repeats on 22q11 (LCR22) by Alu-mediated recombination events during evolution

Low-copy repeats, or segmental duplications, are highly dynamic regions in the genome. The low-copy repeats on chromosome 22q11.2 (LCR22) are a complex mosaic of genes and pseudogenes formed by duplication processes; they mediate chromosome rearrangements associated with velo-cardio-facial syndrome/DiGeorge syndrome, der(22) syndrome, and cat-eye syndrome. The ability to trace the substrates and products of recombination events provides a unique opportunity to identify the mechanisms responsible for shaping LCR22s. We examined the genomic sequence of known LCR22 genes and their duplicated derivatives. We found Alu (SINE) elements at the breakpoints in the substrates and at the junctions in the truncated products of recombination for USP18, GGT, and GGTLA, consistent with Alu-mediated unequal crossing-over events. In addition, we were able to trace a likely interchromosomal Alu-mediated fusion between IGSF3 on 1p13.1 and GGT on 22q11.2. Breakpoints occurred inside Alu elements as well as in the 5' or 3' ends of them. A possible stimulus for the 5' or 3' terminal rearrangements may be the high sequence similarities between different Alu elements, combined with a potential recombinogenic role of retrotransposon target-site duplications flanking the Alu element, containing potentially kinkable DNA sites. Such sites may represent focal points for recombination. Thus, genome shuffling by Alu-mediated rearrangements has contributed to genome architecture during primate evolution.

Eucaryotic genome evolution through the spontaneous duplication of large chromosomal segments

There is growing evidence that duplications have played a major role in eucaryotic genome evolution. Sequencing data revealed the presence of large duplicated regions in the genomes of many eucaryotic organisms, and comparative studies have suggested that duplication of large DNA segments has been a continuing process during evolution. However, little experimental data have been produced regarding this issue. Using a gene dosage assay for growth recovery in Saccharomyces cerevisiae, we demonstrate that a majority of the revertant strains (58%) resulted from the spontaneous duplication of large DNA segments, either intra- or interchromosomally, ranging from 41 to 655 kb in size. These events result in the concomitant duplication of dozens of genes and in some cases in the formation of chimeric open reading frames at the junction of the duplicated blocks. The types of sequences at the breakpoints as well as their superposition with the replication map suggest that spontaneous large segmental duplications result from replication accidents. Aneuploidization events or suppressor mutations that do not involve large-scale rearrangements accounted for the rest of the reversion events (in 26 and 16% of the strains, respectively).

Eukaryotic domain evolution inferred from genome comparisons

Comparative analyses of eukaryotic genomes are providing insights into the mode and tempo of domain family evolution. Gene duplication, the source of family expansion, far exceeds the rate of emergence of domains from non-coding sequence, and the rate of recruitment of domains into novel architectures. Domain families that appear to be restricted to certain lineages are likely to be the result of gene duplication, coupled with rapid sequence diversification. If such families are evidence of past adaptation, then their functions must relate to the underlying mechanism of selection: competition among organisms.

Duplicate, decouple, disperse: the evolutionary transience of human centromeric regions

Human centromeric regions are enriched for segmental duplications, which elsewhere in the genome precipitate both genetic disease and gene formation. Molecular cytogenetic analyses of primate chromosomes have established that centromeres frequently move without altering the surrounding gene order. Recently, the positions of two ancestral centromeres have been mapped to regions of the human genome that are both rich in segmental duplications and are associated with duplication-based clinical phenotypes. This suggests a model for the evolution of euchromatic segmental duplication families involving the localised elevation of recombination rates within the duplication-rich heterochromatin of recently inactivated centromeres, and raises the possibility that the distribution of duplication/deletion syndromes within our genome has been heavily influenced by such events. The relaxation of the heterochromatin environment that must accompany centromere inactivation would also increase the transcriptional activity within previously pericentromeric DNA, increasing the likelihood of chimaeric gene creation through pericentromeric-directed duplication events.

The DNA sequence and analysis of human chromosome 6

Chromosome 6 is a metacentric chromosome that constitutes about 6% of the human genome. The finished sequence comprises 166,880,988 base pairs, representing the largest chromosome sequenced so far. The entire sequence has been subjected to high-quality manual annotation, resulting in the evidence-supported identification of 1,557 genes and 633 pseudogenes. Here we report that at least 96% of the protein-coding genes have been identified, as assessed by multi-species comparative sequence analysis, and provide evidence for the presence of further, otherwise unsupported exons/genes. Among these are genes directly implicated in cancer, schizophrenia, autoimmunity and many other diseases. Chromosome 6 harbours the largest transfer RNA gene cluster in the genome; we show that this cluster co-localizes with a region of high transcriptional activity. Within the essential immune loci of the major histocompatibility complex, we find HLA-B to be the most polymorphic gene on chromosome 6 and in the human genome.

Evolution of primate theta-defensins: a serpentine path to a sweet tooth

Retrocyclins (ancestral human theta-defensins) are cyclic antimicrobial octadecapeptides that interfere with viral uptake and protect human cells from infection by T- and M-tropic strains of HIV-1 in vitro. As are other theta-defensins, retrocyclins are lectins that bind gp120, CD4, and galactosylceramide-all of which are implicated in HIV-1 uptake. Although theta-defensin mRNA transcripts are present in human bone marrow, spleen, thymus, testis, and skeletal muscle, a premature stop codon aborts their translation. We found six theta-defensin (DEFT) genes in the human genome; five on chromosome 8p23 and one on chromosome 1. All six of these pseudogenes, as well as their homologues in chimpanzees and gorillas, contained the same premature stop codon mutation. Whereas we found intact DEFT genes in DNA from several Old World Monkeys, Hylobates syndactylus (a lesser ape) and orangutans, no homologues were present in DNA from six New World Monkeys and five prosimians. We conclude that DEFT genes and theta-defensins arose in Old World Monkeys by mutation of a pre-existing alpha-defensin gene. Although intact DEFT genes survive in some nonhuman primates, our hominid ancestors lost their ability to produce theta-defensins after the orangutan and hominid lineages diverged. It is possible (but may be difficult to prove) that this mutation rendered our species more susceptible to infection by HIV-1.