Biased biological functions of horizontally transferred genes in prokaryotic genomes

Horizontal gene transfer is one of the main mechanisms contributing to microbial genome diversification. To clarify the overall picture of interspecific gene flow among prokaryotes, we developed a new method for detecting horizontally transferred genes and their possible donors by Bayesian inference with training models for nucleotide composition. Our method gives the average posterior probability (horizontal transfer index) for each gene sequence, with a low horizontal transfer index indicating recent horizontal transfer. We found that 14% of open reading frames in 116 prokaryotic complete genomes were subjected to recent horizontal transfer. Based on this data set, we quantitatively determined that the biological functions of horizontally transferred genes, except mobile element genes, are biased to three categories: cell surface, DNA binding and pathogenicity-related functions. Thus, the transferability of genes seems to depend heavily on their functions.     

Restoration of type VII collagen expression and function in dystrophic epidermolysis bullosa

Dystrophic epidermolysis bullosa (DEB) is a family of inherited mechano-bullous disorders caused by mutations in the human type VII collagen gene (COL7A1). Individuals with DEB lack type VII collagen and anchoring fibrils, structures that attach epidermis and dermis. The current lack of treatment for DEB is an impetus to develop gene therapy strategies that efficiently transfer and stably express genes delivered to skin cells in vivo. In this study, we delivered and expressed full-length type VII collagen using a self-inactivating minimal lentivirus-based vector. Transduction of lentiviral vectors containing the COL7A1 transgene into recessive DEB (RDEB) keratinocytes and fibroblasts (in which type VII collagen was absent) resulted in persistent synthesis and secretion of type VII collagen. Unlike RDEB parent cells, the gene-corrected cells had normal morphology, proliferative potential, matrix attachment and motility. We used these gene-corrected cells to regenerate human skin on immune-deficient mice. Human skin regenerated by gene-corrected RDEB cells had restored expression of type VII collagen and formation of anchoring fibrils at the dermal-epidermal junction in vivo. These studies demonstrate that it is possible to restore type VII collagen gene expression in RDEB skin in vivo.     

Epigenome analyses using BAC microarrays identify evolutionary conservation of tissue-specific methylation of SHANK3

CpG islands are present in one-half of all human and mouse genes and typically overlap with promoters or exons. We developed a method for high-resolution analysis of the methylation status of CpG islands genome-wide, using arrays of BAC clones and the methylation-sensitive restriction enzyme NotI. Here we demonstrate the accuracy and specificity of the method. By computationally mapping all NotI sites, methylation events can be defined with single-nucleotide precision throughout the genome. We also demonstrate the unique expandability of the array method using a different methylation-sensitive restriction enzyme, BssHII. We identified and validated new CpG island loci that are methylated in a tissue-specific manner in normal human tissues. The methylation status of the CpG islands is associated with gene expression for several genes, including SHANK3, which encodes a structural protein in neuronal postsynaptic densities. Defects in SHANK3 seem to underlie human 22q13 deletion syndrome. Furthermore, these patterns for SHANK3 are conserved in mice and rats.     

Gene silencing in cancer by histone H3 lysine 27 trimethylation independent of promoter DNA methylation

Epigenetic silencing in cancer cells is mediated by at least two distinct histone modifications, polycomb-based histone H3 lysine 27 trimethylation (H3K27triM) and H3K9 dimethylation. The relationship between DNA hypermethylation and these histone modifications is not completely understood. Using chromatin immunoprecipitation microarrays (ChIP-chip) in prostate cancer cells compared to normal prostate, we found that up to 5% of promoters (16% CpG islands and 84% non-CpG islands) were enriched with H3K27triM. These genes were silenced specifically in prostate cancer, and those CpG islands affected showed low levels of DNA methylation. Downregulation of the EZH2 histone methyltransferase restored expression of the H3K27triM target genes alone or in synergy with histone deacetylase inhibition, without affecting promoter DNA methylation, and with no effect on the expression of genes silenced by DNA hypermethylation. These data establish EZH2-mediated H3K27triM as a mechanism of tumor-suppressor gene silencing in cancer that is potentially independent of promoter DNA methylation.     

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.     

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.     

Combined activation of Ras and Akt in neural progenitors induces glioblastoma formation in mice

Gliomas are the most common primary malignant brain tumours and are classified into four clinical grades, with the most aggressive tumours being grade 4 astrocytomas (also known as glioblastoma multiforme; GBM). Frequent genetic alterations in GBMs (refs 2-5) result in stimulation of common signal transduction pathways involving Ras, Akt and other proteins. It is not known which of these pathways, if any, are sufficient to induce GBM formation. Here we transfer, in a tissue-specific manner, genes encoding activated forms of Ras and Akt to astrocytes and neural progenitors in mice. We found that although neither activated Ras nor Akt alone is sufficient to induce GBM formation, the combination of activated Ras and Akt induces high-grade gliomas with the histological features of human GBMs. These tumours appear to arise after gene transfer to neural progenitors, but not after transfer to differentiated astrocytes. Increased activity of RAS is found in many human GBMs (ref. 11), and we show here that Akt activity is increased in most of these tumours, implying that combined activation of these two pathways accurately models the biology of this disease.     

CTCF-binding sites flank CTG/CAG repeats and form a methylation-sensitive insulator at the DM1 locus

An expansion of a CTG repeat at the DM1 locus causes myotonic dystrophy (DM) by altering the expression of the two adjacent genes, DMPK and SIX5, and through a toxic effect of the repeat-containing RNA. Here we identify two CTCF-binding sites that flank the CTG repeat and form an insulator element between DMPK and SIX5. Methylation of these sites prevents binding of CTCF, indicating that the DM1 locus methylation in congenital DM would disrupt insulator function. Furthermore, CTCF-binding sites are associated with CTG/CAG repeats at several other loci. We suggest a general role for CTG/CAG repeats as components of insulator elements at multiple sites in the human genome.     

Mutations in a novel gene encoding a CRAL-TRIO domain cause human Cayman ataxia and ataxia/dystonia in the jittery mouse

Cayman ataxia is a recessive congenital ataxia restricted to one area of Grand Cayman Island. Comparative mapping suggested that the locus on 19p13.3 associated with Cayman ataxia might be homologous to the locus on mouse chromosome 10 associated with the recessive ataxic mouse mutant jittery. Screening genes in the region of overlap identified mutations in a novel predicted gene in three mouse jittery alleles, including the first mouse mutation caused by an Alu-related (B1 element) insertion. We found two mutations exclusively in all individuals with Cayman ataxia. The gene ATCAY or Atcay encodes a neuron-restricted protein called caytaxin. Caytaxin contains a CRAL-TRIO motif common to proteins that bind small lipophilic molecules. Mutations in another protein containing a CRAL-TRIO domain, alpha-tocopherol transfer protein (TTPA), cause a vitamin E-responsive ataxia. Three-dimensional protein structural modeling predicts that the caytaxin ligand is more polar than vitamin E. Identification of the caytaxin ligand may help develop a therapy for Cayman ataxia.     

Spatial proximity of translocation-prone gene loci in human lymphomas

Cancer cells frequently have disease-specific chromosome rearrangements. It is poorly understood why translocations between chromosomes recur at specific breakpoints in the genome. Here we provide evidence that higher-order spatial genome organization is a contributing factor in the formation of recurrent translocations. We show that MYC, BCL and immunoglobulin loci, which are recurrently translocated in various B-cell lymphomas, are preferentially positioned in close spatial proximity relative to each other in normal B cells. Loci in spatial proximity are non-randomly positioned towards the nuclear interior in normal B cells. This locus proximity is the consequence of higher-order genome structure rather than a property of individual genes. Our results suggest that the formation of specific translocations in human lymphomas, and perhaps other tissues, is determined in part by higher-order spatial organization of the genome.     

Emergence of a new disease as a result of interspecific virulence gene transfer

New diseases of humans, animals and plants emerge regularly. Enhanced virulence on a new host can be facilitated by the acquisition of novel virulence factors. Interspecific gene transfer is known to be a source of such virulence factors in bacterial pathogens (often manifested as pathogenicity islands in the recipient organism) and it has been speculated that interspecific transfer of virulence factors may occur in fungal pathogens. Until now, no direct support has been available for this hypothesis. Here we present evidence that a gene encoding a critical virulence factor was transferred from one species of fungal pathogen to another. This gene transfer probably occurred just before 1941, creating a pathogen population with significantly enhanced virulence and leading to the emergence of a new damaging disease of wheat.     

Genetic dissection and prognostic modeling of overt stroke in sickle cell anemia

Sickle cell anemia (SCA) is a paradigmatic single gene disorder caused by homozygosity with respect to a unique mutation at the beta-globin locus. SCA is phenotypically complex, with different clinical courses ranging from early childhood mortality to a virtually unrecognized condition. Overt stroke is a severe complication affecting 6-8% of individuals with SCA. Modifier genes might interact to determine the susceptibility to stroke, but such genes have not yet been identified. Using Bayesian networks, we analyzed 108 SNPs in 39 candidate genes in 1,398 individuals with SCA. We found that 31 SNPs in 12 genes interact with fetal hemoglobin to modulate the risk of stroke. This network of interactions includes three genes in the TGF-beta pathway and SELP, which is associated with stroke in the general population. We validated this model in a different population by predicting the occurrence of stroke in 114 individuals with 98.2% accuracy.     

Mutations in the pre-replication complex cause Meier-Gorlin syndrome

Meier-Gorlin syndrome (ear, patella and short-stature syndrome) is an autosomal recessive primordial dwarfism syndrome characterized by absent or hypoplastic patellae and markedly small ears1?3. Both pre- and post-natal growth are impaired in this disorder, and although microcephaly is often evident, intellect is usually normal in this syndrome. We report here that individuals with this disorder show marked locus heterogeneity, and we identify mutations in five separate genes: ORC1, ORC4, ORC6, CDT1 and CDC6. All of these genes encode components of the pre-replication complex, implicating defects in replication licensing as the cause of a genetic syndrome with distinct developmental abnormalities.     

Mutations in a gene encoding an ABC transporter cause pseudoxanthoma elasticum

Pseudoxanthoma elasticum (PXE) is a heritable disorder characterized by calcification of elastic fibres in skin, arteries and retina that results in dermal lesions with associated laxity and loss of elasticity, arterial insufficiency and retinal haemorrhages leading to macular degeneration. PXE is usually found as a sporadic disorder, but examples of both autosomal recessive and autosomal dominant forms of PXE have been observed. Partial manifestations of the PXE phenotype have also been described in presumed carriers in PXE families. Linkage of both dominant and recessive forms of PXE to a 5-cM domain on chromosome 16p13.1 has been reported (refs 8,9). We have refined this locus to an 820-kb region containing 6 candidate genes. Here we report the exclusion of five of these genes and the identification of the first mutations responsible for the development of PXE in a gene encoding a protein associated with multidrug resistance (ABCC6).     

Adaptive evolution of bacterial metabolic networks by horizontal gene transfer

Numerous studies have considered the emergence of metabolic pathways, but the modes of recent evolution of metabolic networks are poorly understood. Here, we integrate comparative genomics with flux balance analysis to examine (i) the contribution of different genetic mechanisms to network growth in bacteria, (ii) the selective forces driving network evolution and (iii) the integration of new nodes into the network. Most changes to the metabolic network of Escherichia coli in the past 100 million years are due to horizontal gene transfer, with little contribution from gene duplicates. Networks grow by acquiring genes involved in the transport and catalysis of external nutrients, driven by adaptations to changing environments. Accordingly, horizontally transferred genes are integrated at the periphery of the network, whereas central parts remain evolutionarily stable. Genes encoding physiologically coupled reactions are often transferred together, frequently in operons. Thus, bacterial metabolic networks evolve by direct uptake of peripheral reactions in response to changed environments.