Rare coding variation in paraoxonase-1 is associated with ischemic stroke in the NHLBI Exome Sequencing Project

HDL-associated paraoxonase-1 (PON1) is an enzyme whose activity is associated with cerebrovascular disease. Common PON1 genetic variants have not been consistently associated with cerebrovascular disease. Rare coding variation that likely alters PON1 enzyme function may be more strongly associated with stroke. The National Heart, Lung, and Blood Institute Exome Sequencing Project sequenced the coding regions (exomes) of the genome for heart, lung, and blood-related phenotypes (including ischemic stroke). In this sample of 4,204 unrelated participants, 496 had verified, noncardioembolic ischemic stroke. After filtering, 28 nonsynonymous PON1 variants were identified. Analysis with the sequence kernel association test, adjusted for covariates, identified significant associations between PON1 variants and ischemic stroke (P = 3.01 × 10⁻³). Stratified analyses demonstrated a stronger association of PON1 variants with ischemic stroke in African ancestry (AA) participants (P = 5.03 × 10⁻³). Ethnic differences in the association between PON1 variants with stroke could be due to the effects of PON1_Val109Ile (overall P = 7.88 × 10⁻³; AA P = 6.52 × 10⁻⁴), found at higher frequency in AA participants (1.16% vs. 0.02%) and whose protein is less stable than the common allele. In summary, rare genetic variation in PON1 was associated with ischemic stroke, with stronger associations identified in those of AA. Increased focus on PON1 enzyme function and its role in cerebrovascular disease is warranted.     

Convergent evolution, habitat shifts and variable diversification rates in the ovenbird-woodcreeper family (Furnariidae)

Background  

The Neotropical ovenbird-woodcreeper family (Furnariidae) is an avian group characterized by exceptionally diverse ecomorphological adaptations. For instance, members of the family are known to construct nests of a remarkable variety. This offers a unique opportunity to examine whether changes in nest design, accompanied by expansions into new habitats, facilitates diversification. We present a multi-gene phylogeny and age estimates for the ovenbird-woodcreeper family and use these results to estimate the degree of convergent evolution in both phenotype and habitat utilisation. Furthermore, we discuss whether variation in species richness among ovenbird clades could be explained by differences in clade-specific diversification rates, and whether these rates differ among lineages with different nesting habits. In addition, the systematic positions of some enigmatic ovenbird taxa and the postulated monophyly of some species-rich genera are evaluated.     

TagSNP analyses of the PON gene cluster: effects on PON1 activity, LDL oxidative susceptibility, and vascular disease

Paraoxonase 1 (PON1) activity is consistently predictive of vascular disease, although the genotype at four functional PON1 polymorphisms is not. To address this inconsistency, we investigated the role of all common PON1 genetic variability, as measured by tagging single-nucleotide polymorphisms (tagSNPs), in predicting PON1 activity for phenylacetate hydrolysis, LDL susceptibility to oxidation ex vivo, plasma homocysteine (Hcy) levels, and carotid artery disease (CAAD) status. The biological goal was to establish whether additional common genetic variation beyond consideration of the four known functional SNPs improves prediction of these phenotypes. PON2 and PON3 tagSNPs were secondarily evaluated. Expanded analysis of an additional 26 tagSNPs found evidence of previously undescribed common PON1 polymorphisms that affect PON1 activity independently of the four known functional SNPs. PON1 activity was not significantly correlated with LDL oxidative susceptibility, but genotypes at the PON1(-108) promoter polymorphism and several other PON1 SNPs were. Neither PON1 activity nor PON1 genotype was significantly correlated with plasma Hcy levels. This study revealed previously undetected common functional PON1 polymorphisms that explain 4% of PON1 activity and a high rate of recombination in PON1, but the sum of the common PON1 locus variation does not explain the relationship between PON1 activity and CAAD.     

Novel common and rare genetic determinants of paraoxonase activity: FTO,SERPINA12, and ITGAL

HDL-associated paraoxonase 1 (PON1) activity is associated with cardiovascular and other human diseases. As the role of genetic variants outside of the PON gene cluster on PON1 activity is unknown, we sought to identify common and rare variants in such loci. We typed 33,057 variants on the CVD chip in 1,362 subjects to test for their effects on adjusted-PON1 activity. Three novel genes (FTO, ITGAL, and SERPINA12) and the PON gene cluster had SNPs associated with PON1 arylesterase (AREase) activity. These loci were carried forward for rare-variant analysis using Exome chip genotypes in an overlapping subset of 1,051 subjects using sequence kernel association testing. PON1 (P = 2.24 × 10⁻⁴), PON3 (P = 0.022), FTO (P = 0.019), and SERPINA12 (P = 0.039) had both common and rare variants associated with PON1 AREase. ITGAL variants were associated with PON1 activity when using weighted sequence kernel association testing (SKAT) analysis (P = 2.63 × 10⁻³). When adjusting for the initial common variants, SERPINA12 became marginally significant (P = 0.09), whereas all other findings remained significant (P < 0.05), suggesting independent rare-variant effects. We present novel findings that common and rare variants in FTO, SERPINA12, and ITGAL predict PON1 activity. These results further link PON1 to diabetes and inflammation and may inform the role of HDL in human disease.     

Joint Linkage and Association Analysis with Exome Sequence Data Implicates SLC25A40 in Hypertriglyceridemia

Hypertriglyceridemia (HTG) is a heritable risk factor for cardiovascular disease. Investigating the genetics of HTG may identify new drug targets. There are ∼35 known single-nucleotide variants (SNVs) that explain only ∼10% of variation in triglyceride (TG) level. Because of the genetic heterogeneity of HTG, a family study design is optimal for identification of rare genetic variants with large effect size because the same mutation can be observed in many relatives and cosegregation with TG can be tested. We considered HTG in a five-generation family of European American descent (n = 121), ascertained for familial combined hyperlipidemia. By using Bayesian Markov chain Monte Carlo joint oligogenic linkage and association analysis, we detected linkage to chromosomes 7 and 17. Whole-exome sequence data revealed shared, highly conserved, private missense SNVs in both SLC25A40 on chr7 and PLD2 on chr17. Jointly, these SNVs explained 49% of the genetic variance in TG; however, only the SLC25A40 SNV was significantly associated with TG (p = 0.0001). This SNV, c.374A>G, causes a highly disruptive p.Tyr125Cys substitution just outside the second helical transmembrane region of the SLC25A40 inner mitochondrial membrane transport protein. Whole-gene testing in subjects from the Exome Sequencing Project confirmed the association between TG and SLC25A40 rare, highly conserved, coding variants (p = 0.03). These results suggest a previously undescribed pathway for HTG and illustrate the power of large pedigrees in the search for rare, causal variants.     

Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies

We repurposed existing genotypes in DNA biobanks across the Electronic Medical Records and Genomics network to perform a genome-wide association study for primary hypothyroidism, the most common thyroid disease. Electronic selection algorithms incorporating billing codes, laboratory values, text queries, and medication records identified 1317 cases and 5053 controls of European ancestry within five electronic medical records (EMRs); the algorithms'' positive predictive values were 92.4% and 98.5% for cases and controls, respectively. Four single-nucleotide polymorphisms (SNPs) in linkage disequilibrium at 9q22 near FOXE1 were associated with hypothyroidism at genome-wide significance, the strongest being rs7850258 (odds ratio [OR] 0.74, p = 3.96 × 10⁻⁹). This association was replicated in a set of 263 cases and 1616 controls (OR = 0.60, p = 5.7 × 10⁻⁶). A phenome-wide association study (PheWAS) that was performed on this locus with 13,617 individuals and more than 200,000 patient-years of billing data identified associations with additional phenotypes: thyroiditis (OR = 0.58, p = 1.4 × 10⁻⁵), nodular (OR = 0.76, p = 3.1 × 10⁻⁵) and multinodular (OR = 0.69, p = 3.9 × 10⁻⁵) goiters, and thyrotoxicosis (OR = 0.76, p = 1.5 × 10⁻³), but not Graves disease (OR = 1.03, p = 0.82). Thyroid cancer, previously associated with this locus, was not significantly associated in the PheWAS (OR = 1.29, p = 0.09). The strongest association in the PheWAS was hypothyroidism (OR = 0.76, p = 2.7 × 10⁻¹³), which had an odds ratio that was nearly identical to that of the curated case-control population in the primary analysis, providing further validation of the PheWAS method. Our findings indicate that EMR-linked genomic data could allow discovery of genes associated with many diseases without additional genotyping cost.     

Common ancestry of iron oxide- and iron-sulfide-based biomineralization in magnetotactic bacteria

Magnetosomes are prokaryotic organelles produced by magnetotactic bacteria that consist of nanometer-sized magnetite (Fe₃O₄) or/and greigite (Fe₃S₄) magnetic crystals enveloped by a lipid bilayer membrane. In magnetite-producing magnetotactic bacteria, proteins present in the magnetosome membrane modulate biomineralization of the magnetite crystal. In these microorganisms, genes that encode for magnetosome membrane proteins as well as genes involved in the construction of the magnetite magnetosome chain, the mam and mms genes, are organized within a genomic island. However, partially because there are presently no greigite-producing magnetotactic bacteria in pure culture, little is known regarding the greigite biomineralization process in these organisms including whether similar genes are involved in the process. Here using culture-independent techniques, we now show that mam genes involved in the production of magnetite magnetosomes are also present in greigite-producing magnetotactic bacteria. This finding suggest that the biomineralization of magnetite and greigite did not have evolve independently (that is, magnetotaxis is polyphyletic) as once suggested. Instead, results presented here are consistent with a model in which the ability to biomineralize magnetosomes and the possession of the mam genes was acquired by bacteria from a common ancestor, that is, the magnetotactic trait is monophyletic.     

Charged and hydrophobic surfaces on the a chain of shiga-like toxin 1 recognize the C-terminal domain of ribosomal stalk proteins

Shiga-like toxins are ribosome-inactivating proteins (RIP) produced by pathogenic E. coli strains that are responsible for hemorrhagic colitis and hemolytic uremic syndrome. The catalytic A(1) chain of Shiga-like toxin 1 (SLT-1), a representative RIP, first docks onto a conserved peptide SD[D/E]DMGFGLFD located at the C-terminus of all three eukaryotic ribosomal stalk proteins and halts protein synthesis through the depurination of an adenine base in the sarcin-ricin loop of 28S rRNA. Here, we report that the A(1) chain of SLT-1 rapidly binds to and dissociates from the C-terminal peptide with a monomeric dissociation constant of 13 μM. An alanine scan performed on the conserved peptide revealed that the SLT-1 A(1) chain interacts with the anionic tripeptide DDD and the hydrophobic tetrapeptide motif FGLF within its sequence. Based on these 2 peptide motifs, SLT-1 A(1) variants were generated that displayed decreased affinities for the stalk protein C-terminus and also correlated with reduced ribosome-inactivating activities in relation to the wild-type A(1) chain. The toxin-peptide interaction and subsequent toxicity were shown to be mediated by cationic and hydrophobic docking surfaces on the SLT-1 catalytic domain. These docking surfaces are located on the opposite face of the catalytic cleft and suggest that the docking of the A(1) chain to SDDDMGFGLFD may reorient its catalytic domain to face its RNA substrate. More importantly, both the delineated A(1) chain ribosomal docking surfaces and the ribosomal peptide itself represent a target and a scaffold, respectively, for the design of generic inhibitors to block the action of RIPs.