1
|
King KA, Benton AH, Caudill MT, Stoyanof ST, Kang L, Michalak P, Lahmers KK, Dunman PM, DeHart TG, Ahmad SS, Jutras BL, Poncin K, De Bolle X, Caswell CC. Post-transcriptional control of the essential enzyme MurF by a small regulatory RNA in Brucella abortus. Mol Microbiol 2024; 121:129-141. [PMID: 38082493 DOI: 10.1111/mmi.15207] [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: 09/27/2023] [Revised: 11/22/2023] [Accepted: 11/25/2023] [Indexed: 01/15/2024]
Abstract
Brucella abortus is a facultative, intracellular, zoonotic pathogen that resides inside macrophages during infection. This is a specialized niche where B. abortus encounters various stresses as it navigates through the macrophage. In order to survive this harsh environment, B. abortus utilizes post-transcriptional regulation of gene expression through the use of small regulatory RNAs (sRNAs). Here, we characterize a Brucella sRNAs called MavR (for MurF- and virulence-regulating sRNA), and we demonstrate that MavR is required for the full virulence of B. abortus in macrophages and in a mouse model of chronic infection. Transcriptomic and proteomic studies revealed that a major regulatory target of MavR is MurF. MurF is an essential protein that catalyzes the final cytoplasmic step in peptidoglycan (PG) synthesis; however, we did not detect any differences in the amount or chemical composition of PG in the ΔmavR mutant. A 6-nucleotide regulatory seed region within MavR was identified, and mutation of this seed region resulted in dysregulation of MurF production, as well as significant attenuation of infection in a mouse model. Overall, the present study underscores the importance of sRNA regulation in the physiology and virulence of Brucella.
Collapse
|
2
|
Cereghino C, Roesch F, Carrau L, Hardy A, Ribeiro-Filho HV, Henrion-Lacritick A, Koh C, Marano JM, Bates TA, Rai P, Chuong C, Akter S, Vallet T, Blanc H, Elliott TJ, Brown AM, Michalak P, LeRoith T, Bloom JD, Marques RE, Saleh MC, Vignuzzi M, Weger-Lucarelli J. The E2 glycoprotein holds key residues for Mayaro virus adaptation to the urban Aedes aegypti mosquito. PLoS Pathog 2023; 19:e1010491. [PMID: 37018377 PMCID: PMC10109513 DOI: 10.1371/journal.ppat.1010491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 04/17/2023] [Accepted: 03/13/2023] [Indexed: 04/06/2023] Open
Abstract
Adaptation to mosquito vectors suited for transmission in urban settings is a major driver in the emergence of arboviruses. To better anticipate future emergence events, it is crucial to assess their potential to adapt to new vector hosts. In this work, we used two different experimental evolution approaches to study the adaptation process of an emerging alphavirus, Mayaro virus (MAYV), to Ae. aegypti, an urban mosquito vector of many other arboviruses. We identified E2-T179N as a key mutation increasing MAYV replication in insect cells and enhancing transmission after escaping the midgut of live Ae. aegypti. In contrast, this mutation decreased viral replication and binding in human fibroblasts, a primary cellular target of MAYV in humans. We also showed that MAYV E2-T179N generates reduced viremia and displays less severe tissue pathology in vivo in a mouse model. We found evidence in mouse fibroblasts that MAYV E2-T179N is less dependent on the Mxra8 receptor for replication than WT MAYV. Similarly, exogenous expression of human apolipoprotein receptor 2 and Mxra8 enhanced WT MAYV replication compared to MAYV E2-T179N. When this mutation was introduced in the closely related chikungunya virus, which has caused major outbreaks globally in the past two decades, we observed increased replication in both human and insect cells, suggesting E2 position 179 is an important determinant of alphavirus host-adaptation, although in a virus-specific manner. Collectively, these results indicate that adaptation at the T179 residue in MAYV E2 may result in increased vector competence-but coming at the cost of optimal replication in humans-and may represent a first step towards a future emergence event.
Collapse
|
3
|
Porier DL, Adam A, Kang L, Michalak P, Tupik J, Santos MA, Lee C, Allen IC, Wang T, Auguste AJ. Humoral and T-cell-mediated responses to a pre-clinical Zika vaccine candidate that utilizes a unique insect-specific flavivirus platform. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.01.530296. [PMID: 36909623 PMCID: PMC10002724 DOI: 10.1101/2023.03.01.530296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
Abstract
Vaccination is critical for the control and prevention of viral outbreaks, yet conventional vaccine platforms may involve trade-offs between immunogenicity and safety. Insect-specific viruses have emerged as a novel vaccine platform to overcome this challenge. Detailed studies of humoral and T-cell responses induced by new insect-specific flavivirus (ISFV)-based vaccine platforms are needed to better understand correlates of protection and improve vaccine efficacy. Previously, we used a novel ISFV called Aripo virus (ARPV) to create a Zika virus (ZIKV) vaccine candidate (designated ARPV/ZIKV). ARPV/ZIKV demonstrated exceptional safety and single-dose efficacy, completely protecting mice from a lethal ZIKV challenge. Here, we explore the development of immune responses induced by ARPV/ZIKV immunization and evaluate its correlates of protection. Passive transfer of ARPV/ZIKV-induced immune sera to naïve mice prior to challenge emphasized the importance of neutralizing antibodies as a correlate of protection. Depletion of T-cells in vaccinated mice and adoptive transfer of ARPV/ZIKV-primed T-cells to naïve mice prior to challenge indicated that ARPV/ZIKV-induced CD4 + and CD8 + T-cell responses contribute to the observed protection but may not be essential for protection during ZIKV challenge. However, vaccination of Rag1 KO, Tcra KO, and muMt - mice demonstrated the critical role for ARPV/ZIKV-induced T-cells in developing protective immune responses following vaccination. Overall, both humoral and T-cell-mediated responses induced by ISFV-based vaccines are important for comprehensive immunity, and ISFV platforms continue to be a promising method for future vaccine development.
Collapse
|
4
|
Deeley S, Kang L, Michalak P, Hallerman E, Ford WM. DNA Metabarcoding-Based Evaluation of the Diet of Big Brown Bats (Eptesicus Fuscus) in the Mid-Atlantic Region. Northeast Nat (Steuben) 2023. [DOI: 10.1656/045.029.0405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
5
|
Oakes VJ, Todd SM, Carbonello AA, Michalak P, Lahmers KK. Coinfection of cattle in Virginia with Theileria orientalis Ikeda genotype and Anaplasma marginale. J Vet Diagn Invest 2021; 34:36-41. [PMID: 34763583 DOI: 10.1177/10406387211057627] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Theileria orientalis Ikeda is a newly identified agent of bovine infectious anemia in the United States. Although T. orientalis Ikeda is transmitted by ticks other than the tick that transmits Anaplasma marginale-a bacterial etiology of bovine infectious anemia-the geographic distributions of these 2 infectious organisms overlap, with coinfection reported in some cattle. Only anaplasmosis has an approved effective treatment in the United States. To provide rapid diagnostic information for producers with anemic animals, we developed a duplex real-time PCR (rtPCR) for A. marginale and T. orientalis. With a cutoff of 38 cycles, the duplex assay has a sensitivity of 97.0% and a specificity of 100% for A. marginale; with a cutoff of 45 cycles, the duplex assay has a sensitivity and a specificity of 100% for T. orientalis, compared to existing tests. In addition to providing a tool for improved clinical decision-making for veterinarians and producers, our rtPCR facilitates the study of coinfection of cattle in Virginia. Of 1,359 blood samples analyzed, 174 were positive for T. orientalis, 125 were positive for A. marginale, and 12 samples were positive for both T. orientalis and A. marginale. Hence, coinfection by these 2 agents of bovine infectious anemia does occur within Virginia. It is likely that this pattern of infection will be seen in other regions where T. orientalis and A. marginale infections are endemic, despite the difference in tick vectors.
Collapse
|
6
|
Kang L, Michalak P, Hallerman E, Moncrief ND. A Draft Genome Assembly for the Eastern Fox Squirrel, Sciurus niger. G3-GENES GENOMES GENETICS 2021; 11:6373894. [PMID: 34550334 PMCID: PMC8664420 DOI: 10.1093/g3journal/jkab315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Accepted: 08/27/2021] [Indexed: 11/12/2022]
Abstract
The eastern fox squirrel, Sciurus niger, exhibits marked geographic variation in size and coat color, is a model organism for studies of behavior and ecology, and a potential model for investigating physiological solutions to human porphyrias. We assembled a genome using Illumina HiSeq, PacBio SMRT, and Oxford Nanopore MinION sequencing platforms. Together, the sequencing data resulted in a draft genome of 2.99 Gb, containing 32,830 scaffolds with an average size of 90.9Kb and N50 of 183.8 Kb. Genome completeness was estimated to be 93.78%. A total of 24,443 protein-encoding genes were predicted from the assembly, and 23,079 (94.42%) were annotated. Repeat elements comprised an estimated 38.49% of the genome, with the majority being LINEs (13.92%), SINEs (6.04%), and LTR elements. The topology of the species tree reconstructed using maximum likelihood phylogenetic analysis was congruent with those of previous studies. This genome assembly can prove useful for comparative studies of genome structure and function in this rapidly diversifying lineage of mammals, for studies of population genomics and adaptation, and for biomedical research. Predicted amino acid sequence alignments for genes affecting heme biosynthesis, color vision and hibernation showed point mutations and indels that may affect protein function and ecological adaptation.
Collapse
|
7
|
Aggarwal DD, Rybnikov S, Sapielkin S, Rashkovetsky E, Frenkel Z, Singh M, Michalak P, Korol AB. Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster. Heredity (Edinb) 2021; 127:278-287. [PMID: 34163036 PMCID: PMC8405755 DOI: 10.1038/s41437-021-00449-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Revised: 06/07/2021] [Accepted: 06/07/2021] [Indexed: 02/06/2023] Open
Abstract
Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).
Collapse
|
8
|
Kang L, He G, Sharp AK, Wang X, Brown AM, Michalak P, Weger-Lucarelli J. A selective sweep in the Spike gene has driven SARS-CoV-2 human adaptation. Cell 2021; 184:4392-4400.e4. [PMID: 34289344 PMCID: PMC8260498 DOI: 10.1016/j.cell.2021.07.007] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Revised: 05/07/2021] [Accepted: 07/02/2021] [Indexed: 12/15/2022]
Abstract
The coronavirus disease 2019 (COVID-19) pandemic underscores the need to better understand animal-to-human transmission of coronaviruses and adaptive evolution within new hosts. We scanned more than 182,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes for selective sweep signatures and found a distinct footprint of positive selection located around a non-synonymous change (A1114G; T372A) within the spike protein receptor-binding domain (RBD), predicted to remove glycosylation and increase binding to human ACE2 (hACE2), the cellular receptor. This change is present in all human SARS-CoV-2 sequences but not in closely related viruses from bats and pangolins. As predicted, T372A RBD bound hACE2 with higher affinity in experimental binding assays. We engineered the reversion mutant (A372T) and found that A372 (wild-type [WT]-SARS-CoV-2) enhanced replication in human lung cells relative to its putative ancestral variant (T372), an effect that was 20 times greater than the well-known D614G mutation. Our findings suggest that this mutation likely contributed to SARS-CoV-2 emergence from animal reservoirs or enabled sustained human-to-human transmission.
Collapse
|
9
|
Kinney N, Kang L, Bains H, Lawson E, Husain M, Husain K, Sandhu I, Shin Y, Carter JK, Anandakrishnan R, Michalak P, Garner H. Ethnically biased microsatellites contribute to differential gene expression and glutathione metabolism in Africans and Europeans. PLoS One 2021; 16:e0249148. [PMID: 33765058 PMCID: PMC7993785 DOI: 10.1371/journal.pone.0249148] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 03/11/2021] [Indexed: 12/28/2022] Open
Abstract
Approximately three percent of the human genome is occupied by microsatellites: a type of short tandem repeat (STR). Microsatellites have well established effects on (a) the genetic structure of diverse human populations and (b) expression of nearby genes. These lines of inquiry have uncovered 3,984 ethnically biased microsatellite loci (EBML) and 28,375 expression STRs (eSTRs), respectively. We hypothesize that a combination of EBML, eSTRs, and gene expression data (RNA-seq) can be used to show that microsatellites contribute to differential gene expression and phenotype in human populations. In fact, our previous study demonstrated a degree of mutual overlap between EBML and eSTRs but fell short of quantifying effects on gene expression. The present work aims to narrow the gap. First, we identify 313 overlapping EBML/eSTRs and recapitulate their mutual overlap. The 313 EBML/eSTRs are then characterized across ethnicity and tissue type. We use RNA-seq data to pursue validation of 49 regions that affect whole blood gene expression; 32 out of 54 affected genes are differentially expressed in Africans and Europeans. We quantify the relative contribution of these 32 genes to differential expression; fold change tends to be less than other differentially expressed genes. Repeat length correlates with expression for 15 of the 32 genes; two are conspicuously involved in glutathione metabolism. Finally, we repurpose a mathematical model of glutathione metabolism to investigate how a single polymorphic microsatellite affects phenotype. We conclude with a testable prediction that microsatellite polymorphisms affect GPX7 expression and oxidative stress in Africans and Europeans.
Collapse
|
10
|
Anandakrishnan R, Carpenetti TL, Samuel P, Wasko B, Johnson C, Smith C, Kim J, Michalak P, Kang L, Kinney N, Santo A, Anstrom J, Garner HR, Varghese RT. DNA sequencing of anatomy lab cadavers to provide hands-on precision medicine introduction to medical students. BMC MEDICAL EDUCATION 2020; 20:437. [PMID: 33198737 PMCID: PMC7670733 DOI: 10.1186/s12909-020-02366-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Accepted: 11/09/2020] [Indexed: 05/23/2023]
Abstract
BACKGROUND Medical treatment informed by Precision Medicine is becoming a standard practice for many diseases, and patients are curious about the consequences of genomic variants in their genome. However, most medical students' understanding of Precision Medicine derives from classroom lectures. This format does little to foster an understanding for the potential and limitations of Precision Medicine. To close this gap, we implemented a hands-on Precision Medicine training program utilizing exome sequencing to prepare a clinical genetic report of cadavers studied in the anatomy lab. The program reinforces Precision Medicine related learning objectives for the Genetics curriculum. METHODS Pre-embalmed blood samples and embalmed tissue were obtained from cadavers (donors) used in the anatomy lab. DNA was isolated and sequenced and illustrative genetic reports provided to the students. The reports were used to facilitate discussion with students on the implications of pathogenic genomic variants and the potential correlation of these variants in each "donor" with any anatomical anomalies identified during cadaver dissection. RESULTS In 75% of cases, analysis of whole exome sequencing data identified a variant associated with increased risk for a disease/abnormal condition noted in the donor's cause of death or in the students' anatomical findings. This provided students with real-world examples of the potential relationship between genomic variants and disease risk. Our students also noted that diseases associated with 92% of the pathogenic variants identified were not related to the anatomical findings, demonstrating the limitations of Precision Medicine. CONCLUSION With this study, we have established protocols and classroom procedures incorporating hands-on Precision Medicine training in the medical student curriculum and a template for other medical educators interested in enhancing their Precision Medicine training program. The program engaged students in discovering variants that were associated with the pathophysiology of the cadaver they were studying, which led to more exposure and understanding of the potential risks and benefits of genomic medicine.
Collapse
|
11
|
Budnick JA, Sheehan LM, Ginder MJ, Failor KC, Perkowski JM, Pinto JF, Kohl KA, Kang L, Michalak P, Luo L, Heindl JE, Caswell CC. A central role for the transcriptional regulator VtlR in small RNA-mediated gene regulation in Agrobacterium tumefaciens. Sci Rep 2020; 10:14968. [PMID: 32917931 PMCID: PMC7486931 DOI: 10.1038/s41598-020-72117-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2020] [Accepted: 07/17/2020] [Indexed: 12/28/2022] Open
Abstract
LysR-type transcriptional regulators (LTTRs) are the most common type of transcriptional regulators in prokaryotes and function by altering gene expression in response to environmental stimuli. In the class Alphaproteobacteria, a conserved LTTR named VtlR is critical to the establishment of host-microbe interactions. In the mammalian pathogen Brucella abortus, VtlR is required for full virulence in a mouse model of infection, and VtlR activates the expression of abcR2, which encodes a small regulatory RNA (sRNA). In the plant symbiont Sinorhizobium meliloti, the ortholog of VtlR, named LsrB, is involved in the symbiosis of the bacterium with alfalfa. Agrobacterium tumefaciens is a close relative of both B. abortus and S. meliloti, and this bacterium is the causative agent of crown gall disease in plants. In the present study, we demonstrate that VtlR is involved in the ability of A. tumefaciens to grow appropriately in artificial medium, and an A. tumefaciens vtlR deletion strain is defective in motility, biofilm formation, and tumorigenesis of potato discs. RNA-sequencing analyses revealed that more than 250 genes are dysregulated in the ∆vtlR strain, and importantly, VtlR directly controls the expression of three sRNAs in A. tumefaciens. Taken together, these data support a model in which VtlR indirectly regulates hundreds of genes via manipulation of sRNA pathways in A. tumefaciens, and moreover, while the VtlR/LsrB protein is present and structurally conserved in many members of the Alphaproteobacteria, the VtlR/LsrB regulatory circuitry has diverged in order to accommodate the unique environmental niche of each organism.
Collapse
|
12
|
Pinzari CA, Kang L, Michalak P, Jermiin LS, Price DK, Bonaccorso FJ. Analysis of Genomic Sequence Data Reveals the Origin and Evolutionary Separation of Hawaiian Hoary Bat Populations. Genome Biol Evol 2020; 12:1504-1514. [PMID: 32853363 PMCID: PMC7543519 DOI: 10.1093/gbe/evaa137] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
We examine the genetic history and population status of Hawaiian hoary bats (Lasiurus semotus), the most isolated bats on Earth, and their relationship to northern hoary bats (Lasiurus cinereus), through whole-genome analysis of single-nucleotide polymorphisms mapped to a de novo-assembled reference genome. Profiles of genomic diversity and divergence indicate that Hawaiian hoary bats are distinct from northern hoary bats, and form a monophyletic group, indicating a single ancestral colonization event 1.34 Ma, followed by substantial divergence between islands beginning 0.51 Ma. Phylogenetic analysis indicates Maui is central to the radiation across the archipelago, with the southward expansion to Hawai'i and westward to O'ahu and Kaua'i. Because this endangered species is of conservation concern, a clearer understanding of the population genetic structure of this bat in the Hawaiian Islands is of timely importance.
Collapse
|
13
|
Budnick JA, Sheehan LM, Benton AH, Pitzer JE, Kang L, Michalak P, Roop RM, Caswell CC. Characterizing the transport and utilization of the neurotransmitter GABA in the bacterial pathogen Brucella abortus. PLoS One 2020; 15:e0237371. [PMID: 32845904 PMCID: PMC7449393 DOI: 10.1371/journal.pone.0237371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 08/04/2020] [Indexed: 01/18/2023] Open
Abstract
The neurotransmitter gamma-aminobutyric acid (GABA) is the most abundant inhibitory neurotransmitter in the human brain; however, it is becoming more evident that this non-proteinogenic amino acid plays multiple physiological roles in biology. In the present study, the transport and function of GABA is studied in the highly infectious intracellular bacterium Brucella abortus. The data show that 3H-GABA is imported by B. abortus under nutrient limiting conditions and that the small RNAs AbcR1 and AbcR2 negatively regulate this transport. A specific transport system, gts, is responsible for the transport of GABA as determined by measuring 3H-GABA transport in isogenic deletion strains of known AbcR1/2 regulatory targets; however, this locus is unnecessary for Brucella infection in BALB/c mice. Similar assays revealed that 3H-GABA transport is uninhibited by the 20 standard proteinogenic amino acids, representing preference for the transport of 3H-GABA. Metabolic studies did not show any potential metabolic utilization of GABA by B. abortus as a carbon or nitrogen source, and RNA sequencing analysis revealed limited transcriptional differences between B. abortus 2308 with or without exposure to GABA. While this study provides evidence for GABA transport by B. abortus, questions remain as to why and when this transport is utilized during Brucella pathogenesis.
Collapse
|
14
|
Gemmell NJ, Rutherford K, Prost S, Tollis M, Winter D, Macey JR, Adelson DL, Suh A, Bertozzi T, Grau JH, Organ C, Gardner PP, Muffato M, Patricio M, Billis K, Martin FJ, Flicek P, Petersen B, Kang L, Michalak P, Buckley TR, Wilson M, Cheng Y, Miller H, Schott RK, Jordan MD, Newcomb RD, Arroyo JI, Valenzuela N, Hore TA, Renart J, Peona V, Peart CR, Warmuth VM, Zeng L, Kortschak RD, Raison JM, Zapata VV, Wu Z, Santesmasses D, Mariotti M, Guigó R, Rupp SM, Twort VG, Dussex N, Taylor H, Abe H, Bond DM, Paterson JM, Mulcahy DG, Gonzalez VL, Barbieri CG, DeMeo DP, Pabinger S, Van Stijn T, Clarke S, Ryder O, Edwards SV, Salzberg SL, Anderson L, Nelson N, Stone C. The tuatara genome reveals ancient features of amniote evolution. Nature 2020; 584:403-409. [PMID: 32760000 PMCID: PMC7116210 DOI: 10.1038/s41586-020-2561-9] [Citation(s) in RCA: 78] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Accepted: 06/26/2020] [Indexed: 12/21/2022]
Abstract
The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.
Collapse
|
15
|
Kinney N, Kang L, Eckstrand L, Pulenthiran A, Samuel P, Anandakrishnan R, Varghese RT, Michalak P, Garner HR. Abundance of ethnically biased microsatellites in human gene regions. PLoS One 2019; 14:e0225216. [PMID: 31830051 PMCID: PMC6907796 DOI: 10.1371/journal.pone.0225216] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 10/29/2019] [Indexed: 12/16/2022] Open
Abstract
Microsatellites-a type of short tandem repeat (STR)-have been used for decades as putatively neutral markers to study the genetic structure of diverse human populations. However, recent studies have demonstrated that some microsatellites contribute to gene expression, cis heritability, and phenotype. As a corollary, some microsatellites may contribute to differential gene expression and RNA/protein structure stability in distinct human populations. To test this hypothesis, we investigate genotype frequencies, functional relevance, and adaptive potential of microsatellites in five super-populations (ethnicities) drawn from the 1000 Genomes Project. We discover 3,984 ethnically-biased microsatellite loci (EBML); for each EBML at least one ethnicity has genotype frequencies statistically different from the remaining four. South Asian, East Asian, European, and American EBML show significant overlap; on the contrary, the set of African EBML is mostly unique. We cross-reference the 3,984 EBML with 2,060 previously identified expression STRs (eSTRs); repeats known to affect gene expression (64 total) are over-represented. The most significant pathway enrichments are those associated with the matrisome: a broad collection of genes encoding the extracellular matrix and its associated proteins. At least 14 of the EBML have established links to human disease. Analysis of the 3,984 EBML with respect to known selective sweep regions in the genome shows that allelic variation in some of them is likely associated with adaptive evolution.
Collapse
|
16
|
van Hoek ML, Prickett MD, Settlage RE, Kang L, Michalak P, Vliet KA, Bishop BM. The Komodo dragon (Varanus komodoensis) genome and identification of innate immunity genes and clusters. BMC Genomics 2019; 20:684. [PMID: 31470795 PMCID: PMC6716921 DOI: 10.1186/s12864-019-6029-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2018] [Accepted: 08/12/2019] [Indexed: 12/23/2022] Open
Abstract
Background We report the sequencing, assembly and analysis of the genome of the Komodo dragon (Varanus komodoensis), the largest extant lizard, with a focus on antimicrobial host-defense peptides. The Komodo dragon diet includes carrion, and a complex milieu of bacteria, including potentially pathogenic strains, has been detected in the saliva of wild dragons. They appear to be unaffected, suggesting that dragons have robust defenses against infection. While little information is available regarding the molecular biology of reptile immunity, it is believed that innate immunity, which employs antimicrobial host-defense peptides including defensins and cathelicidins, plays a more prominent role in reptile immunity than it does in mammals. . Results High molecular weight genomic DNA was extracted from Komodo dragon blood cells. Subsequent sequencing and assembly of the genome from the collected DNA yielded a genome size of 1.6 Gb with 45x coverage, and the identification of 17,213 predicted genes. Through further analyses of the genome, we identified genes and gene-clusters corresponding to antimicrobial host-defense peptide genes. Multiple β-defensin-related gene clusters were identified, as well as a cluster of potential Komodo dragon ovodefensin genes located in close proximity to a cluster of Komodo dragon β-defensin genes. In addition to these defensins, multiple cathelicidin-like genes were also identified in the genome. Overall, 66 β-defensin genes, six ovodefensin genes and three cathelicidin genes were identified in the Komodo dragon genome. Conclusions Genes with important roles in host-defense and innate immunity were identified in this newly sequenced Komodo dragon genome, suggesting that these organisms have a robust innate immune system. Specifically, multiple Komodo antimicrobial peptide genes were identified. Importantly, many of the antimicrobial peptide genes were found in gene clusters. We found that these innate immunity genes are conserved among reptiles, and the organization is similar to that seen in other avian and reptilian species. Having the genome of this important squamate will allow researchers to learn more about reptilian gene families and will be a valuable resource for researchers studying the evolution and biology of the endangered Komodo dragon. Electronic supplementary material The online version of this article (10.1186/s12864-019-6029-y) contains supplementary material, which is available to authorized users.
Collapse
|
17
|
Aggarwal DD, Rybnikov S, Cohen I, Frenkel Z, Rashkovetsky E, Michalak P, Korol AB. Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity. Genetica 2019; 147:291-302. [PMID: 31240599 DOI: 10.1007/s10709-019-00070-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2018] [Accepted: 06/18/2019] [Indexed: 12/31/2022]
Abstract
Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.
Collapse
|
18
|
Kinney N, Titus-Glover K, Wren JD, Varghese RT, Michalak P, Liao H, Anandakrishnan R, Pulenthiran A, Kang L, Garner HR. CAGm: a repository of germline microsatellite variations in the 1000 genomes project. Nucleic Acids Res 2019; 47:D39-D45. [PMID: 30329086 PMCID: PMC6323891 DOI: 10.1093/nar/gky969] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Revised: 10/04/2018] [Accepted: 10/05/2018] [Indexed: 12/14/2022] Open
Abstract
The human genome harbors an abundance of repetitive DNA; however, its function continues to be debated. Microsatellites-a class of short tandem repeat-are established as an important source of genetic variation. Array length variants are common among microsatellites and affect gene expression; but, efforts to understand the role and diversity of microsatellite variation has been hampered by several challenges. Without adequate depth, both long-read and short-read sequencing may not detect the variants present in a sample; additionally, large sample sizes are needed to reveal the degree of population-level polymorphism. To address these challenges we present the Comparative Analysis of Germline Microsatellites (CAGm): a database of germline microsatellites from 2529 individuals in the 1000 genomes project. A key novelty of CAGm is the ability to aggregate microsatellite variation by population, ethnicity (super population) and gender. The database provides advanced searching for microsatellites embedded in genes and functional elements. All data can be downloaded as Microsoft Excel spreadsheets. Two use-case scenarios are presented to demonstrate its utility: a mononucleotide (A) microsatellite at the BAT-26 locus and a dinucleotide (CA) microsatellite in the coding region of FGFRL1. CAGm is freely available at http://www.cagmdb.org/.
Collapse
|
19
|
Michalak P, Kang L, Schou MF, Garner HR, Loeschcke V. Genomic signatures of experimental adaptive radiation in Drosophila. Mol Ecol 2018; 28:600-614. [PMID: 30375065 DOI: 10.1111/mec.14917] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2018] [Revised: 10/03/2018] [Accepted: 10/17/2018] [Indexed: 12/12/2022]
Abstract
Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large-scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat-shock proteins. The SNP patterns were highly stressor-specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co-occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat-shock selection regimes.
Collapse
|
20
|
Velmurugan KR, Michalak P, Kang L, Fonville NC, Garner HR. Dysfunctional DNA repair pathway via defective FANCD2 gene engenders multifarious exomic and transcriptomic effects in Fanconi anemia. Mol Genet Genomic Med 2018; 6:1199-1208. [PMID: 30450770 PMCID: PMC6305641 DOI: 10.1002/mgg3.502] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2018] [Revised: 09/20/2018] [Accepted: 10/10/2018] [Indexed: 01/27/2023] Open
Abstract
Background Fanconi anemia (FA) affects only one in 130,000 births, but has severe and diverse clinical consequences. It has been theorized that defects in the FA DNA cross‐link repair complex lead to a spectrum of variants that are responsible for those diverse clinical phenotypes. Methods Using NextGen sequencing, we show that a clinically derived FA cell line had accumulated numerous genetic variants, including high‐impact mutations, such as deletion of start codons, introduction of premature stop codons, missense mutations, and INDELs. Results About 65% of SNPs and 55% of INDELs were found to be commonly present in both the FA dysfunctional and retrovirally corrected cell lines, showing their common origin. The number of INDELs, but not SNPs, is decreased in FANCD2‐corrected samples, suggesting that FANCD2 deficiency preferentially promotes the origin of INDELs. These genetic modifications had a considerable effect on the transcriptome, with statistically significant changes in the expression of 270 genes. These genetic and transcriptomic variants significantly impacted pathways and molecular functions, spanning a diverse spectrum of disease phenotypes/symptoms, consistent with the disease diversity seen in FA patients. Conclusion These results underscore the consequences of defects in the DNA cross‐link repair mechanism and indicate that accumulating diverse mutations from individual parent cells may make it difficult to anticipate the longitudinal clinical behavior of emerging disease states in an individual with FA.
Collapse
|
21
|
Michalak P, Kang L. Unique divergence of the breast cancer 2 ( BRCA2) gene in Neanderthals. Hereditas 2018; 155:34. [PMID: 30410429 PMCID: PMC6215347 DOI: 10.1186/s41065-018-0073-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2018] [Accepted: 10/23/2018] [Indexed: 11/10/2022] Open
Abstract
Unique divergence of the BRCA2, a tumor suppressor gene, in Neanderthals relative to other primates, including modern humans, is highlighted. This divergence with potentially pathogenic consequences raises a question about cancer susceptibility in the archaic species that was replaced by modern humans about 40,000 years ago.
Collapse
|
22
|
Yablonovitch AL, Fu J, Li K, Mahato S, Kang L, Rashkovetsky E, Korol AB, Tang H, Michalak P, Zelhof AC, Nevo E, Li JB. Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates. Nat Commun 2017; 8:1570. [PMID: 29146998 PMCID: PMC5691062 DOI: 10.1038/s41467-017-01658-2] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 10/05/2017] [Indexed: 12/30/2022] Open
Abstract
Determining the mechanisms by which a species adapts to its environment is a key endeavor in the study of evolution. In particular, relatively little is known about how transcriptional processes are fine-tuned to adjust to different environmental conditions. Here we study Drosophila melanogaster from 'Evolution Canyon' in Israel, which consists of two opposing slopes with divergent microclimates. We identify several hundred differentially expressed genes and dozens of differentially edited sites between flies from each slope, correlate these changes with genetic differences, and use CRISPR mutagenesis to validate that an intronic SNP in prominin regulates its editing levels. We also demonstrate that while temperature affects editing levels at more sites than genetic differences, genetically regulated sites tend to be less affected by temperature. This work shows the extent to which gene expression and RNA editing differ between flies from different microclimates, and provides insights into the regulation responsible for these differences.
Collapse
|
23
|
Kang L, George P, Price DK, Sharakhov I, Michalak P. Mapping Genomic Scaffolds to Chromosomes Using Laser Capture Microdissection in Application to Hawaiian Picture-Winged Drosophila. Cytogenet Genome Res 2017; 152:204-212. [PMID: 29130948 DOI: 10.1159/000481790] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/28/2017] [Indexed: 11/19/2022] Open
Abstract
Next-generation sequencing technologies have led to a decreased cost and an increased throughput in genome sequencing. Yet, many genome assemblies based on short sequencing reads have been assembled only to the scaffold level due to the lack of sufficient chromosome mapping information. Traditional ways of mapping scaffolds to chromosomes require a large amount of laboratory work and time to generate genetic and/or physical maps. To address this problem, we conducted a rapid technique which uses laser capture microdissection and enables mapping scaffolds of de novo genome assemblies directly to chromosomes in Hawaiian picture-winged Drosophila. We isolated and sequenced intact chromosome arms from larvae of D. differens. By mapping the reads of each chromosome to the recently assembled scaffolds from 3 Hawaiian picture-winged Drosophila species, at least 67% of the scaffolds were successfully assigned to chromosome arms. Even though the scaffolds are not ordered within a chromosome, the fast-generated chromosome information allows for chromosome-related analyses after genome assembling. We utilize this new information to test the faster-X evolution effect for the first time in these Hawaiian picture-winged Drosophila species.
Collapse
|
24
|
Michalak P, Grzesik J, Rafiński J. TESTS FOR SEXUAL INCOMPATIBILITY BETWEEN TWO NEWT SPECIES, TRITURUS VULGARIS
AND TRITURUS MONTANDONI
: NO-CHOICE MATING DESIGN. Evolution 2017; 51:2045-2050. [DOI: 10.1111/j.1558-5646.1997.tb05128.x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/1997] [Accepted: 08/04/1997] [Indexed: 11/29/2022]
|
25
|
Bishop BM, Juba ML, Russo PS, Devine M, Barksdale SM, Scott S, Settlage R, Michalak P, Gupta K, Vliet K, Schnur JM, van Hoek ML. Discovery of Novel Antimicrobial Peptides from Varanus komodoensis (Komodo Dragon) by Large-Scale Analyses and De-Novo-Assisted Sequencing Using Electron-Transfer Dissociation Mass Spectrometry. J Proteome Res 2017; 16:1470-1482. [PMID: 28164707 DOI: 10.1021/acs.jproteome.6b00857] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Komodo dragons are the largest living lizards and are the apex predators in their environs. They endure numerous strains of pathogenic bacteria in their saliva and recover from wounds inflicted by other dragons, reflecting the inherent robustness of their innate immune defense. We have employed a custom bioprospecting approach combining partial de novo peptide sequencing with transcriptome assembly to identify cationic antimicrobial peptides from Komodo dragon plasma. Through these analyses, we identified 48 novel potential cationic antimicrobial peptides. All but one of the identified peptides were derived from histone proteins. The antimicrobial effectiveness of eight of these peptides was evaluated against Pseudomonas aeruginosa (ATCC 9027) and Staphylococcus aureus (ATCC 25923), with seven peptides exhibiting antimicrobial activity against both microbes and one only showing significant potency against P. aeruginosa. This study demonstrates the power and promise of our bioprospecting approach to cationic antimicrobial peptide discovery, and it reveals the presence of a plethora of novel histone-derived antimicrobial peptides in the plasma of the Komodo dragon. These findings may have broader implications regarding the role that intact histones and histone-derived peptides play in defending the host from infection. Data are available via ProteomeXChange with identifier PXD005043.
Collapse
|