Genome-wide association of echocardiographic dimensions, brachial artery endothelial function and treadmill exercise responses in the Framingham Heart Study

Advancing drug development for atrial fibrillation by prioritising findings from human genetic association studies

BACKGROUND

Drug development for atrial fibrillation (AF) has failed to yield new approved compounds. We sought to identify and prioritise potential druggable targets with support from human genetics, by integrating the available evidence with bioinformatics sources relevant for AF drug development.

METHODS

Genetic hits for AF and related traits were identified through structured search of MEDLINE. Genes derived from each paper were cross-referenced with the OpenTargets platform for drug interactions. Confirmation/validation was demonstrated through structured searches and review of evidence on MEDLINE and ClinialTrials.gov for each drug and its association with AF.

FINDINGS

613 unique drugs were identified, with 21 already included in AF Guidelines. Cardiovascular drugs from classes not currently used for AF (e.g. ranolazine and carperitide) and anti-inflammatory drugs (e.g. dexamethasone and mehylprednisolone) had evidence of potential benefit. Further targets were considered druggable but remain open for drug development.

INTERPRETATION

Our systematic approach, combining evidence from different bioinformatics platforms, identified drug repurposing opportunities and druggable targets for AF.

FUNDING

KK is supported by Barts Charity grant G-002089 and is mentored on the AFGen 2023-24 Fellowship funded by the AFGen NIH/NHLBI grant R01HL092577. RP is supported by the UCL BHF Research Accelerator AA/18/6/34223 and NIHR grant NIHR129463. AFS is supported by the BHF grants PG/18/5033837, PG/22/10989 and UCL BHF Accelerator AA/18/6/34223 as well as the UK Research and Innovation (UKRI) under the UK government's Horizon Europe funding guarantee EP/Z000211/1 and by the UKRI-NIHR grant MR/V033867/1 for the Multimorbidity Mechanism and Therapeutics Research Collaboration. AF is supported by UCL BHF Accelerator AA/18/6/34223. CF is supported by UCL BHF Accelerator AA/18/6/34223.

Transcriptional Variabilities in Human hiPSC-derived Cardiomyocytes: All Genes Are Not Equal and Their Robustness May Foretell Donor's Disease Susceptibility

Human induced pluripotent stem cells (hiPSCs) are frequently used to study disease-associated variations. We characterized transcriptional variability from a hiPSC-derived cardiomyocyte (hiPSC-CM) study of left ventricular hypertrophy (LVH) using donor samples from the HyperGEN study. Multiple hiPSC-CM differentiations over reprogramming events (iPSC generation) across 7 donors were used to assess variabilities from reprogramming, differentiation, and donor LVH status. Variability arising from pathological alterations was assessed using a cardiac stimulant applied to the hiPSC-CMs to trigger hypertrophic responses. We found that for most genes (73.3%~85.5%), technical variability was smaller than biological variability. Further, we identified and characterized lists of "noise" genes showing greater technical variability and "signal" genes showing greater biological variability. Together, they support a "genetic robustness" hypothesis of disease-modeling whereby cellular response to relevant stimuli in hiPSC-derived somatic cells from diseased donors tends to show more transcriptional variability. Our findings suggest that hiPSC-CMs can provide a valid model for cardiac hypertrophy and distinguish between technical and disease-relevant transcriptional changes.

The rs2682826 Polymorphism of the NOS1 Gene Is Associated with the Degree of Disability of Erectile Dysfunction

Erectile dysfunction (ED) is a common male disorder, often associated with cardiovascular disease and ageing. The Sildenafil, a PDE5 inhibitor, can improve the erectile function by prolonging the nitric oxide (NO) downstream effect. NO is a molecule of pivotal importance in erection physiology and is mainly produced by neuronal nitric oxide synthase (nNOS) and endothelial NO synthase (eNOS). While it has been shown that eNOS and nNOS genetic polymorphisms could be associated with Sildenafil responsiveness in ED, no study so far has assessed whether nNOS polymorphisms and PDE5A polymorphism could be associated with increased risk to ED or with intensity of symptoms. A total of 119 ED patients and 114 controls were studied, with evaluation of the clinical disability by the International Index for Erectile Function instrument, plasma assessment of nitrite levels and genomic DNA analysis regarding the rs41279104 and rs2682826 polymorphisms of the NOS1 gene and the rs2389866, rs3733526 and rs13124532 polymorphisms of the PDE5A gene. We have found a significant association of the rs2682826 with lower IIEF scores in the clinical ED group. While this result should be confirmed in other populations, it may be helpful in establishing a genetic panel to better assess disease risk and prognosis on ED therapy.

Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems

Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.

Differentially Methylated DNA Regions and Left Ventricular Hypertrophy in African Americans: A HyperGEN Study

Left ventricular (LV) hypertrophy (LVH) is an independent risk factor for cardiovascular disease, and African Americans experience a disparate high risk of LVH. Genetic studies have identified potential candidate genes and variants related to the condition. Epigenetic modifications may continue to help unravel disease mechanisms. We used methylation and echocardiography data from 636 African Americans selected from the Hypertension Genetic Epidemiology Network (HyperGEN) to identify differentially methylated regions (DMRs) associated with LVH. DNA extracted from whole blood was assayed on Illumina Methyl450 arrays. We fit linear mixed models to examine associations between co-methylated regions and LV traits, and we then conducted single CpG analyses within significant DMRs. We identified associations between DMRs and ejection fraction (XKR6), LV internal diastolic dimension (TRAK1), LV mass index (GSE1, RPS15 A, PSMD7), and relative wall thickness (DNHD1). In single CpG analysis, CpG sites annotated to TRAK1 and DNHD1 were significant. These CpGs were not associated with LV traits in replication cohorts but the direction of effect for DNHD1 was consistent across cohorts. Of note, DNHD1, GSE1, and PSMD7 may contribute to cardiac structural function. Future studies should evaluate relationships between regional DNA methylation patterns and the development of LVH.

Whole genome sequencing reveals a complex introgression history and the basis of adaptation to subarctic climate in wild sheep

To predict species responses to anthropogenic disturbances and climate change, it is reasonable to use species with high sensitivity to such factors. Snow sheep (Ovis nivicola) could represent a good candidate for this; as the only large herbivore species adapted to the cold and alpine habitats of northeastern Siberia, it plays a crucial role in its ecosystem. Despite having an extensive geographical distribution among all ovine species, it is one of the least studied. In this study, we sequenced and analysed six genomes of snow sheep in combination with all other wild sheep species to infer key aspects of their evolutionary history and unveil the genetic basis of their adaptation to subarctic environments. Despite their large census population size, snow sheep genomes showed remarkably low heterozygosity, which could reflect the effect of isolation and historical bottlenecks that we inferred using the pairwise sequential Markovian coalescent and runs of homozygosity. F₄ -statistics indicated instances of introgression involving snow sheep with argali (Ovis ammon) and Dall (Ovis dalli) sheep, suggesting that these species might have been more widespread during the Pleistocene. Furthermore, the introgressed segments, which were identified using mainly minimum relative node depth, covered genes associated with immunity, adipogenesis and morphology-related traits, representing potential targets of adaptive introgression. Genes related to mitochondrial functions and thermogenesis associated with adipose tissue were identified to be under selection. Overall, our data suggest introgression as a mechanism facilitating adaptation in wild sheep species and provide insights into the genetic mechanisms underlying cold adaptation in snow sheep.

The genetics of human performance

Human physiology is likely to have been selected for endurance physical activity. However, modern humans have become largely sedentary, with physical activity becoming a leisure-time pursuit for most. Whereas inactivity is a strong risk factor for disease, regular physical activity reduces the risk of chronic disease and mortality. Although substantial epidemiological evidence supports the beneficial effects of exercise, comparatively little is known about the molecular mechanisms through which these effects operate. Genetic and genomic analyses have identified genetic variation associated with human performance and, together with recent proteomic, metabolomic and multi-omic analyses, are beginning to elucidate the molecular genetic mechanisms underlying the beneficial effects of physical activity on human health.

Genome-wide association study identifies 18 novel loci associated with left atrial volume and function

Aims 

Left atrial (LA) volume and function impose significant impact on cardiovascular pathogenesis if compromised. We aimed at investigating the genetic architecture of LA volume and function using cardiac magnetic resonance imaging data.

Methods and results 

We used the UK Biobank, which is a large prospective population study with available phenotypic and genetic data. On a subset of 35 658 European individuals, we performed genome-wide association studies on five volumetric and functional LA variables, generated using a machine learning algorithm. In total, we identified 18 novel genetic loci, mapped to genes with known roles in cardiomyopathy (e.g. MYO18B, TTN, DSP, ANKRD1) and arrhythmia (e.g. TTN, CASQ2, MYO18B, C9orf3). We observed high genetic correlation between LA volume and function and stroke, which was most pronounced for LA passive emptying fraction (r_g = 0.40, P = 4 × 10⁻⁶). To investigate whether the genetic risk of atrial fibrillation (AF) is associated with LA traits that precede overt AF, we produced a polygenetic risk score for AF. We found that polygenetic risk for AF is associated with increased LA volume and decreased LA function in participants without AF [LA_max 0.25 (mL/m²)/standard deviation (SD), 95% confidence interval (CI) (0.15; 0.36), P = 5.13 × 10⁻⁶; LA_min 0.21 (mL/m²)/SD, 95% CI (0.15; 0.28), P = 1.86 × 10⁻¹⁰; LA active emptying fraction −0.35%/SD, 95% CI (−0.43; −0.26), P = 3.14 × 10⁻¹⁴].

Conclusion 

We report on 18 genetic loci associated with LA volume and function and show evidence for several plausible candidate genes important for LA structure.

Tubulin-folding cofactor E deficiency promotes vascular dysfunction by increased endoplasmic reticulum stress

AIMS

Assessment of endothelial function in humans by measuring flow-mediated dilation (FMD) risk-stratifies individuals with established cardiovascular disease, whereas its predictive value is limited in primary prevention. We therefore aimed to establish and evaluate novel markers of FMD at the population level.

METHODS AND RESULTS

In order to identify novel targets that were negatively correlated with FMD and investigate their contribution to vascular function, we performed a genome-wide association study (GWAS) of 4175 participants of the population based Gutenberg Health Study. Subsequently, conditional knockout mouse models deleting the gene of interest were generated and characterized. GWAS analysis revealed that single-nucleotide polymorphisms (SNPs) in the tubulin-folding cofactor E (TBCE) gene were negatively correlated with endothelial function and TBCE expression. Vascular smooth muscle cell (VSMC)-targeted TBCE deficiency was associated with endothelial dysfunction, aortic wall hypertrophy, and endoplasmic reticulum (ER) stress-mediated VSMC hyperproliferation in mice, paralleled by calnexin up-regulation and exacerbated by the blood pressure hormone angiotensin II. Treating SMMHC-ERT2-Cre+/-TBCEfl/fl mice with the ER stress modulator tauroursodeoxycholic acid amplified Raptor/Beclin-1-dependent autophagy and reversed vascular dysfunction.

CONCLUSION

TBCE and tubulin homeostasis seem to be novel predictors of vascular function and offer a new drug target to ameliorate ER stress-dependent vascular dysfunction.

Genetics of canine myxomatous mitral valve disease

Myxomatous mitral valve disease (MMVD) is the most common heart disease and cause of cardiac death in domestic dogs. MMVD is characterised by slow progressive myxomatous degeneration from the tips of the mitral valves onwards with subsequent mitral valve regurgitation, and left atrial and ventricular dilatation. Although the disease usually has a long asymptomatic period, in dogs with severe disease, mortality is typically secondary to left-sided congestive heart failure. Although it is not uncommon for dogs to survive long enough in the asymptomatic period to die from unrelated causes; a proportion of dogs rapidly advance into congestive heart failure. Heightened prevalence in certain breeds, such as the Cavalier King Charles Spaniel, has indicated that MMVD is under a genetic influence. The genetic characterisation of the factors that underlie the difference in progression of disease is of strong interest to those concerned with dog longevity and welfare. Advanced genomic technologies have the potential to provide information that may impact treatment, prevalence, or severity of MMVD through the elucidation of pathogenic mechanisms and the detection of predisposing genetic loci of major effect. Here we describe briefly the clinical nature of the disorder and consider the physiological mechanisms that might impact its occurrence in the domestic dog. Using results from comparative genomics we suggest possible genetic approaches for identifying genetic risk factors within breeds. The Cavalier King Charles Spaniel breed represents a robust resource for uncovering the genetic basis of MMVD.

Genetic approaches toward understanding the individual variation in cardiac structure, function and responses to exercise training

Cardiovascular disease (CVD) accounts for approximately 30% of all deaths worldwide and its prevalence is constantly increasing despite advancements in medical treatments. Cardiac remodeling and dysfunction are independent risk factors for CVD. Recent studies have demonstrated that cardiac structure and function are genetically influenced, suggesting that understanding the genetic basis for cardiac structure and function could provide new insights into developing novel therapeutic targets for CVD. Regular exercise has long been considered a robust non-therapeutic method of treating or preventing CVD. However, recent studies also indicate that there is inter-individual variation in response to exercise. Nevertheless, the genetic basis for cardiac structure and function as well as their responses to exercise training have yet to be fully elucidated. Therefore, this review summarizes accumulated evidence supporting the genetic contribution to these traits, including findings from population-based studies and unbiased large genomic-scale studies in humans.

Heritability and Sex-Specific Genetic Effects of Self-Reported Physical Activity in a Brazilian Highly Admixed Population

INTRODUCTION

The engagement in sports or habitual physical activity (PA) has shown an extensive protective role against multiple diseases such as cancer, obesity, and many others. Additionally, PA has also a significant impact on life quality, since it aids with managing stress, preserving cognitive function and memory, and preventing fractures in the elderly.

OBJECTIVE

Considering there has been multiple evidence showing that genetic variation underpins variation of PA-related traits, we aimed to estimate the heritability (h2) of these phenotypes in a sample from the Brazilian population and assess whether males and females differ in relation to those estimates.

METHODS

2,027 participants from a highly admixed population from Baependi, MG, Brazil, had information regarding their PA and sedentary behavior (SB) phenotypes collected through a questionnaire (IPAQ-SF). After data cleaning and transformation procedures, we obtained four variables to be evaluated: total PA (TPA MET), walking time, (WK MET), moderate-vigorous PA (MVPA MET), and SB. A model selection procedure was performed using a single-step covariate inclusion approach. We tested for BMI, waist, hip and neck circumferences, smoking, and depression separately, and performed correlation tests among covariates. Linear mixed models, selection procedure, and the variance components approach to estimate h2 were implemented using SOLAR-Eclipse 8.3.1.

RESULTS

We obtained estimates of 0.221, 0.109, 0.226, and 0 for TPA MET, WK MET, MVPA MET, and SB, respectively. We found evidence for gene-sex interactions, with males having higher sex-specific heritabilities than females for TPA MET and MVPA MET. In addition, we found higher estimates of the genetic variance component in males than females for most phenotypes.

DISCUSSION/CONCLUSION

The heritability estimates presented in this work show a moderate heritable set of genetic factors affecting PA in a sample from the Brazilian population. The evaluation of the genetic variance component suggests segregating genetic factors in male individuals are more heterogeneous, which can explain why men globally tend to need to practice more intense PA than women to achieve similar health benefits. Hence, these findings have significant implications for the understanding of the genetic architecture of PA and might aid to promote health in the future.

Genetics and the heart rate response to exercise

The acute heart rate response to exercise, i.e., heart rate increase during and heart rate recovery after exercise, has often been associated with all-cause and cardiovascular mortality. The long-term response of heart rate to exercise results in favourable changes in chronotropic function, including decreased resting and submaximal heart rate as well as increased heart rate recovery. Both the acute and long-term heart rate response to exercise have been shown to be heritable. Advances in genetic analysis enable researchers to investigate this hereditary component to gain insights in possible molecular mechanisms underlying interindividual differences in the heart rate response to exercise. In this review, we comprehensively searched candidate gene, linkage, and genome-wide association studies that investigated the heart rate response to exercise. A total of ten genes were associated with the acute heart rate response to exercise in candidate gene studies. Only one gene (CHRM2), related to heart rate recovery, was replicated in recent genome-wide association studies (GWASs). Additional 17 candidate causal genes were identified for heart rate increase and 26 for heart rate recovery in these GWASs. Nine of these genes were associated with both acute increase and recovery of the heart rate during exercise. These genes can be broadly categorized into four categories: (1) development of the nervous system (CCDC141, PAX2, SOX5, and CAV2); (2) prolongation of neuronal life span (SYT10); (3) cardiac development (RNF220 and MCTP2); (4) cardiac rhythm (SCN10A and RGS6). Additional 10 genes were linked to long-term modification of the heart rate response to exercise, nine with heart rate increase and one with heart rate recovery. Follow-up will be essential to get functional insights in how candidate causal genes affect the heart rate response to exercise. Future work will be required to translate these findings to preventive and therapeutic applications.

A genome-wide association study of energy intake and expenditure

Excessive energy intake or insufficient energy expenditure, which result in energy imbalance, contribute to the development of obesity. Obesity-related genes, such as FTO, are associated with energy traits. No genome-wide association studies (GWAS) have been conducted to detect the genetic associations with energy-related traits, including energy intake and energy expenditure, among European-ancestry populations. In this study, we conducted a genome-wide study using pooled GWAS including 12,030 European-ancestry women and 6,743 European-ancestry men to identify genetic variants associated with these two energy traits. We observed a statistically significant genome-wide SNP heritability for energy intake of 6.05% (95%CI = (1.76, 10.34), P = 0.006); the SNP heritability for expenditure was not statistically significantly greater than zero. We discovered three SNPs on chromosome 12q13 near gene ANKRD33 that were genome-wide significantly associated with increased total energy intake among all men. We also identified signals on region 2q22 that were associated with energy expenditure among lean people. Body mass index related SNPs were found to be significantly associated with energy intake and expenditure through SNP set analyses. Larger GWAS studies of total energy traits are warranted to explore the genetic basis of energy intake, including possible differences between men and women, and the association between total energy intake and other downstream phenotypes, such as diabetes and chronic diseases.

Thirty loci identified for heart rate response to exercise and recovery implicate autonomic nervous system

Impaired capacity to increase heart rate (HR) during exercise (ΔHR^ex), and a reduced rate of recovery post-exercise (ΔHR^rec) are associated with higher cardiovascular mortality rates. Currently, the genetic basis of both phenotypes remains to be elucidated. We conduct genome-wide association studies (GWASs) for ΔHR^ex and ΔHR^rec in ~40,000 individuals, followed by replication in ~27,000 independent samples, all from UK Biobank. Six and seven single-nucleotide polymorphisms for ΔHR^ex and ΔHR^rec, respectively, formally replicate. In a full data set GWAS, eight further loci for ΔHR^ex and nine for ΔHR^rec are genome-wide significant (P ≤ 5 × 10⁻⁸). In total, 30 loci are discovered, 8 being common across traits. Processes of neural development and modulation of adrenergic activity by the autonomic nervous system are enriched in these results. Our findings reinforce current understanding of HR response to exercise and recovery and could guide future studies evaluating its contribution to cardiovascular risk prediction.

Genome-wide association studies have identified multiple loci for resting heart rate (HR) but the genetic factors associated with HR increase during and HR recovery after exercise are less well studied. Here, the authors examine both traits in a two-stage GWAS design in up to 67,257 individuals from UK Biobank.

Genome-wide association study for milking speed in French Holstein cows

Using a combination of data from the BovineSNP50 BeadChip SNP array (Illumina, San Diego, CA) and a EuroGenomics (Amsterdam, the Netherlands) custom single nucleotide polymorphism (SNP) chip with SNP pre-selected from whole genome sequence data, we carried out an association study of milking speed in 32,491 French Holstein dairy cows. Milking speed was measured by a score given by the farmer. Phenotypes were yield deviations as obtained from the French evaluation system. They were analyzed with a linear mixed model for association studies. We identified SNP on 22 chromosomes significantly associated with milking speed. As clinical mastitis and somatic cell score have an unfavorable genetic correlation with milking speed, we tested whether the most significant SNP on these 22 chromosomes associated with milking speed were also associated with clinical mastitis or somatic cell score. Nine hundred seventy-one genome-wide significant SNP were associated with milking speed. Of these, 86 were associated with clinical mastitis and 198 with somatic cell score. The most significant association signals for milking speed were observed on chromosomes 7, 8, 10, 14, and 18. The most significant signal was located on chromosome 14 (ZFAT gene). Eleven novel milking speed quantitative trait loci (QTL) were observed on chromosomes 7, 10, 11, 14, 18, 25, and 26. Twelve candidate SNP for milking speed mapped directly within genes. Of these, 10 were QTL lead SNP, which mapped within the genes HMHA1, POLR2E, GNB5, KLHL29, ZFAT, KCNB2, CEACAM18, CCL24, and LHPP. Limited pleiotropy was observed between milking speed QTL and clinical mastitis.

Sequencing of Linkage Region on Chromosome 12p11 Identifies PKP2 as a Candidate Gene for Left Ventricular Mass in Dominican Families

Increased left ventricular mass (LVM) is an intermediate phenotype for cardiovascular disease (CVD) and a predictor of stroke. Using families from the Dominican Republic, we have previously shown LVM to be heritable and found evidence for linkage to chromosome 12p11. Our current study aimed to further characterize the QTL by sequencing the 1 LOD unit down region in 10 families from the Dominican Republic with evidence for linkage to LVM. Within this region, we tested 5477 common variants [CVs; minor allele frequency (MAF) ≥5%] using the Quantitative Transmission-Disequilibrium Test (QTDT). Gene-based analyses were performed to test rare variants (RVs; MAF < 5%) in 181 genes using the family-based sequence kernel association test. A sample of 618 unrelated Dominicans from the Northern Manhattan Study (NOMAS) and 12 Dominican families with Exome Array data were used for replication analyses. The most strongly associated CV with evidence for replication was rs1046116 (Discovery families P = 9.0 × 10-4; NOMAS P = 0.03; replication families P = 0.46), a missense variant in PKP2 In nonsynonymous RV analyses, PKP2 was one of the most strongly associated genes (P = 0.05) with suggestive evidence for replication in NOMAS (P = 0.05). PKP2 encodes the plakophilin 2 protein and is a desmosomal gene implicated in arrythmogenic right ventricular cardiomyopathy and recently in arrhythmogenic left ventricular cardiomyopathy, which makes PKP2 an excellent candidate gene for LVM. In conclusion, sequencing of our previously reported QTL identified common and rare variants within PKP2 to be associated with LVM. Future studies are necessary to elucidate the role these variants play in influencing LVM.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500–6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east–west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

The Scandinavian peninsula was the last part of Europe to be colonized after the Last Glacial Maximum. The migration routes, cultural networks, and the genetic makeup of the first Scandinavians remain elusive and several hypotheses exist based on archaeology, climate modeling, and genetics. By analyzing the genomes of early Scandinavian hunter-gatherers, we show that their migrations followed two routes: one from the south and another from the northeast along the ice-free Norwegian Atlantic coast. These groups met and mixed in Scandinavia, creating a population more diverse than contemporaneous central and western European hunter-gatherers. As northern Europe is associated with cold and low light conditions, we investigated genomic patterns of adaptation to these conditions and genes known to be involved in skin pigmentation. We demonstrate that Mesolithic Scandinavians had higher levels of light pigmentation variants compared to the respective source populations of the migrations, suggesting adaptation to low light levels and a surprising signal of genetic continuity in TMEM131, a gene that may be involved in long-term adaptation to the cold.

Natural Selection on Genes Related to Cardiovascular Health in High-Altitude Adapted Andeans

The increase in red blood cell mass (polycythemia) due to the reduced oxygen availability (hypoxia) of residence at high altitude or other conditions is generally thought to be beneficial in terms of increasing tissue oxygen supply. However, the extreme polycythemia and accompanying increased mortality due to heart failure in chronic mountain sickness most likely reduces fitness. Tibetan highlanders have adapted to high altitude, possibly in part via the selection of genetic variants associated with reduced polycythemic response to hypoxia. In contrast, high-altitude-adapted Quechua- and Aymara-speaking inhabitants of the Andean Altiplano are not protected from high-altitude polycythemia in the same way, yet they exhibit other adaptive features for which the genetic underpinnings remain obscure. Here, we used whole-genome sequencing to scan high-altitude Andeans for signals of selection. The genes showing the strongest evidence of selection-including BRINP3, NOS2, and TBX5-are associated with cardiovascular development and function but are not in the response-to-hypoxia pathway. Using association mapping, we demonstrated that the haplotypes under selection are associated with phenotypic variations related to cardiovascular health. We hypothesize that selection in response to hypoxia in Andeans could have vascular effects and could serve to mitigate the deleterious effects of polycythemia rather than reduce polycythemia itself.

Functional polymorphism of the renalase gene is associated with cardiac hypertrophy in female patients with aortic stenosis

Renalase decreases circulating catecholamines concentration and is important in maintaining primary cellular metabolism. Renalase acts through the plasma membrane calcium ATPase 4b in the heart, which affects pressure overload but not exercise induced heart hypertrophy. The aim of this study was to test the association between a functional polymorphism Glu37Asp (rs2296545) of the renalase gene and left ventricular hypertrophy in a large cohort of patients with aortic stenosis. The study group consisted of 657 patients with aortic stenosis referred for aortic valve replacement. Preoperative echocardiographic assessment was performed to obtain cardiac phenotypes. Generalized-linear models were implemented to analyze data using crude or full model adjusted for selected clinical factors. In females, the Asp37 variant of the Glu37Asp polymorphism was associated with higher left ventricular mass (p = 0.0021 and p = 0.055 crude and full model respectively), intraventricular septal thickness (p = 0.0003 and p = 0.0143) and posterior wall thickness (p = 0.0005 and p = 0.0219) all indexed to body surface area, as well as relative wall thickness (p = 0.001 and p = 0.0097). No significant associations were found among the male patients. In conclusion, we have found the association of the renalase Glu37Asp polymorphism with left ventricle hypertrophy in large group of females with aortic stenosis. The Glu37Asp polymorphism causes not only amino-acid substitution in FAD binding domain but may also change binding affinity of the hypoxia- and hypertrophy-related transcription factors and influence renalase gene expression. Our data suggest that renalase might play a role in hypertrophic response to pressure overload, but the exact mechanism requires further investigation.

Long-term intra-individual reproducibility of heart rate dynamics during exercise and recovery in the UK Biobank cohort

Background

The heart rate (HR) response to exercise provides useful information about the autonomic function and has prognostic value, but its reproducibility over a long period of time, a critical requirement for using it as a clinical biomarker, is undetermined.

Aim

To determine the intra-individual reproducibility of HR dynamics during sub-maximum exercise and one minute recovery.

Methods

1187 individuals from the Cardio physical fitness assessment test of the UK Biobank repeated a standard exercise stress test twice (recall time 34.2 ± 2.8 months) and were prospectively studied.

Results

821 individuals complied with inclusion criteria for reproducibility analysis, including peak workload differences between assessments ≤10 W. Intra-individual correlation between HR profile during the first and the second assessment was very high and higher than inter-individual correlation (0.92±0.08 vs 0.87±0.11, p<0.01). Intra-individual correlation of indices describing HR dynamics was: ρ = 0.81 for maximum HR during exercise; ρ = 0.71 for minimum HR during recovery; ρ = 0.70 for HR changes during both exercise and recovery; Intra-individual correlation was higher for these indices of HR dynamics than for resting HR (ρ = 0.64). Bland-Altman plots demonstrated good agreement between HR indices estimated during the first and second assessment. A small but consistent bias was registered for all repeated measurements. The intra-individual consistency of abnormal values was about 60–70%.

Conclusions

The HR dynamics during exercise and recovery are reproducible over a period of 3 years, with moderate to strong intra-individual reproducibility of abnormal values.

Factors Associated With Left Atrial Remodeling in the General Population

BACKGROUND

Although contributors to remodeling of the left ventricle (LV) have been well studied in general population cohorts, few data are available describing factors influencing changes in left atrial (LA) structure.

METHODS AND RESULTS

Maximum LA volume was determined by cardiac magnetic resonance imaging among 748 participants in the Dallas Heart Study at 2 visits a mean of 8 years apart. Associations of changes in LA volume (ΔLAV) with traditional risk factors, biomarkers, LV geometry, and remodeling by cardiac magnetic resonance imaging and detailed measurements of global and regional adiposity (by magnetic resonance imaging and dual-energy x ray absorptiometry) were assessed using multivariable linear regression. Greater ΔLAV was independently associated with black and Hispanic race/ethnicity, change in systolic blood pressure, LV mass and ΔLV mass, N-terminal probrain natriuretic peptide and change in N-terminal probrain natriuretic peptide, and body mass index (P<0.05 for each). In subanalyses, the associations of ΔLAV with LV mass parameters were driven by associations with baseline and ΔLV end diastolic volume (P<0.0001 for each) and not wall thickness (P=0.21). Associations of ΔLAV with body mass index were explained exclusively by associations with visceral fat mass (P=0.002), with no association seen between ΔLAV and subcutaneous abdominal fat (P=0.47) or lower body fat (P=0.30).

CONCLUSIONS

Left atrial dilatation in the population is more common in black and Hispanic than in white individuals and is associated with parallel changes in the LV. LA dilatation may be mediated by blood pressure control and the development of visceral adiposity.

Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function

BACKGROUND

Understanding the genetic architecture of cardiac structure and function may help to prevent and treat heart disease. This investigation sought to identify common genetic variations associated with inter-individual variability in cardiac structure and function.

METHODS

A GWAS meta-analysis of echocardiographic traits was performed, including 46,533 individuals from 30 studies (EchoGen consortium). The analysis included 16 traits of left ventricular (LV) structure, and systolic and diastolic function.

RESULTS

The discovery analysis included 21 cohorts for structural and systolic function traits (n = 32,212) and 17 cohorts for diastolic function traits (n = 21,852). Replication was performed in 5 cohorts (n = 14,321) and 6 cohorts (n = 16,308), respectively. Besides 5 previously reported loci, the combined meta-analysis identified 10 additional genome-wide significant SNPs: rs12541595 near MTSS1 and rs10774625 in ATXN2 for LV end-diastolic internal dimension; rs806322 near KCNRG, rs4765663 in CACNA1C, rs6702619 near PALMD, rs7127129 in TMEM16A, rs11207426 near FGGY, rs17608766 in GOSR2, and rs17696696 in CFDP1 for aortic root diameter; and rs12440869 in IQCH for Doppler transmitral A-wave peak velocity. Findings were in part validated in other cohorts and in GWAS of related disease traits. The genetic loci showed associations with putative signaling pathways, and with gene expression in whole blood, monocytes, and myocardial tissue.

CONCLUSION

The additional genetic loci identified in this large meta-analysis of cardiac structure and function provide insights into the underlying genetic architecture of cardiac structure and warrant follow-up in future functional studies.

FUNDING

For detailed information per study, see Acknowledgments.

Exploring Overlaps Between the Genomic and Environmental Determinants of LVH and Stroke: A Multicenter Study in West Africa

BACKGROUND

Whether left ventricular hypertrophy (LVH) is determined by similar genomic and environmental risk factors with stroke, or is simply an intermediate stroke marker, is unknown.

OBJECTIVES

We present a research plan and preliminary findings to explore the overlap in the genomic and environmental determinants of LVH and stroke among Africans participating in the SIREN (Stroke Investigative Research and Education Network) study.

METHODS

SIREN is a transnational, multicenter study involving acute stroke patients and age-, ethnicity-, and sex-matched control subjects recruited from 9 sites in Ghana and Nigeria. Genomic and environmental risk factors and other relevant phenotypes for stroke and LVH are being collected and compared using standard techniques.

RESULTS

This preliminary analysis included only 725 stroke patients (mean age 59.1 ± 13.2 years; 54.3% male). Fifty-five percent of the stroke subjects had LVH with greater proportion among women (51.6% vs. 48.4%; p < 0.001). Those with LVH were younger (57.9 ± 12.8 vs. 60.6 ± 13.4; p = 0.006) and had higher mean systolic and diastolic blood pressure (167.1/99.5 mm Hg vs 151.7/90.6 mm Hg; p < 0.001). Uncontrolled blood pressure at presentation was prevalent in subjects with LVH (76.2% vs. 57.7%; p < 0.001). Significant independent predictors of LVH were age <45 years (adjusted odds ratio [AOR]: 1.91; 95% confidence interval [CI]: 1.14 to 3.19), female sex (AOR: 2.01; 95% CI: 1.44 to 2.81), and diastolic blood pressure > 90 mm Hg (AOR: 2.10; 95% CI: 1.39 to 3.19; p < 0.001).

CONCLUSIONS

The prevalence of LVH was high among stroke patients especially the younger ones, suggesting a genetic component to LVH. Hypertension was a major modifiable risk factor for stroke as well as LVH. It is envisaged that the SIREN project will elucidate polygenic overlap (if present) between LVH and stroke among Africans, thereby defining the role of LVH as a putative intermediate cardiovascular phenotype and therapeutic target to inform interventions to reduce stroke risk in populations of African ancestry.

Genetic modifiers of cardiorespiratory fitness response to lifestyle intervention

PURPOSE

Numerous prospective studies indicate that improved cardiorespiratory fitness reduces type 2 diabetes risk and delays disease progression. We hypothesized that genetic variants modify fitness response to an intensive lifestyle intervention (ILI) in the Action for Health in Diabetes (Look AHEAD) randomized clinical trial, aimed to detect whether ILI will reduce cardiovascular events in overweight/obese subjects with type 2 diabetes compared with a standard of care.

METHODS

Polymorphisms in established fitness genes and in all loci assayed on the Illumina CARe iSelect chip were examined as predictors of change in MET level, estimated using a treadmill test, in response to a 1-yr intervention in 3899 participants.

RESULTS

We identified a significant signal in previously reported fitness-related gene RUNX1 that was associated with 1-yr METs response in ILI (0.19 ± 0.04 MET less improvement per minor allele copy; P = 1.9 × 10(-5)) and genotype-intervention interaction (P = 4.8 × 10(-3)). In the chipwide analysis, FKBP7 rs17225700 showed a significant association with ILI response among subjects not receiving beta-blocker medications (0.47 ± 0.09 METs less improvement; P = 5.3 × 10(-5)) and genotype-treatment interaction (P = 5.3 × 10(-7)). The Gene Relationships Among Implicated Loci pathway-based analysis identified connections between associated genes, including those influencing vascular tone, muscle contraction, cardiac energy substrate dynamics, and muscle protein synthesis.

CONCLUSIONS

This is the first study to identify genetic variants associated with fitness responses to a randomized lifestyle intervention in overweight/obese diabetic individuals. RUNX1 and FKBP7, involved in erythropoesis and muscle protein synthesis, respectively, are related to change in cardiorespiratory fitness in response to exercise.

Genetic Regulation of Endothelial Vasomotor Function

The endothelium plays an important role in the regulation of vasomotor tone and the maintenance of vascular integrity. Endothelial dysfunction, i.e., impaired endothelial dependent dilation, is a fundamental component of the pathogenesis of cardiovascular disease. Although endothelial dysfunction is associated with a number of cardiovascular disease risk factors, those risk factors are not the only determinants of endothelial dysfunction. Despite knowing many molecules involved in endothelial signaling pathways, the genetic contribution to endothelial function has yet to be fully elucidated. This mini-review summarizes current evidence supporting the genetic contribution to endothelial vasomotor function. Findings from population-based studies, association studies for candidate genes, and unbiased large genomic scale studies in humans and rodent models are discussed. A brief synopsis of the current studies addressing the genetic regulation of endothelial responses to exercise training is also included.

A new statistical framework for genetic pleiotropic analysis of high dimensional phenotype data

BACKGROUND

The widely used genetic pleiotropic analyses of multiple phenotypes are often designed for examining the relationship between common variants and a few phenotypes. They are not suited for both high dimensional phenotypes and high dimensional genotype (next-generation sequencing) data. To overcome limitations of the traditional genetic pleiotropic analysis of multiple phenotypes, we develop sparse structural equation models (SEMs) as a general framework for a new paradigm of genetic analysis of multiple phenotypes. To incorporate both common and rare variants into the analysis, we extend the traditional multivariate SEMs to sparse functional SEMs. To deal with high dimensional phenotype and genotype data, we employ functional data analysis and the alternative direction methods of multiplier (ADMM) techniques to reduce data dimension and improve computational efficiency.

RESULTS

Using large scale simulations we showed that the proposed methods have higher power to detect true causal genetic pleiotropic structure than other existing methods. Simulations also demonstrate that the gene-based pleiotropic analysis has higher power than the single variant-based pleiotropic analysis. The proposed method is applied to exome sequence data from the NHLBI's Exome Sequencing Project (ESP) with 11 phenotypes, which identifies a network with 137 genes connected to 11 phenotypes and 341 edges. Among them, 114 genes showed pleiotropic genetic effects and 45 genes were reported to be associated with phenotypes in the analysis or other cardiovascular disease (CVD) related phenotypes in the literature.

CONCLUSIONS

Our proposed sparse functional SEMs can incorporate both common and rare variants into the analysis and the ADMM algorithm can efficiently solve the penalized SEMs. Using this model we can jointly infer genetic architecture and casual phenotype network structure, and decompose the genetic effect into direct, indirect and total effect. Using large scale simulations we showed that the proposed methods have higher power to detect true causal genetic pleiotropic structure than other existing methods.

Genetic Research and Women's Heart Disease: a Primer

PURPOSE OF REVIEW

This review provides a brief synopsis of sexual dimorphism in atherosclerosis with an emphasis on genetic studies aimed to better understand the atherosclerotic process and clinical outcomes in women. Such studies are warranted because development of atherosclerosis, impact of several traditional risk factors, and burden of coronary heart disease (CHD) differ between women and men.

RECENT FINDINGS

While most candidate gene studies pool women and men and adjust for sex, some sex-specific studies provide evidence of association between candidate genes and prevalent and incident CHD in women. So far, most genome-wide association studies (GWAS) also failed to consider sex-specific associations. The few GWAS focused on women tended to have small sample sizes and insufficient power to reject the null hypothesis of no association even if associations exist. Few studies consider that sex can modify the effect of gene variants on CHD. Sufficiently large-scale genetic studies in women of different race/ethnic groups, taking into account possible gene-gene and gene-environment interactions as well as hormone-mediated epigenetic mechanisms, are needed. Using the same disease definition for women and men might not be appropriate. Accurate phenotyping and inclusion of relevant outcomes in women, together with targeting the entire spectrum of atherosclerosis, could help address the contribution of genes to sexual dimorphism in atherosclerosis. Discovered genetic loci should be taken forward for replication and functional studies to elucidate the plausible underlying biological mechanisms. A better understanding of the etiology of atherosclerosis in women would facilitate future prevention efforts and interventions.

Strain survey and genetic analysis of vasoreactivity in mouse aorta

Understanding the genetic influence on vascular reactivity is important for identifying genes underlying impaired vascular function. The purpose of this study was to characterize the genetic contribution to intrinsic vascular function and to identify loci associated with phenotypic variation in vascular reactivity in mice. Concentration response curves to phenylephrine (PE), potassium chloride (KCl), acetylcholine (ACh), and sodium nitroprusside (SNP) were generated in aortic rings from male mice (12 wk old) from 27 inbred mouse strains. Significant strain-dependent differences were found for both maximal responses and sensitivity for each agent, except for SNP Max (%). Strain differences for maximal responses to ACh, PE, and KCl varied by two- to fivefold. On the basis of these large strain differences, we performed genome-wide association mapping (GWAS) to identify loci associated with variation in responses to these agents. GWAS for responses to ACh identified four significant and 19 suggestive loci. Several suggestive loci for responses to SNP, PE, and KCl (including one significant locus for KCl EC₅₀) were also identified. These results demonstrate that intrinsic endothelial function, and more generally vascular function, is genetically determined and associated with multiple genomic loci. Furthermore, these results are supported by the finding that several genes residing in significant and suggestive loci for responses to ACh were previously identified in rat and/or human quantitative trait loci/GWAS for cardiovascular disease. This study represents the first step toward the unbiased comprehensive discovery of genetic determinants that regulate intrinsic vascular function, particularly endothelial function.

Pulmonary artery enlargement and cystic fibrosis pulmonary exacerbations: a cohort study

Background

Acute pulmonary exacerbations are associated with progressive lung function decline and increased mortality in cystic fibrosis (CF). The role of pulmonary vascular disease in pulmonary exacerbations is unknown. We investigated the association between pulmonary artery enlargement (PA:A>1), a marker of pulmonary vascular disease, and exacerbations.

Methods

We analyzed clinical, computed tomography (CT), and prospective exacerbation data in a derivation cohort of 74 adult CF patients, measuring the PA:A at the level of the PA bifurcation. We then replicated our findings in a validation cohort of 190 adult CF patients. Patients were separated into groups based on the presence or absence of a PA:A>1 and were followed for 1-year in the derivation cohort and 2-years in the validation cohort. The primary endpoint was developing ≥1 acute pulmonary exacerbation during follow-up. Linear and logistic regression models were used to determine associations between clinical factors, the PA:A ratio, and pulmonary exacerbations. We used Cox regression to determine time to first exacerbation in the validation cohort.

Findings

We found that PA:A>1 was present in n=37/74 (50%) of the derivation and n=89/190 (47%) of the validation cohort. In the derivation cohort, n=50/74 (68%) had ≥1 exacerbation at 1 year and n=133/190 (70%) in the validation cohort had ≥1 exacerbation after 2 years. PA:A>1 was associated with younger age in both cohorts and with elevated sweat chloride (100.5±10.9 versus 90.4±19.9mmol/L, difference between groups 10.1mmol/L [95%CI 2.5–17.7], P=0.017) in the derivation group. PA:A>1 was associated with exacerbations in the derivation (OR 3.49, 95%CI 1.18–10.3, P=0.023) and validation (OR 2.41, 95%CI 1.06–5.52, P=0.037) cohorts when adjusted for confounders. Time to first exacerbation was shorter in PA:A>1 versus PA:A<1 [HR 1.66 (95%CI 1.18–2.34), P=0.004] in unadjusted analysis, but not when adjusted for sex, BMI, prior exacerbation, positive Pseudomonas status, and FEV1/FVC [HR 1.14 (95%CI 0.80–1.62), P=0.82]).

Interpretation

PA enlargement is prevalent in adult CF patients and is associated with acute pulmonary exacerbation risk in two well-characterized cohorts. PA:A may be a predictive marker in CF.

Cohort Profile: The Framingham Heart Study (FHS): overview of milestones in cardiovascular epidemiology

The Framingham Heart Study (FHS) has conducted seminal research defining cardiovascular disease (CVD) risk factors and fundamentally shaping public health guidelines for CVD prevention over the past five decades. The success of the Original Cohort, initiated in 1948, paved the way for further epidemiological research in preventive cardiology. Due to the keen observations suggesting the role of shared familial factors in the development of CVD, in 1971 the FHS began enroling the second generation cohort, comprising the children of the Original Cohort and the spouses of the children. In 2002, the third generation cohort, comprising the grandchildren of the Original Cohort, was initiated to additionally explore genetic contributions to CVD in greater depth. Additionally, because of the predominance of White individuals of European descent in the three generations of FHS participants noted above, the Heart Study enrolled the OMNI1 and OMNI2 cohorts in 1994 and 2003, respectively, aimed to reflect the current greater racial and ethnic diversity of the town of Framingham. All FHS cohorts have been examined approximately every 2-4 years since the initiation of the study. At these periodic Heart Study examinations, we obtain a medical history and perform a cardiovascular-focused physical examination, 12-lead electrocardiography, blood and urine samples testing and other cardiovascular imaging studies reflecting subclinical disease burden.The FHS has continually evolved along the cutting edge of cardiovascular science and epidemiological research since its inception. Participant studies now additionally include study of cardiovascular imaging, serum and urine biomarkers, genetics/genomics, proteomics, metabolomics and social networks. Numerous ancillary studies have been established, expanding the phenotypes to encompass multiple organ systems including the lungs, brain, bone and fat depots, among others. Whereas the FHS was originally conceived and designed to study the epidemiology of cardiovascular disease, it has evolved over the years with staggering expanded breadth and depth that have far greater implications in the study of the epidemiology of a wide spectrum of human diseases. The FHS welcomes research collaborations using existing or new collection of data. Detailed information regarding the procedures for research application submission and review are available at [http://www.framinghamheartstudy.org/researchers/index.php].

Genetic influence on left ventricular structure and function: a Korean twin and family study

Genetic factors have been suggested to be one of the determinants of the variation of left ventricular (LV) structure and function. However, the heritability range of LV structure varies across studies and the influence of genetics on LV function is not well established, especially in Asian populations. Study subjects were 1,642 healthy Korean adults from 426 families, consisting of 298 pairs of monozygotic twins, 62 pairs of dizygotic twins, one set of triplets, 567 siblings, and 354 parents. LV structure and function were measured by M-mode and 2D echocardiography, and conventional and tissue Doppler imaging (TDI). Pairwise intra-class correlations for various familial relationships and heritability were estimated for LV structure and function. The heritability of LV mass, LV ejection fraction (LVEF), left atrial volume index, the ratio between early and late diastolic velocity of mitral inflow (E/A ratio), and the ratio between early diastolic velocity of mitral inflow and early diastolic mitral annular velocities (E/Ea ratio) was 0.44, 0.27, 0.44, 0.25, and 0.33, respectively. Bivariate genetic analysis showed that LV structural and functional traits had significant genetic correlations with cardiovascular risk factors. Additive genetic correlation (ρG) of LV mass with body mass index, systolic blood pressure, and high density lipoprotein cholesterol were 0.49, 0.42, and -0.15 respectively. LVEF (ρG = 0.33) and left atrial volume index (ρG = 0.24) also had a significant genetic correlation with systolic blood pressure. These findings support the theory that genetic factors have significant influence on these traits and necessitate further work to identify the specific genes involved.

Two for the Price of One: G Protein-Dependent and -Independent Functions of RGS6 In Vivo

Regulator of G protein signaling 6 (RGS6) is unique among the members of the RGS protein family as it remains the only protein with the demonstrated capacity to control G protein-dependent and -independent signaling cascades in vivo. RGS6 inhibits signaling mediated by γ-aminobutyric acid B receptors, serotonin 1A receptors, μ opioid receptors, and muscarinic acetylcholine 2 receptors. RGS6 deletion triggers distinct behavioral phenotypes resulting from potentiated signaling by these G protein-coupled receptors namely ataxia, a reduction in anxiety and depression, enhanced analgesia, and increased parasympathetic tone, respectively. In addition, RGS6 possesses potent proapoptotic and growth suppressive actions. In heart, RGS6-dependent reactive oxygen species (ROS) production promotes doxorubicin (Dox)-induced cardiomyopathy, while in cancer cells RGS6/ROS signaling is necessary for activation of the ataxia telangiectasia mutated/p53/apoptosis pathway required for the chemotherapeutic efficacy of Dox. Further, by facilitating Tip60 (trans-acting regulator protein of HIV type 1-interacting protein 60 kDa)-dependent DNA methyltransferase 1 degradation, RGS6 suppresses cellular transformation in response to oncogenic Ras. The culmination of these G protein-independent actions results in potent tumor suppressor actions of RGS6 in the murine mammary epithelium. This work summarizes evidence from human genetic studies and model animals implicating RGS6 in normal physiology, disease, and the pharmacological actions of multiple drugs. Though efforts by multiple laboratories have contributed to the ever-growing RGS6 oeuvre, the pleiotropic nature of this gene will likely lead to additional work detailing the importance of RGS6 in neuropsychiatric disorders, cardiovascular disease, and cancer.

Telomere length in the two extremes of abnormal fetal growth and the programming effect of maternal arterial hypertension

We tested the hypothesis that leukocyte telomere length (LTL) is associated with birth weight in both extremes of abnormal fetal growth: small (SGA) and large for gestational age newborns (LGA). Clinical and laboratory variables of the mothers and the neonates were explored; 45 newborns with appropriate weight for gestational age (AGA), 12 SGA and 12 LGA were included. Whether the differences might be explained by variation in OBFC1 (rs9419958) and CTC1 (rs3027234) genes associated with LTL was determined. A significant association between birth weight and LTL was observed; LTL was significantly shorter in LGA newborns (1.01 ± 0.12) compared with SGA (1.73 ± 0.19) p < 0.005, mean ± SE. Maternal (Spearman R = −0.6, p = 0.03) and neonatal LTL (R = −0.25, p = 0.03) were significantly and inversely correlated with maternal history of arterial hypertension in previous gestations. Neonatal LTL was not significantly associated with either rs9419950 or rs3027234, suggesting that the association between neonatal LTL and birth weight is not influenced by genetic variation in genes that modify the interindividual LTL. In conclusion, telomere biology seems to be modulated by abnormal fetal growth; modifications in telomere length might be programmed by an adverse environment in utero.

Entropy-based method for assessing the influence of genetic markers and covariates on hypertension: application to Genetic Analysis Workshop 18 data

Many complex diseases are related to genetics, and it is of great interest to evaluate the association between single-nucleotide polymorphisms (SNPs) and disease outcome. The association of genetics with outcome can be modified by covariates such as age, sex, smoking status, and membership to the same pedigree. In this paper, we propose a block entropy method to separate two classes of SNPs, for which the association with hypertension is either sensitive or insensitive to the covariates. We also propose a consistency entropy method to further reduce the number of SNPs that might be associated with the outcome. Based on the data provided by the organizers of Genetic Analysis Workshop 18, we calculated the block entropies for six different blocking strategies. Using block entropy and consistency entropy, we identified 230 SNPs on chromosome 9 that are most likely to be associated with the outcome and whose associations with hypertension are sensitive to the covariates.

Canonical correlation analysis for gene-based pleiotropy discovery

Genome-wide association studies have identified a wealth of genetic variants involved in complex traits and multifactorial diseases. There is now considerable interest in testing variants for association with multiple phenotypes (pleiotropy) and for testing multiple variants for association with a single phenotype (gene-based association tests). Such approaches can increase statistical power by combining evidence for association over multiple phenotypes or genetic variants respectively. Canonical Correlation Analysis (CCA) measures the correlation between two sets of multidimensional variables, and thus offers the potential to combine these two approaches. To apply CCA, we must restrict the number of attributes relative to the number of samples. Hence we consider modules of genetic variation that can comprise a gene, a pathway or another biologically relevant grouping, and/or a set of phenotypes. In order to do this, we use an attribute selection strategy based on a binary genetic algorithm. Applied to a UK-based prospective cohort study of 4286 women (the British Women's Heart and Health Study), we find improved statistical power in the detection of previously reported genetic associations, and identify a number of novel pleiotropic associations between genetic variants and phenotypes. New discoveries include gene-based association of NSF with triglyceride levels and several genes (ACSM3, ERI2, IL18RAP, IL23RAP and NRG1) with left ventricular hypertrophy phenotypes. In multiple-phenotype analyses we find association of NRG1 with left ventricular hypertrophy phenotypes, fibrinogen and urea and pleiotropic relationships of F7 and F10 with Factor VII, Factor IX and cholesterol levels.

Pleiotropy appears when a variation in one gene affects to several non-related phenotypes. The study of this phenomenon can be useful in gene function discovery, but also in the study of the evolution of a gene. In this paper, we present a methodology, based on Canonical Correlation Analysis, which studies gene-centered multiple association of the variation of SNPs in one or a set of genes with one or a set of phenotypes. The resulting methodology can be applied in gene-centered association analysis, multiple association analysis or pleiotropic pattern discovery. We apply this methodology with a genotype dataset and a set of cardiovascular related phenotypes, and discover new gene association between gene NRG1 and phenotypes related with left ventricular hypertrophy, and pleiotropic effects of this gene with other phenotypes as coagulation factors and urea or pleiotropic effects between coagulation related genes F7 and F10 with coagulation factors and cholesterol levels. This methodology could be also used to find multiple associations in other omics datasets.

Association of parental blood pressure with retinal microcirculatory abnormalities indicative of endothelial dysfunction in children

OBJECTIVE

Microcirculatory abnormalities precede the onset of hypertension and may explain its familial nature. We examined the relationship between parental blood pressure (BP) and offspring retinal microvasculature in Pakistani trios [father, mother, and child (aged 9-14 years)].

METHODS

This is a substudy of a population-based trial of BP reduction. Data were available on 358 normotensive, and 410 offspring of at least one hypertensive parent. Retinal vessel characteristics were measured from digital images. Multivariable linear regression models were built to assess the associations between maternal and paternal BP and offspring retinal microvasculature.

RESULTS

Optimality deviation was greatest in offspring of two hypertensive parents, compared with those with one or no hypertensive parent (P=0.030 for trend). Paternal SBP and DBP were each significantly associated with optimality deviation in offspring (P=0.023 and P=0.006, respectively). This relationship persisted after accounting for offspring cardiovascular risk factors [increase in optimality deviation (95% confidence interval, CI) 0.0053 (0.0001-0.0106, P=0.047) and 0.0109 (0.0025-0.0193, P=0.011), for each 10 mmHg increase in paternal SBP and DBP, respectively]. Maternal DBP was inversely associated with offspring arteriovenous ratio -0.0102 (-0.0198 to -0.0007, P=0.035).

CONCLUSION

Microvascular endothelial dysfunction in children is associated with increasing levels of parental hypertension. The association with paternal BP is independent of other cardiovascular risk factors, including the child's BP. Higher maternal DBP is associated with evidence of arteriolar narrowing in offspring. These early microcirculatory changes may help explain familial predisposition to hypertension in people of Pakistani origin at an early age.

VIDEO ABSTRACT

:

WDR72 models of structure and function: a stage-specific regulator of enamel mineralization

Amelogenesis Imperfecta (AI) is a clinical diagnosis that encompasses a group of genetic mutations, each affecting processes involved in tooth enamel formation and thus, result in various enamel defects. The hypomaturation enamel phenotype has been described for mutations involved in the later stage of enamel formation, including Klk4, Mmp20, C4orf26, and Wdr72. Using a candidate gene approach we discovered a novel Wdr72 human mutation in association with AI to be a 5-base pair deletion (c.806_810delGGCAG; p.G255VfsX294). To gain insight into the function of WDR72, we used computer modeling of the full-length human WDR72 protein structure and found that the predicted N-terminal sequence forms two beta-propeller folds with an alpha-solenoid tail at the C-terminus. This domain iteration is characteristic of vesicle coat proteins, such as beta'-COP, suggesting a role for WDR72 in the formation of membrane deformation complexes to regulate intracellular trafficking. Our Wdr72 knockout mouse model (Wdr72(-/-)), containing a LacZ reporter knock-in, exhibited hypomineralized enamel similar to the AI phenotype observed in humans with Wdr72 mutations. MicroCT scans of Wdr72(-/-) mandibles affirmed the hypomineralized enamel phenotype occurring at the onset of the maturation stage. H&E staining revealed a shortened height phenotype in the Wdr72(-/-) ameloblasts with retained proteins in the enamel matrix during maturation stage. H(+)/Cl(-) exchange transporter 5 (CLC5), an early endosome acidifier, was co-localized with WDR72 in maturation-stage ameloblasts and decreased in Wdr72(-/-) maturation-stage ameloblasts. There were no obvious differences in RAB4A and LAMP1 immunostaining of Wdr72(-/-) mice as compared to wildtype controls. Moreover, Wdr72(-/-) ameloblasts had reduced amelogenin immunoreactivity, suggesting defects in amelogenin fragment resorption from the matrix. These data demonstrate that WDR72 has a major role in enamel mineralization, most notably during the maturation stage, and suggest a function involving endocytic vesicle trafficking, possibly in the removal of amelogenin proteins.

Pulse pressure relation to aortic and left ventricular structure in the Age, Gene/Environment Susceptibility (AGES)-Reykjavik Study

High pulse pressure, a major cardiovascular risk factor, has been attributed to medial elastic fiber degeneration and aortic dilation, which transfers hemodynamic load to stiffer collagen. However, recent studies suggest higher pulse pressure is instead associated with smaller aortic diameter. Thus, we sought to elucidate relations of pulse pressure with aortic stiffness and aortic and cardiac dimensions. We used magnetic resonance imaging to examine relationships of pulse pressure with lumen area and wall stiffness and thickness in the thoracic aorta and left ventricular structure in 526 participants (72-94 years of age, 295 women) in the community-based Age, Gene/Environment Susceptibility-Reykjavik Study. In a multivariable model that adjusted for age, sex, height, weight, and standard vascular risk factors, central pulse pressure had a negative relationship with aortic lumen area (all effects expressed as mm Hg/SD; B=-8.1±1.2; P<0.001) and positive relationships with left ventricular end-diastolic volume (B=3.8±1.0; P<0.001), carotid-femoral pulse wave velocity (B=3.6±1.0; P<0.001), and aortic wall area (B=3.0±1.2; P=0.015). Higher pulse pressure in older people is associated with smaller aortic lumen area and greater aortic wall stiffness and thickness and left ventricular volume. Relationships of larger ventricular volume and smaller aortic lumen with higher pulse pressure suggest mismatch in hemodynamic load accommodation by the heart and aorta in older people.

Genomic structure of nucleotide diversity among Lyon rat models of metabolic syndrome

Background

The metabolic syndrome (MetS), a complex disorder involving hypertension, obesity, dyslipidemia and insulin resistance, is a major risk factor for heart disease, stroke, and diabetes. The Lyon Hypertensive (LH), Lyon Normotensive (LN) and Lyon Low-pressure (LL) rats are inbred strains simultaneously derived from a common outbred Sprague Dawley colony by selection for high, normal, and low blood pressure, respectively. Further studies found that LH is a MetS susceptible strain, while LN is resistant and LL has an intermediate phenotype. Whole genome sequencing determined that, while the strains are phenotypically divergent, they are nearly 98% similar at the nucleotide level. Using the sequence of the three strains, we applied an approach that harnesses the distribution of Observed Strain Differences (OSD), or nucleotide diversity, to distinguish genomic regions of identity-by-descent (IBD) from those with divergent ancestry between the three strains. This information was then used to fine-map QTL identified in a cross between LH and LN rats in order to identify candidate genes causing the phenotypes.

Results

We identified haplotypes that, in total, contain at least 95% of the identifiable polymorphisms between the Lyon strains that are likely of differing ancestral origin. By intersecting the identified haplotype blocks with Quantitative Trait Loci (QTL) previously identified in a cross between LH and LN strains, the candidate QTL regions have been narrowed by 78%. Because the genome sequence has been determined, we were further able to identify putative functional variants in genes that are candidates for causing the QTL.

Conclusions

Whole genome sequence analysis between the LH, LN, and LL strains identified the haplotype structure of these three strains and identified candidate genes with sequence variants predicted to affect gene function. This approach, merged with additional integrative genetics approaches, will likely lead to novel mechanisms underlying complex disease and provide new drug targets and therapies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-197) contains supplementary material, which is available to authorized users.

Association of TERC and OBFC1 haplotypes with mean leukocyte telomere length and risk for coronary heart disease

Objective

To replicate the associations of leukocyte telomere length (LTL) with variants at four loci and to investigate their associations with coronary heart disease (CHD) and type II diabetes (T2D), in order to examine possible causal effects of telomere maintenance machinery on disease aetiology.

Methods

Four SNPs at three loci BICD1 (rs2630578 GγC), 18q12.2 (rs2162440 GγT), and OBFC1 (rs10786775 CγG, rs11591710 AγC) were genotyped in four studies comprised of 2353 subjects out of which 1148 had CHD and 566 T2D. Three SNPs (rs12696304 CγG, rs10936601G>T and rs16847897 GγC) at the TERC locus were genotyped in these four studies, in addition to an offspring study of 765 healthy students. For all samples, LTL had been measured using a real-time PCR-based method.

Results

Only one SNP was associated with a significant effect on LTL, with the minor allele G of OBFC1 rs10786775 SNP being associated with longer LTL (β=0.029, P=0.04). No SNPs were significantly associated with CHD or T2D. For OBFC1 the haplotype carrying both rare alleles (rs10786775G and rs11591710C, haplotype frequency 0.089) was associated with lower CHD prevalence (OR: 0.77; 95% CI: 0.61–0.97; P= 0.03). The TERC haplotype GTC (rs12696304G, rs10936601T and rs16847897C, haplotype frequency 0.210) was associated with lower risk for both CHD (OR: 0.86; 95% CI: 0.75-0.99; P=0.04) and T2D (OR: 0.74; 95% CI: 0.61–0.91; P= 0.004), with no effect on LTL. Only the last association remained after adjusting for multiple testing.

Conclusion

Of reported associations, only that between the OBFC1 rs10786775 SNP and LTL was confirmed, although our study has a limited power to detect modest effects. A 2-SNP OBFC1 haplotype was associated with higher risk of CHD, and a 3-SNP TERC haplotype was associated with both higher risk of CHD and T2D. Further work is required to confirm these results and explore the mechanisms of these effects.

RGS6, but not RGS4, is the dominant regulator of G protein signaling (RGS) modulator of the parasympathetic regulation of mouse heart rate

Parasympathetic activity decreases heart rate (HR) by inhibiting pacemaker cells in the sinoatrial node (SAN). Dysregulation of parasympathetic influence has been linked to sinus node dysfunction and arrhythmia. RGS (regulator of G protein signaling) proteins are negative modulators of the parasympathetic regulation of HR and the prototypical M2 muscarinic receptor (M2R)-dependent signaling pathway in the SAN that involves the muscarinic-gated atrial K(+) channel IKACh. Both RGS4 and RGS6-Gβ5 have been implicated in these processes. Here, we used Rgs4(-/-), Rgs6(-/-), and Rgs4(-/-):Rgs6(-/-) mice to compare the relative influence of RGS4 and RGS6 on parasympathetic regulation of HR and M2R-IKACh-dependent signaling in the SAN. In retrogradely perfused hearts, ablation of RGS6, but not RGS4, correlated with decreased resting HR, increased heart rate variability, and enhanced sensitivity to the negative chronotropic effects of the muscarinic agonist carbachol. Similarly, loss of RGS6, but not RGS4, correlated with enhanced sensitivity of the M2R-IKACh signaling pathway in SAN cells to carbachol and a significant slowing of M2R-IKACh deactivation rate. Surprisingly, concurrent genetic ablation of RGS4 partially rescued some deficits observed in Rgs6(-/-) mice. These findings, together with those from an acute pharmacologic approach in SAN cells from Rgs6(-/-) and Gβ5(-/-) mice, suggest that the partial rescue of phenotypes in Rgs4(-/-):Rgs6(-/-) mice is attributable to another R7 RGS protein whose influence on M2R-IKACh signaling is masked by RGS4. Thus, RGS6-Gβ5, but not RGS4, is the primary RGS modulator of parasympathetic HR regulation and SAN M2R-IKACh signaling in mice.

Genome-wide association studies of maximum number of drinks

Maximum number of drinks (MaxDrinks) defined as "Maximum number of alcoholic drinks consumed in a 24-h period" is an intermediate phenotype that is closely related to alcohol dependence (AD). Family, twin and adoption studies have shown that the heritability of MaxDrinks is approximately 0.5. We conducted the first genome-wide association (GWA) study and meta-analysis of MaxDrinks as a continuous phenotype. 1059 individuals were from the Collaborative Study on the Genetics of Alcoholism (COGA) sample and 1628 individuals were from the Study of Addiction - Genetics and Environment (SAGE) sample. Family sample with 3137 individuals was from the Australian twin-family study of alcohol use disorder (OZALC). Two population-based Caucasian samples (COGA and SAGE) with 1 million single-nucleotide polymorphisms (SNPs) were used for gene discovery and one family-based Caucasian sample was used for replication. Through meta-analysis we identified 162 SNPs associated with MaxDirnks (p < 10(-4)). The most significant association with MaxDrinks was observed with SNP rs11128951 (p = 4.27 × 10(-8)) near SGOL1 gene at 3p24.3. Furthermore, several SNPs (rs17144687 near DTWD2, rs12108602 near NDST4, and rs2128158 in KCNB2) showed significant associations with MaxDrinks (p < 5 × 10(-7)) in the meta-analysis. Especially, 8 SNPs in DDC gene showed significant associations with MaxDrinks (p < 5 × 10(-7)) in the SAGE sample. Several flanking SNPs in above genes/regions were confirmed in the OZALC family sample. In conclusions, we identified several genes/regions associated with MaxDrinks. These findings can improve the understanding about the pathogenesis of alcohol consumption phenotypes and alcohol-related disorders.

BET1L and TNRC6B associate with uterine fibroid risk among European Americans

Uterine fibroid (UFs) affect 77 % of women by menopause and account for $9.4 billion in healthcare costs each year. Although UFs are heritable, genetic risk is poorly understood. The first genome-wide association study (GWAS) of UFs was recently performed in a Japanese population, with reported genome-wide significance for single nucleotide polymorphisms (SNPs) across three chromosomal regions. We tested these SNPs for association with UFs in US cohorts. Women were enrolled in the Right from the Start (RFTS) cohort and the BioVU DNA repository. UF status in both cohorts was determined by pelvic imaging. We tested 65 candidate and haplotype-tagging SNPs for association with UFs presence using logistic regression in RFTS and the top three GWAS-associated SNPs in BioVU. We also combined association results from both cohorts using meta-analysis. 1,086 European American (EA) cases and 1,549 controls were examined. Two SNP associations replicated [blocked early in transport 1 homolog (BET1L) rs2280543, RFTS-BioVU meta-odds ratio (OR) = 0.67 95 % confidence interval (CI) 0.38-0.96, Q = 0.70, I = 0, p = 6.9 × 10⁻³; trinucleotide repeat containing 6B (TNRC6B) rs12484776, RFTS-BioVU meta-OR = 1.21, 95 % CI 1.07-1.35, Q = 0.24, I = 28.37, p = 8.7 × 10⁻³). Meta-analyses combining evidence from RFTS, BioVU, and prior GWAS showed little heterogeneity in effect sizes across studies, with meta-p values between 7.45 × 10⁻⁸ and 3.89 × 10⁻⁹, which were stronger than prior GWAS and supported associations observed for all previously identified loci. These data suggest common variants increase risk for UF in both EA and Japanese populations. However, further research is needed to assess the role of these genes across other racial groups.

Genome-wide association study of cardiac structure and systolic function in African Americans: the Candidate Gene Association Resource (CARe) study

BACKGROUND

Using data from 4 community-based cohorts of African Americans, we tested the association between genome-wide markers (single-nucleotide polymorphisms) and cardiac phenotypes in the Candidate-gene Association Resource study.

METHODS AND RESULTS

Among 6765 African Americans, we related age, sex, height, and weight-adjusted residuals for 9 cardiac phenotypes (assessed by echocardiogram or magnetic resonance imaging) to 2.5 million single-nucleotide polymorphisms genotyped using Genome-wide Affymetrix Human SNP Array 6.0 (Affy6.0) and the remainder imputed. Within the cohort, genome-wide association analysis was conducted, followed by meta-analysis across cohorts using inverse variance weights (genome-wide significance threshold=4.0 ×10(-7)). Supplementary pathway analysis was performed. We attempted replication in 3 smaller cohorts of African ancestry and tested lookups in 1 consortium of European ancestry (EchoGEN). Across the 9 phenotypes, variants in 4 genetic loci reached genome-wide significance: rs4552931 in UBE2V2 (P=1.43×10(-7)) for left ventricular mass, rs7213314 in WIPI1 (P=1.68×10(-7)) for left ventricular internal diastolic diameter, rs1571099 in PPAPDC1A (P=2.57×10(-8)) for interventricular septal wall thickness, and rs9530176 in KLF5 (P=4.02×10(-7)) for ejection fraction. Associated variants were enriched in 3 signaling pathways involved in cardiac remodeling. None of the 4 loci replicated in cohorts of African ancestry was confirmed in lookups in EchoGEN.

CONCLUSIONS

In the largest genome-wide association study of cardiac structure and function to date in African Americans, we identified 4 genetic loci related to left ventricular mass, interventricular septal wall thickness, left ventricular internal diastolic diameter, and ejection fraction, which reached genome-wide significance. Replication results suggest that these loci may be unique to individuals of African ancestry. Additional large-scale studies are warranted for these complex phenotypes.

Endothelial nitric oxide gene haplotype reduces the effect of a single bout of exercise on the vascular reactivity in healthy subjects

Polymorphisms in the endothelial nitric oxide synthase (eNOS) gene reduce shear stress-induced nitric oxide production. Thus, we investigated the individual and combined impact of 3 variants in the eNOS gene (-786T>C, intron 4b4a, and 894G>T) on vascular reactivity before and after exercise. Sedentary, healthy subjects were studied (105 women/26 men, age 32 ± 1 years [mean ± standard error of the mean]). Genotypes were determined by polymerase chain reaction restriction fragment length polymorphism, and haplotypes were determined by a Bayesian-based algorithm. Vascular reactivity was evaluated by the percentage of change in forearm vascular conductance provoked by 5 minutes of circulatory occlusion before (baseline) and 10, 60, and 120 minutes after a maximal cardiopulmonary exercise test. Vascular reactivity increased 10 minutes after exercise in the entire sample (baseline: 218 ± 11% vs 10 minutes: 284 ± 15%, P < 0.001), remained increased at 60 minutes (239 ± 12%, P = 0.02 vs baseline), and returned to baseline at 120 minutes (210 ± 10%, P = 0.83 vs baseline). Genotype analysis showed that subjects with the 894G>T polymorphism had lower vascular reactivity than wild counterparts (group effect, P = 0.05). Furthermore, subjects with haplotype 2 (H2), containing the -786T>C and 894G>T polymorphisms, had lower vascular reactivity than wild counterparts (haplotype 1 [H1]) (group effect, P = 0.05), whereas subjects with haplotype 4 (H4), containing only the 894G>T polymorphism, had vascular reactivity similar to that of wild counterparts (H1) (group effect, P = 0.35). Altogether, these results indicate that the 894G>T polymorphism reduced exercise-mediated increase in vascular reactivity, particularly when it occurred concomitantly with the -786T>C polymorphism.

Genetic epidemiology of left ventricular hypertrophy

Left ventricular (LV) hypertrophy is a strong independent predictor of increased cardiovascular morbidity and mortality in clinical and population-based samples. Clinical and hemodynamic stimuli to LV hypertrophy induce not only an increase in cardiac mass and wall thickness but also a fundamental reconfiguration of the protein, cellular and molecular components of the myocardium. Several studies have indicated that LV mass is influenced by genetic factors. The substantial heritability (h(2)) for LV mass in population-based samples of varying ethnicity indicates robust genetic influences on LV hypertrophy. Genome-wide linkage and association studies in diverse populations have been performed to identify genes influencing LV mass, and although several chromosomal regions have been found to be significantly associated with LV mass, the specific genes and functional variants contained in these chromosomal regions have yet to be identified. In addition, multiple studies have tried to link single-nucleotide polymorphisms (SNPs) in regulatory and pathway genes with common forms of LV hypertrophy, but there is little evidence that these genetic variations are functional. Up to this point in time, the results obtained in genetic studies are of limited clinical value. Much of the heritability remains unexplained, the identity of the underlying gene pathways, genes, and functional variants remains unknown, and the promise of genetically-based risk prediction and personalized medicine remain unfulfilled. However, molecular biological technologies continue to improve rapidly, and the long-term potential of sophisticated genetic investigations using these modern genomic technologies, coupled with smart study designs, remains intact. Ultimately, genetic investigations offer much promise for future prevention, early intervention and treatment of this major public health issue.