Search

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability. 1 , 2 , 3 , 4 , 5 , 6 Moreover, small population size may promote gene flow with congeners and outbreeding depression. 7 , 8 , 9 , 10 , 11 , 12 , 13 Here, we examine the connection between habitat availability, effective population size (N e ), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa, 14 , 15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat. 16 , 17 , 18 However, we found long-term low N e , unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low N e and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits., Competing Interests: Declaration of interests M.H. and L.D. are on the advisory board of Colossal Biosciences and hold stock options. S. Bocklandt, K.M.P., G.G., and S. Barish are employed by Colossal Biosciences or FormBio. B.L. is CEO of Colossal Biosciences., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

The importance of COVID-19 surveillance from wastewater continues to grow since case-based surveillance in the general population has been scaled back world-wide. In Berlin, Germany, quantitative and genomic wastewater monitoring for SARS-CoV-2 is performed in three wastewater treatment plants (WWTP) covering 84% of the population since December 2021. The SARS-CoV-2 Omicron sublineage JN.1 (B.2.86.1.1), was first identified from wastewater on 22 October 2023 and rapidly became the dominant sublineage. This change was accompanied by a parallel and still ongoing increase in the notification-based 7-day-hospitalization incidence of COVID-19 and COVID-19 ICU utilization, indicating increasing COVID-19 activity in the (hospital-prone) population and a higher strain on the healthcare system. In retrospect, unique mutations of JN.1 could be identified in wastewater as early as September 2023 but were of unknown relevance at the time. The timely detection of new sublineages in wastewater therefore depends on the availability of new sequences from GISAID and updates to Pango lineage definitions and Nextclade. We show that genomic wastewater surveillance provides timely public health evidence on a regional level, complementing the existing indicators.

The use of RNA sequencing from wastewater samples is a valuable way for estimating infection dynamics and circulating lineages of SARS-CoV-2. This approach is independent from testing individuals and can therefore become the key tool to monitor this and potentially other viruses. However, it is equally important to develop easily accessible and scalable tools which can highlight critical changes in infection rates and dynamics over time across different locations given sequencing data from wastewater. Here, we provide an analysis of lineage dynamics in Berlin and New York City using wastewater sequencing and present PiGx SARS-CoV-2, a highly reproducible computational analysis pipeline with comprehensive reports. This end-to-end pipeline includes all steps from raw data to shareable reports, additional taxonomic analysis, deconvolution and geospatial time series analyses. Using simulated datasets (in silico generated and spiked-in samples) we could demonstrate the accuracy of our pipeline calculating proportions of Variants of Concern (VOC) from environmental as well as pre-mixed samples (spiked-in). By applying our pipeline on a dataset of wastewater samples from Berlin between February 2021 and January 2022, we could reconstruct the emergence of B.1.1.7(alpha) in February/March 2021 and the replacement dynamics from B.1.617.2 (delta) to BA.1 and BA.2 (omicron) during the winter of 2021/2022. Using data from very-short-reads generated in an industrial scale setting, we could see even higher accuracy in our deconvolution. Lastly, using a targeted sequencing dataset from New York City (receptor-binding-domain (RBD) only), we could reproduce the results recovering the proportions of the so-called cryptic lineages shown in the original study. Overall our study provides an in-depth analysis reconstructing virus lineage dynamics from wastewater. While applying our tool on a wide range of different datasets (from different types of wastewater sample locations and sequenced with different methods), we show that PiGx SARS-CoV-2 can be used to identify new mutations and detect any emerging new lineages in a highly automated and scalable way. Our approach can support efforts to establish continuous monitoring and early-warning projects for detecting SARS-CoV-2 or any other pathogen., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022. Published by Elsevier B.V.)

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone., (© 2022. The Author(s), under exclusive licence to The Genetics Society.)

The core microbiome, which refers to a set of consistent microbial features across populations, is of major interest in microbiome research and has been addressed by numerous studies. Understanding the core microbiome can help identify elements that lead to dysbiosis, and lead to treatments for microbiome-related health states. However, defining the core microbiome is a complex task at several levels. In this review, we consider the current state of core human microbiome research. We consider the knowledge that has been gained, the factors limiting our ability to achieve a reliable description of the core human microbiome, and the fields most likely to improve that ability. DNA sequencing technologies and the methods for analyzing metagenomics and amplicon data will most likely facilitate higher accuracy and resolution in describing the microbiome. However, more effort should be invested in characterizing the microbiome's interactions with its human host, including the immune system and nutrition. Other components of this holobiontic system should also be emphasized, such as fungi, protists, lower eukaryotes, viruses, and phages. Most importantly, a collaborative effort of experts in microbiology, nutrition, immunology, medicine, systems biology, bioinformatics, and machine learning is probably required to identify the traits of the core human microbiome.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

During the Miocene, Hyaenidae was a highly diverse family of Carnivora that has since been severely reduced to four species: the bone-cracking spotted, striped, and brown hyenas, and the specialized insectivorous aardwolf. Previous studies investigated the evolutionary histories of the spotted and brown hyenas, but little is known about the remaining two species. Moreover, the genomic underpinnings of scavenging and insectivory, defining traits of the extant species, remain elusive. Here, we generated an aardwolf genome and analyzed it together with the remaining three species to reveal their evolutionary relationships, genomic underpinnings of their scavenging and insectivorous lifestyles, and their respective genetic diversities and demographic histories. High levels of phylogenetic discordance suggest gene flow between the aardwolf lineage and the ancestral brown/striped hyena lineage. Genes related to immunity and digestion in the bone-cracking hyenas and craniofacial development in the aardwolf showed the strongest signals of selection, suggesting putative key adaptations to carrion and termite feeding, respectively. A family-wide expansion in olfactory receptor genes suggests that an acute sense of smell was a key early adaptation. Finally, we report very low levels of genetic diversity within the brown and striped hyenas despite no signs of inbreeding, putatively linked to their similarly slow decline in effective population size over the last ∼2 million years. High levels of genetic diversity and more stable population sizes through time are seen in the spotted hyena and aardwolf. Taken together, our findings highlight how ecological specialization can impact the evolutionary history, demographics, and adaptive genetic changes of an evolutionary lineage., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Since the 19th century, the addax ( Addax nasomaculatus ) has lost approximately 99% of its former range. Along with its close relatives, the blue antelope ( Hippotragus leucophaeus ) and the scimitar-horned oryx ( Oryx dammah ), the addax may be the third large African mammal species to go extinct in the wild in recent times. Despite this, the evolutionary history of this critically endangered species remains virtually unknown. To gain insight into the population history of the addax, we used hybridization capture to generate ten complete mitochondrial genomes from historical samples and assembled a nuclear genome. We found that both mitochondrial and nuclear diversity are low compared to other African bovids. Analysis of mitochondrial genomes revealed a most recent common ancestor ~32 kya (95% CI 11-58 kya) and weak phylogeographic structure, indicating that the addax likely existed as a highly mobile, panmictic population across its Sahelo-Saharan range in the past. PSMC analysis revealed a continuous decline in effective population size since ~2 Ma, with short intermediate increases at ~500 and ~44 kya. Our results suggest that the addax went through a major bottleneck in the Late Pleistocene, remaining at low population size prior to the human disturbances of the last few centuries.

Here, we present the first complete mitochondrial genome of the pentastomid Linguatula arctica collected from the nasal passages of a reindeer ( Rangifer tarandus ) in Norway. The full length mitochondrial genome of L. arctica , which measures 14,789 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of L. arctica with an overall base composition of 32.6% A, 27.5% T, 32.8% C, and 7,1% G., and a GC content of 39.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes., Competing Interests: No potential conflict of interest was reported by the author(s)., (© 2020 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

Background: In December 2016, three cases of serogroup B invasive meningococcal disease, including two children from the same middle school (11 to 15 years old pupils), occurred in the department (administrative district) Côtes-d'Armor (Brittany, France). They were infected by a rare strain (B:P1.7-2,4:F5-9:cc162), covered by the 4CMenB vaccine (Bexsero®). Four months later, two cases due to the same strain occurred in a high school in the same area (15 to 19 years old students). In accordance with French recommendations, vaccination was proposed to students of both schools and to all individuals aged 11-19 years living or studying in the hyperendemic area. We describe these vaccination campaigns, from the alert to the impact evaluation., Methods: The target population included 8884 people: 579 in the middle school, 2007 in the high school and 6298 in the community. In both schools, vaccination sessions were organized directly on site. In the community, teenagers were vaccinated by general practitioners. The vaccination campaign took place from May to October 2017. An active pharmacovigilance follow-up was set up to document adverse effects of the vaccine., Results: Considering the whole target population, the vaccination coverage was estimated at 43% for 1 dose and 34% for 2 doses. Higher vaccination coverage was observed in the schools (79% in the middle school and 42% in the high school for 2 doses) than in the community (27% for 2 doses). The reported adverse effects were consistent with the safety profile of the vaccine and no severe adverse effect was reported., Conclusions: This vaccination campaign was the third one implemented with Bexsero® in France and constitutes a reproducible approach for future targeted vaccination campaigns. No additional cases of the same strain have occurred since the end of the campaigns in the area.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana 1,2 -is an iconic species that is endemic to New Zealand 2,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 amniotes 2,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.

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions., (Copyright © 2019 Elsevier Inc. All rights reserved.)

In the anuran family Dendrobatidae, aposematic species obtain their toxic or unpalatable alkaloids from dietary sources, a process known as sequestering. To understand how toxicity evolved in this family, it is paramount to elucidate the pathways of alkaloid processing (absorption, metabolism, and sequestering). Here, we used an exploratory skin gene expression experiment in which captive-bred dendrobatids were fed alkaloids. Most of these experiments were performed with Dendrobates tinctorius , but some trials were performed with D. auratus , D. leucomelas and Allobates femoralis to explore whether other dendrobatids would show similar patterns of gene expression. We found a consistent pattern of up-regulation of genes related to muscle and mitochondrial processes, probably due to the lack of mutations related to alkaloid resistance in these species. Considering conserved pathways of drug metabolism in vertebrates, we hypothesize alkaloid degradation is a physiological mechanism of resistance, which was evidenced by a strong upregulation of the immune system in D. tinctorius , and of complement C2 across the four species sampled. Probably related to this strong immune response, we found several skin keratins downregulated, which might be linked to a reduction of the cornified layer of the epidermis. Although not conclusive, our results offer candidate genes and testable hypotheses to elucidate alkaloid processing in poison frogs.

In September 2016, a cluster of seven kayakers with clinical symptoms of leptospirosis with onset since July 2016 was reported to French health authorities. Human and animal investigations were undertaken to describe the outbreak, identify the likely place and source of infection and implement necessary control measures. We identified 103 patients with clinical symptoms of leptospirosis between 1 June and 31 October 2016 who lived in the Ille-et-Vilaine district in Brittany. Of these, 14 (including the original seven) reported contacts with the river Vilaine during the incubation period and were defined as outbreak cases: eight were confirmed by serology tests or PCR and six were probable without a laboratory confirmation for leptospirosis. All 14 cases were kayakers. Three distinct contamination sites were identified on a 30 km stretch of the river Vilaine. Nine cases reported having skin wounds while kayaking. None were vaccinated against leptospirosis. The outbreak was attributed to Leptospira kirschneri serogroup Grippotyphosa. Animal investigations did not allow identifying the possible reservoir. Leptospirosis outbreaks associated with freshwater sports are rare in temperate climates. The prevention of such outbreaks requires control of potential animal reservoirs in zones such as the Vilaine valley and that kayakers adopt the recommended individual prevention measures.

Background: Contiguous genome assemblies are a highly valued biological resource because of the higher number of completely annotated genes and genomic elements that are usable compared to fragmented draft genomes. Nonetheless, contiguity is difficult to obtain if only low coverage data and/or only distantly related reference genome assemblies are available., Findings: In order to improve genome contiguity, we have developed Cross-Species Scaffolding-a new pipeline that imports long-range distance information directly into the de novo assembly process by constructing mate-pair libraries in silico., Conclusions: We show how genome assembly metrics and gene prediction dramatically improve with our pipeline by assembling two primate genomes solely based on ∼30x coverage of shotgun sequencing data.

We present the first complete mitochondrial genome of the pentastomid Armillifer grandis (Arthropoda: Pentastomida) collected from the lungs of a rhinoceros viper ( Bitis nasicornis ) in the Democratic Republic of Congo. The full length mitochondrial genome of Armillifer grandis , which measures 16,073 bp in length, contains 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. A clear A + T bias is observed in the mitogenome of Armillifer grandis with an overall base composition of 34.6% A, 29.4% T, 29% C, and 6.9% G, and a GC content of 35.9%. The gene arrangement is identical to that of previously described pentastomid mitogenomes., Competing Interests: The authors thank the ‘Innovationsfond’ of the Museum für Naturkunde Berlin for financial support. The authors declare no conflict of interest. The authors alone are responsible for the content and writing of this article., (© 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

Here we present the first complete mitochondrial genome of the lesser sac-winged bat Saccopteryx leptura (Chiroptera: Emballonuridae) from Costa Rica, assembled from next-generation sequencing data. The mitogenome of Saccopteryx leptura measures 16,577 bp in length, and contains 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. A slight A + T bias was observed in the mitogenome of Saccopteryx leptura with an overall base composition of 31.5% A, 28.3% T, 25.8% C, and 14.2% G, and a GC content of 40.1%. The gene arrangement was identical to that of previously described bat mitogenomes., Competing Interests: The authors declare no conflict of interest., (© 2017 Museum für Naturkunde Berlin.)

Here, we present the first complete mitochondrial genome of the Red Tailed Loach Yasuhikotakia modesta (Teleostei: Botiidae) from Thailand, assembled from next-generation transcriptome sequencing data. The assembled transcript corresponds to the full length mitochondrial genome of Y. modesta , which measured 16,865 bp in length, and contained 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. A slight A + T bias was observed in the mitogenome of Y. modesta with an overall base composition of 32.2% A, 25.8% T, 26.4% C, and 15.4% G, and a GC content of 41.8%. The gene arrangement was identical to that of previously described loach mitogenomes., Competing Interests: This study was funded by the Grant Agency of the Czech Republic, grant number 13 - 37277 S. The authors declare no conflict of interest. The authors alone are responsible for the content and writing of this article., (© 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

The mitochondrial (mt) genome of Telmatobufo australis is a circular molecule of 17,989 bp in length, comprising 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. Gene order and content are identical to those previously reported from other neobatrachian mt genomes. Two protein-coding genes (COI and ATP6) presumably used GTG as start codons while COIII possessed an incomplete stop codon.

Here, we present the first complete mitochondrial genome within the gastropod family Pachychilidae, using the viviparous freshwater snail Tylomelania sarasinorum . This species is a representative member of the lacustrine Tylomelania radiations of the Malili-Lakes-System (Sulawesi, Indonesia). The mitochondrial genome was 16,632 bp long and contained 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. A pronounced A + T bias was observed with an overall base composition of 29.5% A, 35.7% T, 18.3% G and 16.6% C. Tylomelania sarasinorum exhibited a novel mitochondrial gene arrangement, differing from all Caenogastropoda mitochondrial genomes published to date., Competing Interests: The authors report no conflicts of interest. The authors alone are responsible for the content and writing of this article., (© 2016 The Author(s). Published by Taylor & Francis.)

Background: LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology., Results: Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active., Conclusions: The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple "junk DNA" or "harmful genomic parasites" but essential components of the transcriptomic machinery in vertebrates.

We present 5 cases of intrarenal arteriovenous fistulas after renal biopsy. Of these cases 3 followed percutaneous needle biopsy and 2 occurred after open renal biopsy. A nephrectomy was necessary in 1 patient to control the severe hypertension secondary to a hydronephrotic kidney and 1 patient required a partial nephrectomy to control marked hematuria.