HAMAP in 2015: updates to the protein family classification and annotation system.

HAMAP (High-quality Automated and Manual Annotation of Proteins-available at http://hamap.expasy.org/) is a system for the automatic classification and annotation of protein sequences. HAMAP provides annotation of the same quality and detail as UniProtKB/Swiss-Prot, using manually curated profiles for protein sequence family classification and expert curated rules for functional annotation of family members. HAMAP data and tools are made available through our website and as part of the UniRule pipeline of UniProt, providing annotation for millions of unreviewed sequences of UniProtKB/TrEMBL. Here we report on the growth of HAMAP and updates to the HAMAP system since our last report in the NAR Database Issue of 2013. We continue to augment HAMAP with new family profiles and annotation rules as new protein families are characterized and annotated in UniProtKB/Swiss-Prot; the latest version of HAMAP (as of 3 September 2014) contains 1983 family classification profiles and 1998 annotation rules (up from 1780 and 1720). We demonstrate how the complex logic of HAMAP rules allows for precise annotation of individual functional variants within large homologous protein families. We also describe improvements to our web-based tool HAMAP-Scan which simplify the classification and annotation of sequences, and the incorporation of an improved sequence-profile search algorithm.

Genome-wide association study of co-occurring anxiety in major depression.

OBJECTIVES: Co-morbidity between depression and anxiety disorders is common. In this study we define a quantitative measure of anxiety by summating four anxiety items from the SCAN interview in a large collection of major depression (MDD) cases to identify genes contributing to this complex phenotype. METHODS: A total of 1522 MDD cases dichotomised according to those with at least one anxiety item scored (n = 1080) and those without anxiety (n = 442) were analysed, and also compared to 1588 healthy controls at a genome-wide level, to identify genes that may contribute to anxiety in MDD. RESULTS: For the quantitative trait, suggestive evidence of association was detected for two SNPs, and for the dichotomous anxiety present/absent ratings for three SNPs at genome-wide level. In the genome-wide analysis of MDD cases with co-morbid anxiety and healthy controls, two SNPs attained P values of < 5 × 10⁻⁶. Analysing candidate genes, P values ≤ 0.0005 were found with three SNPs for the quantitative trait and three SNPs for the dichotomous trait. CONCLUSIONS: This study provides an initial genome-wide assessment of possible genetic contribution to anxiety in MDD. Although suggestive evidence of association was found for several SNPs, our findings suggest that there are no common variants strongly associated with anxious depression.

Use of green fluorescent protein-based reporters to monitor balanced production of antifungal compounds in the biocontrol agent Pseudomonas fluorescens CHA0.

AIMS: To develop reporter constructs based on stable and unstable variants of the green fluorescent protein (GFP) for monitoring balanced production of antifungal compounds that are crucial for the capacity of the root-colonizing Pseudomonas fluorescens strain CHA0 to control plant diseases caused by soil-borne pathogenic fungi. METHODS AND RESULTS: Pseudomonas fluorescens CHA0 produces the three antifungal metabolites 2,4-diacetylphloroglucinol (DAPG), pyoluteorin (PLT) and pyrrolnitrin (PRN). The gfp[mut3] and gfp[AAV] reporter genes were fused to the promoter regions of the DAPG, PLT and PRN biosynthetic genes. The reporter fusions were then used to follow the kinetics of expression of the three antifungal metabolites in a microplate assay. DAPG and PLT were found to display an inverse relationship in which each metabolite activates its own biosynthesis while repressing the synthesis of the other metabolite. PRN appears not to be involved in this balance. However, the microbial and plant phenolic metabolite salicylate was found to interfere with the expression of both DAPG and PLT. CONCLUSIONS: The results obtained provide evidence that P. fluorescens CHA0 may keep the antifungal compounds DAPG and PLT at a fine-tuned balance that can be affected by certain microbial and plant phenolics. SIGNIFICANCE AND IMPACT OF THE STUDY: To our knowledge, the present study is the first to use stable and unstable GFP variants to study antibiotic gene expression in a biocontrol pseudomonad. The developed reporter fusions will be a highly valuable tool to study in situ expression of this bacterial biocontrol trait on plant roots, i.e. at the site of pathogen suppression.

CUBN is a gene locus for albuminuria.

Identification of genetic risk factors for albuminuria may alter strategies for early prevention of CKD progression, particularly among patients with diabetes. Little is known about the influence of common genetic variants on albuminuria in both general and diabetic populations. We performed a meta-analysis of data from 63,153 individuals of European ancestry with genotype information from genome-wide association studies (CKDGen Consortium) and from a large candidate gene study (CARe Consortium) to identify susceptibility loci for the quantitative trait urinary albumin-to-creatinine ratio (UACR) and the clinical diagnosis microalbuminuria. We identified an association between a missense variant (I2984V) in the CUBN gene, which encodes cubilin, and both UACR (P = 1.1 × 10(-11)) and microalbuminuria (P = 0.001). We observed similar associations among 6981 African Americans in the CARe Consortium. The associations between this variant and both UACR and microalbuminuria were significant in individuals of European ancestry regardless of diabetes status. Finally, this variant associated with a 41% increased risk for the development of persistent microalbuminuria during 20 years of follow-up among 1304 participants with type 1 diabetes in the prospective DCCT/EDIC Study. In summary, we identified a missense CUBN variant that associates with levels of albuminuria in both the general population and in individuals with diabetes.

No evidence for a causal link between uric acid and type 2 diabetes: a Mendelian randomisation approach.

AIMS/HYPOTHESIS: Epidemiological and experimental evidence suggests that uric acid has a role in the aetiology of type 2 diabetes. Using a Mendelian randomisation approach, we investigated whether there is evidence for a causal role of serum uric acid for development of type 2 diabetes. METHODS: We examined the associations of serum-uric-acid-raising alleles of eight common variants recently identified in genome-wide association studies and summarised this in a genetic score with type 2 diabetes in case-control studies including 7,504 diabetes patients and 8,560 non-diabetic controls. We compared the observed effect size to that expected based on: (1) the association between the genetic score and uric acid levels in non-diabetic controls; and (2) the meta-analysed uric acid level to diabetes association. RESULTS: The genetic score showed a linear association with uric acid levels, with a difference of 12.2 μmol/l (95% CI 9.3, 15.1) by score tertile. No significant associations were observed between the genetic score and potential confounders. No association was observed between the genetic score and type 2 diabetes with an OR of 0.99 (95% CI 0.94, 1.04) per score tertile, significantly different (p = 0.046) from that expected (1.04 [95% CI 1.03, 1.05]) based on the observed uric acid difference by score tertile and the uric acid to diabetes association of 1.21 (95% CI 1.14, 1.29) per 60 μmol/l. CONCLUSIONS/INTERPRETATION: Our results do not support a causal role of serum uric acid for the development of type 2 diabetes and limit the expectation that uric-acid-lowering drugs will be effective in the prevention of type 2 diabetes.

Genetik der essentiellen arteriellen Hypertonie [Essential hypertension - which genes can be held responsible for it?].

Hypertension is a common heritable cardiovascular risk factor. Some rare monogenic forms of hypertension have been described, but the majority of patients suffer from essential hypertension, for whom the underlying genetic mechanisms are not clear. Essential hypertension is a complex trait, involving multiple genes and environmental factors. Recently, progress in the identification of common genetic variants associated with essential hypertension has been made due to large-scale international collaborative projects. In this article we review the new research methods used as well as selected recent findings in this field.

Antiretroviral activity of ancestral TRIM5alpha.

The antiretroviral protein TRIM5alpha is known to have evolved different restriction capacities against various retroviruses, driven by positive Darwinian selection. However, how these different specificities have evolved in the primate lineages is not fully understood. Here we used ancestral protein resurrection to estimate the evolution of antiviral restriction specificities of TRIM5alpha on the primate lineage leading to humans. We used TRIM5alpha coding sequences from 24 primates for the reconstruction of ancestral TRIM5alpha sequences using maximum-likelihood and Bayesian approaches. Ancestral sequences were transduced into HeLa and CRFK cells. Stable cell lines were generated and used to test restriction of a panel of extant retroviruses (human immunodeficiency virus type 1 [HIV-1] and HIV-2, simian immunodeficiency virus [SIV] variants SIV(mac) and SIV(agm), and murine leukemia virus [MLV] variants N-MLV and B-MLV). The resurrected TRIM5alpha variant from the common ancestor of Old World primates (Old World monkeys and apes, approximately 25 million years before present) was effective against present day HIV-1. In contrast to the HIV-1 restriction pattern, we show that the restriction efficacy against other retroviruses, such as a murine oncoretrovirus (N-MLV), is higher for more recent resurrected hominoid variants. Ancestral TRIM5alpha variants have generally limited efficacy against HIV-2, SIV(agm), and SIV(mac). Our study sheds new light on the evolution of the intrinsic antiviral defense machinery and illustrates the utility of functional evolutionary reconstruction for characterizing recently emerged protein differences.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Covariation between plumage colour polymorphism and diet in the Barn Owl Tyto alba

Several hypotheses might explain the evolution and maintenance of colour morphs within animal populations. The 'alternative foraging strategy' hypothesis states that alternative colour morphs exploit different ecological niches. This hypothesis predicts that morphs differ in diet, either because foraging success on alternative prey species is morph-dependent or because differently coloured individuals exploit alternative habitats. I examined this prediction in the Barn Owl Tyto alba, a bird that varies in plumage coloration continuously from dark reddish-brown to white. On the European continent, Owls are light-coloured (subspecies T. a. alba) in the south and reddish-brown (T. a. guttata) in the north; in central Europe the two subspecies interbreed, generating many colour variants. If plumage coloration indicates alternative foraging strategies, in sympatry dark- and light-coloured owls should consume prey species that are typical of the diets of T. a. guttata and T. a. alba in allopatry, respectively. In line with this prediction, both in allopatry and in sympatry in Switzerland T. a. guttata fed primarily upon Common Voles Microtus arvalis and T. a. alba upon Wood Mice Apodemus spp. Statistical analyses suggest that morph-dependent diet did not arise from a non-random habitat distribution of owls with respect to plumage coloration. This suggests that foraging success upon alternative prey is morph-dependent.

Identification of an RP1 Prevalent Founder Mutation and Related Phenotype in Spanish Patients with Early-Onset Autosomal Recessive Retinitis.

OBJECTIVE: To identify the genetic causes underlying early-onset autosomal recessive retinitis pigmentosa (arRP) in the Spanish population and describe the associated phenotype. DESIGN: Case series. PARTICIPANTS: A total of 244 unrelated families affected by early-onset arRP. METHODS: Homozygosity mapping or exome sequencing analysis was performed in 3 families segregating arRP. A mutational screening was performed in 241 additional unrelated families for the p.Ser452Stop mutation. Haplotype analysis also was conducted. Individuals who were homozygotes, double heterozygotes, or carriers of mutations in RP1 underwent an ophthalmic evaluation to establish a genotype-phenotype correlation. MAIN OUTCOME MEASURES: DNA sequence variants, homozygous regions, haplotypes, best-corrected visual acuity, visual field assessments, electroretinogram responses, and optical coherence tomography images. RESULTS: Four novel mutations in RP1 were identified. The new mutation p.Ser542Stop was present in 11 of 244 (4.5%) of the studied families. All chromosomes harboring this mutation shared the same haplotype. All patients presented a common phenotype with an early age of onset and a prompt macular degeneration, whereas the heterozygote carriers did not show any signs of retinitis pigmentosa (RP). CONCLUSIONS: p.Ser542Stop is a single founder mutation and the most prevalent described mutation in the Spanish population. It causes early-onset RP with a rapid macular degeneration and is responsible for 4.5% of all cases. Our data suggest that the implication of RP1 in arRP may be underestimated. FINANCIAL DISCLOSURE(S): The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Model structure of human APOBEC3G.

BACKGROUND: APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G) has antiretroviral activity associated with the hypermutation of viral DNA through cytosine deamination. APOBEC3G has two cytosine deaminase (CDA) domains; the catalytically inactive amino-terminal domain of APOBEC3G (N-CDA) carries the Vif interaction domain. There is no 3-D structure of APOBEC3G solved by X-ray or nuclear magnetic resonance. METHODOLOGY/PRINCIPAL FINDINGS: We predicted the structure of human APOBEC3G based on the crystal structure of APOBEC2. To assess the model structure, we evaluated 48 mutants of APOBEC3G N-CDA that identify novel variants altering DeltaVif HIV-1 infectivity and packaging of APOBEC3G. Results indicated that the key residue D128 is exposed at the surface of the model, with a negative local electrostatic potential. Mutation D128K changes the sign of that local potential. In addition, two novel functionally relevant residues that result in defective APOBEC3G encapsidation, R122 and W127, cluster at the surface. CONCLUSIONS/SIGNIFICANCE: The structure model identifies a cluster of residues important for packaging of APOBEC3G into virions, and may serve to guide functional analysis of APOBEC3G.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.

Genetic testing in patients with obesity.

The obesity epidemic is associated with the recent availability of highly palatable and inexpensive caloric food as well as important changes in lifestyle. Genetic factors, however, play a key role in regulating energy balance and numerous twin studies have estimated the BMI heritability between 40 and 70%. While common variants, identified through genome-wide association studies (GWAS) point toward new pathways, their effect size are too low to be of any use in the clinic. This review therefore concentrates on genes and genomic regions associated with very high risks of human obesity. Although there are no consensus guidelines, we review how the knowledge on these "causal factors" can be translated into the clinic for diagnostic purposes. We propose genetic workups guided by clinical manifestations in patients with severe early-onset obesity. While etiological diagnoses are unequivocal in a minority of patients, new genomic tools such as Comparative Genomic Hybridization (CGH) array, have allowed the identification of novel "causal" loci and next-generation sequencing brings the promise of accelerated pace for discoveries relevant to clinical practice.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa.

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes. Copyright © 2011 John Wiley & Sons, Ltd and Eating Disorders Association.

In vivo analysis of efavirenz metabolism in individuals with impaired CYP2A6 function.

INTRODUCTION: The antiretroviral drug efavirenz (EFV) is extensively metabolized into three primary metabolites: 8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV. There is a wide interindividual variability in EFV plasma exposure, explained to a great extent by cytochrome P450 2B6 (CYP2B6), the main isoenzyme responsible for EFV metabolism and involved in the major metabolic pathway (8-hydroxylation) and to a lesser extent in 7-hydroxylation. When CYP2B6 function is impaired, the relevance of CYP2A6, the main isoenzyme responsible for 7-hydroxylation may increase. We hypothesize that genetic variability in this gene may contribute to the particularly high, unexplained variability in EFV exposure in individuals with limited CYP2B6 function. METHODS: This study characterized CYP2A6 variation (14 alleles) in individuals (N=169) previously characterized for functional variants in CYP2B6 (18 alleles). Plasma concentrations of EFV and its primary metabolites (8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV) were measured in different genetic backgrounds in vivo. RESULTS: The accessory metabolic pathway CYP2A6 has a critical role in limiting drug accumulation in individuals characterized as CYP2B6 slow metabolizers. CONCLUSION: Dual CYP2B6 and CYP2A6 slow metabolism occurs at significant frequency in various human populations, leading to extremely high EFV exposure.