1
|
Cen S, Rasmussen DA. Exploring the Accuracy and Limits of Algorithms for Localizing Recombination Breakpoints. Mol Biol Evol 2024; 41:msae133. [PMID: 38917277 PMCID: PMC11229816 DOI: 10.1093/molbev/msae133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Revised: 06/04/2024] [Accepted: 06/11/2024] [Indexed: 06/27/2024] Open
Abstract
Phylogenetic methods are widely used to reconstruct the evolutionary relationships among species and individuals. However, recombination can obscure ancestral relationships as individuals may inherit different regions of their genome from different ancestors. It is, therefore, often necessary to detect recombination events, locate recombination breakpoints, and select recombination-free alignments prior to reconstructing phylogenetic trees. While many earlier studies have examined the power of different methods to detect recombination, very few have examined the ability of these methods to accurately locate recombination breakpoints. In this study, we simulated genome sequences based on ancestral recombination graphs and explored the accuracy of three popular recombination detection methods: MaxChi, 3SEQ, and Genetic Algorithm Recombination Detection. The accuracy of inferred breakpoint locations was evaluated along with the key factors contributing to variation in accuracy across datasets. While many different genomic features contribute to the variation in performance across methods, the number of informative sites consistent with the pattern of inheritance between parent and recombinant child sequences always has the greatest contribution to accuracy. While partitioning sequence alignments based on identified recombination breakpoints can greatly decrease phylogenetic error, the quality of phylogenetic reconstructions depends very little on how breakpoints are chosen to partition the alignment. Our work sheds light on how different features of recombinant genomes affect the performance of recombination detection methods and suggests best practices for reconstructing phylogenies based on recombination-free alignments.
Collapse
Affiliation(s)
- Shi Cen
- Bioinformatics Research Center, North Carolina State University, Raleigh, NC, USA
| | - David A Rasmussen
- Bioinformatics Research Center, North Carolina State University, Raleigh, NC, USA
- Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
2
|
Li L, Qin R, Liu Y, Tseng YS, Zhang W, Yu L, Mietzsch M, Zou X, Liu H, Lu G, Hu H, Mckenna R, Yang J, Wei Y, Agbandje-Mckenna M, Hu J, Yang L. Dissecting positive selection events and immunological drives during the evolution of adeno-associated virus lineages. PLoS Pathog 2024; 20:e1012260. [PMID: 38885242 PMCID: PMC11182496 DOI: 10.1371/journal.ppat.1012260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Accepted: 05/14/2024] [Indexed: 06/20/2024] Open
Abstract
Adeno-associated virus (AAV) serotypes from primates are being developed and clinically used as vectors for human gene therapy. However, the evolutionary mechanism of AAV variants is far from being understood, except that genetic recombination plays an important role. Furthermore, little is known about the interaction between AAV and its natural hosts, human and nonhuman primates. In this study, natural AAV capsid genes were subjected to systemic evolutionary analysis with a focus on selection drives during the diversification of AAV lineages. A number of positively selected sites were identified from these AAV lineages with functional relevance implied by their localization on the AAV structures. The selection drives of the two AAV2 capsid sites were further investigated in a series of biological experiments. These observations did not support the evolution of the site 410 of the AAV2 capsid driven by selection pressure from the human CD4+ T-cell response. However, positive selection on site 548 of the AAV2 capsid was directly related to host humoral immunity because of the profound effects of mutations at this site on the immune evasion of AAV variants from human neutralizing antibodies at both the individual and population levels. Overall, this work provides a novel interpretation of the genetic diversity and evolution of AAV lineages in their natural hosts, which may contribute to their further engineering and application in human gene therapy.
Collapse
Affiliation(s)
- Lirong Li
- Department of Cardiology and Laboratory of Gene Therapy for Heart Diseases, State Key Laboratory of Biotherapy, and Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Runkuan Qin
- Department of Cardiology and Laboratory of Gene Therapy for Heart Diseases, State Key Laboratory of Biotherapy, and Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Yunbo Liu
- Department of Cardiology and Laboratory of Gene Therapy for Heart Diseases, State Key Laboratory of Biotherapy, and Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Yu-Shan Tseng
- Department of Biochemistry and Molecular Biology, Center for Structural Biology, The McKnight Brain Institute, University of Florida, Gainesville, Florida, United States of America
| | - Weihan Zhang
- General Surgery Department, Gastric Cancer Center and Laboratory of Gastric Cancer, State Key Laboratory of Biotherapy and Collaborative Innovation Center of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Lin Yu
- State Key Laboratory of Biotherapy and Cancer Center, Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Mario Mietzsch
- Department of Biochemistry and Molecular Biology, Center for Structural Biology, The McKnight Brain Institute, University of Florida, Gainesville, Florida, United States of America
| | - Xinkai Zou
- Center for Immunology and Hematology, Department of Biotherapy, Cancer Center and State Key Laboratory of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Haizhou Liu
- Computational Virology Group, Center for Bacteria and Viruses Resources and Bioinformation, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan, Hubei, China
| | - Guangwen Lu
- Department of Emergency Medicine, State Key Laboratory of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Hongbo Hu
- Center for Immunology and Hematology, Department of Biotherapy, Cancer Center and State Key Laboratory of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Robert Mckenna
- Department of Biochemistry and Molecular Biology, Center for Structural Biology, The McKnight Brain Institute, University of Florida, Gainesville, Florida, United States of America
| | - Jinliang Yang
- State Key Laboratory of Biotherapy and Cancer Center, Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Yuquan Wei
- Laboratory of Aging Research and Cancer Drug Target, State Key Laboratory of Biotherapy and Cancer Center, National Clinical Research Center for Geriatrics, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Mavis Agbandje-Mckenna
- Department of Biochemistry and Molecular Biology, Center for Structural Biology, The McKnight Brain Institute, University of Florida, Gainesville, Florida, United States of America
| | - Jiankun Hu
- General Surgery Department, Gastric Cancer Center and Laboratory of Gastric Cancer, State Key Laboratory of Biotherapy and Collaborative Innovation Center of Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| | - Lin Yang
- Department of Cardiology and Laboratory of Gene Therapy for Heart Diseases, State Key Laboratory of Biotherapy, and Collaborative Innovation Center for Biotherapy, West China Hospital, Sichuan University, Chengdu, Sichuan, China
| |
Collapse
|
3
|
Molteni C, Forni D, Cagliani R, Mozzi A, Clerici M, Sironi M. Evolution of the orthopoxvirus core genome. Virus Res 2023; 323:198975. [PMID: 36280003 PMCID: PMC9586335 DOI: 10.1016/j.virusres.2022.198975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 10/17/2022] [Accepted: 10/18/2022] [Indexed: 11/07/2022]
Abstract
Orthopoxviruses comprise several relevant pathogens, including the causative agent of smallpox and monkeypox virus. Analysis of orthopoxvirus genome evolution mainly focused on gene gains/losses. We instead analyzed core genes, which are conserved in all orthopoxviruses. We show that, despite their strong constraint, some genes involved in viral morphogenesis and transcription/replication were targets of pervasive positive selection, which was relatively uncommon in immunomodulatory genes. However at least three of the positively selected genes, E3L, A24R, and H3L, might have evolved in response to immune selection. Episodic positive selection was particularly common on the internal branches of the orthopox phylogeny and on the monkeypox virus lineage. The latter showed evidence of episodic positive selection at the D14L gene, which encodes a modulator of complement activation (MOPICE). Notably, two genes (B1R and A33R) targeted by episodic selection on more than one branch are involved in forms of intra-genomic conflict. Finally, we found that, in orthopoxvirus proteomes, intrinsically disordered regions (IDRs) tend to be less constrained and are common targets of positive selection. Extension of our analysis to all poxviruses showed no evidence that the IDR fraction differs with host range. Conversely, we found a strong effect of base composition, which was however not sufficient to explain IDR fraction. We thus suggest that, in poxviruses, the IDR fraction is maintained by modulating GC content to accommodate disorder-promoting codons. Overall, our data provide novel insight in orthopoxvirus evolution and provide a list of genes and sites that are expected to modulate viral phenotypes.
Collapse
Affiliation(s)
- Cristian Molteni
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy.
| | - Diego Forni
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| | - Rachele Cagliani
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| | - Alessandra Mozzi
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| | - Mario Clerici
- University of Milan, Milan, Italy; Don C. Gnocchi Foundation ONLUS, IRCCS, Milan, Italy
| | - Manuela Sironi
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| |
Collapse
|
4
|
Wei X, Meng B, Peng H, Li Y, Liu M, Si H, Wu R, Chen H, Bai Y, Li Y, Feng Q, Wang C, Zhao X. Hemorrhagic fever with renal syndrome caused by destruction of residential area of rodent in a construction site: epidemiological investigation. BMC Infect Dis 2022; 22:761. [PMID: 36175847 PMCID: PMC9521858 DOI: 10.1186/s12879-022-07744-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 09/20/2022] [Indexed: 01/18/2023] Open
Abstract
BACKGROUND An outbreak of hemorrhagic fever with renal syndrome (HFRS), caused by a Hantavirus, affected nine adult males in the southwest area of Xi'an in November 2020 was analyzed in this study. METHODS Clinical and epidemiological data of HFRS patients in this outbreak were retrospectively analyzed. The whole genome of a hantavirus named 201120HV03xa (hv03xa for short) isolated from Apodemus agrarius captured in the construction site was sequenced and analyzed. In addition, nine HFRS patients were monitored for the IgG antibody against the HV N protein at 6 and 12 months, respectively. RESULTS In this study, inhalation of aerosolized excreta and contaminated food may be the main source of infection. Genome analysis and phylogenetic analysis showed that hv03xa is a reassortment strain of HTNV, having an S segment related to A16 of HTN 4, an M segment related to Q37 and Q10 of HTN 4, and an L segment related to prototype strain 76-118 of HTN 7. Potential recombination was detected in the S segment of hv03xa strain. The anti-HV-IgG level of all the patients persist for at least one year after infection. CONCLUSIONS This report documented an HFRS outbreak in Xi'an, China, which provided the basic data for epidemiological surveillance of endemic HTNV infection and facilitated to predict disease risk and implement prevention measures.
Collapse
Affiliation(s)
- Xiao Wei
- Centers for Disease Control and Prevention of PLA, Beijing, China
| | - Biao Meng
- Centers for Disease Control and Prevention of PLA, Beijing, China.,Department of Epidemiology and Biostatistics, School of Public Health, Anhui Medical University, Hefei, China
| | - Hong Peng
- Centers for Disease Control and Prevention of PLA, Beijing, China
| | - Yan Li
- Centers for Disease Control and Prevention of PLA, Beijing, China
| | - Min Liu
- PLA 63750 Military Hospital, Xi'an, Shaanxi, China
| | - Hairui Si
- PLA 63750 Military Hospital, Xi'an, Shaanxi, China
| | - Rui Wu
- Xi'an Center for Disease Control and Prevention, Xi'an, Shaanxi, China
| | - Hailong Chen
- Xi'an Center for Disease Control and Prevention, Xi'an, Shaanxi, China
| | - Ying Bai
- PLA 63750 Military Hospital, Xi'an, Shaanxi, China
| | - Yan Li
- PLA 63750 Military Hospital, Xi'an, Shaanxi, China
| | - Qunling Feng
- PLA 63750 Military Hospital, Xi'an, Shaanxi, China.
| | - Changjun Wang
- Centers for Disease Control and Prevention of PLA, Beijing, China. .,Department of Epidemiology and Biostatistics, School of Public Health, Anhui Medical University, Hefei, China.
| | - Xiangna Zhao
- Centers for Disease Control and Prevention of PLA, Beijing, China. .,Department of Epidemiology and Biostatistics, School of Public Health, Anhui Medical University, Hefei, China.
| |
Collapse
|
5
|
Shikov AE, Malovichko YV, Nizhnikov AA, Antonets KS. Current Methods for Recombination Detection in Bacteria. Int J Mol Sci 2022; 23:ijms23116257. [PMID: 35682936 PMCID: PMC9181119 DOI: 10.3390/ijms23116257] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 05/30/2022] [Accepted: 05/30/2022] [Indexed: 02/05/2023] Open
Abstract
The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.
Collapse
Affiliation(s)
- Anton E. Shikov
- Laboratory for Proteomics of Supra-Organismal Systems, All-Russia Research Institute for Agricultural Microbiology (ARRIAM), 196608 St. Petersburg, Russia; (A.E.S.); (Y.V.M.); (A.A.N.)
- Faculty of Biology, St. Petersburg State University (SPbSU), 199034 St. Petersburg, Russia
| | - Yury V. Malovichko
- Laboratory for Proteomics of Supra-Organismal Systems, All-Russia Research Institute for Agricultural Microbiology (ARRIAM), 196608 St. Petersburg, Russia; (A.E.S.); (Y.V.M.); (A.A.N.)
- Faculty of Biology, St. Petersburg State University (SPbSU), 199034 St. Petersburg, Russia
| | - Anton A. Nizhnikov
- Laboratory for Proteomics of Supra-Organismal Systems, All-Russia Research Institute for Agricultural Microbiology (ARRIAM), 196608 St. Petersburg, Russia; (A.E.S.); (Y.V.M.); (A.A.N.)
- Faculty of Biology, St. Petersburg State University (SPbSU), 199034 St. Petersburg, Russia
| | - Kirill S. Antonets
- Laboratory for Proteomics of Supra-Organismal Systems, All-Russia Research Institute for Agricultural Microbiology (ARRIAM), 196608 St. Petersburg, Russia; (A.E.S.); (Y.V.M.); (A.A.N.)
- Faculty of Biology, St. Petersburg State University (SPbSU), 199034 St. Petersburg, Russia
- Correspondence:
| |
Collapse
|
6
|
Alex CE, Canuti M, Schlesinger MS, Jackson KA, Needle D, Jardine C, Nituch L, Bourque L, Lang AS, Pesavento PA. Natural disease and evolution of an amdoparvovirus endemic in striped skunks (
Mephitis mephitis
). Transbound Emerg Dis 2022; 69:e1758-e1767. [DOI: 10.1111/tbed.14511] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2022] [Revised: 03/07/2022] [Accepted: 03/08/2022] [Indexed: 11/30/2022]
Affiliation(s)
- Charles E. Alex
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine University of California‐Davis Davis CA USA
| | - Marta Canuti
- Department of Biology Memorial University of Newfoundland St. John's NL Canada
| | - Maya S. Schlesinger
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine University of California‐Davis Davis CA USA
| | - Kenneth A. Jackson
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine University of California‐Davis Davis CA USA
| | - David Needle
- New Hampshire Veterinary Diagnostic Laboratory, College of Life Sciences and Agriculture University of New Hampshire Durham NH USA
| | - Claire Jardine
- Department of Pathobiology, Canadian Wildlife Health Cooperative University of Guelph Guelph ON Canada
| | - Larissa Nituch
- Ontario Ministry of Northern Development Mines, Natural Resources and Forestry Peterborough ON Canada
| | - Laura Bourque
- Canadian Wildlife Health Cooperative – Atlantic Region University of Prince Edward Island 550 University Ave Charlottetown PE C1A4P3 Canada
| | - Andrew S. Lang
- Department of Biology Memorial University of Newfoundland St. John's NL Canada
| | - Patricia A. Pesavento
- Department of Pathology, Microbiology, and Immunology, School of Veterinary Medicine University of California‐Davis Davis CA USA
| |
Collapse
|
7
|
Baker EP, Sayegh R, Kohler KM, Borman W, Goodfellow CK, Brush ER, Barber MF. Evolution of host-microbe cell adherence by receptor domain shuffling. eLife 2022; 11:73330. [PMID: 35076392 PMCID: PMC8860441 DOI: 10.7554/elife.73330] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Accepted: 01/22/2022] [Indexed: 11/29/2022] Open
Abstract
Stable adherence to epithelial surfaces is required for colonization by diverse host-associated microbes. Successful attachment of pathogenic microbes to host cells via adhesin molecules is also the first step in many devastating infections. Despite the primacy of epithelial adherence in establishing host-microbe associations, the evolutionary processes that shape this crucial interface remain enigmatic. Carcinoembryonic antigen-related cell adhesion molecules (CEACAMs) encompass a multifunctional family of vertebrate cell surface proteins which are recurrent targets of bacterial adhesins at epithelial barriers. Here, we show that multiple members of the primate CEACAM family exhibit evidence of repeated natural selection at protein surfaces targeted by bacteria, consistent with pathogen-driven evolution. Divergence of CEACAM proteins between even closely related great apes is sufficient to control molecular interactions with a range of bacterial adhesins. Phylogenetic analyses further reveal that repeated gene conversion of CEACAM extracellular domains during primate divergence plays a key role in limiting bacterial adhesin host tropism. Moreover, we demonstrate that gene conversion has continued to shape CEACAM diversity within human populations, with abundant human CEACAM1 variants mediating evasion of adhesins from pathogenic Neisseria. Together this work reveals a mechanism by which gene conversion shapes first contact between microbes and animal hosts. Trillions of bacteria live in and on the human body. Most of them are harmless but some can cause serious infections. To grow in or on the body, bacteria often attach to proteins on the surface of cells that make up the lining of tissues like the gut or the throat. In some cases, bacteria use these proteins to invade the cells causing an infection. Genetic mutations in the genes encoding these proteins that protect against infection are more likely to be passed on to future generations. This may lead to rapid spread of these beneficial genes in a population. A family of proteins called CEACAMs are frequent targets of infection-causing bacteria. These proteins have been shown to play a role in cancer progression. But they also play many helpful roles in the body, including helping transmit messages between cells, aiding cell growth, and helping the immune system recognize pathogens. Scientists are not sure if these multi-tasking CEACAM proteins can evolve to evade bacteria without affecting their other roles. Baker et al. show that CEACAM proteins targeted by bacteria have undergone rapid evolution in primates. In the experiments, human genes encoding CEACAMs were compared with equivalent genes from 19 different primates. Baker et al. found the changes in human and primate CEACAMs often occur through a process called gene conversion. Gene conversion occurs when DNA sections are copied and pasted from one gene to another. Using laboratory experiments, they showed that some of these changes enabled CEACAM proteins to prevent certain harmful bacteria from binding. The experiments suggest that some versions of CEACAM genes may protect humans or other primates against bacterial infections. Studies in natural populations are needed to test if this is the case. Learning more about how CEACAM proteins evolve and what they do may help scientists better understand the role they play in cancer and help improve cancer care. Studying CEACAM evolution may also help scientists understand how bacteria and other pathogens drive protein evolution in the body.
Collapse
Affiliation(s)
- EmilyClare P Baker
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Ryan Sayegh
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Kristin M Kohler
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Wyatt Borman
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Claire K Goodfellow
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Eden R Brush
- Institute of Ecology and Evolution, University of Oregon, Eugene, United States
| | - Matthew F Barber
- Department of Biology, University of Oregon, Eugene, United States
| |
Collapse
|
8
|
Gendreau KL, Hornsby AD, Hague MTJ, McGlothlin JW. Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts. Mol Biol Evol 2021; 38:4077-4094. [PMID: 34129031 PMCID: PMC8476164 DOI: 10.1093/molbev/msab182] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Nav) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Nav gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Nav paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Nav gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Nav1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Nav1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.
Collapse
Affiliation(s)
- Kerry L Gendreau
- Department of Biological Sciences, Virginia Tech, Blacksburg, United States
| | - Angela D Hornsby
- Department of Biological Sciences, Virginia Tech, Blacksburg, United States.,Philip L. Wright Zoological Museum, Division of Biological Sciences, University of Montana, Missoula, United States
| | - Michael T J Hague
- Division of Biological Sciences, University of Montana, Missoula, MT, United States
| | - Joel W McGlothlin
- Department of Biological Sciences, Virginia Tech, Blacksburg, United States
| |
Collapse
|
9
|
Techera C, Marandino A, Tomás G, Grecco S, Hernández M, Hernández D, Panzera Y, Pérez R. Origin, spreading and genetic variability of chicken anaemia virus. Avian Pathol 2021; 50:311-320. [PMID: 33928817 DOI: 10.1080/03079457.2021.1919289] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Chicken anaemia virus (CAV) is a widespread pathogen that causes immunosuppression in chickens. The virus-induced immunosuppression often results in secondary infections and a sub-optimal response to vaccinations, leading to high mortality rates and significant economic losses in the poultry industry. The small circular ssDNA genome (2.3 kb) has three partially overlapping genes: vp1, vp2 and vp3. VP1 capsid protein is highly variable and contains the neutralizing epitopes. Here, we analysed CAV strains from Uruguay using the full-length vp1 gene and performed a global comparative analysis to provide new evidence about the origin, dispersion and genetic variability of the virus. The phylogenetic analysis classified CAV in three or four major clades. Two clades (II and III) grouped most of the strains circulating worldwide including the Uruguayan strains. The phylodynamic analyses indicated that CAV emerged in the early 1900s and diverged to originate clade II and III. This early period of viral emergence was characterised by local diversification promoted by the extremely high substitution rate inferred for the virus (3.8 × 10-4 substitutions/site/year). Later, the virus underwent a global spreading by intra- and inter-continental migrations that correlates with a significant rise in the effective population size. In South America, CAV was introduced in three different migratory events and spread across the continent. Our findings suggest that the current CAV distribution is the consequence of its continuous expansion capability that homogenizes the populations and prevents the detection of clear temporal and geographic patterns of evolution in most strains.RESEARCH HIGHLIGHTS Current strains of chicken anaemia virus emerged in Asia in the early 1900s.Chicken anaemia virus has a high substitution rate.The phylogenetic analysis classified chicken anaemia virus in four major clades.Evolution in South America was characterized by long migration and local spreading.
Collapse
Affiliation(s)
- Claudia Techera
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Ana Marandino
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Gonzalo Tomás
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Sofía Grecco
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Martín Hernández
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Diego Hernández
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Yanina Panzera
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - Ruben Pérez
- Sección Genética Evolutiva, Departamento de Biología Animal, Instituto de Biología, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| |
Collapse
|
10
|
Forni D, Cagliani R, Clerici M, Pozzoli U, Sironi M. You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations. Mol Biol Evol 2020; 37:442-454. [PMID: 31593241 DOI: 10.1093/molbev/msz227] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.
Collapse
Affiliation(s)
- Diego Forni
- Scientific Institute, IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Lecco, Italy
| | - Rachele Cagliani
- Scientific Institute, IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Lecco, Italy
| | - Mario Clerici
- Department of Physiopathology and Transplantation, University of Milan, Milan, Italy.,IRCCS Fondazione Don Carlo Gnocchi, Milan, Italy
| | - Uberto Pozzoli
- Scientific Institute, IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Lecco, Italy
| | - Manuela Sironi
- Scientific Institute, IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Lecco, Italy
| |
Collapse
|
11
|
Ribeiro JMC, Mans BJ. TickSialoFam (TSFam): A Database That Helps to Classify Tick Salivary Proteins, a Review on Tick Salivary Protein Function and Evolution, With Considerations on the Tick Sialome Switching Phenomenon. Front Cell Infect Microbiol 2020; 10:374. [PMID: 32850476 PMCID: PMC7396615 DOI: 10.3389/fcimb.2020.00374] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2020] [Accepted: 06/17/2020] [Indexed: 01/09/2023] Open
Abstract
Tick saliva contains a complex mixture of peptides and non-peptides that counteract their hosts' hemostasis, immunity, and tissue-repair reactions. Recent transcriptomic studies have revealed over one thousand different transcripts coding for secreted polypeptides in a single tick species. Not only do these gene products belong to many expanded families, such as the lipocalins, metalloproteases, Antigen-5, cystatins, and apyrases, but also families that are found exclusively in ticks, such as the evasins, Isac, DAP36, and many others. Phylogenetic analysis of the deduced protein sequences indicate that the salivary genes exhibit an increased rate of evolution due to a lower evolutionary constraint and/or positive selection, allowing for a large diversity of tick salivary proteins. Thus, for each new tick species that has its salivary transcriptome sequenced and assembled, a formidable task of annotation of these transcripts awaits. Currently, as of November 2019, there are over 287 thousand coding sequences deposited at the National Center for Biotechnology Information (NCBI) that are derived from tick salivary gland mRNA. Here, from these 287 thousand sequences we identified 45,264 potential secretory proteins which possess a signal peptide and no transmembrane domains on the mature peptide. By using the psiblast tools, position-specific matrices were constructed and assembled into the TickSialoFam (TSF) database. The TSF is a rpsblastable database that can help with the annotation of tick sialotranscriptomes. The TSA database identified 136 tick salivary secreted protein families, as well as 80 families of endosomal-related products, mostly having a protein modification function. As the number of sequences increases, and new annotation details become available, new releases of the TSF database may become available.
Collapse
Affiliation(s)
- José M. C. Ribeiro
- Section of Vector Biology, Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, Rockville, MD, United States
| | - Ben J. Mans
- Epidemiology, Parasites and Vectors, Agricultural Research Council - Onderstepoort Veterinary Research, Pretoria, South Africa
- The Department of Veterinary Tropical Diseases, University of Pretoria, Pretoria, South Africa
- Department of Life and Consumer Sciences, University of South Africa, Pretoria, South Africa
| |
Collapse
|
12
|
Cagliani R, Forni D, Clerici M, Sironi M. Coding potential and sequence conservation of SARS-CoV-2 and related animal viruses. INFECTION GENETICS AND EVOLUTION 2020; 83:104353. [PMID: 32387562 PMCID: PMC7199688 DOI: 10.1016/j.meegid.2020.104353] [Citation(s) in RCA: 61] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Revised: 04/14/2020] [Accepted: 05/02/2020] [Indexed: 12/13/2022]
Abstract
In December 2019, a novel human-infecting coronavirus (SARS-CoV-2) was recognized in China. In a few months, SARS-CoV-2 has caused thousands of disease cases and deaths in several countries. Phylogenetic analyses indicated that SARS-CoV-2 clusters with SARS-CoV in the Sarbecovirus subgenus and viruses related to SARS-CoV-2 were identified from bats and pangolins. Coronaviruses have long and complex genomes with high plasticity in terms of gene content. To date, the coding potential of SARS-CoV-2 remains partially unknown. We thus used available sequences of bat and pangolin viruses to determine the selective events that shaped the genome structure of SARS-CoV-2 and to assess its coding potential. By searching for signals of significantly reduced variability at synonymous sites (dS), we identified six genomic regions, one of these corresponding to the programmed −1 ribosomal frameshift. The most prominent signal of dS reduction was observed within the E gene. A genome-wide analysis of conserved RNA structures indicated that this region harbors a putative functional RNA element that is shared with the SARS-CoV lineage. Additional signals of reduced dS indicated the presence of internal ORFs. Whereas the presence ORF9a (internal to N) was previously proposed by homology with a well characterized protein of SARS-CoV, ORF3h (for hypothetical, within ORF3a) was not previously described. The predicted product of ORF3h has 90% identity with the corresponding predicted product of SARS-CoV and displays features suggestive of a viroporin. Finally, analysis of the putative ORF10 revealed high dN/dS (3.82) in SARS-CoV-2 and related coronaviruses. In the SARS-CoV lineage, the ORF is predicted to encode a truncated protein and is neutrally evolving. These data suggest that ORF10 encodes a functional protein in SARS-CoV-2 and that positive selection is driving its evolution. Experimental analyses will be necessary to validate and characterize the coding and non-coding functional elements we identified. We analyzed the coding region of SARS-CoV-2 and related bat/pangolin viruses. We identified six regions of significantly low variability at sysnonymous sites. One of these corresponds to a conserved RNA structure shared with the SARS-CoV lineage. The dS reduction within ORF3a corresponds to a potential ORF encoding a viroporin. In SARS-CoV-2 and related viruses, the putative 3′ terminal ORF10 has high dN/dS.
Collapse
Affiliation(s)
- Rachele Cagliani
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| | - Diego Forni
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy
| | - Mario Clerici
- Department of Physiopathology and Transplantation, University of Milan, Milan, Italy; Don C. Gnocchi Foundation ONLUS, IRCCS, Milan, Italy
| | - Manuela Sironi
- Scientific Institute IRCCS E. MEDEA, Bioinformatics, Bosisio Parini, Italy.
| |
Collapse
|
13
|
Forni D, Cagliani R, Clerici M, Pozzoli U, Sironi M. Evolutionary analysis of exogenous and integrated HHV-6A/HHV-6B populations. Virus Evol 2020; 6:veaa035. [PMID: 32551136 PMCID: PMC7293831 DOI: 10.1093/ve/veaa035] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6A/HHV-6B, iciHHV-6A/iciHHV-6B). We reconstructed the population structures of HHV-6A and HHV-6B, showing that HHV-6A diverged less than HHV-6B genomes from the projected common ancestral population. Thus, HHV-6B genomes experienced stronger drift, as also supported by calculation of nucleotide diversity and Tajima's D. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, ancestry proportions were similar in exogenous HHV-6B viruses and iciHHV-6B. In line with previous indications, this suggests the distinct exogenous viral populations that originated iciHHV-6B in subjects with European and Asian ancestry are still causing infections in the corresponding geographic areas. Notably, for both iciHHV-6A and iciHHV-6B, we found that European and American sequences tend to have high proportions of ancestry from viral populations that experienced considerable drift, suggesting that they underwent one or more bottlenecks followed by population expansion. Finally, analysis of HHV-6B exogenous viruses sampled in Japan indicated that proportions of ancestry components of most of these viruses are different from the majority of those sampled in the USA. More generally, we show that, in both viral species, both integrated and exogenous viral genomes have different ancestry components, partially depending on geographic location. It would be extremely important to determine whether such differences account for the diversity of HHV-6A/HHV-6B-associated clinical symptoms and epidemiology. Also, the sequencing of additional exogenous and integrated viral genomes will be instrumental to confirm and expand our conclusions, which are based on a relatively small number of genomes, sequenced with variable quality, and with unequal sampling in terms of geographic origin.
Collapse
Affiliation(s)
- Diego Forni
- Bioinformatics, Scientific Institute IRCCS E. MEDEA, 23842 Bosisio Parini, Lecco, Italy
| | - Rachele Cagliani
- Bioinformatics, Scientific Institute IRCCS E. MEDEA, 23842 Bosisio Parini, Lecco, Italy
| | - Mario Clerici
- Department of Physiopathology and Transplantation, University of Milan, 20090 Milan, Italy.,IRCCS Fondazione Don Carlo Gnocchi, 20148 Milan, Italy
| | - Uberto Pozzoli
- Bioinformatics, Scientific Institute IRCCS E. MEDEA, 23842 Bosisio Parini, Lecco, Italy
| | - Manuela Sironi
- Bioinformatics, Scientific Institute IRCCS E. MEDEA, 23842 Bosisio Parini, Lecco, Italy
| |
Collapse
|
14
|
Population structure of Lassa Mammarenavirus in West Africa. Viruses 2020; 12:v12040437. [PMID: 32294960 PMCID: PMC7232344 DOI: 10.3390/v12040437] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 04/07/2020] [Accepted: 04/10/2020] [Indexed: 01/01/2023] Open
Abstract
Lassa mammarenavirus (LASV) is the etiologic agent of Lassa fever. In endemic regions in West Africa, LASV genetic diversity tends to cluster by geographic area. Seven LASV lineages are recognized, but the role of viral genetic determinants on disease presentation in humans is uncertain. We investigated the geographic structure and distribution of LASV in West Africa. We found strong spatial clustering of LASV populations, with two major east–west and north–south diversity gradients. Analysis of ancestry components indicated that known LASV lineages diverged from an ancestral population that most likely circulated in Nigeria, although alternative locations, such as Togo, cannot be excluded. Extant sequences carrying the largest contribution of this ancestral population include the prototype Pinneo strain, the Togo isolates, and a few viruses isolated in Nigeria. The LASV populations that experienced the strongest drift circulate in Mali and the Ivory Coast. By focusing on sequences form a single LASV sublineage (IIg), we identified an ancestry component possibly associated with protection from a fatal disease outcome. Although the same ancestry component tends to associate with lower viral loads in plasma, the small sample size requires that these results are treated with extreme caution.
Collapse
|
15
|
Forni D, Cagliani R, Clerici M, Pozzoli U, Sironi M. A complex evolutionary relationship between HHV-6A and HHV-6B. Virus Evol 2019; 5:vez043. [PMID: 31649826 PMCID: PMC6800887 DOI: 10.1093/ve/vez043] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6, iciHHV6). We reconstructed the population structure of HHV-6 and we show that HHV-6A genomes diverged less than HHV-6B genomes from the ancestral common HHV-6A/B population. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, exogenous viral and iciHHV-6B populations were similar in terms of ancestry components, with no evident geographic structuring. Most HHV-6B genomes sampled to date derive from viral populations that experienced considerable drift. However, a population of HHV-6 exogenous viruses, currently classified as HHV-6B and sampled in New York state, formed a separate cluster (NY cluster) and harbored a considerable portion of HHV-6A-like ancestry. Recombination detection methods identified these viruses as interspecies recombinants, but phylogenetic reconstruction indicated that the recombination signals are due to shared ancestry. In analogy to iciHHV-6A, NY cluster viruses have high nucleotide diversity and constant population size. We propose that HHV-6A sequences and the NY cluster population diverged from an ancestral HHV-6A-like population. A relatively recent bottleneck of the NY (or a related) population with subsequent expansion originated most HHV-6B genomes currently sampled. Our findings indicate that the distinction between HHV-6A and -6B is not as clear-cut as previously thought. More generally, epidemiological and clinical surveys would benefit from taking HHV-6 genetic diversity into account.
Collapse
Affiliation(s)
- Diego Forni
- Bioinformatics, Scientific Institute IRCCS E. Medea, Bosisio Parini, Lecco, Italy
- Corresponding author: E-mail:
| | - Rachele Cagliani
- Bioinformatics, Scientific Institute IRCCS E. Medea, Bosisio Parini, Lecco, Italy
| | - Mario Clerici
- Department of Physiopathology and Transplantation, University of Milan, Milan, Italy
- IRCCS Fondazione Don Carlo Gnocchi, Milan, Italy
| | - Uberto Pozzoli
- Bioinformatics, Scientific Institute IRCCS E. Medea, Bosisio Parini, Lecco, Italy
| | - Manuela Sironi
- Bioinformatics, Scientific Institute IRCCS E. Medea, Bosisio Parini, Lecco, Italy
| |
Collapse
|
16
|
Abduriyim S, Nishita Y, Kosintsev PA, Raichev E, Väinölä R, Kryukov AP, Abramov AV, Kaneko Y, Masuda R. Evolution of MHC class I genes in Eurasian badgers, genus Meles (Carnivora, Mustelidae). Heredity (Edinb) 2019; 122:205-218. [PMID: 29959426 PMCID: PMC6327056 DOI: 10.1038/s41437-018-0100-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Revised: 05/30/2018] [Accepted: 05/30/2018] [Indexed: 11/09/2022] Open
Abstract
Because of their role in immune defense against pathogens, major histocompatibility complex (MHC) genes are useful in evolutionary studies on how wild vertebrates adapt to their environments. We investigated the molecular evolution of MHC class I (MHCI) genes in four closely related species of Eurasian badgers, genus Meles. All four species of badgers showed similarly high variation in MHCI sequences compared to other Carnivora. We identified 7-21 putatively functional MHCI sequences in each of the badger species, and 2-7 sequences per individual, indicating the existence of 1-4 loci. MHCI exon 2 and 3 sequences encoding domains α1 and α2 exhibited different clade topologies in phylogenetic networks. Non-synonymous nucleotide substitutions at codons for antigen-binding sites exceeded synonymous substitutions for domain α1 but not for domain α2, suggesting that the domains α1 and α2 likely had different evolutionary histories in these species. Positive selection and recombination seem to have shaped the variation in domain α2, whereas positive selection was dominant in shaping the variation in domain α1. In the separate phylogenetic analyses for exon 2, exon 3, and intron 2, each showed three clades of Meles alleles, with rampant trans-species polymorphism, indicative of the long-term maintenance of ancestral MHCI polymorphism by balancing selection.
Collapse
Affiliation(s)
- Shamshidin Abduriyim
- Department of Natural History Sciences, Graduate School of Science, Hokkaido University, Sapporo, 060-0810, Japan
| | - Yoshinori Nishita
- Department of Natural History Sciences, Graduate School of Science, Hokkaido University, Sapporo, 060-0810, Japan
- Department of Biological Sciences, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan
| | - Pavel A Kosintsev
- Institute of Plant and Animal Ecology, Ural Branch, Russian Academy of Sciences, Ekaterinburg, 620144, Russia
| | - Evgeniy Raichev
- Agricultural Faculty, Trakia University, 6000, Stara Zagora, Bulgaria
| | - Risto Väinölä
- Finnish Museum of Natural History, University of Helsinki, P.O. Box 17, FI-00014, Helsinki, Finland
| | - Alexey P Kryukov
- Federal Scientific Center of the East Asia Terrestrial Biodiversity, Far Eastern Branch of the Russian Academy of Sciences, Vladivostok, 690022, Russia
| | - Alexei V Abramov
- Zoological Institute, Russian Academy of Sciences, Saint Petersburg, 199034, Russia
| | - Yayoi Kaneko
- Carnivore Ecology and Conservation Research Group, Tokyo University of Agriculture and Technology, Tokyo, 183-8509, Japan
| | - Ryuichi Masuda
- Department of Natural History Sciences, Graduate School of Science, Hokkaido University, Sapporo, 060-0810, Japan.
- Department of Biological Sciences, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan.
| |
Collapse
|
17
|
Pontremoli C, Forni D, Cagliani R, Sironi M. Analysis of Reptarenavirus genomes indicates different selective forces acting on the S and L segments and recent expansion of common genotypes. INFECTION GENETICS AND EVOLUTION 2018; 64:212-218. [DOI: 10.1016/j.meegid.2018.06.031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Revised: 06/28/2018] [Accepted: 06/28/2018] [Indexed: 01/20/2023]
|
18
|
Gao Y, Suárez NM, Wilkie GS, Dong C, Bergmann S, Lee PYA, Davison AJ, Vanderplasschen AFC, Boutier M. Genomic and biologic comparisons of cyprinid herpesvirus 3 strains. Vet Res 2018; 49:40. [PMID: 29716648 PMCID: PMC5930815 DOI: 10.1186/s13567-018-0532-z] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2018] [Accepted: 03/30/2018] [Indexed: 11/21/2022] Open
Abstract
Cyprinid herpesvirus 3 (CyHV-3) is the archetypal fish alloherpesvirus and the etiologic agent of a lethal disease in common and koi carp. To date, the genome sequences of only four CyHV-3 isolates have been published, but no comparisons of the biologic properties of these strains have been reported. We have sequenced the genomes of a further seven strains from various geographical sources, and have compared their growth in vitro and virulence in vivo. The major findings were: (i) the existence of the two genetic lineages previously described as European and Asian was confirmed, but inconsistencies between the geographic origin and genotype of some strains were revealed; (ii) potential inter-lineage recombination was detected in one strain, which also suggested the existence of a third, as yet unidentified lineage; (iii) analysis of genetic disruptions led to the identification of non-essential genes and their potential role in virulence; (iv) comparison of the in vitro and in vivo properties of strains belonging to the two lineages revealed that inter-lineage polymorphisms do not contribute to the differences in viral fitness observed; and (v) a negative correlation was observed among strains between viral growth in vitro and virulence in vivo. This study illustrates the importance of coupling genomic and biologic comparisons of viral strains in order to enhance understanding of viral evolution and pathogenesis.
Collapse
Affiliation(s)
- Yuan Gao
- Immunology-Vaccinology, Department of Infectious and Parasitic Diseases (B43b), Fundamental and Applied Research for Animals & Health (FARAH), Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| | - Nicolás M Suárez
- MRC-University of Glasgow Centre for Virus Research, Glasgow, UK
| | - Gavin S Wilkie
- MRC-University of Glasgow Centre for Virus Research, Glasgow, UK
| | - Chuanfu Dong
- MOE Key Laboratory of Aquatic Food Safety/State Key Laboratory for Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou, China
| | - Sven Bergmann
- Friedrich-Loeffler Institut, Federal Research Institute for Animal Health, Institute of Infectology, Greifswald-Insel Riems, Germany
| | - Pei-Yu Alison Lee
- Department of Research and Development, GeneReach, Biotechnology Corporation, Taichung, China
| | - Andrew J Davison
- MRC-University of Glasgow Centre for Virus Research, Glasgow, UK
| | - Alain F C Vanderplasschen
- Immunology-Vaccinology, Department of Infectious and Parasitic Diseases (B43b), Fundamental and Applied Research for Animals & Health (FARAH), Faculty of Veterinary Medicine, University of Liège, Liège, Belgium.
| | - Maxime Boutier
- Immunology-Vaccinology, Department of Infectious and Parasitic Diseases (B43b), Fundamental and Applied Research for Animals & Health (FARAH), Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| |
Collapse
|
19
|
Forni D, Cagliani R, Pontremoli C, Pozzoli U, Vertemara J, De Gioia L, Clerici M, Sironi M. Evolutionary Analysis Provides Insight Into the Origin and Adaptation of HCV. Front Microbiol 2018; 9:854. [PMID: 29765366 PMCID: PMC5938362 DOI: 10.3389/fmicb.2018.00854] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2018] [Accepted: 04/13/2018] [Indexed: 12/12/2022] Open
Abstract
Hepatitis C virus (HCV) belongs to the Hepacivirus genus and is genetically heterogeneous, with seven major genotypes further divided into several recognized subtypes. HCV origin was previously dated in a range between ∼200 and 1000 years ago. Hepaciviruses have been identified in several domestic and wild mammals, the largest viral diversity being observed in bats and rodents. The closest relatives of HCV were found in horses/donkeys (equine hepaciviruses, EHV). However, the origin of HCV as a human pathogen is still an unsolved puzzle. Using a selection-informed evolutionary model, we show that the common ancestor of extant HCV genotypes existed at least 3000 years ago (CI: 3192–5221 years ago), with the oldest genotypes being endemic to Asia. EHV originated around 1100 CE (CI: 291–1640 CE). These time estimates exclude that EHV transmission was mainly sustained by widespread veterinary practices and suggest that HCV originated from a single zoonotic event with subsequent diversification in human populations. We also describe a number of biologically important sites in the major HCV genotypes that have been positively selected and indicate that drug resistance-associated variants are significantly enriched at positively selected sites. HCV exploits several cell-surface molecules for cell entry, but only two of these (CD81 and OCLN) determine the species-specificity of infection. Herein evolutionary analyses do not support a long-standing association between primates and hepaciviruses, and signals of positive selection at CD81 were only observed in Chiroptera. No evidence of selection was detected for OCLN in any mammalian order. These results shed light on the origin of HCV and provide a catalog of candidate genetic modulators of HCV phenotypic diversity.
Collapse
Affiliation(s)
- Diego Forni
- Bioinformatics Laboratory, Scientific Institute IRCCS E.Medea, Bosisio Parini, Italy
| | - Rachele Cagliani
- Bioinformatics Laboratory, Scientific Institute IRCCS E.Medea, Bosisio Parini, Italy
| | - Chiara Pontremoli
- Bioinformatics Laboratory, Scientific Institute IRCCS E.Medea, Bosisio Parini, Italy
| | - Uberto Pozzoli
- Bioinformatics Laboratory, Scientific Institute IRCCS E.Medea, Bosisio Parini, Italy
| | - Jacopo Vertemara
- Department of Biotechnology and Biosciences, University of Milan-Bicocca, Milan, Italy
| | - Luca De Gioia
- Department of Biotechnology and Biosciences, University of Milan-Bicocca, Milan, Italy
| | - Mario Clerici
- Department of Pathophysiology and Transplantation, University of Milan, Milan, Italy.,Don C. Gnocchi Foundation Onlus, IRCCS, Milan, Italy
| | - Manuela Sironi
- Bioinformatics Laboratory, Scientific Institute IRCCS E.Medea, Bosisio Parini, Italy
| |
Collapse
|
20
|
Lima K, Leal É, Cavalcanti AMS, Salustiano DM, de Medeiros LB, da Silva SP, Lacerda HR. Increase in human immunodeficiency virus 1 diversity and detection of various subtypes and recombinants in north-eastern Brazil. J Med Microbiol 2017; 66:526-535. [PMID: 28425872 DOI: 10.1099/jmm.0.000447] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
PURPOSE Diverse human immunodeficiency virus 1 (HIV-1) subtypes and circulating recombinant forms are found in Brazil. The majority of HIV-1 molecular epidemiological studies in Brazil have been conducted in the southern and south-eastern regions of the country, although several recent studies in the north-eastern region have addressed this issue. The objective of this study was to molecularly characterize HIV-1 circulating in Pernambuco, north-eastern Brazil. METHODOLOGY A total of 64 samples were collected from 2002 to 2003, and another 103 were collected from 2007 to 2009. The protease and partial reverse transcriptase regions of the HIV-1 polymerase-encoding (pol) gene were sequenced, and subtyping, recombination and phylogenetic analyses were performed.Results/Key findings. Subtype B (60.9 %) was found to be predominant, followed by HIV-1 F (31.4 %). Several BF recombinants (4.2 %), and BC and AG recombinants were also identified. The intra-subtype genetic diversity was estimated to be 0.065 (sd±0.004) for HIV-1 B and 0.055 (sd±0.004) for HIV-1 F, reflecting a greater accumulation of mutations in subtype B (P<0.01). More codons were found to be under positive selective pressure in samples collected from 2007 to 2009, from individuals with a T-cell count≥200 cells mm-3 and from women. Coalescence data indicated that the subtype F population has been continuously expanding. CONCLUSIONS HIV-1 shows high genetic diversity in the state of Pernambuco. Thus, additional molecular evaluations of circulating strains will provide a better understanding of the epidemic and may lead to more effective preventive strategies.
Collapse
Affiliation(s)
- Kledoaldo Lima
- Tropical Medicine, Federal University of Pernambuco, Recife, PE, Brazil
| | - Élcio Leal
- Institute of Biotechnology, Federal University of Pará, Belém, PA, Brazil
| | | | | | | | | | | |
Collapse
|
21
|
Cicero MF, Pena NM, Santana LC, Arnold R, Azevedo RG, Leal ÉDS, Diaz RS, Komninakis SV. Is Hepatitis Delta infections important in Brazil? BMC Infect Dis 2016; 16:525. [PMID: 27686363 PMCID: PMC5041555 DOI: 10.1186/s12879-016-1856-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2015] [Accepted: 09/20/2016] [Indexed: 12/18/2022] Open
Abstract
Background The Hepatitis Delta Virus (HDV) can increase the incidence of fulminant hepatitis. For this infection occurs, the host must also be infected with Hepatitis B Virus. Previous studies demonstrated the endemicity and near exclusivity of this infection in the Amazon region, and as a consequence of the difficulty in accessing this area we used dried blood spots (DBS) in sample collection. The aims of this study were to investigate the presence of recombination, to analyze the epidemiology, ancestry and evolutionary pressures on HDV in Brazil. Methods Blood samples from 50 individuals were collected using dried-blood spots (DBS 903, Whatman), and sent via regular mail to Retrovirology Laboratory from Federal University of São Paulo, where the samples were processed. In the analysis the following software were used: PhyML, RDP, BEAST, jModelTest and CODEML. Results Our results confirm the prevalence of HDV-3 in the Amazon region of Brazil, with the absence of inter-genotypic recombination. It was identified a positive selection in probable epitopes of HDV on B lymphocytes that might indicate that the virus is changing to escape the humoral response of the host. The analysis of the time of the most common ancestor demonstrated the exponential growth of this virus in late 1970s that lasted until 1995, after which it remained constant. It was also observed a probable founder effect in two cities, which demonstrate the need to focus on prevention methods against HBV/HDV infection. Conclusion We confirmed the prevalence of HDV-3 in the Amazon region of Brazil, without inter-genotypic recombination. The analysis of the time of the most common ancestor showed that this infection remain constant in the studied area. Taking into account the probable founder effect established in the cities of Rio Branco and Porto Velho, a focus on preventive methods is recommended against these infections.
Collapse
Affiliation(s)
- Maira Ferreira Cicero
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Nathalia Mantovani Pena
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Luiz Claudio Santana
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Rafael Arnold
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Rafael Gonçalves Azevedo
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Élcio de Souza Leal
- Institute of Biological Sciences, Federal University of Pará, 01 Augusto Corrêa Street, Belém, Pará, Brazil
| | - Ricardo Sobhie Diaz
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil
| | - Shirley Vasconcelos Komninakis
- Retrovirology Laboratory, Federal University of São Paulo, 781 Pedro de Toledo Street, São Paulo, São Paulo, Brazil. .,School of Medicine of ABC (FMABC), Clinical Immunology Laboratory, 821 Príncipe de Gales Avenue, Santo André, São Paulo, Brazil.
| |
Collapse
|
22
|
Bielawski JP, Baker JL, Mingrone J. Inference of Episodic Changes in Natural Selection Acting on Protein Coding Sequences via CODEML. ACTA ACUST UNITED AC 2016; 54:6.15.1-6.15.32. [PMID: 27322407 DOI: 10.1002/cpbi.2] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This unit provides protocols for using the CODEML program from the PAML package to make inferences about episodic natural selection in protein-coding sequences. The protocols cover inference tasks such as maximum likelihood estimation of selection intensity, testing the hypothesis of episodic positive selection, and identifying sites with a history of episodic evolution. We provide protocols for using the rich set of models implemented in CODEML to assess robustness, and for using bootstrapping to assess if the requirements for reliable statistical inference have been met. An example dataset is used to illustrate how the protocols are used with real protein-coding sequences. The workflow of this design, through automation, is readily extendable to a larger-scale evolutionary survey. © 2016 by John Wiley & Sons, Inc.
Collapse
Affiliation(s)
- Joseph P Bielawski
- Department of Biology, Dalhousie University, Halifax, Nova Scotia, Canada.,Department of Mathematics & Statistics, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Jennifer L Baker
- Center for Research on Genomics and Global Health, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland
| | - Joseph Mingrone
- Department of Mathematics & Statistics, Dalhousie University, Halifax, Nova Scotia, Canada
| |
Collapse
|
23
|
Analysis of full-length genomes of porcine teschovirus (PTV) and the effect of purifying selection on phylogenetic trees. Arch Virol 2016; 161:1199-208. [DOI: 10.1007/s00705-015-2744-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2015] [Accepted: 12/21/2015] [Indexed: 10/22/2022]
|
24
|
Levy Karin E, Rabin A, Ashkenazy H, Shkedy D, Avram O, Cartwright RA, Pupko T. Inferring Indel Parameters using a Simulation-based Approach. Genome Biol Evol 2015; 7:3226-38. [PMID: 26537226 PMCID: PMC4700945 DOI: 10.1093/gbe/evv212] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.
Collapse
Affiliation(s)
- Eli Levy Karin
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
| | - Avigayel Rabin
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
| | - Haim Ashkenazy
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
| | - Dafna Shkedy
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
| | - Oren Avram
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel The Blavatnik School of Computer Science, Tel-Aviv University, Tel-Aviv, Israel
| | - Reed A Cartwright
- The Biodesign Institute, Arizona State University, Tempe School of Life Sciences, Arizona State University, Tempe
| | - Tal Pupko
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
| |
Collapse
|
25
|
Woegerbauer M, Kuffner M, Domingues S, Nielsen KM. Involvement of aph(3')-IIa in the formation of mosaic aminoglycoside resistance genes in natural environments. Front Microbiol 2015; 6:442. [PMID: 26042098 PMCID: PMC4437187 DOI: 10.3389/fmicb.2015.00442] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2015] [Accepted: 04/24/2015] [Indexed: 11/13/2022] Open
Abstract
Intragenic recombination leading to mosaic gene formation is known to alter resistance profiles for particular genes and bacterial species. Few studies have examined to what extent aminoglycoside resistance genes undergo intragenic recombination. We screened the GenBank database for mosaic gene formation in homologs of the aph(3')-IIa (nptII) gene. APH(3')-IIa inactivates important aminoglycoside antibiotics. The gene is widely used as a selectable marker in biotechnology and enters the environment via laboratory discharges and the release of transgenic organisms. Such releases may provide opportunities for recombination in competent environmental bacteria. The retrieved GenBank sequences were grouped in three datasets comprising river water samples, duck pathogens and full-length variants from various bacterial genomes and plasmids. Analysis for recombination in these datasets was performed with the Recombination Detection Program (RDP4), and the Genetic Algorithm for Recombination Detection (GARD). From a total of 89 homologous sequences, 83% showed 99-100% sequence identity with aph(3')-IIa originally described as part of transposon Tn5. Fifty one were unique sequence variants eligible for recombination analysis. Only a single recombination event was identified with high confidence and indicated the involvement of aph(3')-IIa in the formation of a mosaic gene located on a plasmid of environmental origin in the multi-resistant isolate Pseudomonas aeruginosa PA96. The available data suggest that aph(3')-IIa is not an archetypical mosaic gene as the divergence between the described sequence variants and the number of detectable recombination events is low. This is in contrast to the numerous mosaic alleles reported for certain penicillin or tetracycline resistance determinants.
Collapse
Affiliation(s)
- Markus Woegerbauer
- Integrative Risk Assessment - Data - Statistics, GMO Risk Assessment, Austrian Agency for Health and Food Safety Vienna, Austria
| | - Melanie Kuffner
- Integrative Risk Assessment - Data - Statistics, GMO Risk Assessment, Austrian Agency for Health and Food Safety Vienna, Austria
| | - Sara Domingues
- Faculty of Pharmacy and Center for Neuroscience and Cell Biology, University of Coimbra Coimbra, Portugal
| | - Kaare M Nielsen
- Department of Pharmacy, University of Tromsø Tromsø, Norway ; Genøk-Center for Biosafety Tromsø Tromsø, Norway
| |
Collapse
|
26
|
Pérez R, Calleros L, Marandino A, Sarute N, Iraola G, Grecco S, Blanc H, Vignuzzi M, Isakov O, Shomron N, Carrau L, Hernández M, Francia L, Sosa K, Tomás G, Panzera Y. Phylogenetic and genome-wide deep-sequencing analyses of canine parvovirus reveal co-infection with field variants and emergence of a recent recombinant strain. PLoS One 2014; 9:e111779. [PMID: 25365348 PMCID: PMC4218814 DOI: 10.1371/journal.pone.0111779] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2014] [Accepted: 09/30/2014] [Indexed: 11/28/2022] Open
Abstract
Canine parvovirus (CPV), a fast-evolving single-stranded DNA virus, comprises three antigenic variants (2a, 2b, and 2c) with different frequencies and genetic variability among countries. The contribution of co-infection and recombination to the genetic variability of CPV is far from being fully elucidated. Here we took advantage of a natural CPV population, recently formed by the convergence of divergent CPV-2c and CPV-2a strains, to study co-infection and recombination. Complete sequences of the viral coding region of CPV-2a and CPV-2c strains from 40 samples were generated and analyzed using phylogenetic tools. Two samples showed co-infection and were further analyzed by deep sequencing. The sequence profile of one of the samples revealed the presence of CPV-2c and CPV-2a strains that differed at 29 nucleotides. The other sample included a minor CPV-2a strain (13.3% of the viral population) and a major recombinant strain (86.7%). The recombinant strain arose from inter-genotypic recombination between CPV-2c and CPV-2a strains within the VP1/VP2 gene boundary. Our findings highlight the importance of deep-sequencing analysis to provide a better understanding of CPV molecular diversity.
Collapse
Affiliation(s)
- Ruben Pérez
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
- * E-mail:
| | - Lucía Calleros
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Ana Marandino
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Nicolás Sarute
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Gregorio Iraola
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Sofia Grecco
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Hervé Blanc
- Institut Pasteur, Viral Populations and Pathogenesis Unit, Centre National de la Recherche Scientifique, Paris, France
| | - Marco Vignuzzi
- Institut Pasteur, Viral Populations and Pathogenesis Unit, Centre National de la Recherche Scientifique, Paris, France
| | - Ofer Isakov
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Noam Shomron
- Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Lucía Carrau
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Martín Hernández
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Lourdes Francia
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Katia Sosa
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Gonzalo Tomás
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Yanina Panzera
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| |
Collapse
|
27
|
RIG-I-Like Receptors Evolved Adaptively in Mammals, with Parallel Evolution at LGP2 and RIG-I. J Mol Biol 2014; 426:1351-65. [DOI: 10.1016/j.jmb.2013.10.040] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2013] [Revised: 10/11/2013] [Accepted: 10/30/2013] [Indexed: 01/18/2023]
|
28
|
Hester SE, Park J, Goodfield LL, Feaga HA, Preston A, Harvill ET. Horizontally acquired divergent O-antigen contributes to escape from cross-immunity in the classical bordetellae. BMC Evol Biol 2013; 13:209. [PMID: 24067113 PMCID: PMC3849452 DOI: 10.1186/1471-2148-13-209] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2013] [Accepted: 09/13/2013] [Indexed: 11/10/2022] Open
Abstract
Background Horizontal gene transfer (HGT) allows for rapid spread of genetic material between species, increasing genetic and phenotypic diversity. Although HGT contributes to adaptation and is widespread in many bacteria, others show little HGT. This study builds on previous work to analyze the evolutionary mechanisms contributing to variation within the locus encoding a prominent antigen of the classical bordetellae. Results We observed amongst classical bordetellae discrete regions of the lipopolysaccharide O-antigen locus with higher sequence diversity than the genome average. Regions of this locus had less than 50% sequence similarity, low dN/dS ratios and lower GC content compared to the genome average. Additionally, phylogenetic tree topologies based on genome-wide SNPs were incongruent with those based on genes within these variable regions, suggesting portions of the O-antigen locus may have been horizontally transferred. Furthermore, several predicted recombination breakpoints correspond with the ends of these variable regions. To examine the evolutionary forces that might have selected for this rare example of HGT in bordetellae, we compared in vitro and in vivo phenotypes associated with different O-antigen types. Antibodies against O1- and O2-serotypes were poorly cross-reactive, and did not efficiently kill or mediate clearance of alternative O-type bacteria, while a distinct and poorly immunogenic O-antigen offered no protection against colonization. Conclusions This study suggests that O-antigen variation was introduced to the classical bordetellae via HGT through recombination. Additionally, genetic variation may be maintained within the O-antigen locus because it can provide escape from immunity to different O-antigen types, potentially allowing for the circulation of different Bordetella strains within the same host population.
Collapse
Affiliation(s)
- Sara E Hester
- Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, W-210 Millennium Science Complex, University Park, PA, 16802, USA.
| | | | | | | | | | | |
Collapse
|
29
|
Privman E, Wurm Y, Keller L. Duplication and concerted evolution in a master sex determiner under balancing selection. Proc Biol Sci 2013; 280:20122968. [PMID: 23466984 DOI: 10.1098/rspb.2012.2968] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.
Collapse
Affiliation(s)
- Eyal Privman
- Department of Ecology and Evolution, University of Lausanne, Lausanne 1015, Switzerland
| | | | | |
Collapse
|
30
|
Inference of functional divergence among proteins when the evolutionary process is non-stationary. J Mol Evol 2013; 76:205-15. [PMID: 23443835 DOI: 10.1007/s00239-013-9549-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2012] [Accepted: 02/09/2013] [Indexed: 10/27/2022]
Abstract
Functional shifts during protein evolution are expected to yield shifts in substitution rate, and statistical methods can test for this at both codon and amino acid levels. Although methods based on models of sequence evolution serve as powerful tools for studying evolutionary processes, violating underlying assumptions can lead to false biological conclusions. It is not unusual for functional shifts to be accompanied by changes in other aspects of the evolutionary process, such as codon or amino acid frequencies. However, models used to test for functional divergence assume these frequencies remain constant over time. We employed simulation to investigate the impact of non-stationary evolution on functional divergence inference. We investigated three likelihood ratio tests based on codon models and found varying degrees of sensitivity. Joint effects of shifts in frequencies and selection pressures can be large, leading to false signals for positive selection. Amino acid-based tests (FunDi and Bivar) were also compromised when several aspects of the substitution process were not adequately modeled. We applied the same tests to a core genome "scan" for functional divergence between light-adapted ecotypes of the cyanobacteria Prochlorococcus, and carried out gene-specific simulations for ten genes. Results of those simulations illustrated how the inference of functional divergence at the genomic level can be seriously impacted by model misspecification. Although computationally costly, simulations motivated by data in hand are warranted when several aspects of the substitution process are either misspecified or not included in the models upon which the statistical tests were built.
Collapse
|
31
|
Leal É, Villanova FE, Lin W, Hu F, Liu Q, Liu Y, Cui S. Interclade recombination in porcine parvovirus strains. J Gen Virol 2012; 93:2692-2704. [DOI: 10.1099/vir.0.045765-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/30/2023] Open
Abstract
A detailed analysis of the Ns1/Vp1Vp2 genome region of the porcine parvovirus (PPV) strains isolated from vaccinated animals was performed. We found many inconsistencies in the phylogenetic trees of these viral isolates, such as low statistical support and strains with long branches in the phylogenetic trees. Thus, we used distance-based and phylogenetic methods to distinguish de facto recombinants from spurious recombination signals. We found a mosaic virus in which the Ns1 gene was acquired from one PPV clade and the Vp1Vp2 gene was acquired from a distinct phylogenetic clade. We also described the interclade mosaic structure of the Vp1Vp2 gene of a reference strain. If recombination is an adaptive mechanism over the course of PPV evolution, we would likely observe increasing numbers of chimeric strains over time. However, when the PPV sequences isolated from 1964 to 2011 were analysed, only two chimeric strains were detected. Thus, PPV recombination is an independent event, resulting from close contact between animals housed in high-density conditions.
Collapse
Affiliation(s)
- Élcio Leal
- Federal University of Pará, Belém, Brazil
| | | | - Wencheng Lin
- Division of Swine Infectious Diseases, State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute of CAAS, Heilongjiang, PR China
| | - Feng Hu
- Division of Swine Infectious Diseases, State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute of CAAS, Heilongjiang, PR China
| | - Qinfang Liu
- Division of Swine Infectious Diseases, State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute of CAAS, Heilongjiang, PR China
| | - Yebing Liu
- China Institute of Veterinary Drug Control, Beijing 100081, PR China
| | - Shangjin Cui
- Division of Swine Infectious Diseases, State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute of CAAS, Heilongjiang, PR China
| |
Collapse
|