Characterizing redescriptions using persistent homology to isolate genetic pathways contributing to pathogenesis

Background: Complex diseases may have multiple pathways leading to disease. E.g. coronary artery disease evolves from arterial damage to their epithelial layers, but has multiple causal pathways. More challenging, those pathways are highly correlated within metabolic syndrome. The challenge is to identify specific clusters of phenotype characteristics (composite phenotypes) that may reflect these different etiologies. Further, GWAS seeking to identify SNPs satisfying multiple composite phenotype descriptions allows for lower false positive rates at lower α thresholds, allowing for the possibility of reducing false negatives. This may provide a window into the missing heritability problem. Methods: We identify significant phenotype patterns, and identify fuzzy redescriptions among those patterns using Jaccard distances. Further, we construct Vietoris-Rips complexes from the Jaccard distances and compute the persistent homology associated with those. The patterns comprising these topological features are identified as composite phenotpyes, whose genetic associations are explored with logistic regression applied to pathways and to GWAS. Results: We identified several phenotypes that tended to be dominated by metabolic syndrome descriptions, and which were distinct among the combinations of metabolic syndrome conditions. Among SNPs marking the RAAS complex, various SNPs associated specifically with different groups of composite phenotypes, as well as distinguishing between the composite phenotypes and simple phenotypes. Each of these showed different genetic associations, namely rs6693954, rs762551, rs1378942, and rs1133323. GWAS identified SNPs that associated with composite phenotypes included rs12365545, rs6847235, and rs701319. Eighteen GWAS identified SNPs appeared in combinations supported in composite combinations with greater power than for any individual phenotype. Conclusions: We do find systematic associations among metabolic syndrome variates that show distinctive genetic association profiles. Further, the systematic characterization involves composite phenotype descriptions that allow for combined power of individual phenotype GWAS tests, yielding more significance for lower individual thresholds, permitting the exploration of SNPs that would otherwise show as false negatives

T2DM GWAS in the Lebanese population confirms the role of TCF7L2 and CDKAL1 in disease susceptibility

Genome-wide association studies (GWAS) of multiple populations with distinctive genetic and lifestyle backgrounds are crucial to the understanding of Type 2 Diabetes Mellitus (T2DM) pathophysiology. We report a GWAS on the genetic basis of T2DM in a 3,286 Lebanese participants. More than 5,000,000 SNPs were directly genotyped or imputed using the 1000 Genomes Project reference panels. We identify genome-wide significant variants in two loci CDKAL1 and TCF7L2, independent of sex, age and BMI, with leading variants rs7766070 (OR = 1.39, P = 4.77 × 10(−9)) and rs34872471 (OR = 1.35, P = 1.01 × 10(−8)) respectively. The current study is the first GWAS to find genomic regions implicated in T2DM in the Lebanese population. The results support a central role of CDKAL1 and TCF7L2 in T2DM susceptibility in Southwest Asian populations and provide a plausible component for understanding molecular mechanisms involved in the disease

Multivariate epidemiologic analysis of type 2 diabetes mellitus risks in the Lebanese population

Background: The burden of diabetes in Lebanon requires well-targeted interventions for screening type 2 diabetes mellitus (T2DM) and prediabetes and prevention of risk factors. Newly recruited 998 Lebanese individuals, in addition to 7,292 already available, were studied to investigate the prevalence of diabetes, prediabetes and their associated risk factors. Methods: Participants had fasting blood sugar and glycohemoglobin tests in addition to a lipid profile. Clinical and demographic information were obtained from a detailed questionnaire. The relationship between T2DM, its risk factors, and its complications were tested. Comparisons of these risk factors among diabetics, healthy, and coronary artery disease (CAD) patients were performed. Results: The prevalence of T2DM significantly increased with increasing BMI (p < 0.0001). Exercise activity level negatively correlated with the disease (p = 0.002), whereas the prevalence of T2DM (p < 0.0001) and CAD family history (p = 0.006) positively correlated with the affection status. The mean levels of triglycerides and LDL-C were significantly higher in diabetics (1.87; 1.35) compared to individuals with prediabetes (1.63; 1.26) and unaffected controls (1.49; 1.19). People with T2DM showed a significant decrease in HDL-C levels. A strong correlation of overall hyperlipidemia with the diabetes affection status was shown (p < 0.0001). Other comorbid factors such as hypertension (p < 0.0001) and self-reported obesity (p < 0.0001) were highly associated with T2DM and prediabetes. Reproductive health of women showed a strong correlation between giving birth to a baby with a high weight and the occurrence of T2DM and prediabetes later in life (p < 0.0001). Retinopathy and peripheral neuropathy were significantly correlated with diabetes and prediabetes (p < 0.0001). Conclusions: The present study shows an alarming prevalence of diabetes and prediabetes in the studied subgroups representative of the Lebanese population. Electronic supplementary material The online version of this article (doi:10.1186/1758-5996-6-89) contains supplementary material, which is available to authorized users

Y-Chromosome and mtDNA Genetics Reveal Significant Contrasts in Affinities of Modern Middle Eastern Populations with European and African Populations

The Middle East was a funnel of human expansion out of Africa, a staging area for the Neolithic Agricultural Revolution, and the home to some of the earliest world empires. Post LGM expansions into the region and subsequent population movements created a striking genetic mosaic with distinct sex-based genetic differentiation. While prior studies have examined the mtDNA and Y-chromosome contrast in focal populations in the Middle East, none have undertaken a broad-spectrum survey including North and sub-Saharan Africa, Europe, and Middle Eastern populations. In this study 5,174 mtDNA and 4,658 Y-chromosome samples were investigated using PCA, MDS, mean-linkage clustering, AMOVA, and Fisher exact tests of FST's, RST's, and haplogroup frequencies. Geographic differentiation in affinities of Middle Eastern populations with Africa and Europe showed distinct contrasts between mtDNA and Y-chromosome data. Specifically, Lebanon's mtDNA shows a very strong association to Europe, while Yemen shows very strong affinity with Egypt and North and East Africa. Previous Y-chromosome results showed a Levantine coastal-inland contrast marked by J1 and J2, and a very strong North African component was evident throughout the Middle East. Neither of these patterns were observed in the mtDNA. While J2 has penetrated into Europe, the pattern of Y-chromosome diversity in Lebanon does not show the widespread affinities with Europe indicated by the mtDNA data. Lastly, while each population shows evidence of connections with expansions that now define the Middle East, Africa, and Europe, many of the populations in the Middle East show distinctive mtDNA and Y-haplogroup characteristics that indicate long standing settlement with relatively little impact from and movement into other populations

From cheek swabs to consensus sequences : an A to Z protocol for high-throughput DNA sequencing of complete human mitochondrial genomes

Background: Next-generation DNA sequencing (NGS) technologies have made huge impacts in many fields of biological research, but especially in evolutionary biology. One area where NGS has shown potential is for high-throughput sequencing of complete mtDNA genomes (of humans and other animals). Despite the increasing use of NGS technologies and a better appreciation of their importance in answering biological questions, there remain significant obstacles to the successful implementation of NGS-based projects, especially for new users. Results: Here we present an ‘A to Z’ protocol for obtaining complete human mitochondrial (mtDNA) genomes – from DNA extraction to consensus sequence. Although designed for use on humans, this protocol could also be used to sequence small, organellar genomes from other species, and also nuclear loci. This protocol includes DNA extraction, PCR amplification, fragmentation of PCR products, barcoding of fragments, sequencing using the 454 GS FLX platform, and a complete bioinformatics pipeline (primer removal, reference-based mapping, output of coverage plots and SNP calling). Conclusions: All steps in this protocol are designed to be straightforward to implement, especially for researchers who are undertaking next-generation sequencing for the first time. The molecular steps are scalable to large numbers (hundreds) of individuals and all steps post-DNA extraction can be carried out in 96-well plate format. Also, the protocol has been assembled so that individual ‘modules’ can be swapped out to suit available resources

Genome-Wide Association Study in a Lebanese Cohort Confirms PHACTR1 as a Major Determinant of Coronary Artery Stenosis

The manifestation of coronary artery disease (CAD) follows a well-choreographed series of events that includes damage of arterial endothelial cells and deposition of lipids in the sub-endothelial layers. Genome-wide association studies (GWAS) of multiple populations with distinctive genetic and lifestyle backgrounds are a crucial step in understanding global CAD pathophysiology. In this study, we report a GWAS on the genetic basis of arterial stenosis as measured by cardiac catheterization in a Lebanese population. The locus of the phosphatase and actin regulator 1 gene (PHACTR1) showed association with coronary stenosis in a discovery experiment with genome wide data in 1,949 individuals (rs9349379, OR = 1.37, p = 1.57×10−5). The association was replicated in an additional 2,547 individuals (OR = 1.31, p = 8.85×10−6), leading to genome-wide significant association in a combined analysis (OR = 1.34, p = 8.02×10−10). Results from this GWAS support a central role of PHACTR1 in CAD susceptibility irrespective of lifestyle and ethnic divergences. This association provides a plausible component for understanding molecular mechanisms involved in the formation of stenosis in cardiac vessels and a potential drug target against CAD

Genomic Ancestry of North Africans Supports Back-to-Africa Migrations

North African populations are distinct from sub-Saharan Africans based on cultural, linguistic, and phenotypic attributes; however, the time and the extent of genetic divergence between populations north and south of the Sahara remain poorly understood. Here, we interrogate the multilayered history of North Africa by characterizing the effect of hypothesized migrations from the Near East, Europe, and sub-Saharan Africa on current genetic diversity. We present dense, genome-wide SNP genotyping array data (730,000 sites) from seven North African populations, spanning from Egypt to Morocco, and one Spanish population. We identify a gradient of likely autochthonous Maghrebi ancestry that increases from east to west across northern Africa; this ancestry is likely derived from “back-to-Africa” gene flow more than 12,000 years ago (ya), prior to the Holocene. The indigenous North African ancestry is more frequent in populations with historical Berber ethnicity. In most North African populations we also see substantial shared ancestry with the Near East, and to a lesser extent sub-Saharan Africa and Europe. To estimate the time of migration from sub-Saharan populations into North Africa, we implement a maximum likelihood dating method based on the distribution of migrant tracts. In order to first identify migrant tracts, we assign local ancestry to haplotypes using a novel, principal component-based analysis of three ancestral populations. We estimate that a migration of western African origin into Morocco began about 40 generations ago (approximately 1,200 ya); a migration of individuals with Nilotic ancestry into Egypt occurred about 25 generations ago (approximately 750 ya). Our genomic data reveal an extraordinarily complex history of migrations, involving at least five ancestral populations, into North Africa

Synchronization in periodically driven and coupled stochastic systems-A discrete state approach

Wir untersuchen das Verhalten von stochastischen bistabilen und erregbaren Systemen auf der Basis einer Modellierung mit diskreten Zuständen. In Ergänzung zum bekannten Markovschen Zwei-Zustandsmodell bistabiler stochastischer Dynamik stellen wir ein nicht Markovsches Drei-Zustandsmodell für erregbare Systeme vor. Seine relative Einfachheit, verglichen mit stochastischen Modellen erregbarer Dynamik mit kontinuierlichem Phasenraum, ermöglicht eine teilweise analytische Auswertung in verschiedenen Zusammenhängen. Zunächst untersuchen wir den gemeinsamen Einfluß eines periodischen Treibens und Rauschens. Dieser wird entweder mit Hilfe spektraler Größen oder durch Synchronisation des Systems mit dem treibenden Signal charakterisiert. Wir leiten analytische Ausdrücke für die spektrale Leistungsverstärkung und das Signal-zu-Rauschen Verhältnis für periodisch getriebene Renewal-Prozesse her und wenden diese auf das diskrete Modell für erregbare Dynamik an. Stochastische Synchronization des Systems mit dem treibenden Signal wird auf der Basis der Diffusionseigenschaften der Übergangsereignisse zwischen den diskreten Zuständen untersucht. Wir leiten allgemeine Formeln her, um die mittlere Häufigkeit dieser Ereignisse sowie deren effektiven Diffusionskoeffizienten zu berechnen. Über die konkrete Anwendung auf die untersuchten diskreten Modelle hinaus stellen diese Ergebnisse ein neues Werkzeug für die Untersuchung periodischer Renewal-Prozesse dar. Schließlich betrachten wir noch das Verhalten global gekoppelter bistabiler und erregbarer Systeme. Im Gegensatz zu bistabilen System können erregbare Systeme synchronisiert werden und zeigen kohärente Oszillationen. Alle Untersuchungen des nicht Markovschen Drei-Zustandsmodells werden mit dem prototypischen Modell für erregbare Dynamik, dem FitzHugh-Nagumo System, verglichen und zeigen eine gute Übereinstimmung.We investigate the behavior of stochastic bistable and excitable dynamics based on a discrete state modeling. In addition to the well known Markovian two state model for bistable dynamics we introduce a non Markovian three state model for excitable systems. Its relative simplicity compared to stochastic models of excitable dynamics with continuous phase space allows to obtain analytical results in different contexts. First, we study the joint influence of periodic signals and noise, both based on a characterization in terms of spectral quantities and in terms of synchronization with the periodic driving. We present expressions for the spectral power amplification and signal to noise ratio for renewal processes driven by periodic signals and apply these results to the discrete model for excitable systems. Stochastic synchronization of the system to the driving signal is investigated based on diffusion properties of the transition events between the discrete states. We derive general results for the mean frequency and effective diffusion coefficient which, beyond the application to the discrete models considered in this work, provide a new tool in the study of periodically driven renewal processes. Finally the behavior of globally coupled excitable and bistable units is investigated based on the discrete state description. In contrast to the bistable systems, the excitable system exhibits synchronization and thus coherent oscillations. All investigations of the non Markovian three state model are compared with the prototypical continuous model for excitable dynamics, the FitzHugh-Nagumo system, revealing a good agreement between both models

Genetically determined height and coronary artery disease.

BACKGROUND: The nature and underlying mechanisms of an inverse association between adult height and the risk of coronary artery disease (CAD) are unclear. METHODS: We used a genetic approach to investigate the association between height and CAD, using 180 height-associated genetic variants. We tested the association between a change in genetically determined height of 1 SD (6.5 cm) with the risk of CAD in 65,066 cases and 128,383 controls. Using individual-level genotype data from 18,249 persons, we also examined the risk of CAD associated with the presence of various numbers of height-associated alleles. To identify putative mechanisms, we analyzed whether genetically determined height was associated with known cardiovascular risk factors and performed a pathway analysis of the height-associated genes. RESULTS: We observed a relative increase of 13.5% (95% confidence interval [CI], 5.4 to 22.1; P<0.001) in the risk of CAD per 1-SD decrease in genetically determined height. There was a graded relationship between the presence of an increased number of height-raising variants and a reduced risk of CAD (odds ratio for height quartile 4 versus quartile 1, 0.74; 95% CI, 0.68 to 0.84; P<0.001). Of the 12 risk factors that we studied, we observed significant associations only with levels of low-density lipoprotein cholesterol and triglycerides (accounting for approximately 30% of the association). We identified several overlapping pathways involving genes associated with both development and atherosclerosis. CONCLUSIONS: There is a primary association between a genetically determined shorter height and an increased risk of CAD, a link that is partly explained by the association between shorter height and an adverse lipid profile. Shared biologic processes that determine achieved height and the development of atherosclerosis may explain some of the association. (Funded by the British Heart Foundation and others.).Supported by the British Heart Foundation, the United Kingdom National Institute for Health Research, the European Union project CVgenes@target, and a grant from the Leducq Foundation.This is the final published version. It first appeared at http://www.nejm.org/doi/full/10.1056/NEJMoa1404881

Large Scale Association Analysis Identifies Three Susceptibility Loci for Coronary Artery Disease

Genome wide association studies (GWAS) and their replications that have associated DNA variants with myocardial infarction (MI) and/or coronary artery disease (CAD) are predominantly based on populations of European or Eastern Asian descent. Replication of the most significantly associated polymorphisms in multiple populations with distinctive genetic backgrounds and lifestyles is crucial to the understanding of the pathophysiology of a multifactorial disease like CAD. We have used our Lebanese cohort to perform a replication study of nine previously identified CAD/MI susceptibility loci (LTA, CDKN2A-CDKN2B, CELSR2-PSRC1-SORT1, CXCL12, MTHFD1L, WDR12, PCSK9, SH2B3, and SLC22A3), and 88 genes in related phenotypes. The study was conducted on 2,002 patients with detailed demographic, clinical characteristics, and cardiac catheterization results. One marker, rs6922269, in MTHFD1L was significantly protective against MI (OR = 0.68, p = 0.0035), while the variant rs4977574 in CDKN2A-CDKN2B was significantly associated with MI (OR = 1.33, p = 0.0086). Associations were detected after adjustment for family history of CAD, gender, hypertension, hyperlipidemia, diabetes, and smoking. The parallel study of 88 previously published genes in related phenotypes encompassed 20,225 markers, three quarters of which with imputed genotypes The study was based on our genome-wide genotype data set, with imputation across the whole genome to HapMap II release 22 using HapMap CEU population as a reference. Analysis was conducted on both the genotyped and imputed variants in the 88 regions covering selected genes. This approach replicated HNRNPA3P1-CXCL12 association with CAD and identified new significant associations of CDKAL1, ST6GAL1, and PTPRD with CAD. Our study provides evidence for the importance of the multifactorial aspect of CAD/MI and describes genes predisposing to their etiology