The role of DNA methylation in setting up chromatin structure during development

DNA methylation inhibits gene expression in animal cells, probably by affecting chromatin structure. Biochemical studies suggest that this process may be mediated by methyl-specific binding proteins that recruit enzymatic machinery capable of locally altering histone modification. To test whether DNA methylation actually has a role in the assembly of chromatin during normal development, we used cell transfection and a transgene construct genetically programmed to be either methylated or unmethylated in all cell types of the mouse. Chromatin immunoprecipitation (ChIP) analysis shows that the presence of DNA methylation brings about the deacetylation of histone H4 and methylation of Lys9 of histone H3 (H3 Lys9) and prevents methylation of Lys4 of histone H3 (H3 Lys4), thus generating a structure identical to that of methylated sequences in the genome. These results indicate that the methylation pattern established in early embryogenesis is profoundly important in setting up the structural profile of the genome.     

Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodeling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.     

Dynamic assembly of silent chromatin during thymocyte maturation

Considerable knowledge has been gained from temporal analyses of molecular events culminating in gene activation, but technical hurdles have hindered comparable studies of gene silencing. Here we describe the temporal assembly of silent chromatin at the mouse terminal transferase gene (Dntt), which is silenced and repositioned to pericentromeric heterochromatin during thymocyte maturation. Silencing was nucleated at the Dntt promoter by the ordered deacetylation of histone H3 at Lys9 (H3-Lys9), loss of methylation at H3-Lys4 and acquisition of methylation at H3-Lys9, followed by bidirectional spreading of each event. Deacetylation at H3-Lys9 coincided with pericentromeric repositioning, and neither of these early events required de novo protein synthesis. CpG methylation increased primarily in mature T cells that had left the thymus. A transformed thymocyte line supported reversible inactivation of Dntt without repositioning. In these cells, histone modification changes were nucleated at the promoter but did not spread. These results provide a foundation for elucidating the mechanisms of silent chromatin assembly during development.     

Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide. We sequenced and analyzed the whole genomes of 27 HCCs, 25 of which were associated with hepatitis B or C virus infections, including two sets of multicentric tumors. Although no common somatic mutations were identified in the multicentric tumor pairs, their whole-genome substitution patterns were similar, suggesting that these tumors developed from independent mutations, although their shared etiological backgrounds may have strongly influenced their somatic mutation patterns. Statistical and functional analyses yielded a list of recurrently mutated genes. Multiple chromatin regulators, including ARID1A, ARID1B, ARID2, MLL and MLL3, were mutated in ～50% of the tumors. Hepatitis B virus genome integration in the TERT locus was frequently observed in a high clonal proportion. Our whole-genome sequencing analysis of HCCs identified the influence of etiological background on somatic mutation patterns and subsequent carcinogenesis, as well as recurrent mutations in chromatin regulators in HCCs.     

Genomic analysis of Bartonella identifies type IV secretion systems as host adaptability factors

The bacterial genus Bartonella comprises 21 pathogens causing characteristic intraerythrocytic infections. Bartonella bacilliformis is a severe pathogen representing an ancestral lineage, whereas the other species are benign pathogens that evolved by radial speciation. Here, we have used comparative and functional genomics to infer pathogenicity genes specific to the radiating lineage, and we suggest that these genes may have facilitated adaptation to the host environment. We determined the complete genome sequence of Bartonella tribocorum by shotgun sequencing and functionally identified 97 pathogenicity genes by signature-tagged mutagenesis. Eighty-one pathogenicity genes belong to the core genome (1,097 genes) of the radiating lineage inferred from genome comparison of B. tribocorum, Bartonella henselae and Bartonella quintana. Sixty-six pathogenicity genes are present in B. bacilliformis, and one has been lost by deletion. The 14 pathogenicity genes specific for the radiating lineage encode two laterally acquired type IV secretion systems, suggesting that these systems have a role in host adaptability.     

Nonsense-mediated decay microarray analysis identifies mutations of EPHB2 in human prostate cancer

The identification of tumor-suppressor genes in solid tumors by classical cancer genetics methods is difficult and slow. We combined nonsense-mediated RNA decay microarrays and array-based comparative genomic hybridization for the genome-wide identification of genes with biallelic inactivation involving nonsense mutations and loss of the wild-type allele. This approach enabled us to identify previously unknown mutations in the receptor tyrosine kinase gene EPHB2. The DU 145 prostate cancer cell line, originating from a brain metastasis, carries a truncating mutation of EPHB2 and a deletion of the remaining allele. Additional frameshift, splice site, missense and nonsense mutations are present in clinical prostate cancer samples. Transfection of DU 145 cells, which lack functional EphB2, with wild-type EPHB2 suppresses clonogenic growth. Taken together with studies indicating that EphB2 may have an essential role in cell migration and maintenance of normal tissue architecture, our findings suggest that mutational inactivation of EPHB2 may be important in the progression and metastasis of prostate cancer.     

Whole-genome analysis of diverse Chlamydia trachomatis strains identifies phylogenetic relationships masked by current clinical typing

Chlamydia trachomatis is responsible for both trachoma and sexually transmitted infections, causing substantial morbidity and economic cost globally. Despite this, our knowledge of its population and evolutionary genetics is limited. Here we present a detailed phylogeny based on whole-genome sequencing of representative strains of C. trachomatis from both trachoma and lymphogranuloma venereum (LGV) biovars from temporally and geographically diverse sources. Our analysis shows that predicting phylogenetic structure using ompA, which is traditionally used to classify Chlamydia, is misleading because extensive recombination in this region masks any true relationships present. We show that in many instances, ompA is a chimera that can be exchanged in part or as a whole both within and between biovars. We also provide evidence for exchange of, and recombination within, the cryptic plasmid, which is another key diagnostic target. We used our phylogenetic framework to show how genetic exchange has manifested itself in ocular, urogenital and LGV C. trachomatis strains, including the epidemic LGV serotype L2b.     

Exon capture analysis of G protein-coupled receptors identifies activating mutations in GRM3 in melanoma

G protein-coupled receptors (GPCRs), the largest human gene family, are important regulators of signaling pathways. However, knowledge of their genetic alterations is limited. In this study, we used exon capture and massively parallel sequencing methods to analyze the mutational status of 734 GPCRs in melanoma. This investigation revealed that one family member, GRM3, was frequently mutated and that one of its mutations clustered within one position. Biochemical analysis of GRM3 alterations revealed that mutant GRM3 selectively regulated the phosphorylation of MEK, leading to increased anchorage-independent growth and migration. Melanoma cells expressing mutant GRM3 had reduced cell growth and cellular migration after short hairpin RNA-mediated knockdown of GRM3 or treatment with a selective MEK inhibitor, AZD-6244, which is currently being used in phase 2 clinical trials. Our study yields the most comprehensive map of genetic alterations in the GPCR gene family.     

Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits

We present an approximate conditional and joint association analysis that can use summary-level statistics from a meta-analysis of genome-wide association studies (GWAS) and estimated linkage disequilibrium (LD) from a reference sample with individual-level genotype data. Using this method, we analyzed meta-analysis summary data from the GIANT Consortium for height and body mass index (BMI), with the LD structure estimated from genotype data in two independent cohorts. We identified 36 loci with multiple associated variants for height (38 leading and 49 additional SNPs, 87 in total) via a genome-wide SNP selection procedure. The 49 new SNPs explain approximately 1.3% of variance, nearly doubling the heritability explained at the 36 loci. We did not find any locus showing multiple associated SNPs for BMI. The method we present is computationally fast and is also applicable to case-control data, which we demonstrate in an example from meta-analysis of type 2 diabetes by the DIAGRAM Consortium.     

Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4

We conducted a combined genome-wide association study (GWAS) of 7,481 individuals with bipolar disorder (cases) and 9,250 controls as part of the Psychiatric GWAS Consortium. Our replication study tested 34 SNPs in 4,496 independent cases with bipolar disorder and 42,422 independent controls and found that 18 of 34 SNPs had P < 0.05, with 31 of 34 SNPs having signals with the same direction of effect (P = 3.8 × 10(-7)). An analysis of all 11,974 bipolar disorder cases and 51,792 controls confirmed genome-wide significant evidence of association for CACNA1C and identified a new intronic variant in ODZ4. We identified a pathway comprised of subunits of calcium channels enriched in bipolar disorder association intervals. Finally, a combined GWAS analysis of schizophrenia and bipolar disorder yielded strong association evidence for SNPs in CACNA1C and in the region of NEK4-ITIH1-ITIH3-ITIH4. Our replication results imply that increasing sample sizes in bipolar disorder will confirm many additional loci.     

Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy. Here, we performed high-resolution copy-number analysis on 125 HCC tumors and whole-exome sequencing on 24 of these tumors. We identified 135 homozygous deletions and 994 somatic mutations of genes with predicted functional consequences. We found new recurrent alterations in four genes (ARID1A, RPS6KA3, NFE2L2 and IRF2) not previously described in HCC. Functional analyses showed tumor suppressor properties for IRF2, whose inactivation, exclusively found in hepatitis B virus (HBV)-related tumors, led to impaired TP53 function. In contrast, inactivation of chromatin remodelers was frequent and predominant in alcohol-related tumors. Moreover, association of mutations in specific genes (RPS6KA3-AXIN1 and NFE2L2-CTNNB1) suggested that Wnt/β-catenin signaling might cooperate in liver carcinogenesis with both oxidative stress metabolism and Ras/mitogen-activated protein kinase (MAPK) pathways. This study provides insight into the somatic mutational landscape in HCC and identifies interactions between mutations in oncogene and tumor suppressor gene mutations related to specific risk factors.     

Targets and dynamics of promoter DNA methylation during early mouse development

DNA methylation is extensively reprogrammed during the early phases of mammalian development, yet genomic targets of this process are largely unknown. We optimized methylated DNA immunoprecipitation for low numbers of cells and profiled DNA methylation during early development of the mouse embryonic lineage in vivo. We observed a major epigenetic switch during implantation at the transition from the blastocyst to the postimplantation epiblast. During this period, DNA methylation is primarily targeted to repress the germline expression program. DNA methylation in the epiblast is also targeted to promoters of lineage-specific genes such as hematopoietic genes, which are subsequently demethylated during terminal differentiation. De novo methylation during early embryogenesis is catalyzed by Dnmt3b, and absence of DNA methylation leads to ectopic gene activation in the embryo. Finally, we identify nonimprinted genes that inherit promoter DNA methylation from parental gametes, suggesting that escape of post-fertilization DNA methylation reprogramming is prevalent in the mouse genome.     

Conditional inactivation of Fgf4 reveals complexity of signalling during limb bud development

Development of the vertebrate limb bud depends on reciprocal interactions between the zone of polarizing activity (ZPA) and the apical ectodermal ridge (AER). Sonic hedgehog (SHH) and fibroblast growth factors (FGFs) are key signalling molecules produced in the ZPA and AER, respectively. Experiments in chicks suggested that SHH expression in the ZPA is maintained by FGF4 expression in the AER, and vice versa, providing a molecular mechanism for coordinating the activities of these two signalling centres. This SHH/FGF4 feedback loop model is supported by genetic evidence showing that Fgf4 expression is not maintained in Shh-/- mouse limbs. We report here that Shh expression is maintained and limb formation is normal when Fgf4 is inactivated in mouse limbs, thus contradicting the model. We also found that maintenance of Fgf9 and Fgf17 expression is dependent on Shh, whereas Fgf8 expression is not. We discuss a model in which no individual Fgf expressed in the AER (AER-Fgf) is solely necessary to maintain Shh expression, but, instead, the combined activities of two or more AER-Fgfs function in a positive feedback loop with Shh to control limb development.