ToppGene Suite for gene list enrichment analysis and candidate gene prioritization

ToppGene Suite (http://toppgene.cchmc.org; this web site is free and open to all users and does not require a login to access) is a one-stop portal for (i) gene list functional enrichment, (ii) candidate gene prioritization using either functional annotations or network analysis and (iii) identification and prioritization of novel disease candidate genes in the interactome. Functional annotation-based disease candidate gene prioritization uses a fuzzy-based similarity measure to compute the similarity between any two genes based on semantic annotations. The similarity scores from individual features are combined into an overall score using statistical meta-analysis. A P-value of each annotation of a test gene is derived by random sampling of the whole genome. The protein–protein interaction network (PPIN)-based disease candidate gene prioritization uses social and Web networks analysis algorithms (extended versions of the PageRank and HITS algorithms, and the K-Step Markov method). We demonstrate the utility of ToppGene Suite using 20 recently reported GWAS-based gene–disease associations (including novel disease genes) representing five diseases. ToppGene ranked 19 of 20 (95%) candidate genes within the top 20%, while ToppNet ranked 12 of 16 (75%) candidate genes among the top 20%

Genetic Affinities of the Central Indian Tribal Populations

Background: The central Indian state Madhya Pradesh is often called as ‘heart of India ’ and has always been an important region functioning as a trinexus belt for three major language families (Indo-European, Dravidian and Austroasiatic). There are less detailed genetic studies on the populations inhabited in this region. Therefore, this study is an attempt for extensive characterization of genetic ancestries of three tribal populations, namely; Bharia, Bhil and Sahariya, inhabiting this region using haploid and diploid DNA markers. Methodology/Principal Findings: Mitochondrial DNA analysis showed high diversity, including some of the older sublineages of M haplogroup and prominent R lineages in all the three tribes. Y-chromosomal biallelic markers revealed high frequency of Austroasiatic-specific M95-O2a haplogroup in Bharia and Sahariya, M82-H1a in Bhil and M17-R1a in Bhil and Sahariya. The results obtained by haploid as well as diploid genetic markers revealed strong genetic affinity of Bharia (a Dravidian speaking tribe) with the Austroasiatic (Munda) group. The gene flow from Austroasiatic group is further confirmed by their Y-STRs haplotype sharing analysis, where we determined their founder haplotype from the North Munda speaking tribe, while, autosomal analysis was largely in concordant with the haploid DNA results. Conclusions/Significance: Bhil exhibited largely Indo-European specific ancestry, while Sahariya and Bharia showed admixed genetic package of Indo-European and Austroasiatic populations. Hence, in a landscape like India, linguistic labe

Targeted Genome-Wide Enrichment of Functional Regions

Only a small fraction of large genomes such as that of the human contains the functional regions such as the exons, promoters, and polyA sites. A platform technique for selective enrichment of functional genomic regions will enable several next-generation sequencing applications that include the discovery of causal mutations for disease and drug response. Here, we describe a powerful platform technique, termed “functional genomic fingerprinting” (FGF), for the multiplexed genomewide isolation and analysis of targeted regions such as the exome, promoterome, or exon splice enhancers. The technique employs a fixed part of a uniquely designed Fixed-Randomized primer, while the randomized part contains all the possible sequence permutations. The Fixed-Randomized primers bind with full sequence complementarity at multiple sites where the fixed sequence (such as the splice signals) occurs within the genome, and multiplex amplify many regions bounded by the fixed sequences (e.g., exons). Notably, validation of this technique using cardiac myosin binding protein-C (MYBPC3) gene as an example strongly supports the application and efficacy of this method. Further, assisted by genomewide computational analyses of such sequences, the FGF technique may provide a unique platform for high-throughput sample production and analysis of targeted genomic regions by the next-generation sequencing techniques, with powerful applications in discovering disease and drug response genes

Cardiac Myosin Binding Protein C and MAP-Kinase Activating Death Domain-Containing Gene Polymorphisms and Diastolic Heart Failure

OBJECTIVE: Myosin binding protein C (MYBPC3) plays a role in ventricular relaxation. The aim of the study was to investigate the association between cardiac myosin binding protein C (MYBPC3) gene polymorphisms and diastolic heart failure (DHF) in a human case-control study. METHODS: A total of 352 participants of 1752 consecutive patients from the National Taiwan University Hospital and its affiliated hospital were enrolled. 176 patients diagnosed with DHF confirmed by echocardiography were recruited. Controls were matched 1-to-1 by age, sex, hypertension, diabetes, renal function and medication use. We genotyped 12 single nucleotide polymorphisms (SNPs) according to HapMap Han Chinese Beijing databank across a 40 kb genetic region containing the MYBPC3 gene and the neighboring DNA sequences to capture 100% of haplotype variance in all SNPs with minor allele frequencies ≥ 5%. We also analyzed associations of these tagging SNPs and haplotypes with DHF and linkage disequilibrium (LD) structure of the MYBPC3 gene. RESULTS: In a single locus analysis, SNP rs2290149 was associated with DHF (allele-specific p = 0.004; permuted p = 0.031). The SNP with a minor allele frequency of 9.4%, had an odds ratio 2.14 (95% CI 1.25-3.66; p = 0.004) for the additive model and 2.06 for the autosomal dominant model (GG+GA : AA, 95% CI 1.17-3.63; p = 0.013), corresponding to a population attributable risk fraction of 12.02%. The haplotypes in a LD block of rs2290149 (C-C-G-C) was also significantly associated with DHF (odds ratio 2.10 (1.53-2.89); permuted p = 0.029). CONCLUSIONS: We identified a SNP (rs2290149) among the tagging SNP set that was significantly associated with early DHF in a Chinese population

Genome-Wide Association Studies of the PR Interval in African Americans

The PR interval on the electrocardiogram reflects atrial and atrioventricular nodal conduction time. The PR interval is heritable, provides important information about arrhythmia risk, and has been suggested to differ among human races. Genome-wide association (GWA) studies have identified common genetic determinants of the PR interval in individuals of European and Asian ancestry, but there is a general paucity of GWA studies in individuals of African ancestry. We performed GWA studies in African American individuals from four cohorts (n = 6,247) to identify genetic variants associated with PR interval duration. Genotyping was performed using the Affymetrix 6.0 microarray. Imputation was performed for 2.8 million single nucleotide polymorphisms (SNPs) using combined YRI and CEU HapMap phase II panels. We observed a strong signal (rs3922844) within the gene encoding the cardiac sodium channel (SCN5A) with genome-wide significant association (p<2.5×10−8) in two of the four cohorts and in the meta-analysis. The signal explained 2% of PR interval variability in African Americans (beta  = 5.1 msec per minor allele, 95% CI  = 4.1–6.1, p = 3×10−23). This SNP was also associated with PR interval (beta = 2.4 msec per minor allele, 95% CI = 1.8–3.0, p = 3×10−16) in individuals of European ancestry (n = 14,042), but with a smaller effect size (p for heterogeneity <0.001) and variability explained (0.5%). Further meta-analysis of the four cohorts identified genome-wide significant associations with SNPs in SCN10A (rs6798015), MEIS1 (rs10865355), and TBX5 (rs7312625) that were highly correlated with SNPs identified in European and Asian GWA studies. African ancestry was associated with increased PR duration (13.3 msec, p = 0.009) in one but not the other three cohorts. Our findings demonstrate the relevance of common variants to African Americans at four loci previously associated with PR interval in European and Asian samples and identify an association signal at one of these loci that is more strongly associated with PR interval in African Americans than in Europeans

Evolutionary Analyses of Entire Genomes Do Not Support the Association of mtDNA Mutations with Ras/MAPK Pathway Syndromes

BACKGROUND: There are several known autosomal genes responsible for Ras/MAPK pathway syndromes, including Noonan syndrome (NS) and related disorders (such as LEOPARD, neurofibromatosis type 1), although mutations of these genes do not explain all cases. Due to the important role played by the mitochondrion in the energetic metabolism of cardiac muscle, it was recently proposed that variation in the mitochondrial DNA (mtDNA) genome could be a risk factor in the Noonan phenotype and in hypertrophic cardiomyopathy (HCM), which is a common clinical feature in Ras/MAPK pathway syndromes. In order to test these hypotheses, we sequenced entire mtDNA genomes in the largest series of patients suffering from Ras/MAPK pathway syndromes analyzed to date (n = 45), most of them classified as NS patients (n = 42). METHODS/PRINCIPAL FINDINGS: The results indicate that the observed mtDNA lineages were mostly of European ancestry, reproducing in a nutshell the expected haplogroup (hg) patterns of a typical Iberian dataset (including hgs H, T, J, and U). Three new branches of the mtDNA phylogeny (H1j1, U5b1e, and L2a5) are described for the first time, but none of these are likely to be related to NS or Ras/MAPK pathway syndromes when observed under an evolutionary perspective. Patterns of variation in tRNA and protein genes, as well as redundant, private and heteroplasmic variants, in the mtDNA genomes of patients were as expected when compared with the patterns inferred from a worldwide mtDNA phylogeny based on more than 8700 entire genomes. Moreover, most of the mtDNA variants found in patients had already been reported in healthy individuals and constitute common polymorphisms in human population groups. CONCLUSIONS/SIGNIFICANCE: As a whole, the observed mtDNA genome variation in the NS patients was difficult to reconcile with previous findings that indicated a pathogenic role of mtDNA variants in NS

Myosin binding protein C: implications for signal-transduction

Myosin binding protein C (MYBPC) is a crucial component of the sarcomere and an important regulator of muscle function. While mutations in different myosin binding protein C (MYBPC) genes are well known causes of various human diseases, such as hypertrophic (HCM) and dilated (DCM) forms of cardiomyopathy as well as skeletal muscular disorders, the underlying molecular mechanisms remain not well understood. A variety of MYBPC3 (cardiac isoform) mutations have been studied in great detail and several corresponding genetically altered mouse models have been generated. Most MYBPC3 mutations may cause haploinsufficiency and with it they may cause a primary increase in calcium sensitivity which is potentially able to explain major features observed in HCM patients such as the hypercontractile phenotype and the well known secondary effects such as myofibrillar disarray, fibrosis, myocardial hypertrophy and remodelling including arrhythmogenesis. However the presence of poison peptides in some cases cannot be fully excluded and most probably other mechanisms are also at play. Here we shall discuss MYBPC interacting proteins and possible pathways linked to cardiomyopathy and heart failure