Regional loss of imprinting and growth deficiency in mice with a targeted deletion of KvDMR1

Genomic imprinting is an epigenetic modification that results in expression from only one of the two parental copies of a gene. Differences in methylation between the two parental chromosomes are often observed at or near imprinted genes. Beckwith-Wiedemann syndrome (BWS), which predisposes to cancer and excessive growth, results from a disruption of imprinted gene expression in chromosome band 11p15.5. One third of individuals with BWS lose maternal-specific methylation at KvDMR1, a putative imprinting control region within intron 10 of the KCNQ1 gene, and it has been proposed that this epimutation results in aberrant imprinting and, consequently, BWS1, 2. Here we show that paternal inheritance of a deletion of KvDMR1 results in the de-repression in cis of six genes, including Cdkn1c, which encodes cyclin-dependent kinase inhibitor 1C. Furthermore, fetuses and adult mice that inherited the deletion from their fathers were 20-25% smaller than their wildtype littermates. By contrast, maternal inheritance of this deletion had no effect on imprinted gene expression or growth. Thus, the unmethylated paternal KvDMR1 allele regulates imprinted expression by silencing genes on the paternal chromosome. These findings support the hypothesis that loss of methylation in BWS patients activates the repressive function of KvDMR1 on the maternal chromosome, resulting in abnormal silencing of CDKN1C and the development of BWS.     

Whole-genome sequencing and comprehensive variant analysis of a Japanese individual using massively parallel sequencing

We report the analysis of a Japanese male using high-throughput sequencing to × 40 coverage. More than 99% of the sequence reads were mapped to the reference human genome. Using a Bayesian decision method, we identified 3,132,608 single nucleotide variations (SNVs). Comparison with six previously reported genomes revealed an excess of singleton nonsense and nonsynonymous SNVs, as well as singleton SNVs in conserved non-coding regions. We also identified 5,319 deletions smaller than 10 kb with high accuracy, in addition to copy number variations and rearrangements. De novo assembly of the unmapped sequence reads generated around 3 Mb of novel sequence, which showed high similarity to non-reference human genomes and the human herpesvirus 4 genome. Our analysis suggests that considerable variation remains undiscovered in the human genome and that whole-genome sequencing is an invaluable tool for obtaining a complete understanding of human genetic variation.     

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.     

Closing gaps in the human genome with fosmid resources generated from multiple individuals

The human genome sequence has been finished to very high standards; however, more than 340 gaps remained when the finished genome was published by the International Human Genome Sequencing Consortium in 2004. Using fosmid resources generated from multiple individuals, we targeted gaps in the euchromatic part of the human genome. Here we report 2,488,842 bp of previously unknown euchromatic sequence, 363,114 bp of which close 26 of 250 euchromatic gaps, or 10%, including two remaining euchromatic gaps on chromosome 19. Eight (30.7%) of the closed gaps were found to be polymorphic. These sequences allow complete annotation of several human genes as well as the assignment of mRNAs. The gap sequences are 2.3-fold enriched in segmentally duplicated sequences compared to the whole genome. Our analysis confirms that not all gaps within 'finished' genomes are recalcitrant to subcloning and suggests that the paired-end-sequenced fosmid libraries could prove to be a rich resource for completion of the human euchromatic genome.     

Evidence that positive selection drives Y-chromosome degeneration in Drosophila miranda

Why does the Y chromosome harbor so few functional loci? Evolutionary theory predicts that Y chromosomes degenerate because they lack genetic recombination. Both positive and negative selection models have been invoked to explain this degeneration, as both can result in the recurrent fixation of linked deleterious mutations on a nonrecombining Y chromosome. To distinguish between these models, I investigated patterns of nucleotide variability along 37 kb of the recently formed neo-Y chromosome in Drosophila miranda. Levels of nucleotide variability on this chromosome are 30 times lower than in highly recombining portions of the genome. Both positive and negative selection models can result in reduced variability levels, but their effects on the frequency spectrum of mutations differ. Using coalescent simulations, I show that the patterns of nucleotide variability on the neo-Y chromosome are unlikely under deleterious mutation models (including background selection and Muller's ratchet) but are expected under recent positive selection. These results implicate positive selection as an important force driving the degeneration of Y chromosomes; adaptation at a few loci, possibly increasing male fitness, occurs at the cost of most other genes on this chromosome.     

Candidate lung tumor susceptibility genes identified through whole-genome association analyses in inbred mice

We performed a whole-genome association analysis of lung tumor susceptibility using dense SNP maps ( approximately 1 SNP per 20 kb) in inbred mice. We reproduced the pulmonary adenoma susceptibility 1 (Pas1) locus identified in previous linkage studies and further narrowed this quantitative trait locus (QTL) to a region of less than 0.5 Mb in which at least two genes, Kras2 (Kirsten rat sarcoma oncogene 2) and Casc1 (cancer susceptibility candidate 1; also known as Las1), are strong candidates. Casc1 knockout mouse tumor bioassays showed that Casc1-deficient mice were susceptible to chemical induction of lung tumors. We also found three more genetic loci for lung adenoma development. Analysis of one of these candidate loci identified a previously uncharacterized gene Lasc1, bearing a nonsynonymous substitution (D102E). We found that the Lasc1 Glu102 allele preferentially promotes lung tumor cell growth. Our findings demonstrate the prospects for using dense SNP maps in laboratory mice to refine previous QTL regions and identify genetic determinants of complex traits.     

An azoospermic man with a de novo point mutation in the Y-chromosomal gene USP9Y

In humans, deletion of any one of three Y-chromosomal regions- AZFa, AZFb or AZFc-disrupts spermatogenesis, causing infertility in otherwise healthy men. Although candidate genes have been identified in all three regions, no case of spermatogenic failure has been traced to a point mutation in a Y-linked gene, or to a deletion of a single Y-linked gene. We sequenced the AZFa region of the Y chromosome and identified two functional genes previously described: USP9Y (also known as DFFRY) and DBY (refs 7,8). Screening of the two genes in 576 infertile and 96 fertile men revealed several sequence variants, most of which appear to be heritable and of little functional consequence. We found one de novo mutation in USP9Y: a 4-bp deletion in a splice-donor site, causing an exon to be skipped and protein truncation. This mutation was present in a man with nonobstructive azoospermia (that is, no sperm was detected in semen), but absent in his fertile brother, suggesting that the USP9Y mutation caused spermatogenic failure. We also identified a single-gene deletion associated with spermatogenic failure, again involving USP9Y, by re-analysing a published study.     

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.     

Evidence in the Legionella pneumophila genome for exploitation of host cell functions and high genome plasticity

Legionella pneumophila, the causative agent of Legionnaires' disease, replicates as an intracellular parasite of amoebae and persists in the environment as a free-living microbe. Here we have analyzed the complete genome sequences of L. pneumophila Paris (3,503,610 bp, 3,077 genes), an endemic strain that is predominant in France, and Lens (3,345,687 bp, 2,932 genes), an epidemic strain responsible for a major outbreak of disease in France. The L. pneumophila genomes show marked plasticity, with three different plasmids and with about 13% of the sequence differing between the two strains. Only strain Paris contains a type V secretion system, and its Lvh type IV secretion system is encoded by a 36-kb region that is either carried on a multicopy plasmid or integrated into the chromosome. Genetic mobility may enhance the versatility of L. pneumophila. Numerous genes encode eukaryotic-like proteins or motifs that are predicted to modulate host cell functions to the pathogen's advantage. The genome thus reflects the history and lifestyle of L. pneumophila, a human pathogen of macrophages that coevolved with fresh-water amoebae.     

Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes

Bacterial pathogens evolve during the infection of their human host(1-8), but separating adaptive and neutral mutations remains challenging(9-11). Here we identify bacterial genes under adaptive evolution by tracking recurrent patterns of mutations in the same pathogenic strain during the infection of multiple individuals. We conducted a retrospective study of a Burkholderia dolosa outbreak among subjects with cystic fibrosis, sequencing the genomes of 112 isolates collected from 14 individuals over 16 years. We find that 17 bacterial genes acquired nonsynonymous mutations in multiple individuals, which indicates parallel adaptive evolution. Mutations in these genes affect important pathogenic phenotypes, including antibiotic resistance and bacterial membrane composition and implicate oxygen-dependent regulation as paramount in lung infections. Several genes have not previously been implicated in pathogenesis and may represent new therapeutic targets. The identification of parallel molecular evolution as a pathogen spreads among multiple individuals points to the key selection forces it experiences within human hosts.     

Linkage of an important gene locus for tuberous sclerosis to a chromosome 16 marker for polycystic kidney disease.

Tuberous sclerosis complex (TSC) is an autosomal dominant disorder of unknown aetiology that affects numerous body systems including skin, brain and kidneys. Some TSC has been linked to chromosome 9, additional TSC genes on chromosomes 11 and 12 have been proposed, but the majority of TSC families remain unlinked. Using TSC families in which data had excluded linkage to chromosome 9, we failed to detect linkage with loci on chromosomes 11, 12 and others. One marker examined was D16S283, the closest locus on the proximal side of the polycystic kidney disease type 1 (PKD1) gene. Linkage between TSC and D16S283 demonstrated a lod score of 9.50 at theta = 0.02 with one family independently presenting a lod score of 4.44 at theta = 0.05. These data reveal an important TSC locus near the region of PKD1 on chromosome 16p13.     

Evidence of en bloc duplication in vertebrate genomes

It has been 30 years since it was first proposed that the vertebrate genome evolved through several rounds of genome-wide duplications (polyploidizations). Despite rapid advances in genetics, including sequencing of the complete genomes of several divergent species, this hypothesis has not been tested rigorously and is still a matter of debate. If polyploidizations occurred during chordate evolution, there should be a network of paralogous regions in the present-day jawed vertebrate (Gnathostomata) genomes. Here we present an investigation of the major histocompatibility complex (MHC) paralogous regions, which we accomplished by characterizing the corresponding region in amphioxus by identifying nine anchor genes and sequencing both the anchor genes and the regions that flank them (a total of 400 kb). Phylogenetic analysis of 31 genes (including the anchor genes) in these regions shows that duplications occurred after the divergence of cephalochordates and vertebrates but before the Gnathostomata radiation. The distribution of human and amphioxus orthologs in their respective genomes and the relationship between these distributions support the en bloc duplication events. Our analysis represents the first step towards demonstrating that the human ancestral genome has undergone polyploidization. Moreover, reconstruction of the pre-duplicated region indicates that one of the duplicated regions retains the ancestral organization.     

Susceptibility to testicular germ-cell tumours in a 129.MOLF-Chr 19 chromosome substitution strain.

The identification of genes that control susceptibility to testicular germ-cell tumours (TGCTs), the most common cancer affecting young men, has been difficult. In laboratory mice, TGCTs arise from primordial germ cells in only the 129 inbred strains, and susceptibility is under multigenic control. The spontaneously arising mutation Ter (ref. 5) on mouse chromosome 18 (Refs 6,7) increases TGCT frequency on a 129/Sv background. We originally used Ter in genetic crosses to identify loci that control tumorigenesis. A genome scan of tumour-bearing progeny from backcrosses between the 129/Sv-Ter/+ and MOLF/Ei strains provided modest evidence that MOLF-derived alleles on chromosome 19 enhance development of bilateral TGCTs (ref. 9). To obtain independent evidence for linkage to the MOLF chromosome, we made an autosomal chromosome substitution strain (CSS; or 'consomic strain') in which chromosome 19 of 129/Sv+/+ was replaced by its MOLF-derived homologue. The unusually high frequency of TGCTs in this CSS (even in the absence of the Ter mutation) provides evidence confirming the genome survey results, identifies linkage for a naturally occurring strain variant allele that confers susceptibility to TGCTs and illustrates the power of CSSs in complex trait analysis.     

Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors

Aberrant methylation of CpG islands and genomic deletion are two predominant mechanisms of gene inactivation in tumorigenesis, but the extent to which they interact is largely unknown. The lack of an integrated approach to study these mechanisms has limited the understanding of tumor genomes and cancer genes. Restriction landmark genomic scanning (RLGS; ref. 1) is useful for global analysis of aberrant methylation of CpG islands, but has not been amenable to alignment with deletion maps because the identity of most RLGS fragments is unknown. Here, we determined the nucleotide sequence and exact chromosomal position of RLGS fragments throughout the genome using the whole chromosome of origin of the fragments and in silico restriction digestion of the human genome sequence. To study the interaction of these gene-inactivation mechanisms in primary brain tumors, we integrated RLGS-based methylation analysis with high-resolution deletion maps from microarray-based comparative genomic hybridization (array CGH; ref. 3). Certain subsets of gene-associated CpG islands were preferentially affected by convergent methylation and deletion, including genes that exhibit tumor-suppressor activity, such as CISH1 (encoding SOCS1; ref. 4), as well as genes such as COE3 that have been missed by traditional non-integrated approaches. Our results show that most aberrant methylation events are focal and independent of deletions, and the rare convergence of these mechanisms can pinpoint biallelic gene inactivation without the use of positional cloning.     

Mutations in ATP2A2, encoding a Ca2+ pump, cause Darier disease

Darier disease (DD) is an autosomal-dominant skin disorder characterized by loss of adhesion between epidermal cells (acantholysis) and abnormal keratinization. Recently we constructed a 2.4-Mb, P1-derived artificial chromosome contig spanning the DD candidate region on chromosome 12q23-24.1. After screening several genes that mapped to this region, we identified mutations in the ATP2A2 gene, which encodes the sarco/endoplasmic reticulum Ca2(+)-ATPase type 2 isoform (SERCA2) and is highly expressed in keratinocytes. Thirteen mutations were identified, including frameshift deletions, in-frame deletions or insertions, splice-site mutations and non-conservative missense mutations in functional domains. Our results demonstrate that mutations in ATP2A2 cause DD and disclose a role for this pump in a Ca(2+)-signalling pathway regulating cell-to-cell adhesion and differentiation of the epidermis.     

Global patterns of human mitochondrial DNA and Y-chromosome structure are not influenced by higher migration rates of females versus males

Global-scale patterns of human population structure may be influenced by the rate of migration among populations that is nearly eight times higher for females than for males. This difference is attributed mainly to the widespread practice of patrilocality, in which women move into their mates' residences after marriage. Here we directly test this hypothesis by comparing global patterns of DNA sequence variation on the Y chromosome and mitochondrial DNA (mtDNA) in the same panel of 389 individuals from ten populations (four from Africa and two each from Europe, Asia and Oceania). We introduce a new strategy to assay Y-chromosome variation that identifies a high density of single-nucleotide polymorphisms, allows complete sequencing of all individuals rather than relying on predetermined markers and provides direct sequence comparisons with mtDNA. We found the overall proportion of between-group variation (Phi(ST)) to be 0.334 for the Y chromosome and 0.382 for mtDNA. Genetic differentiation between populations was similar for the Y chromosome and mtDNA at all geographic scales that we tested. Although patrilocality may be important at the local scale, patterns of genetic structure on the continental and global scales are not shaped by the higher rate of migration among females than among males.