A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Nephron number varies widely between 0.3 and 1.3 million per kidney in humans. During fetal life, the rate of nephrogenesis is influenced by local retinoic acid (RA) level such that even moderate maternal vitamin A deficiency limits the final nephron number in rodents. Inactivation of genes in the RA pathway causes renal agenesis in mice; however, the impact of retinoids on human kidney development is unknown. To resolve this, we tested for associations between variants of genes involved in RA metabolism (ALDH1A2, CYP26A1, and CYP26B1) and kidney size among normal newborns. Homozygosity for a common (1 in 5) variant, rs7169289(G), within an Sp1 transcription factor motif of the ALDH1A2 gene, showed a significant 22% increase in newborn kidney volume when adjusted for body surface area. Infants bearing this allele had higher umbilical cord blood RA levels compared to those with homozygous wild-type ALDH1A2 rs7169289(A) alleles. Furthermore, the effect of the rs7169289(G) variant was evident in subgroups with or without a previously reported hypomorphic RET 1476(A) proto-oncogene allele that is critical in determining final nephron number. As maternal vitamin A deficiency is widespread in developing countries and may compromise availability of retinol for fetal RA synthesis, our study suggests that the ALDH1A2 rs7169289(G) variant might be protective for such individuals.

Familial clustering of colorectal cancer occurs in 15-20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer.

Parent-of-origin-dependent expression of alleles, imprinting, has been suggested to impact a substantial proportion of mammalian genes. Its discovery requires allele-specific detection of expressed transcripts, but in some cases detected allelic expression bias has been interpreted as imprinting without demonstrating compatible transmission patterns and excluding heritable variation. Therefore, we utilized a genome-wide tool exploiting high density genotyping arrays in parallel measurements of genotypes in RNA and DNA to determine allelic expression across the transcriptome in lymphoblastoid cell lines (LCLs) and skin fibroblasts derived from families.

Recurrent mutations affecting the histone H3.3 residues Lys27 or indirectly Lys36 are frequent drivers of pediatric high-grade gliomas (over 30% of HGGs). To identify additional driver mutations in HGGs, we investigated a cohort of 60 pediatric HGGs using whole-exome sequencing (WES) and compared them to 543 exomes from non-cancer control samples. We identified mutations in SETD2, a H3K36 trimethyltransferase, in 15% of pediatric HGGs, a result that was genome-wide significant (FDR = 0.029). Most SETD2 alterations were truncating mutations. Sequencing the gene in this cohort and another validation cohort (123 gliomas from all ages and grades) showed SETD2 mutations to be specific to high-grade tumors affecting 15% of pediatric HGGs (11/73) and 8% of adult HGGs (5/65) while no SETD2 mutations were identified in low-grade diffuse gliomas (0/45). Furthermore, SETD2 mutations were mutually exclusive with H3F3A mutations in HGGs (P = 0.0492) while they partly overlapped with IDH1 mutations (4/14), and SETD2-mutant tumors were found exclusively in the cerebral hemispheres (P = 0.0055). SETD2 is the only H3K36 trimethyltransferase in humans, and SETD2-mutant tumors showed a substantial decrease in H3K36me3 levels (P < 0.001), indicating that the mutations are loss-of-function. These data suggest that loss-of-function SETD2 mutations occur in older children and young adults and are specific to HGG of the cerebral cortex, similar to the H3.3 G34R/V and IDH mutations. Taken together, our results suggest that mutations disrupting the histone code at H3K36, including H3.3 G34R/V, IDH1 and/or SETD2 mutations, are central to the genesis of hemispheric HGGs in older children and young adults.

New technologies have enabled genome-wide association studies to be conducted with hundreds of thousands of genotyped SNPs. Several different first-generation genome-wide panels of SNPs have been commercialized. The total amount of common genetic variation is still unknown; however, the coverage of commercial panels can be evaluated against reference population samples genotyped by the International HapMap project. Less information is available about coverage in samples from other populations.

Coronary artery disease (CAD) is a major health concern in both developed and developing countries. With a heritability estimated at ~50%, there is a strong rationale to better define the genetic contribution to CAD. This project involves the analysis of 884 individuals from 142 families (with average sibships of 5.7) as well as 558 case and control subjects from the Saguenay Lac St-Jean region of northeastern Quebec, with the use of 1,536 single-nucleotide polymorphisms (SNPs) in 103 candidate genes for CAD. By use of clusters of SNPs to generate multiallelic haplotypes at candidate loci for segregation studies within families, suggestive linkage for high-density lipoprotein (HDL) cholesterol is observed on chromosome 1p36.22. Furthermore, several associations that remain significant after Bonferroni correction are observed with lipoprotein-related traits as well as plasma concentrations of adiponectin. Of note, HDL cholesterol levels are associated with an amino acid substitution (lysine/asparagine) at codon 198 (rs5370) of endothelin-1 (EDN1) in a sex-specific manner, as well as with a SNP (rs2292318) located 7.7 kb upstream of lecithin cholesterol acyl-transferase (LCAT). Whereas the other observed associations are described in the current literature, these two are new. Using an independent validation sample of 806 individuals, we confirm the EDN1 association (P<.005), whereas the LCAT association was nonsignificant (P=.12).

To date, epimutations reported in man have been somatic and erased in germlines. Here, we identify a cause of the autosomal recessive cblC class of inborn errors of vitamin B12 metabolism that we name "epi-cblC". The subjects are compound heterozygotes for a genetic mutation and for a promoter epimutation, detected in blood, fibroblasts, and sperm, at the MMACHC locus; 5-azacytidine restores the expression of MMACHC in fibroblasts. MMACHC is flanked by CCDC163P and PRDX1, which are in the opposite orientation. The epimutation is present in three generations and results from PRDX1 mutations that force antisense transcription of MMACHC thereby possibly generating a H3K36me3 mark. The silencing of PRDX1 transcription leads to partial hypomethylation of the epiallele and restores the expression of MMACHC. This example of epi-cblC demonstrates the need to search for compound epigenetic-genetic heterozygosity in patients with typical disease manifestation and genetic heterozygosity in disease-causing genes located in other gene trios.

We recently reported that atypical teratoid rhabdoid tumors (ATRTs) comprise at least two transcriptional subtypes with different clinical outcomes; however, the mechanisms underlying therapeutic heterogeneity remained unclear. In this study, we analyzed 191 primary ATRTs and 10 ATRT cell lines to define the genomic and epigenomic landscape of ATRTs and identify subgroup-specific therapeutic targets. We found ATRTs segregated into three epigenetic subgroups with distinct genomic profiles, SMARCB1 genotypes, and chromatin landscape that correlated with differential cellular responses to a panel of signaling and epigenetic inhibitors. Significantly, we discovered that differential methylation of a PDGFRB-associated enhancer confers specific sensitivity of group 2 ATRT cells to dasatinib and nilotinib, and suggest that these are promising therapies for this highly lethal ATRT subtype.

Multiple sclerosis (MS) is a chronic disease of the central nervous system responsible for a large portion of neurological disabilities in young adults. Similar to what occurs in numerous complex diseases, both unknown environmental factors and genetic predisposition are required to generate MS. We ascertained a set of 63 Finnish MS families, originating from a high-risk region of the country, to identify a susceptibility gene within the previously established 3.4-Mb region on 17q24. Initial single nucleotide polymorphism (SNP)-based association implicated PRKCA (protein kinase C alpha) gene, and this association was replicated in an independent set of 148 Finnish MS families (p = 0.0004; remaining significant after correction for multiple testing). Further, a dense set of 211 SNPs evenly covering the PRKCA gene and the flanking regions was selected from the dbSNP database and analyzed in two large, independent MS cohorts: in 211 Finnish and 554 Canadian MS families. A multipoint SNP analysis indicated linkage to PRKCA and its telomeric flanking region in both populations, and SNP haplotype and genotype combination analyses revealed an allelic variant of PRKCA, which covers the region between introns 3 and 8, to be over-represented in Finnish MS cases (odds ratio = 1.34, 95% confidence interval 1.07-1.68). A second allelic variant, covering the same region of the PRKCA gene, showed somewhat stronger evidence for association in the Canadian families (odds ratio = 1.64, 95% confidence interval 1.39-1.94). Initial functional relevance for disease predisposition was suggested by the expression analysis: The transcript levels of PRKCA showed correlation with the copy number of the Finnish and Canadian "risk" haplotypes in CD4-negative mononuclear cells of five Finnish multiplex families and in lymphoblast cell lines of 11 Centre d'Etude du Polymorphisme Humain (CEPH) individuals of European origin.

Left-sided congenital heart disease (CHD) encompasses a spectrum of malformations that range from bicuspid aortic valve to hypoplastic left heart syndrome. It contributes significantly to infant mortality and has serious implications in adult cardiology. Although left-sided CHD is known to be highly heritable, the underlying genetic determinants are largely unidentified. In this study, we sought to determine the impact of structural genomic variation on left-sided CHD and compared multiplex families (464 individuals with 174 affecteds (37.5%) in 59 multiplex families and 8 trios) to 1,582 well-phenotyped controls. 73 unique inherited or de novo CNVs in 54 individuals were identified in the left-sided CHD cohort. After stringent filtering, our gene inventory reveals 25 new candidates for LS-CHD pathogenesis, such as SMC1A, MFAP4, and CTHRC1, and overlaps with several known syndromic loci. Conservative estimation examining the overlap of the prioritized gene content with CNVs present only in affected individuals in our cohort implies a strong effect for unique CNVs in at least 10% of left-sided CHD cases. Enrichment testing of gene content in all identified CNVs showed a significant association with angiogenesis. In this first family-based CNV study of left-sided CHD, we found that both co-segregating and de novo events associate with disease in a complex fashion at structural genomic level. Often viewed as an anatomically circumscript disease, a subset of left-sided CHD may in fact reflect more general genetic perturbations of angiogenesis and/or vascular biology.

Adaptor protein (AP) complexes regulate clathrin-coated vesicle assembly, protein cargo sorting, and vesicular trafficking between organelles in eukaryotic cells. Because disruption of the various subunits of the AP complexes is embryonic lethal in the majority of cases, characterization of their function in vivo is still lacking. Here, we describe the first mutation in the human AP1S1 gene, encoding the small subunit sigma1A of the AP-1 complex. This founder splice mutation, which leads to a premature stop codon, was found in four families with a unique syndrome characterized by mental retardation, enteropathy, deafness, peripheral neuropathy, ichthyosis, and keratodermia (MEDNIK). To validate the pathogenic effect of the mutation, we knocked down Ap1s1 expression in zebrafish using selective antisens morpholino oligonucleotides (AMO). The knockdown phenotype consisted of perturbation in skin formation, reduced pigmentation, and severe motility deficits due to impaired neural network development. Both neural and skin defects were rescued by co-injection of AMO with wild-type (WT) human AP1S1 mRNA, but not by co-injecting the truncated form of AP1S1, consistent with a loss-of-function effect of this mutation. Together, these results confirm AP1S1 as the gene responsible for MEDNIK syndrome and demonstrate a critical role of AP1S1 in development of the skin and spinal cord.

To determine if 16 single nucleotide polymorphisms (SNPs) associated with type 2 diabetes (T2DM) in Europeans are also associated with T2DM in South Asians and Latinos and if they can add to the prediction of incident T2DM in a high-risk population.

Diffuse intrinsic pontine glioma (DIPG) is a fatal brain cancer that arises in the brainstem of children, with no effective treatment and near 100% fatality. The failure of most therapies can be attributed to the delicate location of these tumors and to the selection of therapies on the basis of assumptions that DIPGs are molecularly similar to adult disease. Recent studies have unraveled the unique genetic makeup of this brain cancer, with nearly 80% found to harbor a p.Lys27Met histone H3.3 or p.Lys27Met histone H3.1 alteration. However, DIPGs are still thought of as one disease, with limited understanding of the genetic drivers of these tumors. To understand what drives DIPGs, we integrated whole-genome sequencing with methylation, expression and copy number profiling, discovering that DIPGs comprise three molecularly distinct subgroups (H3-K27M, silent and MYCN) and uncovering a new recurrent activating mutation affecting the activin receptor gene ACVR1 in 20% of DIPGs. Mutations in ACVR1 were constitutively activating, leading to SMAD phosphorylation and increased expression of the downstream activin signaling targets ID1 and ID2. Our results highlight distinct molecular subgroups and novel therapeutic targets for this incurable pediatric cancer.

Mycobacterium avium subsp. paratuberculosis (MAP) is an obligate intracellular pathogen that infects many ruminant species. The acquisition of foreign genes via horizontal gene transfer has been postulated to contribute to its pathogenesis, as these genetic elements are absent from its putative ancestor, M. avium subsp. hominissuis (MAH), an environmental organism with lesser pathogenicity. In this study, high-throughput sequencing of MAP transposon libraries were analyzed to qualitatively and quantitatively determine the contribution of individual genes to bacterial survival during infection.

The pufferfish Fugu rubripes (Fugu) with its compact genome is increasingly recognized as an important vertebrate model for comparative genomic studies. In particular, large regions of conserved synteny between human and Fugu genomes indicate its utility to identify disease-causing genes. The human chromosome 12p12 is frequently deleted in various hematological malignancies and solid tumors, but the actual tumor suppressor gene remains unidentified.

Clinical laboratories are adopting array genomic hybridization as a standard clinical test. A number of whole genome array genomic hybridization platforms are available, but little is known about their comparative performance in a clinical context.

One third of humans are infected lifelong with the brain-dwelling, protozoan parasite, Toxoplasma gondii. Approximately fifteen million of these have congenital toxoplasmosis. Although neurobehavioral disease is associated with seropositivity, causality is unproven. To better understand what this parasite does to human brains, we performed a comprehensive systems analysis of the infected brain: We identified susceptibility genes for congenital toxoplasmosis in our cohort of infected humans and found these genes are expressed in human brain. Transcriptomic and quantitative proteomic analyses of infected human, primary, neuronal stem and monocytic cells revealed effects on neurodevelopment and plasticity in neural, immune, and endocrine networks. These findings were supported by identification of protein and miRNA biomarkers in sera of ill children reflecting brain damage and T. gondii infection. These data were deconvoluted using three systems biology approaches: "Orbital-deconvolution" elucidated upstream, regulatory pathways interconnecting human susceptibility genes, biomarkers, proteomes, and transcriptomes. "Cluster-deconvolution" revealed visual protein-protein interaction clusters involved in processes affecting brain functions and circuitry, including lipid metabolism, leukocyte migration and olfaction. Finally, "disease-deconvolution" identified associations between the parasite-brain interactions and epilepsy, movement disorders, Alzheimer's disease, and cancer. This "reconstruction-deconvolution" logic provides templates of progenitor cells' potentiating effects, and components affecting human brain parasitism and diseases.

Huntington disease (HD) is an autosomal-dominant disorder that results from >or=36 CAG repeats in the HD gene (HTT). Approximately 10% of patients inherit a chromosome that underwent CAG expansion from an unaffected parent with <36 CAG repeats. This study is a comprehensive analysis of genetic diversity in HTT and reveals that HD patients of European origin (n = 65) have a significant enrichment (95%) of a specific set of 22 tagging single nucleotide polymorphisms (SNPs) that constitute a single haplogroup. The disease association of many SNPs is much stronger than any previously reported polymorphism and was confirmed in a replication cohort (n = 203). Importantly, the same haplogroup is also significantly enriched (83%) in individuals with 27-35 CAG repeats (intermediate alleles, n = 66), who are unaffected by the disease, but have increased CAG tract sizes relative to the general population (n = 116). These data support a stepwise model for CAG expansion into the affected range (>or=36 CAG) and identifies specific haplogroup variants in the general population associated with this instability. The specific variants at risk for CAG expansion are not present in the general population in China, Japan, and Nigeria where the prevalence of HD is much lower. The current data argue that cis-elements have a major predisposing influence on CAG instability in HTT. The strong association between specific SNP alleles and CAG expansion also provides an opportunity of personalized therapeutics in HD where the clinical development of only a small number of allele-specific targets may be sufficient to treat up to 88% of the HD patient population.

The Haplotype Map (HapMap) project recently generated genotype data for more than 1 million single-nucleotide polymorphisms (SNPs) in four population samples. The main application of the data is in the selection of tag single-nucleotide polymorphisms (tSNPs) to use in association studies. The usefulness of this selection process needs to be verified in populations outside those used for the HapMap project. In addition, it is not known how well the data represent the general population, as only 90-120 chromosomes were used for each population and since the genotyped SNPs were selected so as to have high frequencies. In this study, we analyzed more than 1,000 individuals from Estonia. The population of this northern European country has been influenced by many different waves of migrations from Europe and Russia. We genotyped 1,536 randomly selected SNPs from two 500-kbp ENCODE regions on Chromosome 2. We observed that the tSNPs selected from the CEPH (Centre d'Etude du Polymorphisme Humain) from Utah (CEU) HapMap samples (derived from US residents with northern and western European ancestry) captured most of the variation in the Estonia sample. (Between 90% and 95% of the SNPs with a minor allele frequency of more than 5% have an r2 of at least 0.8 with one of the CEU tSNPs.) Using the reverse approach, tags selected from the Estonia sample could almost equally well describe the CEU sample. Finally, we observed that the sample size, the allelic frequency, and the SNP density in the dataset used to select the tags each have important effects on the tagging performance. Overall, our study supports the use of HapMap data in other Caucasian populations, but the SNP density and the bias towards high-frequency SNPs have to be taken into account when designing association studies.