Adaptation shapes patterns of genome evolution on sexual and asexual chromosomes in Drosophila

What advantage might sexual recombination confer? Population genetics theory predicts that asexual genomes are less efficient at eliminating deleterious mutations and incorporating beneficial alleles. Here, I compare patterns of genome evolution in a 40-kb gene-rich region on homologous neo-sex chromosomes of Drosophila miranda. Genes on the non-recombining neo-Y show various signs of degeneration, including transposable-element insertions, frameshift mutations and a higher rate of amino-acid substitution. In contrast, loci on the recombining neo-X show intact open reading frames and generally low rates of amino-acid substitution. One exceptional gene on the neo-X shows evidence for adaptive protein evolution, affecting patterns of variability at neighboring regions along the chromosome. These findings illustrate the limits to natural selection in an asexual genome. Deleterious mutations, including repetitive DNA, accumulate on a non-recombining chromosome, whereas rapid protein evolution due to positive selection is confined to the recombining homolog.     

Kallmann syndrome gene on the X and Y chromosomes: implications for evolutionary divergence of human sex chromosomes.

The recently identified gene for X-linked Kallmann syndrome (hypogonadotropic hypogonadism and anosmia) has a closely related homologue on the Y chromosome. The X and Y copies of this gene are located in a large region of X/Y homology, on Xp22.3 and Yq11.2, respectively. Comparison of the structure of the X-linked Kallmann syndrome gene and its Y homologue shed light on the evolutionary history of this region of the human sex chromosomes. Our data show that the Y homologue is not functional. Comparative analysis of X/Y sequence identity at several loci on Xp22.3 and Yq11.2 suggests that the homology between these two regions is the result of a complex series of events which occurred in the recent evolution of sex chromosomes.     

Dichotomy of single-nucleotide polymorphism haplotypes in olfactory receptor genes and pseudogenes

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies. Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome. We report here the population sequence diversity of genomic segments within a 450-kb cluster of olfactory receptor (OR) genes on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes, and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.     

Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway

The mammalian Y chromosome has unique characteristics compared with the autosomes or X chromosomes. Here we report the finished sequence of the chimpanzee Y chromosome (PTRY), including 271 kb of the Y-specific pseudoautosomal region 1 and 12.7 Mb of the male-specific region of the Y chromosome. Greater sequence divergence between the human Y chromosome (HSAY) and PTRY (1.78%) than between their respective whole genomes (1.23%) confirmed the accelerated evolutionary rate of the Y chromosome. Each of the 19 PTRY protein-coding genes analyzed had at least one nonsynonymous substitution, and 11 genes had higher nonsynonymous substitution rates than synonymous ones, suggesting relaxation of selective constraint, positive selection or both. We also identified lineage-specific changes, including deletion of a 200-kb fragment from the pericentromeric region of HSAY, expansion of young Alu families in HSAY and accumulation of young L1 elements and long terminal repeat retrotransposons in PTRY. Reconstruction of the common ancestral Y chromosome reflects the dynamic changes in our genomes in the 5-6 million years since speciation.     

Mutations in SPTLC1, encoding serine palmitoyltransferase, long chain base subunit-1, cause hereditary sensory neuropathy type I

Hereditary sensory neuropathy type I (HSN1) is the most common hereditary disorder of peripheral sensory neurons. HSN1 is an autosomal dominant progressive degeneration of dorsal root ganglia and motor neurons with onset in the second or third decades. Initial symptoms are sensory loss in the feet followed by distal muscle wasting and weakness. Loss of pain sensation leads to chronic skin ulcers and distal amputations. The HSN1 locus has been mapped to chromosome 9q22.1-22.3 (refs. 3,4). Here we map the gene SPTLC1, encoding serine palmitoyltransferase, long chain base subunit-1, to this locus. Mutation screening revealed 3 different missense mutations resulting in changes to 2 amino acids in all affected members of 11 HSN1 families. We found two mutations to be located in exon 5 (C133Y and C133W) and one mutation to be located in exon 6 of SPTLC1 (V144D). All families showing definite or probable linkage to chromosome 9 had mutations in these two exons. These mutations are associated with increased de novo glucosyl ceramide synthesis in lymphoblast cell lines in affected individuals. Increased de novo ceramide synthesis triggers apoptosis and is associated with massive cell death during neural tube closure, raising the possibility that neural degeneration in HSN1 is due to ceramide-induced apoptotic cell death.     

High mutation rates have driven extensive structural polymorphism among human Y chromosomes

Although much structural polymorphism in the human genome has been catalogued, the kinetics of underlying change remain largely unexplored. Because human Y chromosomes are clonally inherited, it has been possible to capture their detailed relationships in a robust, worldwide genealogical tree. Examination of structural variation across this tree opens avenues for investigating rates of underlying mutations. We selected one Y chromosome from each of 47 branches of this tree and searched for large-scale variation. Four chromosomal regions showed extensive variation resulting from numerous large-scale mutations. Within the tree encompassed by the studied chromosomes, the distal-Yq heterochromatin changed length > or = 12 times, the TSPY gene array changed length > or = 23 times, the 3.6-Mb IR3/IR3 region changed orientation > or = 12 times and the AZFc region was rearranged > or = 20 times. After determining the total time spanned by all branches of this tree (approximately 1.3 million years or 52,000 generations), we converted these mutation counts to lower bounds on rates: > or = 2.3 x 10(-4), > or = 4.4 x 10(-4), > or = 2.3 x 10(-4) and > or = 3.8 x 10(-4) large-scale mutations per father-to-son Y transmission, respectively. Thus, high mutation rates have driven extensive structural polymorphism among human Y chromosomes. At the same time, we found limited variation in the copy number of Y-linked genes, which raises the possibility of selective constraints.     

A common CFH haplotype, with deletion of CFHR1 and CFHR3, is associated with lower risk of age-related macular degeneration

Age-related macular degeneration (AMD; OMIM #603075) is the most frequent cause of visual impairment in the elderly population, with severe disease affecting nearly 10% of individuals of European descent over the age of 75 years. It is a complex disease in which genetic and environmental factors contribute to susceptibility. Complement factor H (CFH) has recently been identified as a major AMD susceptibility gene, and the Y402H polymorphism has been proposed as the likely causative factor. We genotyped polymorphisms spanning the cluster of CFH and five CFH-related genes on chromosome 1q23 in 173 individuals with severe neovascular AMD and 170 elderly controls with no signs of AMD. Detailed analysis showed a common haplotype associated with decreased risk of AMD that was present on 20% of chromosomes of controls and 8% of chromosomes of individuals with AMD. We found that this haplotype carried a deletion of CFHR1 and CFHR3, and the proteins encoded by these genes were absent in serum of homozygotes. The protective effect of the deletion haplotype cannot be attributed to linkage disequilibrium with Y402H and was replicated in an independent sample.     

Mutation of CDH23, encoding a new member of the cadherin gene family, causes Usher syndrome type 1D

Usher syndrome type I (USH1) is an autosomal recessive disorder characterized by congenital sensorineural hearing loss, vestibular dysfunction and visual impairment due to early onset retinitis pigmentosa (RP). So far, six loci (USH1A-USH1F) have been mapped, but only two USH1 genes have been identified: MYO7A for USH1B and the gene encoding harmonin for USH1C. We identified a Cuban pedigree linked to the locus for Usher syndrome type 1D (MIM 601067) within the q2 region of chromosome 10). Affected individuals present with congenital deafness and a highly variable degree of retinal degeneration. Using a positional candidate approach, we identified a new member of the cadherin gene superfamily, CDH23. It encodes a protein of 3,354 amino acids with a single transmembrane domain and 27 cadherin repeats. In the Cuban family, we detected two different mutations: a severe course of the retinal disease was observed in individuals homozygous for what is probably a truncating splice-site mutation (c.4488G-->C), whereas mild RP is present in individuals carrying the homozygous missense mutation R1746Q. A variable expression of the retinal phenotype was seen in patients with a combination of both mutations. In addition, we identified two mutations, Delta M1281 and IVS51+5G-->A, in a German USH1 patient. Our data show that different mutations in CDH23 result in USH1D with a variable retinal phenotype. In an accompanying paper, it is shown that mutations in the mouse ortholog cause disorganization of inner ear stereocilia and deafness in the waltzer mouse.     

Homozygous Tsix mutant mice reveal a sex-ratio distortion and revert to random X-inactivation

Tsix controls X-chromosome inactivation (XCI) by blocking the accumulation of Xist RNA on the future active X chromosome. Deleting Tsix on one X chromosome (X(Delta)X) skews XCI toward the mutated X chromosome in the female soma. Here I have generated homozygous Tsix-null mice (X(Delta)X(Delta)) to test how deleting the second allele affects the choice of XCI. Homozygosity leads to extremely low fertility and reveals two previously unknown non-mendelian patterns of inheritance. First, the sex ratio is skewed against female births so that one daughter is born for every two to three sons. Second, the pattern of XCI unexpectedly returns to random in surviving X(Delta)X(Delta) mice. Thus, with respect to choice, mutation of Tsix yields a phenotypic abnormality in heterozygotes but not homozygotes. To reconcile the paradox of female loss with apparent reversion to random choice, I propose that deleting both Tsix alleles results in chaotic choice and that randomness in X(Delta)X(Delta) survivors reflects a fortuitous selection of distinct X chromosomes as active and inactive.     

Global diversity and evidence for coevolution of KIR and HLA

The killer immunoglobulin-like receptor (KIR) gene cluster shows extensive genetic diversity, as do the HLA class I loci, which encode ligands for KIR molecules. We genotyped 1,642 individuals from 30 geographically distinct populations to examine population-level evidence for coevolution of these two functionally related but unlinked gene clusters. We observed strong negative correlations between the presence of activating KIR genes and their corresponding HLA ligand groups across populations, especially KIR3DS1 and its putative HLA-B Bw4-80I ligands (r = -0.66, P = 0.038). In contrast, we observed weak positive relationships between the various inhibitory KIR genes and their ligands. We observed a negative correlation between distance from East Africa and frequency of activating KIR genes and their corresponding ligands, suggesting a balance between selection on HLA and KIR loci. Most KIR-HLA genetic association studies indicate a primary influence of activating KIR-HLA genotypes in disease risk; concomitantly, activating receptor-ligand pairs in this study show the strongest signature of coevolution of these two complex genetic systems as compared with inhibitory receptor-ligand pairs.     

Abnormal cerebellar development and axonal decussation due to mutations in AHI1 in Joubert syndrome

Joubert syndrome is a congenital brain malformation of the cerebellar vermis and brainstem with abnormalities of axonal decussation (crossing in the brain) affecting the corticospinal tract and superior cerebellar peduncles. Individuals with Joubert syndrome have motor and behavioral abnormalities, including an inability to walk due to severe clumsiness and 'mirror' movements, and cognitive and behavioral disturbances. Here we identified a locus associated with Joubert syndrome, JBTS3, on chromosome 6q23.2-q23.3 and found three deleterious mutations in AHI1, the first gene to be associated with Joubert syndrome. AHI1 is most highly expressed in brain, particularly in neurons that give rise to the crossing axons of the corticospinal tract and superior cerebellar peduncles. Comparative genetic analysis of AHI1 indicates that it has undergone positive evolutionary selection along the human lineage. Therefore, changes in AHI1 may have been important in the evolution of human-specific motor behaviors.     

Natural selection has driven population differentiation in modern humans

The considerable range of observed phenotypic variation in human populations may reflect, in part, distinctive processes of natural selection and adaptation to variable environmental conditions. Although recent genome-wide studies have identified candidate regions under selection, it is not yet clear how natural selection has shaped population differentiation. Here, we have analyzed the degree of population differentiation at 2.8 million Phase II HapMap single-nucleotide polymorphisms. We find that negative selection has globally reduced population differentiation at amino acid-altering mutations, particularly in disease-related genes. Conversely, positive selection has ensured the regional adaptation of human populations by increasing population differentiation in gene regions, primarily at nonsynonymous and 5'-UTR variants. Our analyses identify a fraction of loci that have contributed, and probably still contribute, to the morphological and disease-related phenotypic diversity of current human populations.     

Fine mapping of regulatory loci for mammalian gene expression using radiation hybrids

We mapped regulatory loci for nearly all protein-coding genes in mammals using comparative genomic hybridization and expression array measurements from a panel of mouse-hamster radiation hybrid cell lines. The large number of breaks in the mouse chromosomes and the dense genotyping of the panel allowed extremely sharp mapping of loci. As the regulatory loci result from extra gene dosage, we call them copy number expression quantitative trait loci, or ceQTLs. The -2log10P support interval for the ceQTLs was <150 kb, containing an average of <2-3 genes. We identified 29,769 trans ceQTLs with -log10P > 4, including 13 hotspots each regulating >100 genes in trans. Further, this work identifies 2,761 trans ceQTLs harboring no known genes, and provides evidence for a mode of gene expression autoregulation specific to the X chromosome.     

Nested chromosomal deletions induced with retroviral vectors in mice

Chromosomal deletions, especially nested deletions, are major genetic tools in diploid organisms that facilitate the functional analysis of large chromosomal regions and allow the rapid localization of mutations to specific genetic intervals. In mice, well-characterized overlapping deletions are only available at a few chromosomal loci, partly due to drawbacks of existing methods. Here we exploit the random integration of a retrovirus to generate high-resolution sets of nested deletions around defined loci in embryonic stem (ES) cells, with sizes extending from a few kilobases to several megabases. This approach expands the application of Cre-loxP-based chromosome engineering because it not only allows the construction of hundreds of overlapping deletions, but also provides molecular entry points to regions based on the retroviral tags. Our approach can be extended to any region of the mouse genome.     

Telomere-associated chromosome fragmentation: applications in genome manipulation and analysis.

Telomere-associated chromosome fragmentation (TACF) is a new approach for chromosome mapping based on the non-targeted introduction of cloned telomeres into mammalian cells. TACF has been used to generate a panel of somatic cell hybrids with nested terminal deletions of the long arm of the human X chromosome, extending from Xq26 to the centromere. This panel has been characterized using a series of X chromosome loci. Recovery of the end clones by plasmid rescue produces a telomeric marker for each cell line and partial sequencing will allow the generation of sequence tagged sites (STSs). TACF provides a powerful and widely applicable method for genome analysis, a general way of manipulating mammalian chromosomes and a first step towards constructing artificial mammalian chromosomes.     

A comprehensive linkage analysis for myocardial infarction and its related risk factors

Coronary artery disease and myocardial infarction (MI) are leading causes of death in the western world. Numerous studies have shown that risk factors such as diabetes mellitus, arterial hypertension and hypercholesterolemia contribute to the development of the disease. Although each risk factor by itself is partly under genetic control, a positive family history is an independent predictor, which suggests that there are additional susceptibility genes. We have scanned the whole genome in 513 families to identify chromosomal regions linked to myocardial infarction and related risk factors that are known to be under genetic control. Here we show, by using variance component analysis and incorporating risk factors, that risk of myocardial infarction maps to a single region on chromosome 14 with a significant lod score of 3.9 (pointwise P=0.00015, genome-wide P<0.05), providing evidence of a principal MI locus. To characterize this locus we analyzed each risk factor by itself. Serum concentrations of lipoprotein (a) show linkage to both the apolipoprotein (a) locus (lod score 26.99) and a new locus on chromosome 1 (lod score 3.8). There is suggestive linkage for diabetes mellitus on chromosome 6 (lod score 2.96), for hypertension on chromosomes 1 and 6, for high-density and low-density lipoprotein cholesterol on chromosomes 1 and 17, and for triglyceride concentrations on chromosome 9. Although some of these risk factors overlap with previously identified loci, none overlaps with the newly identified susceptibility locus for myocardial infarction and coronary artery disease.     

Genetic inroads in familial ALS

Amyotrophic lateral sclerosis (ALS) is a common neurodegenerative disease causing cell death of motor neurons and progressive muscle weakness. The disease is familial in ten percent of cases, of which one-fifth are due to mutations in the gene encoding Cu/Zn superoxide dismutase (SOD1). Two papers in this issue of Nature Genetics describe homozygous mutations in a new gene on chromosome 2q33 in 4 families of Arabian origin with a rare form of juvenile onset ALS (ALS2). The predicted protein structure has domains homologous to GTPase regulatory proteins, and both the types of mutation and the pattern of inheritance suggest that motor neuron degeneration is the result of a loss of function. Further work will determine the relevance of this breakthrough to other, more common forms of ALS.     

Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility

A genome-wide association scan in individuals with Crohn's disease by the Wellcome Trust Case Control Consortium detected strong association at four novel loci. We tested 37 SNPs from these and other loci for association in an independent case-control sample. We obtained replication for the autophagy-inducing IRGM gene on chromosome 5q33.1 (replication P = 6.6 x 10(-4), combined P = 2.1 x 10(-10)) and for nine other loci, including NKX2-3, PTPN2 and gene deserts on chromosomes 1q and 5p13.     

Quantitative trait loci mapped to single-nucleotide resolution in yeast

Identifying the genetic variation underlying quantitative trait loci remains problematic. Consequently, our molecular understanding of genetically complex, quantitative traits is limited. To address this issue directly, we mapped three quantitative trait loci that control yeast sporulation efficiency to single-nucleotide resolution in a noncoding regulatory region (RME1) and to two missense mutations (TAO3 and MKT1). For each quantitative trait locus, the responsible polymorphism is rare among a diverse set of 13 yeast strains, suggestive of genetic heterogeneity in the control of yeast sporulation. Additionally, under optimal conditions, we reconstituted approximately 92% of the sporulation efficiency difference between the two genetically distinct parents by engineering three nucleotide changes in the appropriate yeast genome. Our results provide the highest resolution to date of the molecular basis of a quantitative trait, showing that the interaction of a few genetic variants can have a profound phenotypic effect.     

Meiotic pairing and imprinted X chromatin assembly in Caenorhabditis elegans

The genetic imprinting of individual loci or whole chromosomes, as in imprinted X-chromosome inactivation in mammals, is established and reset during gametogenesis; defects in this process in the parent can result in disease in the offspring. We describe a sperm-specific chromatin-based imprinting of the X chromosome in the nematode Caenorhabditis elegans that is restricted to histone H3 modifications. The epigenetic imprint is established during spermatogenesis and its stability in the offspring is affected by the presence of a pairing partner during meiosis in the parental germ line. We observed that DNA lacking a pairing partner during meiosis, the normal situation for the X chromosome in males, is targeted for methylation of histone H3 at Lys9 (H3-Lys9) and can be silenced. Targeting unpaired DNA for silencing during meiosis, a potential hallmark of genome defense, could therefore have a conserved role in imprinted X-chromosome inactivation and, ultimately, in sex chromosome evolution.