Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Coevolution of restriction factors and retroviruses

Les virus exploitent la machinerie cellulaire de l'hôte pour se répliquer. Ils doivent s'adapter pour infecter la cellule hôte de manière optimale tout en échappant à la vigilance du système de défense de l'hôte. Ainsi l'hôte et les virus se livrent à de constantes batailles évolutives. Mon travail de thèse a porté sur l'étude des signatures évolutives de facteurs de l'hôte agissant comme des 'facteurs de restriction' en bloquant la réplication rétrovirale chez les primates. Plus spécifiquement, mon travail a visé à utiliser des données évolutives pour renseigner les analyses fonctionnelles et la biologie. Nous avons étudié le facteur anti-VIH-1 nommé TRIM5a (i) chez les prosimiens pour mieux comprendre son rôle dans le contrôle d'un lentivirus endogène, (ii) dans son activité contre d'autres anciennes infections représentées par des rétrovirus endogènes humains et (iii) en tant que protéine capable de générer des mutants de la capside. Premièrement nous nous sommes intéressés à TRIM5a chez deux espèces de lémuriens dont Microcebus murinus qui porte le lentivirus endogène PSIV dans son génome depuis plusieurs millions d'années,. Nous avons observé que TRIM5a chez M. murinus a un spectre d'activité antivirale réduit à l'opposé de TRIM5a chez le Lemur catta - non porteur du PSIV endogène - qui bloque une large variété de rétrovirus dont le PSIV. De ce fait TRIM5a aurait pu contribuer à protéger certaines espèces de lémuriens vis-à-vis d'anciennes infections par le PSIV. A l'inverse du PSIV, des virus dérivés des rétrovirus endogènes humains HERV-K and HERV-H se sont révélés largement résistants à l'inhibition par TRIM5a. Ces données illustrent une absence de protection par TRIM5a face à d'autres anciennes infections rétrovirales. Puis, pour évaluer l'impact de la protéine TRIM5a humaine sur le VIH-1, nous avons testé l'effet de mutations des résidues sous sélection positive dans la capside du VIH-1 sur l'inhibition par TRIM5a. Nos résultats montrent que TRIM5a ne jouerait pas un rôle significatif dans l'évolution de la capside du VIH-1. Enfin notre travail a porté sur le facteur anti-VIH-1 SAMHD1 récemment découvert, que nous avons séquencé chez 25 espèces de primates. L'analyse évolutive des sites sous sélection positive et des expériences fonctionnelles ont permis d'identifier le domaine de SAMHD1 interagissant avec la protéine lentivirale Vpx. De même que d'autres protéines virales contrecarrent les facteurs de restriction en les menant à la dégradation, nous avons observé que Vpx induit la dégradation de SAMHD1 de manière spécifique à l'espèce. Ces découvertes contribuent à comprendre comment les facteurs de restriction et les virus co-évoluent pour se neutraliser l'un l'autre. - Viruses hijack the host cellular machinery to replicate. They adapt to infect optimally host cells while escaping host defense systems. Viruses and the host coevolve in an evolutionary struggle. My thesis work has been devoted to study the evolutionary signatures of host factors acting as restriction factors that block retroviral replication in primates. Specifically, my work aimed at using evolutionary data to inform functional analyses and biology. We studied the anti-HIV-1 factor TRIM5a (i) in prosimians to better understand its possible role in the control of an endogenous lentivirus, (ii) in its activity against other ancient infections - as represented by HERVs, and (iii) as a protein capable of generating escape mutants in the viral capsid. First, my work focused on two lemur species, one of which was the gray mouse lemur that carries the endogenous lentivirus PSIV integrated in its genome for several million years. TRIM5a from gray mouse lemur exhibited a limited antiviral spectrum as opposed to TRIM5a from ring-tailed lemur - not a host of PSIV - that is able to block diverse retroviruses notably PSIV. These results support the possible contribution of TRIM5a in protecting lemur species from ancient infection by PSIV. In contrast, chimeric viruses derived from two human endogenous retroviruses were broadly resistant to TRIM5a-mediated restriction, suggesting TRIM5a lack of activity against other types of ancient infections. To evaluate the recent impact of human TRIM5a on HIV-1 evolution, we tested whether variants at positively selected sites in the HIV-1 capsid affected the ability of human TRIM5a alleles to restrict HIV-1. Our results indicate that TRIM5a does not play a significant role in the evolution of HIV1 capsid. At last, our work concentrated on the newly discovered anti-HIV-1 restriction factor SAMHD1. We determined its coding sequence in a panel of 25 species of primates. Evolutionary analyses of positively selected sites in SAMHD1 domains and functional assays identified the domain of SAMHD1 interacting with the lentiviral protein Vpx. Similar to other viral countermeasures targeting cellular restriction factors, Vpx was responsible of the degradation of SAMHD1 orthologs in a species-specific manner. These findings contributed to understanding how restriction factors and viruses evolve to counteract each other.

The specificity of TRIM5 alpha-mediated restriction is influenced by its coiled-coil domain.

Retroviruses are both powerful evolutionary forces and dangerous threats to genome integrity. As such, they have imposed strong selective pressure on their hosts, notably triggering the emergence of restriction factors, such as TRIM5 alpha, that act as potent barriers to their cross-species transmission. TRIM5 alpha orthologues from different primates have distinct retroviral restriction patterns, largely dictated by the sequence of their C-terminal PRYSPRY domain, which binds the capsid protein of incoming virions. Here, by combining genetic and functional analyses of human and squirrel monkey TRIM5 alpha, we demonstrate that the coiled-coil domain of this protein, thus far essentially known for mediating oligomerization, also conditions the spectrum of antiretroviral activity. Furthermore, we identify three coiled-coil residues responsible for this effect, one of which has been under positive selection during primate evolution, notably in New World monkeys. These results indicate that the PRYSPRY and coiled-coil domains cooperate to determine the specificity of TRIM5 alpha-mediated capture of retroviral capsids, shedding new light on this complex event.

Role of retroviral restriction factors in the interferon-α-mediated suppression of HIV-1 in vivo.

The antiviral potency of the cytokine IFN-α has been long appreciated but remains poorly understood. A number of studies have suggested that induction of the apolipoprotein B mRNA editing enzyme, catalytic polypeptide 3 (APOBEC3) and bone marrow stromal cell antigen 2 (BST-2/tetherin/CD317) retroviral restriction factors underlies the IFN-α-mediated suppression of HIV-1 replication in vitro. We sought to characterize the as-yet-undefined relationship between IFN-α treatment, retroviral restriction factors, and HIV-1 in vivo. APOBEC3G, APOBEC3F, and BST-2 expression levels were measured in HIV/hepatitis C virus (HCV)-coinfected, antiretroviral therapy-naïve individuals before, during, and after pegylated IFN-α/ribavirin (IFN-α/riba) combination therapy. IFN-α/riba therapy decreased HIV-1 viral load by -0.921 (±0.858) log(10) copies/mL in HIV/HCV-coinfected patients. APOBEC3G/3F and BST-2 mRNA expression was significantly elevated during IFN-α/riba treatment in patient-derived CD4+ T cells (P < 0.04 and P < 0.008, paired Wilcoxon), and extent of BST-2 induction was correlated with reduction in HIV-1 viral load during treatment (P < 0.05, Pearson's r). APOBEC3 induction during treatment was correlated with degree of viral hypermutation (P < 0.03, Spearman's ρ), and evolution of the HIV-1 accessory protein viral protein U (Vpu) during IFN-α/riba treatment was suggestive of increased BST-2-mediated selection pressure. These data suggest that host restriction factors play a critical role in the antiretroviral capacity of IFN-α in vivo, and warrant investigation into therapeutic strategies that specifically enhance the expression of these intrinsic immune factors in HIV-1-infected individuals.

The Peroxisomal Enzyme L-PBE Is Required to Prevent the Dietary Toxicity of Medium-Chain Fatty Acids.

Specific metabolic pathways are activated by different nutrients to adapt the organism to available resources. Although essential, these mechanisms are incompletely defined. Here, we report that medium-chain fatty acids contained in coconut oil, a major source of dietary fat, induce the liver ω-oxidation genes Cyp4a10 and Cyp4a14 to increase the production of dicarboxylic fatty acids. Furthermore, these activate all ω- and β-oxidation pathways through peroxisome proliferator activated receptor (PPAR) α and PPARγ, an activation loop normally kept under control by dicarboxylic fatty acid degradation by the peroxisomal enzyme L-PBE. Indeed, L-pbe(-/-) mice fed coconut oil overaccumulate dicarboxylic fatty acids, which activate all fatty acid oxidation pathways and lead to liver inflammation, fibrosis, and death. Thus, the correct homeostasis of dicarboxylic fatty acids is a means to regulate the efficient utilization of ingested medium-chain fatty acids, and its deregulation exemplifies the intricate relationship between impaired metabolism and inflammation.

Evolutionary and Functional Analyses of the Interaction between the Myeloid Restriction Factor SAMHD1 and the Lentiviral Vpx Protein.

SAMHD1 has recently been identified as an HIV-1 restriction factor operating in myeloid cells. As a countermeasure, the Vpx accessory protein from HIV-2 and certain lineages of SIV have evolved to antagonize SAMHD1 by inducing its ubiquitin-proteasome-dependent degradation. Here, we show that SAMHD1 experienced strong positive selection episodes during primate evolution that occurred in the Catarrhini ancestral branch prior to the separation between hominoids (gibbons and great apes) and Old World monkeys. The identification of SAMHD1 residues under positive selection led to mapping the Vpx-interaction domain of SAMHD1 to its C-terminal region. Importantly, we found that while SAMHD1 restriction activity toward HIV-1 is evolutionarily maintained, antagonism of SAMHD1 by Vpx is species-specific. The distinct evolutionary signature of SAMHD1 sheds light on the development of its antiviral specificity.

Analysis of the contribution of the β-oxidation auxiliary enzymes in the degradation of the dietary conjugated linoleic acid 9-cis-11-trans-octadecadienoic acid in the peroxisomes of Saccharomyces cerevisiae

Beta-oxidation of the conjugated linoleic acid 9-cis,11-trans-octadecadienoic acid (rumenic acid) was analyzed in vivo in Saccharomyces cerevisiae by monitoring polyhydroxyalkanoate production in the peroxisome. Polyhydroxyalkanoate is synthesized by the polymerization of the beta-oxidation intermediates 3-hydroxyacyl-CoAs via a bacterial polyhydroxyalkanoate synthase targeted to the peroxisome. The amount of polyhydroxyalkanaote synthesized from the degradation of rumenic acid was found to be similar to the amount synthesized from the degradation of 10-trans,12-cis-octadecadienoic acid, oleic acid or 10-cis-heptadecenoic acid. Furthermore, the degradation of 10-cis-heptadecenoic acid was found to be unaffected by the presence of rumenic acid in the media. Efficient degradation of rumenic acid was found to be independent of the Delta(3,5),Delta(2,4)-dienoyl-CoA isomerase but instead relied on the presence of Delta(3),Delta(2)-enoyl-CoA isomerase activity. The presence of the unsaturated monomer 3-hydroxydodecenoic acid in polyhydroxyalkanoate derived from rumenic acid degradation was found to be dependent on the presence of a Delta(3),Delta(2)-enoyl-CoA isomerase activity. Together, these data indicate that rumenic acid is mainly degraded in vivo in S. cerevisiae through a pathway requiring only the participation of the auxiliary enzymes Delta(3),Delta(2)-enoyl-CoA isomerase, along with the enzyme of the core beta-oxidation cycle.

Interactions of dietary whole-grain intake with fasting glucose- and insulin-related genetic loci in individuals of European descent: a meta-analysis of 14 cohort studies.

Whole-grain foods are touted for multiple health benefits, including enhancing insulin sensitivity and reducing type 2 diabetes risk. Recent genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs) associated with fasting glucose and insulin concentrations in individuals free of diabetes. We tested the hypothesis that whole-grain food intake and genetic variation interact to influence concentrations of fasting glucose and insulin. Via meta-analysis of data from 14 cohorts comprising ∼ 48,000 participants of European descent, we studied interactions of whole-grain intake with loci previously associated in GWAS with fasting glucose (16 loci) and/or insulin (2 loci) concentrations. For tests of interaction, we considered a P value <0.0028 (0.05 of 18 tests) as statistically significant. Greater whole-grain food intake was associated with lower fasting glucose and insulin concentrations independent of demographics, other dietary and lifestyle factors, and BMI (β [95% CI] per 1-serving-greater whole-grain intake: -0.009 mmol/l glucose [-0.013 to -0.005], P < 0.0001 and -0.011 pmol/l [ln] insulin [-0.015 to -0.007], P = 0.0003). No interactions met our multiple testing-adjusted statistical significance threshold. The strongest SNP interaction with whole-grain intake was rs780094 (GCKR) for fasting insulin (P = 0.006), where greater whole-grain intake was associated with a smaller reduction in fasting insulin concentrations in those with the insulin-raising allele. Our results support the favorable association of whole-grain intake with fasting glucose and insulin and suggest a potential interaction between variation in GCKR and whole-grain intake in influencing fasting insulin concentrations.

Sodium and potassium balance depends on αENaC expression in connecting tubule.

Mutations in α, β, or γ subunits of the epithelial sodium channel (ENaC) can downregulate ENaC activity and cause a severe salt-losing syndrome with hyperkalemia and metabolic acidosis, designated pseudohypoaldosteronism type 1 in humans. In contrast, mice with selective inactivation of αENaC in the collecting duct (CD) maintain sodium and potassium balance, suggesting that the late distal convoluted tubule (DCT2) and/or the connecting tubule (CNT) participates in sodium homeostasis. To investigate the relative importance of ENaC-mediated sodium absorption in the CNT, we used Cre-lox technology to generate mice lacking αENaC in the aquaporin 2-expressing CNT and CD. Western blot analysis of microdissected cortical CD (CCD) and CNT revealed absence of αENaC in the CCD and weak αENaC expression in the CNT. These mice exhibited a significantly higher urinary sodium excretion, a lower urine osmolality, and an increased urine volume compared with control mice. Furthermore, serum sodium was lower and potassium levels were higher in the genetically modified mice. With dietary sodium restriction, these mice experienced significant weight loss, increased urinary sodium excretion, and hyperkalemia. Plasma aldosterone levels were significantly elevated under both standard and sodium-restricted diets. In summary, αENaC expression within the CNT/CD is crucial for sodium and potassium homeostasis and causes signs and symptoms of pseudohypoaldosteronism type 1 if missing.