Search

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant ‘islands’. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which ‘islands’ of habitat are most in need of management and protection.

Many initiatives have addressed the global need to upskill biologists in bioinformatics tools and techniques. Australia is not unique in its requirement for such training, but due to its large size and relatively small and geographically dispersed population, Australia faces specific challenges. A combined training approach was implemented by the authors to overcome these challenges. The “hybrid” method combines guidance from experienced trainers with the benefits of both webinar-style delivery and concurrent face-to-face hands-on practical exercises in classrooms. Since 2017, the hybrid method has been used to conduct 9 hands-on bioinformatics training sessions at international scale in which over 800 researchers have been trained in diverse topics on a range of software platforms. The method has become a key tool to ensure scalable and more equitable delivery of short-course bioinformatics training across Australia and can be easily adapted to other locations, topics, or settings.

Natural history museum collections worldwide represent a tremendous resource of information on past and present biodiversity. Fish, reptiles, amphibians and many invertebrate collections have often been preserved in ethanol for decades or centuries and our knowledge on the genomic and metagenomic research potential of such material is limited. Here, we use ancient DNA protocols, combined with shotgun sequencing to test the molecular preservation in liver, skin and bone tissue from five old (1842 to 1964) museum specimens of the common garter snake (Thamnophis sirtalis). When mapping reads to a T. sirtalis reference genome, we find that the DNA molecules are highly damaged with short average sequence lengths (38–64 bp) and high C-T deamination, ranging from 9% to 21% at the first position. Despite this, the samples displayed relatively high endogenous DNA content, ranging from 26% to 56%, revealing that genome-scale analyses are indeed possible from all specimens and tissues included here. Of the three tested types of tissue, bone shows marginally but significantly higher DNA quality in these metrics. Though at least one of the snakes had been exposed to formalin, neither the concentration nor the quality of the obtained DNA was affected. Lastly, we demonstrate that these specimens display a diverse and tissue-specific microbial genetic profile, thus offering authentic metagenomic data despite being submerged in ethanol for many years. Our results emphasize that historical museum collections continue to offer an invaluable source of information in the era of genomics.

The Aeolian wall lizard Podarcis raffonei is an island endemic that survives only on three tiny islets, and on the Capo Grosso peninsula of the Vulcano island, thus is among the European vertebrates with the smallest range and one of the most threatened by extinction. This species is declining due to competition and hybridization with the non-native lizard Podarcis siculus, but a regular monitoring program is lacking. Here we assessed the size and status of the Capo Grosso population of P. raffonei on Vulcano. In September 2015 we captured 30 individuals showing the typical brown phenotype of P. raffonei, while one single male showed a green phenotype, apparently intermediate between P. raffonei and the non-native Podarcis siculus. In May 2017, only 47% of 131 individuals showed the typical brown phenotype (P. raffonei-like) and 53% showed the green phenotype (P. siculus-like). Based on N-mixture models and removal sampling the estimated size of the Capo Grosso population was of 800-1300 individuals in 2017, being similar to 2015; available data suggest that the total range of the species could be as small as 2 ha. The frequency of individuals with the typical P. raffonei phenotype dramatically dropped between two samplings with a parallel increase of individuals displaying the green phenotype. Observation on outdoor captive-bred individuals demonstrates plasticity for colouration in P. raffonei individuals from Capo Grosso, with several individuals showing the typical brown pattern in September 2017 and a green pattern in March 2021. Non-exclusive hypotheses, including hybridization with P. siculus and plasticity in colour pattern of P. raffonei, are discussed to explain the phenotypic shifts of the P. raffonei population of Capo Grosso. While genomic evidence is required to reach conclusions and investigate eventual hybridization, it is urgent to undertake a programme for the monitoring and management of this lizard.

Squamate reptiles exhibit high variation in their traits and geographical distribution and are therefore fascinating taxa for evolutionary and ecological research. However, high-quality genomic recourses are very limited for this group of species, which inhibits some research efforts. To address this gap, we assembled a high-quality genome of the common lizard Zootoca vivipara (Lacertidae) using a combination of high coverage Illumina (shotgun and mate-pair) and PacBio sequence data, with RNAseq data and genetic linkage maps. The 1.46 Gbp genome assembly has scaffold N50 of 11.52 Mbp with N50 contig size of 220.4 Kbp and only 2.96% gaps. A BUSCO analysis indicates that 97.7% of the single-copy Tetrapoda orthologs were recovered in the assembly. In total 19,829 gene models were annotated in the genome using a combination of three ab initio and homology-based methods. To improve the chromosome-level assembly, we generated a high-density linkage map from wild-caught families and developed a novel analytical pipeline to accommodate multiple paternity and unknown father genotypes. We successfully anchored and oriented almost 90% of the genome on 19 linkage groups. This annotated and oriented chromosome-level reference genome represents a valuable resource to facilitate evolutionary studies in squamate reptiles.

Targeted sequencing of genomic regions is a cost- and time-efficient approach for screening patient cohorts. We present a fast and efficient workflow to analyze highly imbalanced, targeted next-generation sequencing data generated using molecular inversion probe (MIP) capture. Our Snakemake pipeline performs sample demultiplexing, overlap paired-end merging, alignment, MIP-arm trimming, variant calling, coverage analysis and report generation. Further, we support the analysis of probes specifically designed to capture certain structural variants and can assign sex using Y-chromosome-unique probes. In a user-friendly HTML report, we summarize all these results including covered, incomplete or missing regions, called variants and their predicted effects. We developed and tested our pipeline using the hemophilia A & B MIP design from the “My Life, Our Future” initiative. HemoMIPs is available as an open-source tool on GitHub at: https://github.com/kircherlab/hemoMIPs, Author summary Next generation sequencing techniques enable researchers to identify causal variants for patients in large cohorts. Targeted sequencing approaches capture genomic regions of interest to allow high throughput and cost efficient patient-specific data generation. HemoMIPs is an open-source software that analyses targeted sequencing datasets generated using molecular inversion probes (MIPs) and provides HTML reports of pathogenic and benign variants, patient sex, existence of known structural variants as well as performance statistics on the sequencing run.

A geographically isolated maize landrace cultivated on nitrogen-depleted fields without synthetic fertilizer in the Sierra Mixe region of Oaxaca, Mexico utilizes nitrogen derived from the atmosphere and develops an extensive network of mucilage-secreting aerial roots that harbors a diazotrophic microbiota. Targeting these diazotrophs, we selected nearly 600 microbes from a collection isolated from these plants and confirmed their ability to incorporate heavy nitrogen (15N2) metabolites in vitro. Sequencing their genomes and conducting comparative bioinformatic analyses showed that these genomes had substantial phylogenetic diversity. We examined each diazotroph genome for the presence of nif genes essential to nitrogen fixation (nifHDKENB) and carbohydrate utilization genes relevant to the mucilage polysaccharide digestion. These analyses identified diazotrophs that possessed canonical nif gene operons, as well as many other operon configurations with concomitant fixation and release of >700 different 15N labeled metabolites. We further demonstrated that many diazotrophs possessed alternative nif gene operons and confirmed their genomic potential to derive chemical energy from mucilage polysaccharide to fuel nitrogen fixation. These results confirm that some diazotrophic bacteria associated with Sierra Mixe maize were capable of incorporating atmospheric nitrogen into their small molecule extracellular metabolites through multiple nif gene configurations while others were able to fix nitrogen without the canonical (nifHDKENB) genes.Data SummaryGenetic resources, including biological materials and nucleic acid sequences, were accessed under an Access and Benefit Sharing (ABS) Agreement between the Sierra Mixe community and the Mars Corporation, and with authorization from the Mexican government. An internationally recognized certificate of compliance has been issued by the Mexican government under the Nagoya Protocol for such activities (ABSCH-IRCC-MX-207343-3). Any party seeking access to the nucleic acid sequences underlying the analysis reported here is subject to the full terms and obligations of the ABS agreement and the authorization from the government of Mexico. Individuals wishing to access nucleic acid sequence data for scientific research activities should contact Mars Incorporated Chief Science Officer at CSO@effem.com.

ZMAT2 is among the least-studied of mammalian proteins and genes, even though it is the ortholog of Snu23, a protein involved in pre-mRNA splicing in yeast. Here we have used data from genomic and gene expression repositories to examine the Zmat2 gene and locus in 8 terrestrial vertebrates, 10 ray-finned fish, and 1 lobe-finned fish representing > 500 million years of evolutionary diversification. The analyses revealed that vertebrate Zmat2 genes are similar to their mammalian counterparts, as in 16/19 species studied they contain 6 exons, and in 18/19 encode a single conserved protein. However, unlike in mammals, no Zmat2 pseudogenes were identified in these vertebrates, although an expressed Zmat2 paralog was characterized in flycatcher that resembled a DNA copy of a processed and retro-transposed mRNA, and thus could be a proto-pseudogene captured during its evolutionary journey from active to inert. The Zmat2 locus in terrestrial vertebrates, and in spotted gar and coelacanth, also shares additional genes with its mammalian counterparts, including Histidyl-tRNA synthetase (Hars), Hars2, and others, but these are absent from the Zmat2 locus in teleost fish, in which Stem-loop-binding protein 2 (Slbp2) and Lymphocyte cytosolic protein 2a (Lcp2a) are present instead. Taken together, these observations argue that a recognizable Zmat2 was present in the earliest vertebrate ancestors, and postulate that during chromosomal tetraploidization and subsequent re-diploidization during modern teleost evolution, the duplicated Zmat2 gene was retained and the original lost. This study also highlights how information from genomic resources can be leveraged to reveal new biologically significant insights.

Mycobacterium chelonae is a rapidly growing nontuberculous mycobacterium that is a common cause of nosocomial infections. Here we describe investigation of a possible nosocomial transmission of M. chelonae at the Hospital of the University of Pennsylvania (HUP). M. chelonae strains with similar high-level antibiotic resistance patterns were isolated from two patients who developed post-operative infections at HUP in 2017, suggesting a possible point source infection. The isolates, along with other clinical isolates from other patients, were sequenced using the Illumina and Oxford Nanopore technologies. The resulting short and long reads were hybrid assembled into draft genomes. The genomes were compared by quantifying single nucleotide variants in the core genome and assessed using a control dataset to quantify error rates in comparisons of identical genomes. We show that all M. chelonae isolates tested were highly dissimilar, as indicated by high pairwise SNV values, consistent with environmental acquisition and not a nosocomial point source. Our control dataset determined a threshold for evaluating identity between strains while controlling for sequencing error. Finally, antibiotic resistance genes were predicted for our isolates, and several single nucleotide variants were identified that have the potential to modulated drug resistance.

In this review, we summarize the current knowledge concerning the eukaryotic protozoan parasite Leishmania tarentolae, with a main focus on its potential for biotechnological applications. We will also discuss the genus, subgenus, and species-level classification of this parasite, its life cycle and geographical distribution, and similarities and differences to human-pathogenic species, as these aspects are relevant for the evaluation of biosafety aspects of L. tarentolae as host for recombinant DNA/protein applications. Studies indicate that strain LEM-125 but not strain TARII/UC of L. tarentolae might also be capable of infecting mammals, at least transiently. This could raise the question of whether the current biosafety level of this strain should be reevaluated. In addition, we will summarize the current state of biotechnological research involving L. tarentolae and explain why this eukaryotic parasite is an advantageous and promising human recombinant protein expression host. This summary includes overall biotechnological applications, insights into its protein expression machinery (especially on glycoprotein and antibody fragment expression), available expression vectors, cell culture conditions, and its potential as an immunotherapy agent for human leishmaniasis treatment. Furthermore, we will highlight useful online tools and, finally, discuss possible future applications such as the humanization of the glycosylation profile of L. tarentolae or the expression of mammalian recombinant proteins in amastigote-like cells of this species or in amastigotes of avirulent human-pathogenic Leishmania species., Author summary Recent studies indicate that the classification of Leishmania species into one of the four currently recognized subgenera (L. [Leishmania], L. [Sauroleishmania], L. [Viannia], and L. [Mundinia]), or even the related genera, can be a challenging task. A monotonous morphology of this group is coupled with a high genetic variability in nature, a multitude of clinical manifestations, and a propensity for rapid evolution in culture. As some geographic regions have not been adequately sampled, the known biodiversity of this group may be an underestimate of its true dimensions. Therefore, the taxonomic classification of the genus Leishmania and its close relatives is not definitively settled. In this review, we will focus on biotechnological applications of L. tarentolae, a protozoan parasite of geckos. As this species belongs to the genus Leishmania and subgenus L. (Sauroleishmania), and because many Leishmania are human-pathogenic, the taxonomic status and phylogenetic position of L. tarentolae are relevant for biosafety and also discussed herein. The subgenus of L. (Sauroleishmania) mostly contains species that are infectious to reptiles, but some species and strains, such as L. adleri and L. tarentolae strain LEM-125, were also shown to be (transiently) infectious to humans. However, most strains of L. tarentolae are nonpathogenic to humans and can be easily handled as laboratory culture. Therefore, this species has been a successful model system for representing other Leishmania in basic research. It has also been of great interest for the scientific community, as it represents a promising host for the expression of human recombinant proteins (including glycoproteins and its future application for the expression of full-length antibodies) and immunotherapy agent for human leishmaniasis treatment. These two topics will also be discussed.

Mycobacterium chelonae is a member of the Mycobacterium chelonae-abscessus complex and a cause of opportunistic disease in fish, reptiles, birds, and mammals including humans. Isolates in the complex are often difficult to identify and have differing antimicrobial susceptibilities. Thirty-one previously identified rapidly-growing, non-tuberculous Mycobacterium sp. isolates cultured from biofilms, fish, reptiles, mammals, including humans, and three ATCC reference strains were evaluated with nine M. chelonae-abscessus complex whole genome sequences from GenBank by phylogenomic analysis, targeted gene comparisons, and in-vitro antimicrobial susceptibility patterns to assess strain variation among isolates from different sources. Results revealed minimal genetic variation among the M. chelonae strains. However, the core genomic alignment and SNP pattern of the complete 16S rRNA sequence clearly separated the turtle type strain ATCC 35752T from the clinical isolates and human reference strain "M. chelonae chemovar niacinogenes" ATCC 19237, providing evidence of two distinct subspecies. Concatenation of the partial rpoB (752 bp) and complete hsp65 (1,626 bp) sequence produced the same species/subspecies delineations as the core phylogeny. Partial rpoB and hsp65 sequences identified all the clinical isolates to the appropriate species level when respective cut-offs of 98% and 98.4% identity to the M. chelonae type strain ATCC 35752T were employed. The human strain, ATCC19237, was the most representative strain for the evaluated human, veterinary, and environmental strains. Additionally, two isolates were identified as Mycobacterium saopaulense, its first identification in a non-fish or non-human host.

Mycoplasma agassizii is a common cause of upper respiratory tract disease in Mojave desert tortoises (Gopherus agassizii). So far, only two strains of this bacterium have been sequenced, and very little is known about its patterns of genetic diversity. Understanding genetic variability of this pathogen is essential to implement conservation programs for their threatened, long-lived hosts. We used next generation sequencing to explore the genomic diversity of 86 cultured samples of M. agassizii collected from mostly healthy Mojave and Sonoran desert tortoises in 2011 and 2012. All samples with enough sequencing coverage exhibited a higher similarity to M. agassizii strain PS6(T) (collected in Las Vegas Valley, Nevada) than to strain 723 (collected in Sanibel Island, Florida). All eight genomes with a sequencing coverage over 2x were subjected to multiple analyses to detect single-nucleotide polymorphisms (SNPs). Strikingly, even though we detected 1373 SNPs between strains PS6(T) and 723, we did not detect any SNP between PS6(T) and our eight samples. Our whole genome analyses reveal that M. agassizii strain PS6(T) may be present across a wide geographic extent in healthy Mojave and Sonoran desert tortoises. US Fish and Wildlife ServiceUS Fish & Wildlife Service [F19AC00070]; Nevada INBRE (National Institute of General Medical Sciences, National Institutes of Health) [P20GM103440]; National Science FoundationNational Science Foundation (NSF) [MCB 1818288] This work was supported by grant F19AC00070 from the US Fish and Wildlife Service, and partially supported by Nevada INBRE (funded by grant P20GM103440 from the National Institute of General Medical Sciences, National Institutes of Health). AL and DAP were also supported by a grant from the National Science Foundation (MCB 1818288). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Comparative genomics studies may be used to acquire new knowledge regarding genome architecture, which defines the rules for combining sets of genes in the genome of living organisms. Hundreds of thousands of prokaryotic genomes have been sequenced and assembled. However, computational tools capable of simultaneously comparing large numbers of genomes are lacking. We developed the Genome Complexity Browser, a tool that allows the visualization of gene contexts, in a graph-based format, and the quantification of variability for different segments of a genome. The graph-based visualization allows the inspection of changes in gene contents and neighborhoods across hundreds of genomes, simultaneously, which may facilitate the identification of conserved and variable segments of operons or the estimation of the overall variability associated with a particular genome locus. We introduced a measure called complexity, to quantify genome variability. Intraspecies and interspecies comparisons revealed that regions with high complexity values tended to be located in areas that are conserved across different strains and species., Author summary The comparison of genomes among different bacteria and archaea species has revealed that many species frequently exchange genes. Occasionally, such horizontal gene transfer events result in the acquisition of pathogenic properties or antibiotic resistance in the recipient organism. Previously, the probabilities of gene insertions were found to vary, with unequal distributions along a chromosome. At some loci, referred to as hotspots, changes occur with much higher frequencies compared with other regions of the chromosome. We developed a computational method and a software tool, called Genome Complexity Browser, that allows the identification of genome variability hotspots and the visualization of changes. We compared the localization of various hotspots and revealed that some demonstrate conserved localizations, even across species, whereas others are transient. Our tool allows users to visually inspect the patterns of gene changes in graph-based format, which presents the visualization in a format that is both compact and informative.

Maximum lifespan for most animal species is difficult to define. This is challenging for wildlife management as it is critical for estimating important aspects of population biology such as mortality rate, population viability, and period of reproductive potential. Recently, it has been shown cytosine-phosphate-guanine (CpG) density is predictive of maximum lifespan in vertebrates. This has made it possible to predict lifespan in long-lived species, which are generally the most intractable. In this study, we use gene promoter CpG density to predict the lifespan of five marine turtle species. Marine turtles are a particularly difficult group for lifespan estimation because of their migratory behaviour, longevity and high juvenile mortality rates, which all restrict individual tracking over their lifespan. Sanger sequencing was used to determine the CpG density in selected promoters. We predicted the lifespans for marine turtle species ranged from 50.4 years (flatback turtle, Natator depressus) to 90.4 years (leatherback turtle, Dermochelys coriacea). These lifespan predictions have broad applications in marine turtle research such as better understanding life cycles and determining population viability.

The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.

We generated a phylogeny for Caucasian rock lizards (Darevskia), and included six other families of true lizards (Lacertini), based on complete mitochondrial genome analysis. Next-generation sequencing (NGS) of genomic DNA was used to obtain 16 new mitogenomes of Darevskia. These, along with 35 sequences downloaded from GenBank: genera Darevskia, Zootoca, Podarcis, Phoenicolacerta, Takydromus, Lacerta, and Eremias-were used in the analysis. All four analytical methods (Bayesian Inference, BI; Maximum Likelihood, ML; Maximum Parsimony, MP; and Neighbor-Joining, NJ) showed almost congruent intra-generic topologies for Darevskia and other lizard genera. However, ML and NJ methods on one side, and BI and MP methods on the other harvested conflicting phylogenies. The ML/NJ topology supports earlier published separation of Darevskia into three mitochondrial clades (Murphy, Fu, Macculloch, Darevsky, and Kupinova, 2000), but BI and MP topologies support that the basal branching occurred between D. parvula from the western Lesser Caucasus and the rest of Darevskia. All topologies altered the phylogenetic position of some individual species, including D. daghestanica, D. derjugini, and D. chlorogaster. Reanalysis after excluding four saturated genes from the data set, and excluding genus Eremias gives fully convergent topologies. The most basal branching for true lizards was between Far Eastern Takydromus and the Western Eurasian genera (BI). Comparing phylogenetic performance of individual genes relative to whole mitogenome data, concatenated 16S RNA (the least saturated gene in our analyses) and Cytochrome b genes generate a robust phylogeny that is fully congruent with that based on the complete mitogenome.

Background Reptile-associated Salmonella bacteria are a major, but often neglected cause of both gastrointestinal and bloodstream infection in humans globally. The diversity of Salmonella enterica has not yet been determined in venomous snakes, however other ectothermic animals have been reported to carry a broad range of Salmonella bacteria. We investigated the prevalence and diversity of Salmonella in a collection of venomous snakes and non-venomous reptiles. Methodology/Principle findings We used a combination of selective enrichment techniques to establish a unique dataset of reptilian isolates to study Salmonella enterica species-level evolution and ecology and used whole-genome sequencing to investigate the relatedness of phylogenetic groups. We observed that 91% of venomous snakes carried Salmonella, and found that a diverse range of serovars (n = 58) were carried by reptiles. The Salmonella serovars belonged to four of the six Salmonella enterica subspecies: diarizonae, enterica, houtanae and salamae. Subspecies enterica isolates were distributed among two distinct phylogenetic clusters, previously described as clade A (52%) and clade B (48%). We identified metabolic differences between S. diarizonae, S. enterica clade A and clade B involving growth on lactose, tartaric acid, dulcitol, myo-inositol and allantoin. Significance We present the first whole genome-based comparative study of the Salmonella bacteria that colonise venomous and non-venomous reptiles and shed new light on Salmonella evolution. Venomous snakes examined in this study carried a broad range of Salmonella, including serovars which have been associated with disease in humans such as S. Enteritidis. The findings raise the possibility that venomous snakes could be a reservoir for Salmonella serovars associated with human salmonellosis., Author summary Salmonella enterica is a remarkable bacterial species that causes Neglected Tropical Diseases globally. The burden of disease is greatest in some of the most poverty-afflicted regions of Africa, where salmonellosis frequently causes bloodstream infection with fatal consequences. The bacteria have the ability to colonise the gastrointestinal tract of a wide range of animals including reptiles. Direct or indirect contact between reptiles and humans can cause salmonellosis. In this study, we determined the prevalence and diversity of Salmonella in a collection of African venomous snakes for the first time. We showed that the majority of venomous snakes (91%) in our study carry Salmonella, and used bacterial genomics to assign two thirds of isolates to the S. enterica subspecies enterica which is associated with human salmonellosis. We identified two evolutionary groups of S. enterica subspecies enterica that display distinct metabolic profiles with infection-relevant carbon sources. Our findings could have a broad significance because venomous snakes can move freely around human dwellings in tropical regions of the world such as Africa, and could potentially shed contaminated faecal matter onto surfaces and into water supplies.

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies.

The family Arenaviridae comprises three genera, Mammarenavirus, Reptarenavirus and the most recently added Hartmanivirus. Arenaviruses have a bisegmented genome with ambisense coding strategy. For mammarenaviruses and reptarenaviruses the L segment encodes the Z protein (ZP) and the RNA-dependent RNA polymerase, and the S segment encodes the glycoprotein precursor and the nucleoprotein. Herein we report the full length genome and characterization of Haartman Institute snake virus-1 (HISV-1), the putative type species of hartmaniviruses. The L segment of HISV-1 lacks an open-reading frame for ZP, and our analysis of purified HISV-1 particles by SDS-PAGE and electron microscopy further support the lack of ZP. Since we originally identified HISV-1 in co-infection with a reptarenavirus, one could hypothesize that co-infecting reptarenavirus provides the ZP to complement HISV-1. However, we observed that co-infection does not markedly affect the amount of hartmanivirus or reptarenavirus RNA released from infected cells in vitro, indicating that HISV-1 does not benefit from reptarenavirus ZP. Furthermore, we succeeded in generating a pure HISV-1 isolate showing the virus to replicate without ZP. Immunofluorescence and ultrastructural studies demonstrate that, unlike reptarenaviruses, HISV-1 does not produce the intracellular inclusion bodies typical for the reptarenavirus-induced boid inclusion body disease (BIBD). While we observed HISV-1 to be slightly cytopathic for cultured boid cells, the histological and immunohistological investigation of HISV-positive snakes showed no evidence of a pathological effect. The histological analyses also revealed that hartmaniviruses, unlike reptarenaviruses, have a limited tissue tropism. By nucleic acid sequencing, de novo genome assembly, and phylogenetic analyses we identified additional four hartmanivirus species. Finally, we screened 71 individuals from a collection of snakes with BIBD by RT-PCR and found 44 to carry hartmaniviruses. These findings suggest that harmaniviruses are common in captive snake populations, but their relevance and pathogenic potential needs yet to be revealed., Author summary From the 1930s to 2015 arenaviruses were known as mainly rodent-borne viruses, which occasionally infect humans, causing a severe disease. After isolation of novel arenaviruses from snakes, the family Arenaviridae now comprises three genera, Mammarenavirus, Reptarenavirus, and Hartmanivirus. Here we characterize a pure isolate of the putative hartmanivirus type species, Haartman Institute snake virus-1 (HISV-1) and report three new hartmanivirus species: Old Schoolhouse virus, Dante Muikkunen virus and Veterinary Pathology Zurich virus. The genomes of all these hartmaniviruses lack the matrix protein found in other arenaviruses, and using HISV-1 as a model we demonstrate that hartmaniviruses can nonetheless form infectious particles. Although HISV-1 induces cytopathic changes in cell culture, there was no evidence of a cytopathogenic effect in infected snakes, where hartmaniviruses have a restricted cell tropism. We further identified hartmanivirus infection in over 60% of snakes in a single breeding collection, indicating that hartmaniviruses are relatively common in boid snakes. Several animals of the collection studied were diagnosed with BIBD, and it is not clear whether hartmaniviruses play a role in the development of BIBD in reptarenavirus-infected snakes. The co-infection with viruses from two genera of the same virus family provides an interesting model for future studies.

Tissue samples from Australian carpet pythons (Morelia spilota) with neurological disease were screened for viruses using next-generation sequencing. Coding complete genomes of two bornaviruses were identified with the gene order 3’-N-X-P-G-M-L, representing a transposition of the G and M genes compared to other bornaviruses and most mononegaviruses. Use of these viruses to search available vertebrate genomes enabled recognition of further endogenous bornavirus-like elements (EBLs) in diverse placental mammals, including humans. Codivergence patterns and shared integration sites revealed an ancestral laurasiatherian EBLG integration (77 million years ago [MYA]) and a previously identified afrotherian EBLG integration (83 MYA). The novel python bornaviruses clustered more closely with these EBLs than with other exogenous bornaviruses, suggesting that these viruses diverged from previously known bornaviruses prior to the end-Cretaceous (K-Pg) extinction, 66 MYA. It is possible that EBLs protected mammals from ancient bornaviral disease, providing a selective advantage in the recovery from the K-Pg extinction. A degenerate PCR primer set was developed to detect a highly conserved region of the bornaviral polymerase gene. It was used to detect 15 more genetically distinct bornaviruses from Australian pythons that represent a group that is likely to contain a number of novel species., Author summary Unlike plants and animals, viruses don’t leave fossilised remains that can be used to study their ancient history. Rarely, however, some virus sequences end up as integrated copies in host genomes. These ‘fossilised’ virus sequences provide a window into virus evolution on a geological time scale. In this study, we identified bornaviruses in pythons, some of which had neurological disease. Bornaviruses have been associated with persistent neurological infections in a variety of mammalian and avian hosts. Here, we reveal that bornaviruses are associated with similar disease in non-avian reptiles as well. These new snake viruses turned out to be most closely related to bornavirus-like sequences that integrated into the DNA of mammals prior to the major extinction event at the end of the age of dinosaurs, the end-Cretaceous (K-Pg) extinction, 66 million years ago, allowing us to infer details of their prehistoric evolutionary history and origin. We also developed a test to detect these viruses from samples collected from live snakes, and used it to identify more related snake bornaviruses, indicating that this group will likely continue to expand.

Comparative functional genomic studies require the proper identification of gene orthologs to properly exploit animal biomedical research models. To identify gene orthologs, comprehensive, conserved gene synteny analyses are necessary to unwind gene histories that are convoluted by two rounds of early vertebrate genome duplication, and in the case of the teleosts, a third round, the teleost genome duplication (TGD). Recently, the genome of the spotted gar, a holostean outgroup to the teleosts that did not undergo this third genome duplication, was sequenced and applied as an orthology bridge to facilitate the identification of teleost orthologs to human genes and to enhance the power of teleosts as biomedical models. In this study, we apply the spotted gar orthology bridge to help describe the gene history of the vertebrate TNFAIP8 family. Members of the TNFAIP8 gene family have been linked to regulation of immune function and homeostasis and the development of multiple cancer types. Through a conserved gene synteny analysis, we identified zebrafish orthologs to human TNFAIP8L1 and TNFAIP8L3 genes and two co-orthologs to human TNFAIP8L2, but failed to identify an ortholog to human TNFAIP8. Through the application of the orthology bridge, we determined that teleost orthologs to human TNFAIP8 genes were likely lost in a genome inversion event after their divergence from their common ancestor with spotted gar. These findings demonstrate the value of this enhanced approach to gene history analysis and support the development of teleost models to study complex questions related to an array of biomedical issues, including immunity and cancer.

Nine protein arginine methyltransferases (PRMTs) are conserved in mammals and fish. Among these, PRMT1 is the major type I PRMT for asymmetric dimethylarginine (ADMA) formation and is the most conserved and widely distributed one. Two chicken prmt1 splicing variants were assembled and confirmed by RT-PCR experiments. However, only two scaffolds containing single separate prmt1 exon with high GC contents are present in the current chicken genome assembly. Besides, prmt1 exons are scattered in separate small scaffolds in most avian species. Complete prmt1 gene has only been predicted from two falcon species with few neighboring genes. Crocodilians are considered close to the common ancestor shared by crocodilians and birds. The gene arrangements around prmt1 in American alligator are different from that in birds but are largely conserved in human. Orthologues of genes in a large segment of human chromosomal 19 around PRMT1 are missing or not assigned to the current chicken chromosomes. In comparison, prmt8, the prmt1 paralogue, is on chicken chromosome 1 with the gene arrangements downstream of prmt8 highly conserved in birds, crocodilians, and human. However, the ones upstream vary greatly in birds. Biochemically, we found that though prmt1 transcripts were detected, limited or none PRMT1 protein was present in chicken tissues. Moreover, a much higher level of PRMT8 protein was detected in chicken brain than in mouse brain. While PRMT8 is brain specific in other vertebrate species studied, low level of PRMT8 was present in chicken but not mouse liver and muscle. We also showed that the ADMA level in chicken was similar to that in mouse. This study provides the critical information of chicken PRMT1 and PRMT8 for future analyses of the function of protein arginine methyltransferases in birds.

Here, we report the results of a large-scale PCR survey on the prevalence and diversity of adenoviruses (AdVs) in samples collected randomly from free-living reptiles. On the territories of the Guadarrama Mountains National Park in Central Spain and of the Chafarinas Islands in North Africa, cloacal swabs were taken from 318 specimens of eight native species representing five squamate reptilian families. The healthy-looking animals had been captured temporarily for physiological and ethological examinations, after which they were released. We found 22 AdV-positive samples in representatives of three species, all from Central Spain. Sequence analysis of the PCR products revealed the existence of three hitherto unknown AdVs in 11 Carpetane rock lizards (Iberolacerta cyreni), nine Iberian worm lizards (Blanus cinereus), and two Iberian green lizards (Lacerta schreiberi), respectively. Phylogeny inference showed every novel putative virus to be a member of the genus Atadenovirus. This is the very first description of the occurrence of AdVs in amphisbaenian and lacertid hosts. Unlike all squamate atadenoviruses examined previously, two of the novel putative AdVs had A+T rich DNA, a feature generally deemed to mirror previous host switch events. Our results shed new light on the diversity and evolution of atadenoviruses., The authors gratefully acknowledge the financial support provided by the Hungarian Scientific Research Fund (OTKA grant K100163) and by the projects MICIIN-CGL2011-24150/BOS and MINECO CGL2014-53523-P.

The Forkhead box transcription factor FoxP3 is pivotal to the development and function of regulatory T cells (Tregs), which make a major contribution to peripheral tolerance. FoxP3 is believed to perform a regulatory role in all the vertebrate species in which it has been detected. The prevailing view is that FoxP3 is absent in birds and that avian Tregs rely on alternative developmental and suppressive pathways. Prompted by the automated annotation of foxp3 in the ground tit (Parus humilis) genome, we have questioned this assumption. Our analysis of all available avian genomes has revealed that the foxp3 locus is missing, incomplete or of poor quality in the relevant genomic assemblies for nearly all avian species. Nevertheless, in two species, the peregrine falcon (Falco peregrinus) and the saker falcon (F. cherrug), there is compelling evidence for the existence of exons showing synteny with foxp3 in the ground tit. A broader phylogenomic analysis has shown that FoxP3 sequences from these three species are similar to crocodilian sequences, the closest living relatives of birds. In both birds and crocodilians, we have also identified a highly proline-enriched region at the N terminus of FoxP3, a region previously identified only in mammals.