Activity of ancestral restriction factors against ancient retroviruses

Analysis of TRIM5α and APOBEC3G genes suggests that these two restriction factors underwent strong positive selection throughout primate evolution. This pressure was possibly imposed by ancient exogenous retroviruses, of which endogenous retroviruses are remnants. Our study aims to assess in vitro the activity of these factors against ancient retroviruses by reconstructing their ancestral gag sequences, as well as the ancestral TRIM5α and APOBEC3G for primates. Based on evolutionary genomics approach, we reconstructed ancestors of the two largest families of human endogenous retroviruses (HERV), namely HERV-K and HERV-H, as well as primate ancestral TRIM5α and APOBEC3G variants. The oldest TRIM5α sequence was the catarhinne TRIM5α, common ancestor of Old World monkeys and hominoids, dated from 25 million years ago (mya). From the oldest, to the youngest, ancestral TRIM5α variants showed less restriction of HIV-1 in vitro [1]. Likewise three ancestral APOBEC3Gs sequences common to hominoids (18 mya), Old World monkeys, and catarhinnes (25 mya) were reconstructed. All ancestral APOBEC3G variants inhibited efficiently HIV-1Δvif in vitro, compared to modern APOBEC3Gs. The ability of Vif proteins (HIV-1, HIV-2, SIVmac and SIVagm) to counteract their activity tallied with the residue 128 on ancestral APOBEC3Gs. Moreover we are attempting to reconstruct older ancestral sequences of both restriction factors by using prosimian orthologue sequences. An infectious onemillion- years-old HERV-KCON previously reconstituted was shown to be resistant to modern TRIM5α and APOBEC3G [2]. Our ancestral TRIM5α and APOBEC3G variants were inactive against HERV-KCON. Besides we reconstructed chimeric HERV-K bearing ancestral capsids (up to 7 mya) that resulted in infectious viruses resistant to modern and ancestral TRIM5α. Likewise HERV-K viruses bearing ancestral nucleocapsids will be tested for ancestral and modern APOBEC3G restriction. In silico reconstruction and structural modeling of ancestral HERV-H capsids resulted in structures homologous to that of the gammaretrovirus MLV. Thus we are attempting to construct chimeric MLV virus bearing HERV-H ancestral capsids. These chimeric ancestral HERVs will be tested for infectivity and restriction by ancestral TRIM5α. Similarly chimeric MLV viruses bearing ancestral HERV-H nucleocapsids will be reconstructed and tested for APOBEC3G restriction.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy.

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

Serum uric acid and adiposity: deciphering causality using a bidirectional Mendelian randomization approach.

Although the relationship between serum uric acid (SUA) and adiposity is well established, the direction of the causality is still unclear in the presence of conflicting evidences. We used a bidirectional Mendelian randomization approach to explore the nature and direction of causality between SUA and adiposity in a population-based study of Caucasians aged 35 to 75 years. We used, as instrumental variables, rs6855911 within the SUA gene SLC2A9 in one direction, and combinations of SNPs within the adiposity genes FTO, MC4R and TMEM18 in the other direction. Adiposity markers included weight, body mass index, waist circumference and fat mass. We applied a two-stage least squares regression: a regression of SUA/adiposity markers on our instruments in the first stage and a regression of the response of interest on the fitted values from the first stage regression in the second stage. SUA explained by the SLC2A9 instrument was not associated to fat mass (regression coefficient [95% confidence interval]: 0.05 [-0.10, 0.19] for fat mass) contrasting with the ordinary least square estimate (0.37 [0.34, 0.40]). By contrast, fat mass explained by genetic variants of the FTO, MC4R and TMEM18 genes was positively and significantly associated to SUA (0.31 [0.01, 0.62]), similar to the ordinary least square estimate (0.27 [0.25, 0.29]). Results were similar for the other adiposity markers. Using a bidirectional Mendelian randomization approach in adult Caucasians, our findings suggest that elevated SUA is a consequence rather than a cause of adiposity.

The Growing Importance of CNVs: New Insights for Detection and Clinical Interpretation.

Differences between genomes can be due to single nucleotide variants, translocations, inversions, and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 500 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease. Hence there is a need for better-tailored and more robust tools for the detection and genome-wide analyses of CNVs. While a link between a given CNV and a disease may have often been established, the relative CNV contribution to disease progression and impact on drug response is not necessarily understood. In this review we discuss the progress, challenges, and limitations that occur at different stages of CNV analysis from the detection (using DNA microarrays and next-generation sequencing) and identification of recurrent CNVs to the association with phenotypes. We emphasize the importance of germline CNVs and propose strategies to aid clinicians to better interpret structural variations and assess their clinical implications.

The shared allelic architecture of adiponectin levels and coronary artery disease.

OBJECTIVE: A large body of epidemiologic data strongly suggests an association between excess adiposity and coronary artery disease (CAD). Low adiponectin levels, a hormone secreted only from adipocytes, have been associated with an increased risk of CAD in observational studies. However, these associations cannot clarify whether this relationship is causal or due to a shared set of causal factors or even confounding. Genome-wide association studies have identified common variants that influence adiponectin levels, providing valuable tools to examine the genetic relationship between adiponectin and CAD. METHODS: Using 145 genome wide significant SNPs for adiponectin from the ADIPOGen consortium (n = 49,891), we tested whether adiponectin-decreasing alleles influenced risk of CAD in the CARDIoGRAM consortium (n = 85,274). RESULTS: In single-SNP analysis, 5 variants among 145 SNPs were associated with increased risk of CAD after correcting for multiple testing (P < 4.4 × 10(-4)). Using a multi-SNP genotypic risk score to test whether adiponectin levels and CAD have a shared genetic etiology, we found that adiponectin-decreasing alleles increased risk of CAD (P = 5.4 × 10(-7)). CONCLUSION: These findings demonstrate that adiponectin levels and CAD have a shared allelic architecture and provide rationale to undertake a Mendelian randomization studies to understand if this relationship is causal.

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Improving binding affinity and stability of peptide ligands by substituting glycines with D-amino acids.

Improving the binding affinity and/or stability of peptide ligands often requires testing of large numbers of variants to identify beneficial mutations. Herein we propose a type of mutation that promises a high success rate. In a bicyclic peptide inhibitor of the cancer-related protease urokinase-type plasminogen activator (uPA), we observed a glycine residue that has a positive ϕ dihedral angle when bound to the target. We hypothesized that replacing it with a D-amino acid, which favors positive ϕ angles, could enhance the binding affinity and/or proteolytic resistance. Mutation of this specific glycine to D-serine in the bicyclic peptide indeed improved inhibitory activity (1.75-fold) and stability (fourfold). X-ray-structure analysis of the inhibitors in complex with uPA showed that the peptide backbone conformation was conserved. Analysis of known cyclic peptide ligands showed that glycine is one of the most frequent amino acids, and that glycines with positive ϕ angles are found in many protein-bound peptides. These results suggest that the glycine-to-D-amino acid mutagenesis strategy could be broadly applied.

Splice variant-specific stabilization of JNKs by IB1/JIP1.

Islet-Brain 1 (IB1) (also called JNK-interacting protein 1; JIP1) is a scaffold protein that tethers components of the JNK mitogen-activated protein kinase pathway inducing a modulation of the activity and the target specificity of the JNK kinases. Dysfunctions in IB1 have been associated with diseases such as early type II diabetes. To gain more insight in the functions of IB1, its ability to modulate the expression levels of the various JNK proteins was assessed. Each of the three JNK genes gives rise to several splice variants encoding short or long proteins. The expression levels of the short JNK proteins, but not of the long variants, were systematically higher in rat tissues and in transformed cell lines expressing high IB1 levels compared to tissues and cells with no or low IB1 expression. HEK293 cells bearing a tetracycline-inducible IB1 construct showed a specific increase of the short JNK endogenous splice variants in the presence of tetracycline. The augmented expression level of the short JNK splice variants induced by IB1 resulted from an increased stability towards degradation. Modulation of the stability of specific JNK splice variants represents therefore a newly identified mechanism used by IB1 to regulate the JNK MAPK pathway.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

Publication trends of shared decision making in 15 high impact medical journals: a full-text review with bibliometric analysis.

BACKGROUND: Shared Decision Making (SDM) is increasingly advocated as a model for medical decision making. However, there is still low use of SDM in clinical practice. High impact factor journals might represent an efficient way for its dissemination. We aimed to identify and characterize publication trends of SDM in 15 high impact medical journals. METHODS: We selected the 15 general and internal medicine journals with the highest impact factor publishing original articles, letters and editorials. We retrieved publications from 1996 to 2011 through the full-text search function on each journal website and abstracted bibliometric data. We included publications of any type containing the phrase "shared decision making" or five other variants in their abstract or full text. These were referred to as SDM publications. A polynomial Poisson regression model with logarithmic link function was used to assess the evolution across the period of the number of SDM publications according to publication characteristics. RESULTS: We identified 1285 SDM publications out of 229,179 publications in 15 journals from 1996 to 2011. The absolute number of SDM publications by journal ranged from 2 to 273 over 16 years. SDM publications increased both in absolute and relative numbers per year, from 46 (0.32% relative to all publications from the 15 journals) in 1996 to 165 (1.17%) in 2011. This growth was exponential (P < 0.01). We found fewer research publications (465, 36.2% of all SDM publications) than non-research publications, which included non-systematic reviews, letters, and editorials. The increase of research publications across time was linear. Full-text search retrieved ten times more SDM publications than a similar PubMed search (1285 vs. 119 respectively). CONCLUSION: This review in full-text showed that SDM publications increased exponentially in major medical journals from 1996 to 2011. This growth might reflect an increased dissemination of the SDM concept to the medical community.

The balanced lethal system of crested newts: a ghost of sex chromosomes past?

Balanced lethal systems are more than biological curiosities: as theory predicts, they should quickly be eliminated through the joint forces of recombination and selection. That such systems might become fixed in natural populations poses a challenge to evolutionary theory. Here we address the case of a balanced lethal system fixed in crested newts and related species, which makes 50% of offspring die early in development. All adults are heteromorphic for chromosome pair 1. The two homologues (1A and 1B) have different recessive deleterious alleles fixed on a nonrecombining segment, so that heterozygotes are viable, while homozygotes are lethal. Given such a strong segregation load, how could autosomes stop recombining? We propose a role for a sex-chromosome turnover from pair 1 (putative ancestral sex chromosome) to pair 4 (currently active sex chromosome). Accordingly, 1A and 1B represent two variants (Y(A) and Y(B)) of the Y chromosome from an ancestral male-heterogametic system. We formalize a scenario in which turnovers are driven by sex ratio selection stemming from gene-environment interactions on sex determination. Individual-based simulations show that a balanced lethal system can be fixed with significant likelihood, provided the masculinizing allele on chromosome 4 appears after the elimination of the feminizing allele on chromosome 1. Our study illustrates how strikingly maladaptive traits might evolve through natural selection.

Genetic overlap in kallmann syndrome, combined pituitary hormone deficiency, and septo-optic dysplasia.

Context: Kallmann syndrome (KS), combined pituitary hormone deficiency (CPHD), and septo-optic dysplasia (SOD) all result from development defects of the anterior midline in the human forebrain. Objective: The objective of the study was to investigate whether KS, CPHD, and SOD have shared genetic origins. Design and Participants: A total of 103 patients with either CPHD (n = 35) or SOD (n = 68) were investigated for mutations in genes implicated in the etiology of KS (FGFR1, FGF8, PROKR2, PROK2, and KAL1). Consequences of identified FGFR1, FGF8, and PROKR2 mutations were investigated in vitro. Results: Three patients with SOD had heterozygous mutations in FGFR1; these were either shown to alter receptor signaling (p.S450F, p.P483S) or predicted to affect splicing (c.336C>T, p.T112T). One patient had a synonymous change in FGF8 (c.216G>A, p.T72T) that was shown to affect splicing and ligand signaling activity. Four patients with CPHD/SOD were found to harbor heterozygous rare loss-of-function variants in PROKR2 (p.R85G, p.R85H, p.R268C). Conclusions: Mutations in FGFR1/FGF8/PROKR2 contributed to 7.8% of our patients with CPHD/SOD. These data suggest a significant genetic overlap between conditions affecting the development of anterior midline in the human forebrain.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

Contribution de l'étude CoLaus à l'élucidation des déterminants de l'acide urique sérique [Contribution of the CoLaus study to decipher the determinants of serum uric acid].

Asymptomatic hyperuricemia affects one in five adults in the general population and is associated with elevated cardiovascular risk. It is however not clear whether asymptomatic hyperuricemia is a cause or simply a marker of conditions associated with high cardiovascular risk. Sex, age, obesity, renal function and selected drugs are major determinants of serum uric acid. Moreover, recent genome-wide association studies have identified new genes involved in the control of serum uric acid levels, in particular SLC2A9, which encodes a urate transporter located in the kidney. A genetic score based on several genetic variants associated with serum uric acid is strongly associated with the risk of gout, but not with cardiovascular events so far.

Analysis of the cytolytic T lymphocyte response of melanoma patients to the naturally HLA-A*0201-associated tyrosinase peptide 368-376.

The human tyrosinase gene codes for two distinct antigens that are recognized by HLA-A*0201-restricted CTLs. For one of them, tyrosinase peptide 368-376, the sequence identified by mass spectrometry in melanoma cell eluates differs from the gene-encoded sequence as a result of posttranslational modification of amino acid residue 370 (asparagine to aspartic acid). Here, we used fluorescent tetrameric complexes ("tetramers") of HLA-A*0201 and tyrosinase peptide 368-376 (YMDGTMSQV) to characterize the CD8+ T-cell response to this antigen in lymphoid cell populations from HLA-A2 melanoma patients. Taking advantage of the presence of significant numbers of tetramer-positive CD8+ T cells in tumor-infiltrated lymph node cells from a melanoma patient, we derived polyclonal and monoclonal tyrosinase peptide 368-376-specific CTLs by tetramer-guided flow cytometric sorting. These CTLs efficiently and specifically lysed HLA-A*0201- and tyrosinase-positive melanoma cells. As assessed with tyrosinase peptide variants, the fine antigen specificity of the CTLs was quite diverse at the clonal level. Flow cytometric analysis of PBMCs stained with tetramers showed that tyrosinase peptide 368-376-specific CD8+ T cells were hardly detectable in peripheral blood of melanoma patients. However, significant numbers of such cells were detected after short-term stimulation of CD8+ lymphocytes with tyrosinase peptide 368-376 in 6 of 10 HLA-A2 melanoma patients. Taken together, these findings emphasize the significant contribution of the natural tyrosinase peptide 368-376 to the antigenic specificities recognized by the tumor-reactive CTLs that may develop in HLA-A2 melanoma patients.