Epigenetic characterization of the early embryo with a chromatin immunoprecipitation protocol applicable to small cell populations

Chromatin immunoprecipitation (ChIP) defines the genomic distribution of proteins and their modifications but is limited by the cell numbers required (ideally >10(7)). Here we describe a protocol that uses carrier chromatin and PCR, 'carrier' ChIP (CChIP), to permit analysis of as few as 100 cells. We assayed histone modifications at key regulator genes (such as Nanog, Pou5f1 (also known as Oct4) and Cdx2) by CChIP in mouse embryonic stem (ES) cells and in inner cell mass (ICM) and trophectoderm of cultured blastocysts. Activating and silencing modifications (H4 acetylation and H3K9 methylation) mark active and silent promoters as predicted, and we find close correlation between values derived from CChIP (1,000 ES cells) and conventional ChIP (5 x 10(7) ES cells). Studies on genes silenced in both ICM and ES cells (Cdx2, Cfc1, Hhex and Nkx2-2, also known as Nkx) show that the intensity of silencing marks is relatively diminished in ES cells, indicating a possible relaxation of some components of silencing on adaptation to culture.     

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.     

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.     

Imprinted microRNA genes transcribed antisense to a reciprocally imprinted retrotransposon-like gene

MicroRNAs (miRNAs) are an abundant class of RNAs that are approximately 21-25 nucleotides (nt) long, interact with mRNAs and trigger either translation repression or RNA cleavage (RNA interference, RNAi) depending on the degree of complementarity with their targets. Here we show that the imprinted mouse distal chromosome 12 locus encodes two miRNA genes expressed from the maternally inherited chromosome and antisense to a retrotransposon-like gene (Rtl1) expressed only from the paternal allele.     

Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease

We report results of a nonsynonymous SNP scan for ulcerative colitis and identify a previously unknown susceptibility locus at ECM1. We also show that several risk loci are common to ulcerative colitis and Crohn's disease (IL23R, IL12B, HLA, NKX2-3 and MST1), whereas autophagy genes ATG16L1 and IRGM, along with NOD2 (also known as CARD15), are specific for Crohn's disease. These data provide the first detailed illustration of the genetic relationship between these common inflammatory bowel diseases.     

MicroRNA-responsive 'sensor' transgenes uncover Hox-like and other developmentally regulated patterns of vertebrate microRNA expression

MicroRNAs (miRNAs) are a class of short ( approximately 22-nt) noncoding RNA molecules that downregulate expression of their mRNA targets. Since their discovery as regulators of developmental timing in Caenorhabditis elegans, hundreds of miRNAs have been identified in both animals and plants. Here, we report a technique for visualizing detailed miRNA expression patterns in mouse embryos. We elucidate the tissue-specific expression of several miRNAs during embryogenesis, including two encoded by genes embedded in homeobox (Hox) clusters, miR-10a and miR-196a. These two miRNAs are expressed in patterns that are markedly reminiscent of those of Hox genes. Furthermore, miR-196a negatively regulates Hoxb8, indicating that its restricted expression pattern probably reflects a role in the patterning function of the Hox complex.     

A mouse model for spinal muscular atrophy

The survival motor neuron gene is present in humans in a telomeric copy, SMN1, and several centromeric copies, SMN2. Homozygous mutation of SMN1 is associated with proximal spinal muscular atrophy (SMA), a severe motor neuron disease characterized by early childhood onset of progressive muscle weakness. To understand the functional role of SMN1 in SMA, we produced mouse lines deficient for mouse Smn and transgenic mouse lines that expressed human SMN2. Smn-/- mice died during the peri-implantation stage. In contrast, transgenic mice harbouring SMN2 in the Smn-/- background showed pathological changes in the spinal cord and skeletal muscles similar to those of SMA patients. The severity of the pathological changes in these mice correlated with the amount of SMN protein that contained the region encoded by exon 7. Our results demonstrate that SMN2 can partially compensate for lack of SMN1. The variable phenotypes of Smn-/-SMN2 mice reflect those seen in SMA patients, providing a mouse model for this disease.     

Bayesian method to predict individual SNP genotypes from gene expression data

RNA profiling can be used to capture the expression patterns of many genes that are associated with expression quantitative trait loci (eQTLs). Employing published putative cis eQTLs, we developed a Bayesian approach to predict SNP genotypes that is based only on RNA expression data. We show that predicted genotypes can accurately and uniquely identify individuals in large populations. When inferring genotypes from an expression data set using eQTLs of the same tissue type (but from an independent cohort), we were able to resolve 99% of the identities of individuals in the cohort at P(adjusted) ≤ 1 × 10(-5). When eQTLs derived from one tissue were used to predict genotypes using expression data from a different tissue, the identities of 90% of the study subjects could be resolved at P(adjusted) ≤ 1 × 10(-5). We discuss the implications of deriving genotypic information from RNA data deposited in the public domain.     

Common variants near MBNL1 and NKX2-5 are associated with infantile hypertrophic pyloric stenosis

Infantile hypertrophic pyloric stenosis (IHPS) is a severe condition characterized by hypertrophy of the pyloric sphincter muscle. We conducted a genome-wide association study (GWAS) on 1,001 surgery-confirmed cases and 2,401 controls from Denmark. The six most strongly associated loci were tested in a replication set of 796 cases and 876 controls. Three SNPs reached genome-wide significance. One of these SNPs, rs11712066 (odds ratio (OR) = 1.61; P = 1.5 × 10(-17)) at 3p25.1, is located 150 kb upstream of MBNL1, which encodes a factor that regulates splicing transitions occurring shortly after birth. The second SNP, rs573872 (OR = 1.41; P = 4.3 × 10(-12)), maps to an intergenic region at 3p25.2 approximately 1.3 Mb downstream of MBNL1. The third SNP, rs29784 (OR = 1.42; P = 1.5 × 10(-15)) at 5q35.2, is 64 kb downstream of NKX2-5, which is involved in development of cardiac muscle tissue and embryonic gut development.     

Application of comparative functional genomics to identify best-fit mouse models to study human cancer

Genetically modified mice have been extensively used for analyzing the molecular events that occur during tumor development. In many, if not all, cases, however, it is uncertain to what extent the mouse models reproduce features observed in the corresponding human conditions. This is due largely to lack of precise methods for direct and comprehensive comparison at the molecular level of the mouse and human tumors. Here we use global gene expression patterns of 68 hepatocellular carcinomas (HCCs) from seven different mouse models and 91 human HCCs from predefined subclasses to obtain direct comparison of the molecular features of mouse and human HCCs. Gene expression patterns in HCCs from Myc, E2f1 and Myc E2f1 transgenic mice were most similar to those of the better survival group of human HCCs, whereas the expression patterns in HCCs from Myc Tgfa transgenic mice and in diethylnitrosamine-induced mouse HCCs were most similar to those of the poorer survival group of human HCCs. Gene expression patterns in HCCs from Acox1(-/-) mice and in ciprofibrate-induced HCCs were least similar to those observed in human HCCs. We conclude that our approach can effectively identify appropriate mouse models to study human cancers.     

A functional switch from lung cancer resistance to susceptibility at the Pas1 locus in Kras2LA2 mice

Pulmonary adenoma susceptibility 1 (Pas1) is the major mouse lung cancer susceptibility locus on chromosome 6 (ref. 1). Kras2 is a common target of somatic mutation in chemically induced mouse lung tumors and is a candidate Pas1 gene. M. spretus mice (SPRET/Ei) carry a Pas1 resistance haplotype for chemically induced lung tumors. We demonstrate that the SPRET/Ei Pas1 allele is switched from resistance to susceptibility by fixation of the parental origin of the mutant Kras2 allele. This switch correlates with low expression of endogenous Kras2 in SPRET/Ei. We propose that the Pas1 modifier effect is due to Kras2, and that a sensitive balance between the expression levels of wild-type and mutant alleles determines lung tumor susceptibility. These data demonstrate that cancer predisposition should also be considered in the context of somatic events and could have major implications for the design of human association studies to identify cancer susceptibility genes.