Indigenous Arabs are descendants of the earliest split from ancient Eurasian populations.

An open question in the history of human migration is the identity of the earliest Eurasian populations that have left contemporary descendants. The Arabian Peninsula was the initial site of the out-of-Africa migrations that occurred between 125,000 and 60,000 yr ago, leading to the hypothesis that the first Eurasian populations were established on the Peninsula and that contemporary indigenous Arabs are direct descendants of these ancient peoples. To assess this hypothesis, we sequenced the entire genomes of 104 unrelated natives of the Arabian Peninsula at high coverage, including 56 of indigenous Arab ancestry. The indigenous Arab genomes defined a cluster distinct from other ancestral groups, and these genomes showed clear hallmarks of an ancient out-of-Africa bottleneck. Similar to other Middle Eastern populations, the indigenous Arabs had higher levels of Neanderthal admixture compared to Africans but had lower levels than Europeans and Asians. These levels of Neanderthal admixture are consistent with an early divergence of Arab ancestors after the out-of-Africa bottleneck but before the major Neanderthal admixture events in Europe and other regions of Eurasia. When compared to worldwide populations sampled in the 1000 Genomes Project, although the indigenous Arabs had a signal of admixture with Europeans, they clustered in a basal, outgroup position to all 1000 Genomes non-Africans when considering pairwise similarity across the entire genome. These results place indigenous Arabs as the most distant relatives of all other contemporary non-Africans and identify these people as direct descendants of the first Eurasian populations established by the out-of-Africa migrations.

Point-of-care whole-exome sequencing of idiopathic male infertility.

PURPOSE: Nonobstructive azoospermia (NOA) affects 1% of the male population; however, despite state-of-the-art clinical assessment, for most patients the cause is unknown. We capitalized on an analysis of multiplex families in the Middle East to identify highly penetrant genetic causes. METHODS: We used whole-exome sequencing (WES) in 8 consanguineous families and combined newly discovered genes with previously reported ones to create a NOA gene panel, which was used to identify additional variants in 75 unrelated idiopathic NOA subjects and 74 fertile controls. RESULTS: In five of eight families, we identified rare deleterious recessive variants in CCDC155, NANOS2, SPO11, TEX14, and WNK3 segregating with disease. These genes, which are novel to human NOA, have remarkable testis-specific expression, and murine functional evidence supports roles for them in spermatogenesis. Among 75 unrelated NOA subjects, we identified 4 (~5.3%) with additional recessive variants in these newly discovered genes and 6 with deleterious variants in previously reported NOA genes, yielding an overall genetic etiology for 13.3% subjects versus 0 fertile controls (p = 0.001). CONCLUSION: NOA affects millions of men, many of whom remain idiopathic despite extensive laboratory evaluation. The genetic etiology for a substantial fraction of these patients (>50% familial and >10% sporadic) may be discovered by WES at the point of care.

Reconstructing Native American migrations from whole-genome and whole-exome data.

There is great scientific and popular interest in understanding the genetic history of populations in the Americas. We wish to understand when different regions of the continent were inhabited, where settlers came from, and how current inhabitants relate genetically to earlier populations. Recent studies unraveled parts of the genetic history of the continent using genotyping arrays and uniparental markers. The 1000 Genomes Project provides a unique opportunity for improving our understanding of population genetic history by providing over a hundred sequenced low coverage genomes and exomes from Colombian (CLM), Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore the genomic contributions of African, European, and especially Native American ancestry to these populations. Estimated Native American ancestry is 48% in MXL, 25% in CLM, and 13% in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern American ancestry of the Taíno people of the Caribbean. We present new methods to estimate the allele frequencies in the Native American fraction of the populations, and model their distribution using a demographic model for three ancestral Native American populations. These ancestral populations likely split in close succession: the most likely scenario, based on a peopling of the Americas 16 thousand years ago (kya), supports that the MXL Ancestors split 12.2kya, with a subsequent split of the ancestors to CLM and PUR 11.7kya. The model also features effective populations of 62,000 in Mexico, 8,700 in Colombia, and 1,900 in Puerto Rico. Modeling Identity-by-descent (IBD) and ancestry tract length, we show that post-contact populations also differ markedly in their effective sizes and migration patterns, with Puerto Rico showing the smallest effective size and the earlier migration from Europe. Finally, we compare IBD and ancestry assignments to find evidence for relatedness among European founders to the three populations.

Copy number variations in the genome of the Qatari population.

BACKGROUND: The populations of the Arabian Peninsula remain the least represented in public genetic databases, both in terms of single nucleotide variants and of larger genomic mutations. We present the first high-resolution copy number variation (CNV) map for a Gulf Arab population, using a hybrid approach that integrates array genotyping intensity data and next-generation sequencing reads to call CNVs in the Qatari population. METHODS: CNVs were detected in 97 unrelated Qatari individuals by running two calling algorithms on each of two primary datasets: high-resolution genotyping (Illumina Omni 2.5M) and high depth whole-genome sequencing (Illumina PE 100bp). The four call-sets were integrated to identify high confidence CNV regions, which were subsequently annotated for putative functional effect and compared to public databases of CNVs in other populations. The availability of genome sequence was leveraged to identify tagging SNPs in high LD with common deletions in this population, enabling their imputation from genotyping experiments in the future. RESULTS: Genotyping intensities and genome sequencing data from 97 Qataris were analyzed with four different algorithms and integrated to discover 16,660 high confidence CNV regions (CNVRs) in the total population, affecting ~28 Mb in the median Qatari genome. Up to 40% of all CNVs affected genes, including novel CNVs affecting Mendelian disease genes, segregating at different frequencies in the 3 major Qatari subpopulations, including those with Bedouin, Persian/South Asian, and African ancestry. Consistent with high consanguinity levels in the Bedouin subpopulation, we found an increased burden for homozygous deletions in this group. In comparison to known CNVs in the comprehensive Database of Genomic Variants, we found that 5% of all CNVRs in Qataris were completely novel, with an enrichment of CNVs affecting several known chromosomal disorder loci and genes known to regulate sugar metabolism and type 2 diabetes in the Qatari cohort. Finally, we leveraged the availability of genome sequence to find suitable tagging SNPs for common deletions in this population. CONCLUSION: We combine four independently generated datasets from 97 individuals to study CNVs for the first time at high-resolution in a Gulf Arab population.

Exome sequencing identifies potential risk variants for Mendelian disorders at high prevalence in Qatar.

Exome sequencing of families of related individuals has been highly successful in identifying genetic polymorphisms responsible for Mendelian disorders. Here, we demonstrate the value of the reverse approach, where we use exome sequencing of a sample of unrelated individuals to analyze allele frequencies of known causal mutations for Mendelian diseases. We sequenced the exomes of 100 individuals representing the three major genetic subgroups of the Qatari population (Q1 Bedouin, Q2 Persian-South Asian, Q3 African) and identified 37 variants in 33 genes with effects on 36 clinically significant Mendelian diseases. These include variants not present in 1000 Genomes and variants at high frequency when compared with 1000 Genomes populations. Several of these Mendelian variants were only segregating in one Qatari subpopulation, where the observed subpopulation specificity trends were confirmed in an independent population of 386 Qataris. Premarital genetic screening in Qatar tests for only four out of the 37, such that this study provides a set of Mendelian disease variants with potential impact on the epidemiological profile of the population that could be incorporated into the testing program if further experimental and clinical characterization confirms high penetrance.

Cysteine-Altering NOTCH3 Variants Are Associated with an Increased Risk of Autoimmune Diseases.

Autoimmune conditions have been reported among patients with cysteine-altering NOTCH3 variants and CADASIL. This study aimed to investigate the occurrence of autoimmune illnesses and markers of inflammation in such populations. Cases were identified who had a NOTCH3 cysteine-altering variant from the Geisinger MyCode® Community Health Initiative (MyCode®). We further performed external validation using the UK Biobank cohort. A cohort of 121 individuals with a NOTCH3 cysteine-altering variant from MyCode® was compared to a control group with no non-synonymous variation in NOTCH3 (n = 184). Medical records were evaluated for inflammatory markers and autoimmune conditions, which were grouped by the organ systems involved. A similar analysis was conducted using data from the UK Biobank (n~450,000). An overall increase in inflammatory markers among participants with a NOTCH3 cysteine-altering variant was observed when compared to an age- and sex-matched MyCode® control group (out of participants with laboratory testing: 50.9% versus 26.7%; p = 0.0047; out of total participants: 23.1% versus 10.9%; p = 0.004). Analysis of UK Biobank data indicated any autoimmune diagnosis (1.63 [1.14, 2.09], p= 2.665 × 10-3) and multiple sclerosis (3.42 [1.67, 6.02], p = 9.681 × 10-4) are associated with a NOTCH3 cysteine-altering variant in any domain. Our findings suggest a possible association between NOTCH3 cysteine-altering variants and autoimmune conditions.

RNA-Seq quantification of the human small airway epithelium transcriptome.

BACKGROUND: The small airway epithelium (SAE), the cell population that covers the human airway surface from the 6th generation of airway branching to the alveoli, is the major site of lung disease caused by smoking. The focus of this study is to provide quantitative assessment of the SAE transcriptome in the resting state and in response to chronic cigarette smoking using massive parallel mRNA sequencing (RNA-Seq). RESULTS: The data demonstrate that 48% of SAE expressed genes are ubiquitous, shared with many tissues, with 52% enriched in this cell population. The most highly expressed gene, SCGB1A1, is characteristic of Clara cells, the cell type unique to the human SAE. Among other genes expressed by the SAE are those related to Clara cell differentiation, secretory mucosal defense, and mucociliary differentiation. The high sensitivity of RNA-Seq permitted quantification of gene expression related to infrequent cell populations such as neuroendocrine cells and epithelial stem/progenitor cells. Quantification of the absolute smoking-induced changes in SAE gene expression revealed that, compared to ubiquitous genes, more SAE-enriched genes responded to smoking with up-regulation, and those with the highest basal expression levels showed most dramatic changes. Smoking had no effect on SAE gene splicing, but was associated with a shift in molecular pattern from Clara cell-associated towards the mucus-secreting cell differentiation pathway with multiple features of cancer-associated molecular phenotype. CONCLUSIONS: These observations provide insights into the unique biology of human SAE by providing quantitative assessment of the global transcriptome under physiological conditions and in response to the stress of chronic cigarette smoking.

Inferring genome-wide patterns of admixture in Qataris using fifty-five ancestral populations.

BACKGROUND: Populations of the Arabian Peninsula have a complex genetic structure that reflects waves of migrations including the earliest human migrations from Africa and eastern Asia, migrations along ancient civilization trading routes and colonization history of recent centuries. RESULTS: Here, we present a study of genome-wide admixture in this region, using 156 genotyped individuals from Qatar, a country located at the crossroads of these migration patterns. Since haplotypes of these individuals could have originated from many different populations across the world, we have developed a machine learning method "SupportMix" to infer loci-specific genomic ancestry when simultaneously analyzing many possible ancestral populations. Simulations show that SupportMix is not only more accurate than other popular admixture discovery tools but is the first admixture inference method that can efficiently scale for simultaneous analysis of 50-100 putative ancestral populations while being independent of prior demographic information. CONCLUSIONS: By simultaneously using the 55 world populations from the Human Genome Diversity Panel, SupportMix was able to extract the fine-scale ancestry of the Qatar population, providing many new observations concerning the ancestry of the region. For example, as well as recapitulating the three major sub-populations in Qatar, composed of mainly Arabic, Persian, and African ancestry, SupportMix additionally identifies the specific ancestry of the Persian group to populations sampled in Greater Persia rather than from China and the ancestry of the African group to sub-Saharan origin and not Southern African Bantu origin as previously thought.

The QChip1 knowledgebase and microarray for precision medicine in Qatar.

Risk genes for Mendelian (single-gene) disorders (SGDs) are consistent across populations, but pathogenic risk variants that cause SGDs are typically population-private. The goal was to develop "QChip1," an inexpensive genotyping microarray to comprehensively screen newborns, couples, and patients for SGD risk variants in Qatar, a small nation on the Arabian Peninsula with a high degree of consanguinity. Over 108 variants in 8445 Qatari were identified for inclusion in a genotyping array containing 165,695 probes for 83,542 known and potentially pathogenic variants in 3438 SGDs. QChip1 had a concordance with whole-genome sequencing of 99.1%. Testing of QChip1 with 2707 Qatari genomes identified 32,674 risk variants, an average of 134 pathogenic alleles per Qatari genome. The most common pathogenic variants were those causing homocystinuria (1.12% risk allele frequency), and Stargardt disease (2.07%). The majority (85%) of Qatari SGD pathogenic variants were not present in Western populations such as European American, South Asian American, and African American in New York City and European and Afro-Caribbean in Puerto Rico; and only 50% were observed in a broad collection of data across the Greater Middle East including Kuwait, Iran, and United Arab Emirates. This study demonstrates the feasibility of developing accurate screening tools to identify SGD risk variants in understudied populations, and the need for ancestry-specific SGD screening tools.

Pro-hormone secretogranin II regulates dense core secretory granule biogenesis in catecholaminergic cells.

Processes underlying the formation of dense core secretory granules (DCGs) of neuroendocrine cells are poorly understood. Here, we present evidence that DCG biogenesis is dependent on the secretory protein secretogranin (Sg) II, a member of the granin family of pro-hormone cargo of DCGs in neuroendocrine cells. Depletion of SgII expression in PC12 cells leads to a decrease in both the number and size of DCGs and impairs DCG trafficking of other regulated hormones. Expression of SgII fusion proteins in a secretory-deficient PC12 variant rescues a regulated secretory pathway. SgII-containing dense core vesicles share morphological and physical properties with bona fide DCGs, are competent for regulated exocytosis, and maintain an acidic luminal pH through the V-type H(+)-translocating ATPase. The granulogenic activity of SgII requires a pH gradient along this secretory pathway. We conclude that SgII is a critical factor for the regulation of DCG biogenesis in neuroendocrine cells, mediating the formation of functional DCGs via its pH-dependent aggregation at the trans-Golgi network.

Mitochondrial haplogroup J associated with higher risk of obesity in the Qatari population.

Obesity, a major risk factor for metabolic disorders, is highly prevalent in Qatari population. Maternal transmission of obesity traits can be significant; for example, X haplogroup is known to be associated with lower BMI and body fat mass in Northern Europeans and T haplogroup which is a sister haplogroup of J is known to be associated with obesity in Caucasian subjects from Austria and Southern Italy. We aimed to delineate the mitochondrial haplogroups and variants associated with obesity in Qatari population. Mitochondrial genomes of 864 Qatari individuals were extracted from whole exome sequencing data with an average coverage of 77X. We distributed the participants into 2 sub-cohorts: obese (BMI ≥ 30) and non-obese (BMI < 30); the mean value of BMI from these two groups were 36.5 ± 5.7 and 26.5 ± 2.6, respectively. Mitochondrial haplogroup profiling followed by uni- and multivariant association tests adjusted for covariates were performed. Qatari individuals with mitochondrial haplogroup J had an increased (twofold) risk of obesity (odds ratio [OR] 1.925; 95% CI 1.234-3.002; P = 0.0038; the Bonferroni adjusted P value threshold is 0.0041), whereas the individuals with haplogroup X were at low risk of obesity (OR 0.387; 95% CI 0.175-0.857; P = 0.019). Further, a set of 38 mitochondrial variants were found to be associated (at P ≤ 0.05) with obesity in models adjusted for age, sex and haplogroup.

Large-Scale Identification of Clonal Hematopoiesis and Mutations Recurrent in Blood Cancers.

Clonal hematopoiesis of indeterminate potential (CHIP) is characterized by detectable hematopoietic-associated gene mutations in a person without evidence of hematologic malignancy. We sought to identify additional cancer-presenting mutations useable for CHIP detection by performing a data mining analysis of 48 somatic mutation studies reporting mutations at diagnoses of 7,430 adult and pediatric patients with hematologic malignancies. Following extraction of 20,141 protein-altering mutations, we identified 434 significantly recurrent mutation hotspots, 364 of which occurred at loci confidently assessable for CHIP. We then performed an additional large-scale analysis of whole exome sequencing data from 4,538 persons belonging to three non-cancer cohorts for clonal mutations. We found the combined cohort prevalence of CHIP with mutations identical to those reported at blood cancer mutation hotspots to be 1.8%, and that some of these CHIP mutations occurred in children. Our findings may help to improve CHIP detection and pre-cancer surveillance for both children and adults.

Genome diversity in Ukraine.

BACKGROUND: The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage. RESULTS: The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population. CONCLUSIONS: Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

Discovery of common human genetic variants of GTP cyclohydrolase 1 (GCH1) governing nitric oxide, autonomic activity, and cardiovascular risk.

GTP cyclohydrolase 1 (GCH1) is rate limiting in the provision of the cofactor tetrahydrobiopterin for biosynthesis of catecholamines and NO. We asked whether common genetic variation at GCH1 alters transmitter synthesis and predisposes to disease. Here we undertook a systematic search for polymorphisms in GCH1, then tested variants' contributions to NO and catecholamine release as well as autonomic function in twin pairs. Renal NO and neopterin excretions were significantly heritable, as were baroreceptor coupling (heart rate response to BP fluctuation) and pulse interval (1/heart rate). Common GCH1 variant C+243T in the 3'-untranslated region (3'-UTRs) predicted NO excretion, as well as autonomic traits: baroreceptor coupling, maximum pulse interval, and pulse interval variability, though not catecholamine secretion. In individuals with the most extreme BP values in the population, C+243T affected both diastolic and systolic BP, principally in females. In functional studies, C+243T decreased reporter expression in transfected 3'-UTRs plasmids. We conclude that human NO secretion traits are heritable, displaying joint genetic determination with autonomic activity by functional polymorphism at GCH1. Our results document novel pathophysiological links between a key biosynthetic locus and NO metabolism and suggest new strategies for approaching the mechanism, diagnosis, and treatment of risk predictors for cardiovascular diseases such as hypertension.

Conserved regulatory motifs at phenylethanolamine N-methyltransferase (PNMT) are disrupted by common functional genetic variation: an integrated computational/experimental approach.

The adrenomedullary hormone epinephrine transduces environmental stressors into cardiovascular events (tachycardia and hypertension). Although the epinephrine biosynthetic enzyme PNMT genetic locus displays both linkage and association to such traits, genetic variation underlying these quantitative phenotypes is not established. Using an integrated suite of computational and experimental approaches, we elucidate a functional mechanism for common (minor allele frequencies > 30%) genetic variants at PNMT. Transcription factor binding motif prediction on mammalian PNMT promoter alignments identified two variant regulatory motifs, SP1 and EGR1, disrupted by G-367A (rs3764351), and SOX17 motif created by G-161A (rs876493). Electrophoretic mobility shifts of approximately 30-bp oligonucleotides containing ancestral versus variant alleles validated the computational hypothesis. Queried against chromaffin cell nuclear protein extracts, only the G-367 and -161A alleles shifted. Specific antibodies applied in electrophoretic gel shift experiments confirmed binding of SP1 and EGR1 to G-367 and SOX17 to -161A. The in vitro allele-specific binding was verified in cella through promoter reporter assays: lower activity for -367A haplotypes cotransfected by SP1 (p = 0.002) and EGR1 (p = 0.034); and enhanced inhibition of -161A haplotypes (p = 0.0003) cotransfected with SP1 + SOX17. Finally, we probed cis/trans regulation with endogenous factors by chromatin immunoprecipitation using SP1/EGR1/SOX17 antibodies. We describe the systematic application of complementary computational and experimental techniques to detect and document functional genetic variation in a trait-associated regulatory region. The results provide insight into cis and trans transcriptional mechanisms whereby common variation at PNMT can give rise to quantitative changes in human physiological and disease traits. Thus, PNMT variants in cis may interact with nuclear factors in trans to govern adrenergic activity.

Human tyrosine hydroxylase natural allelic variation: influence on autonomic function and hypertension.

The catecholamine biosynthetic pathway consists of several enzymatic steps in series, beginning with the amino acids phenylalanine and tyrosine, and eventuating in the catecholamines norepinephrine (noradrenaline) and epinephrine (adrenaline). Since the enzyme tyrosine hydroxylase (TH; tyrosine 3-mono-oxygenase; EC 1.14.16.2; chromosome 11p15.5) is generally considered to be rate-limiting in this pathway, probed as to whether common genetic variation at the TH gene occurred, and whether such variants contributed to inter-individual alterations in autonomic function, either biochemical or physiological. We began with sequencing a tetranucleotide (TCAT) repeat in the first intron, and found that the two most common versions, (TCAT)(6) and (TCAT)(10i), predicted heritable autonomic traits in twin pairs. We then conducted systematic polymorphism discovery across the ~8 kbp locus, and discovered numerous variants, principally non-coding. The proximal promoter block contained four common variants, and its haplotypes and SNPs (especially C-824T, rs10770141) predicted catecholamine secretion, environmental stress-induced BP increments, and hypertension. Finally, we found that two of the common promoter variants, C-824T (rs10770141) and A-581G (rs10770140), were functional in that they differentially affected transcriptional activity of the isolated promoter, disrupted recognition motifs for specific transcription factor binding, altered the promoter responses to the co-transfected (exogenous) factors, and bound the endogenous factors in the chromatin fraction of the nucleus. We concluded that common variation in the proximal TH promoter is functional, giving rise to changes in autonomic function and consequently cardiovascular risk.

Genomes of Three Closely Related Caribbean Amazons Provide Insight for Species History and Conservation.

Islands have been used as model systems for studies of speciation and extinction since Darwin published his observations about finches found on the Galapagos. Amazon parrots inhabiting the Greater Antillean Islands represent a fascinating model of species diversification. Unfortunately, many of these birds are threatened as a result of human activity and some, like the Puerto Rican parrot, are now critically endangered. In this study we used a combination of de novo and reference-assisted assembly methods, integrating it with information obtained from related genomes to perform genome reconstruction of three amazon species. First, we used whole genome sequencing data to generate a new de novo genome assembly for the Puerto Rican parrot (Amazona vittata). We then improved the obtained assembly using transcriptome data from Amazona ventralis and used the resulting sequences as a reference to assemble the genomes Hispaniolan (A. ventralis) and Cuban (Amazona leucocephala) parrots. Finally, we, annotated genes and repetitive elements, estimated genome sizes and current levels of heterozygosity, built models of demographic history and provided interpretation of our findings in the context of parrot evolution in the Caribbean.

Tyrosine hydroxylase, the rate-limiting enzyme in catecholamine biosynthesis: discovery of common human genetic variants governing transcription, autonomic activity, and blood pressure in vivo.

BACKGROUND: Tyrosine hydroxylase (TH) is the rate-limiting enzyme in catecholamine biosynthesis. Does common genetic variation at human TH alter autonomic activity and predispose to cardiovascular disease? We undertook systematic polymorphism discovery at the TH locus and then tested variants for contributions to sympathetic function and blood pressure. METHODS AND RESULTS: We resequenced 80 ethnically diverse individuals across the TH locus. One hundred seventy-two twin pairs were evaluated for sympathetic traits, including catecholamine production, reflex control of the circulation, and environmental (cold) stress responses. To evaluate hypertension, we genotyped subjects selected from the most extreme diastolic blood pressure percentiles in the population. Human TH promoter haplotype/reporter plasmids were transfected into chromaffin cells. Forty-nine single-nucleotide polymorphisms were discovered, but coding region polymorphism did not account for common phenotypic variation. A block of linkage disequilibrium spanned 4 common variants in the proximal promoter. Catecholamine secretory traits were significantly heritable (h2), as were stress-induced blood pressure changes. In the TH promoter, significant associations were found for urinary catecholamine excretion and for blood pressure response to stress. TH promoter haplotype 2 (TGGG) showed pleiotropy, increasing both norepinephrine excretion and blood pressure during stress. Coalescent simulations suggest that TH haplotype 2 likely arose approximately 380,000 years ago. In hypertension, 2 independent case-control studies (1266 subjects with 53% women and 927 subjects with 24% women) replicated the effect of C-824T in the determination of blood pressure. CONCLUSIONS: We conclude that human catecholamine secretory traits are heritable, displaying joint genetic determination (pleiotropy) with autonomic activity and finally with blood pressure in the population. Catecholamine secretion is influenced by genetic variation in the adrenergic pathway encoding catecholamine synthesis, especially at the classically rate-limiting step, TH. The results suggest novel pathophysiological links between a key adrenergic locus, catecholamine metabolism, and blood pressure and suggest new strategies to approach the mechanism, diagnosis, and treatment of systemic hypertension.

Whole-exome sequencing identifies common and rare variant metabolic QTLs in a Middle Eastern population.

Metabolomics-genome-wide association studies (mGWAS) have uncovered many metabolic quantitative trait loci (mQTLs) influencing human metabolic individuality, though predominantly in European cohorts. By combining whole-exome sequencing with a high-resolution metabolomics profiling for a highly consanguineous Middle Eastern population, we discover 21 common variant and 12 functional rare variant mQTLs, of which 45% are novel altogether. We fine-map 10 common variant mQTLs to new metabolite ratio associations, and 11 common variant mQTLs to putative protein-altering variants. This is the first work to report common and rare variant mQTLs linked to diseases and/or pharmacological targets in a consanguineous Arab cohort, with wide implications for precision medicine in the Middle East.

Exome sequencing-based identification of novel type 2 diabetes risk allele loci in the Qatari population.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and lifestyle factors. To date, the majority of genetic studies of T2D have been in populations of European and Asian descent. The focus of this study is on genetic variations underlying T2D in Qataris, a population with one of the highest incidences of T2D worldwide. RESULTS: Illumina HiSeq exome sequencing was performed on 864 Qatari subjects (574 T2D cases, 290 controls). Sequence kernel association test (SKAT) gene-based analysis identified an association for low frequency potentially deleterious variants in 6 genes. However, these findings were not replicated by SKAT analysis in an independent cohort of 12,699 exomes, primarly due to the absence of low frequency potentially deleterious variants in 5 of the 6 genes. Interestingly one of the genes identified, catenin beta 1 (CTNNB1, ß-catenin), is the key effector of the Wnt pathway and interacts with the nuclear receptor transcription factor 7-like 2 (TCF7L2), variants which are the most strongly associated with risk of developing T2D worldwide. Single variant analysis did not identify any associated variants, suggesting the SKAT association signal was not driven by individual variants. None of the 6 associated genes were among 634 previously described T2D genes. CONCLUSIONS: The observation that genes not previously linked to T2D in prior studies of European and Asian populations are associated with T2D in Qatar provides new insights into the complexity of T2D pathogenesis and emphasizes the importance of understudied populations when assessing genetic variation in the pathogenesis of common disorders.