Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Coast-to-Coast Spread of SARS-CoV-2 during the Early Epidemic in the United States.

The novel coronavirus SARS-CoV-2 was first detected in the Pacific Northwest region of the United States in January 2020, with subsequent COVID-19 outbreaks detected in all 50 states by early March. To uncover the sources of SARS-CoV-2 introductions and patterns of spread within the United States, we sequenced nine viral genomes from early reported COVID-19 patients in Connecticut. Our phylogenetic analysis places the majority of these genomes with viruses sequenced from Washington state. By coupling our genomic data with domestic and international travel patterns, we show that early SARS-CoV-2 transmission in Connecticut was likely driven by domestic introductions. Moreover, the risk of domestic importation to Connecticut exceeded that of international importation by mid-March regardless of our estimated effects of federal travel restrictions. This study provides evidence of widespread sustained transmission of SARS-CoV-2 within the United States and highlights the critical need for local surveillance.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States.

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

BMP6 and VDR gene polymorphisms are associated with osteonecrosis in a sickle cell anaemia cohort.

The occurrence and severity of osteonecrosis in sickle cell anaemia (SCA) vary due to risk factors, including genetic modifiers. Bone morphogenetic proteins (BMPs), particularly BMP6, and the vitamin D receptor (VDR) play key roles in cartilage and bone metabolism, making them potential contributors to orthopaedic outcomes in SCA. Here, we evaluated the association of polymorphisms in BMP6 (rs3812163, rs270393 and rs449853) and VDR (FokI rs2228570 and Cdx2 rs11568820) genes with osteonecrosis risk in a Brazilian SCA cohort. A total of 177 unrelated SCA patients were selected. The AA genotype of BMP6 rs3812163 was independently associated with a lower osteonecrosis risk (p = 0.015; odds ratio (OR): 0.38; 95% confidence interval (CI): 0.18-0.83) and with the long-term cumulative incidence of osteonecrosis (p = 0.029; hazard ratio: 0.56, 95% CI: 0.34-0.94). The VDR rs2228570 TT genotype was independently associated with a lower osteonecrosis risk (p = 0.039; OR: 0.14; 95% CI: 0.02-0.90). In summary, our results provide evidence that BMP6 rs3812163 and the VDR rs2228570 might be implicated in osteonecrosis pathophysiology in SCA and might help identify individuals at high risk.

Previous fracture and subsequent fracture risk: a meta-analysis to update FRAX.

A large international meta-analysis using primary data from 64 cohorts has quantified the increased risk of fracture associated with a previous history of fracture for future use in FRAX. INTRODUCTION: The aim of this study was to quantify the fracture risk associated with a prior fracture on an international basis and to explore the relationship of this risk with age, sex, time since baseline and bone mineral density (BMD). METHODS: We studied 665,971 men and 1,438,535 women from 64 cohorts in 32 countries followed for a total of 19.5 million person-years. The effect of a prior history of fracture on the risk of any clinical fracture, any osteoporotic fracture, major osteoporotic fracture, and hip fracture alone was examined using an extended Poisson model in each cohort. Covariates examined were age, sex, BMD, and duration of follow-up. The results of the different studies were merged by using the weighted ß-coefficients. RESULTS: A previous fracture history, compared with individuals without a prior fracture, was associated with a significantly increased risk of any clinical fracture (hazard ratio, HR = 1.88; 95% CI = 1.72-2.07). The risk ratio was similar for the outcome of osteoporotic fracture (HR = 1.87; 95% CI = 1.69-2.07), major osteoporotic fracture (HR = 1.83; 95% CI = 1.63-2.06), or for hip fracture (HR = 1.82; 95% CI = 1.62-2.06). There was no significant difference in risk ratio between men and women. Subsequent fracture risk was marginally downward adjusted when account was taken of BMD. Low BMD explained a minority of the risk for any clinical fracture (14%), osteoporotic fracture (17%), and for hip fracture (33%). The risk ratio for all fracture outcomes related to prior fracture decreased significantly with adjustment for age and time since baseline examination. CONCLUSION: A previous history of fracture confers an increased risk of fracture of substantial importance beyond that explained by BMD. The effect is similar in men and women. Its quantitation on an international basis permits the more accurate use of this risk factor in case finding strategies.

Real-time public health communication of local SARS-CoV-2 genomic epidemiology.

Genomic epidemiology can provide a unique, real-time understanding of SARS-CoV-2 transmission patterns. Yet the potential for genomic analyses to guide local policy and community-based behavioral decisions is limited because they are often oriented towards specially trained scientists and conducted on a national or global scale. Here, we propose a new paradigm: Phylogenetic analyses performed on a local level (municipal, county, or state), with results communicated in a clear, timely, and actionable manner to strengthen public health responses. We believe that presenting results rapidly, and tailored to a non-expert audience, can serve as a template for effective public health response to COVID-19 and other emerging viral diseases.

Biomass and carbon balance in a dry tropical forest area in northeast Brazil.

Forest ecosystems play an important role in mitigating the concentration increase of carbon dioxide (CO2) in the atmosphere through carbon sequestration by plants and its storage in biomass and soil. The objective was to determine the aerial biomass carbon stock in a dry tropical forest in Brazil. It was developed between 2012 to 2015, in an area with an advanced regeneration stage (50 years) in the semi-arid region of Pernambuco and it was used 40 permanent plots (400 m²) distant 80 m apart, with 50 m from the border, totaling 1.6 ha of the area to sample the shrubby-arboreal component, where all individuals with circumference at breast height (1.30 m of the soil) equal or greater than 6 cm were identified, measured and labeled in 2012 and remeasured in 2015. It was calculated the biomass and carbon stocks through developed equations available in the literature. The results showed that the total biomass and carbon stock in the first year was 27.97 e 12.92 Mg.ha-1 while in 2015 it was 18.49 and 8.39 Mg.ha-1 respectively. The results showed a biomass and carbon stock reduction of more than 30% in the period evaluated, even this, the area manages to present values within the expected pattern for the region, assuring the importance of sustainable forest management of these native/natural vegetation areas.

Longitudinal Immune Profiling of a Severe Acute Respiratory Syndrome Coronavirus 2 Reinfection in a Solid Organ Transplant Recipient.

BACKGROUND: The underlying immunologic deficiencies enabling severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reinfection are currently unknown. We describe deep longitudinal immune profiling of a transplant recipient hospitalized twice for coronavirus disease 2019 (COVID-19). METHODS: A 66-year-old male renal transplant recipient was hospitalized with COVID-19 March 2020 then readmitted to the hospital with COVID-19 233 days after initial diagnosis. Virologic and immunologic investigations were performed on samples from the primary and secondary infections. RESULTS: Whole viral genome sequencing and phylogenetic analysis revealed that viruses causing both infections were caused by distinct genetic lineages without evidence of immune escape mutations. Longitudinal comparison of cellular and humoral responses during primary SARS-CoV-2 infection revealed that this patient responded to the primary infection with low neutralization titer anti-SARS-CoV-2 antibodies that were likely present at the time of reinfection. CONCLUSIONS: The development of neutralizing antibodies and humoral memory responses in this patient failed to confer protection against reinfection, suggesting that they were below a neutralizing titer threshold or that additional factors may be required for efficient prevention of SARS-CoV-2 reinfection. Development of poorly neutralizing antibodies may have been due to profound and relatively specific reduction in naive CD4 T-cell pools. Seropositivity alone may not be a perfect correlate of protection in immunocompromised patients.

Update of the fracture risk prediction tool FRAX: a systematic review of potential cohorts and analysis plan.

We describe the collection of cohorts together with the analysis plan for an update of the fracture risk prediction tool FRAX with respect to current and novel risk factors. The resource comprises 2,138,428 participants with a follow-up of approximately 20 million person-years and 116,117 documented incident major osteoporotic fractures. INTRODUCTION: The availability of the fracture risk assessment tool FRAX® has substantially enhanced the targeting of treatment to those at high risk of fracture with FRAX now incorporated into more than 100 clinical osteoporosis guidelines worldwide. The aim of this study is to determine whether the current algorithms can be further optimised with respect to current and novel risk factors. METHODS: A computerised literature search was performed in PubMed from inception until May 17, 2019, to identify eligible cohorts for updating the FRAX coefficients. Additionally, we searched the abstracts of conference proceedings of the American Society for Bone and Mineral Research, European Calcified Tissue Society and World Congress of Osteoporosis. Prospective cohort studies with data on baseline clinical risk factors and incident fractures were eligible. RESULTS: Of the 836 records retrieved, 53 were selected for full-text assessment after screening on title and abstract. Twelve cohorts were deemed eligible and of these, 4 novel cohorts were identified. These cohorts, together with 60 previously identified cohorts, will provide the resource for constructing an updated version of FRAX comprising 2,138,428 participants with a follow-up of approximately 20 million person-years and 116,117 documented incident major osteoporotic fractures. For each known and candidate risk factor, multivariate hazard functions for hip fracture, major osteoporotic fracture and death will be tested using extended Poisson regression. Sex- and/or ethnicity-specific differences in the weights of the risk factors will be investigated. After meta-analyses of the cohort-specific beta coefficients for each risk factor, models comprising 10-year probability of hip and major osteoporotic fracture, with or without femoral neck bone mineral density, will be computed. CONCLUSIONS: These assembled cohorts and described models will provide the framework for an updated FRAX tool enabling enhanced assessment of fracture risk (PROSPERO (CRD42021227266)).

COVID-19 one year into the pandemic: from genetics and genomics to therapy, vaccination, and policy.

COVID-19 has engulfed the world and it will accompany us all for some time to come. Here, we review the current state at the milestone of 1 year into the pandemic, as declared by the WHO (World Health Organization). We review several aspects of the on-going pandemic, focusing first on two major topics: viral variants and the human genetic susceptibility to disease severity. We then consider recent and exciting new developments in therapeutics, such as monoclonal antibodies, and in prevention strategies, such as vaccines. We also briefly discuss how advances in basic science and in biotechnology, under the threat of a worldwide emergency, have accelerated to an unprecedented degree of the transition from the laboratory to clinical applications. While every day we acquire more and more tools to deal with the on-going pandemic, we are aware that the path will be arduous and it will require all of us being community-minded. In this respect, we lament past delays in timely full investigations, and we call for bypassing local politics in the interest of humankind on all continents.

Identification and selection of a new Saccharomyces cerevisiae strain isolated from Brazilian ethanol fermentation process for application in beer production.

The fermented beverage industry is always pursuing alternatives to make products that delight consumers with special or unique characteristics. The identification and improvement of new yeast strains emerge as an opportunity; however, wild strains usually have a limitation in maltose fermentation and/or off-flavors production. Here we report the production of a Blond-style ale beer using a bioethanol isolated strain (LBGA-287) with flavor complexity approved in sensorial panels. LBGA-287 also showed an increase in maltose consumption, growth and fermentation rates when compared to the commercial yeast. Using qPCR analysis, genes related to the (i) efficiency of fermentation (ii) production of aromas/off-flavors, and (iii) metabolization of carbohydrates were found as differentially expressed in the isolated strains when compared to industrial yeast. This suggests that LBGA-287 could have an important impact on beer production, improving brewing efficiency, quality and diversity of this beverage, and most importantly satisfying the final consumer.

Comparison of blind and endoscopic-guided orotracheal intubation on laryngeal and tracheal damage in domestic rabbits (Oryctolagus cuniculus).

OBJECTIVE: To compare blind and endoscopic-guided techniques for orotracheal intubation in rabbits and the number of intubation attempts with laryngeal/tracheal damage. STUDY DESIGN: Prospective, randomized experimental study. ANIMALS: A total of 24 healthy, intact female New Zealand White rabbits, weighing 2.2 ± 0.2 kg (mean ± standard deviation). METHODS: Rabbits were randomly assigned to blind (group B) or endoscopic-guided (group E) orotracheal intubation with a 2.0 mm internal diameter uncuffed tube. Intramuscular (IM) alfaxalone (7 mg kg-1), hydromorphone (0.1 mg kg-1) and dexmedetomidine (0.005 mg kg-1) were administered, and additional IM alfaxalone (3-5 mg kg-1) and dexmedetomidine (0.025 mg kg-1) were administered to rabbits with strong jaw tone. An intubation attempt was defined as the advancement of the endotracheal tube from the incisors to the laryngeal entrance. Tracheal intubation was confirmed via capnography and anesthesia was maintained with isoflurane for 2 hours. Following euthanasia, laryngeal and tracheal tissues were submitted for histopathology. Quality of anesthesia for orotracheal intubation, intubation procedure and tissue damage were numerically scored. Data were analyzed using Poisson regression, Spearman's correlation, t test, mixed anova, Mann-Whitney U test, Friedman and Chi square tests as appropriate. RESULTS: Median (range) intubation attempts were 2 (1-8) and 1 (1-3) for groups B and E, respectively. More rabbits in group E (91.6%) required additional alfaxalone and dexmedetomidine than in group B (16.7%). Median (range) cumulative histopathology scores were 6 (3-10) and 6 (2-9) for groups B and E, respectively. Scores were highest in the cranial trachea, but there was no difference between groups and no correlation between laryngeal/tracheal damage and the number of intubation attempts. CONCLUSIONS AND CLINICAL RELEVANCE: Both orotracheal intubation techniques were associated with laryngeal/tracheal damage. Although blind orotracheal intubation was associated with a higher number of attempts, the tissue damage was similar between groups.

Fungal Planet description sheets: 1436-1477.

Novel species of fungi described in this study include those from various countries as follows: Argentina, Colletotrichum araujiae on leaves, stems and fruits of Araujia hortorum. Australia, Agaricus pateritonsus on soil, Curvularia fraserae on dying leaf of Bothriochloa insculpta, Curvularia millisiae from yellowing leaf tips of Cyperus aromaticus, Marasmius brunneolorobustus on well-rotted wood, Nigrospora cooperae from necrotic leaf of Heteropogon contortus, Penicillium tealii from the body of a dead spider, Pseudocercospora robertsiorum from leaf spots of Senna tora, Talaromyces atkinsoniae from gills of Marasmius crinis-equi and Zasmidium pearceae from leaf spots of Smilaxglyciphylla. Brazil, Preussia bezerrensis from air. Chile, Paraconiothyrium kelleni from the rhizosphere of Fragaria chiloensis subsp. chiloensis f. chiloensis. Finland, Inocybe udicola on soil in mixed forest with Betula pendula, Populus tremula, Picea abies and Alnus incana. France, Myrmecridium normannianum on dead culm of unidentified Poaceae. Germany, Vexillomyces fraxinicola from symptomless stem wood of Fraxinus excelsior. India, Diaporthe limoniae on infected fruit of Limonia acidissima, Didymella naikii on leaves of Cajanus cajan, and Fulvifomes mangroviensis on basal trunk of Aegiceras corniculatum. Indonesia, Penicillium ezekielii from Zea mays kernels. Namibia, Neocamarosporium calicoremae and Neocladosporium calicoremae on stems of Calicorema capitata, and Pleiochaeta adenolobi on symptomatic leaves of Adenolobus pechuelii. Netherlands, Chalara pteridii on stems of Pteridium aquilinum, Neomackenziella juncicola (incl. Neomackenziella gen. nov.) and Sporidesmiella junci from dead culms of Juncus effusus. Pakistan, Inocybe longistipitata on soil in a Quercus forest. Poland, Phytophthora viadrina from rhizosphere soil of Quercus robur, and Septoria krystynae on leaf spots of Viscum album. Portugal (Azores), Acrogenospora stellata on dead wood or bark. South Africa, Phyllactinia greyiae on leaves of Greyia sutherlandii and Punctelia anae on bark of Vachellia karroo. Spain, Anteaglonium lusitanicum on decaying wood of Prunus lusitanica subsp. lusitanica, Hawksworthiomyces riparius from fluvial sediments, Lophiostoma carabassense endophytic in roots of Limbarda crithmoides, and Tuber mohedanoi from calcareus soils. Spain (Canary Islands), Mycena laurisilvae on stumps and woody debris. Sweden, Elaphomyces geminus from soil under Quercus robur. Thailand, Lactifluus chiangraiensis on soil under Pinus merkusii, Lactifluus nakhonphanomensis and Xerocomus sisongkhramensis on soil under Dipterocarpus trees. Ukraine, Valsonectria robiniae on dead twigs of Robinia hispida. USA, Spiralomyces americanus (incl. Spiralomyces gen. nov.) from office air. Morphological and culture characteristics are supported by DNA barcodes. Citation: Tan YP, Bishop-Hurley SL, Shivas RG, et al. 2022. Fungal Planet description sheets: 1436-1477. Persoonia 49: 261-350. https://doi.org/10.3767/persoonia.2022.49.08.

Multiple Transmission Chains within COVID-19 Cluster, Connecticut, USA, 20201.

In fall 2020, a coronavirus disease cluster comprising 16 cases occurred in Connecticut, USA. Epidemiologic and genomic evidence supported transmission among persons at a school and fitness center but not a workplace. The multiple transmission chains identified within this cluster highlight the necessity of a combined investigatory approach.

Twenty years of West Nile virus spread and evolution in the Americas visualized by Nextstrain.

It has been 20 years since West Nile virus first emerged in the Americas, and since then, little progress has been made to control outbreaks caused by this virus. After its first detection in New York in 1999, West Nile virus quickly spread across the continent, causing an epidemic of human disease and massive bird die-offs. Now the virus has become endemic to the United States, where an estimated 7 million human infections have occurred, making it the leading mosquito-borne virus infection and the most common cause of viral encephalitis in the country. To bring new attention to one of the most important mosquito-borne viruses in the Americas, we provide an interactive review using Nextstrain: a visualization tool for real-time tracking of pathogen evolution (nextstrain.org/WNV/NA). Nextstrain utilizes a growing database of more than 2,000 West Nile virus genomes and harnesses the power of phylogenetics for students, educators, public health workers, and researchers to visualize key aspects of virus spread and evolution. Using Nextstrain, we use virus genomics to investigate the emergence of West Nile virus in the U S, followed by its rapid spread, evolution in a new environment, establishment of endemic transmission, and subsequent international spread. For each figure, we include a link to Nextstrain to allow the readers to directly interact with and explore the underlying data in new ways. We also provide a brief online narrative that parallels this review to further explain the data and highlight key epidemiological and evolutionary features (nextstrain.org/narratives/twenty-years-of-WNV). Mirroring the dynamic nature of outbreaks, the Nextstrain links provided within this paper are constantly updated as new West Nile virus genomes are shared publicly, helping to stay current with the research. Overall, our review showcases how genomics can track West Nile virus spread and evolution, as well as potentially uncover novel targeted control measures to help alleviate its public health burden.

Spatial multicriteria approach to support water resources management with multiple sources in semi-arid areas in Brazil.

Semi-arid regions often face severe drought events that reduce agricultural and livestock production. In recent years, some international studies have used multicriteria decision analysis (MCDA) approaches combined with geographic information systems (GIS-MCDA) to support decision-makers in assessing the suitability of agricultural land for irrigation in semi-arid regions. Unlike previous studies, which have only considered a single source of water for crop irrigation, this study proposes a GIS-MCDA approach that considers all potentially available local water sources (e.g., groundwater, surface water, and wastewater) as possible alternatives for better multisource water resource management (MWRM) in regions facing water shortages. The geospatial multicriteria evaluation implemented in this study considers a series of technical, environmental, and agricultural productivity criteria using the analytical hierarchy process (AHP) method. Three independent baseline maps were generated, showing the spatial distribution of suitable areas for crop irrigation for each considered water source in the studied area. Surface water, groundwater, and wastewater offered suitable crop irrigation for 83%, 70%, and 26% of the study area, respectively. Overlapping these areas produced a final map showing all the feasible areas for each crop irrigation alternative at the same time. The MWRM approach considering all water sources increased the coverage of suitable areas to be irrigated in the study area by 2.2%, 20.4%, and more than 225% compared to considering surface water, groundwater, and wastewater, respectively, independently. The GIS-MCDA framework proposed in this study provides better support for decision-makers and stakeholders, favouring a reduction in possible conflicts over water scarcity, the diversification of irrigated crops, and an improvement in the quality-quantitative management of water resources in semi-arid regions.

Amazon plant diversity revealed by a taxonomically verified species list.

Recent debates on the number of plant species in the vast lowland rain forests of the Amazon have been based largely on model estimates, neglecting published checklists based on verified voucher data. Here we collate taxonomically verified checklists to present a list of seed plant species from lowland Amazon rain forests. Our list comprises 14,003 species, of which 6,727 are trees. These figures are similar to estimates derived from nonparametric ecological models, but they contrast strongly with predictions of much higher tree diversity derived from parametric models. Based on the known proportion of tree species in neotropical lowland rain forest communities as measured in complete plot censuses, and on overall estimates of seed plant diversity in Brazil and in the neotropics in general, it is more likely that tree diversity in the Amazon is closer to the lower estimates derived from nonparametric models. Much remains unknown about Amazonian plant diversity, but this taxonomically verified dataset provides a valid starting point for macroecological and evolutionary studies aimed at understanding the origin, evolution, and ecology of the exceptional biodiversity of Amazonian forests.

Global gene expression reveals an increase of HMGB1 and APEX1 proteins and their involvement in oxidative stress, apoptosis and inflammation pathways among beta-thalassaemia intermedia and major phenotypes.

Beta-thalassaemia (BT) is classified according to blood transfusion requirement as minor (BTMi), intermedia (BTI) and major (BTM). BTM is the most severe form, requiring regular transfusions while transfusion need is only occasional in BTI. Differential gene expression between patients has not been assessed so far. Here, we evaluated the global gene expression profiles during differentiation of human erythroid cells of two patients carrying the same mutation [CD39, (C → T)], though displaying different phenotypes (BTI and BTM). Considering the role of reactive oxygen species (ROS) in the pathophysiology of thalassaemia, we focused on differentially expressed genes involved in metabolic pathways triggered by ROS, such as inflammation and apoptosis, and, from these, we selected the Apurinic/Apyrimidinic Endodeoxyribonuclease 1 (APEX1) and High Mobility Group Box1 (HMGB1) genes, whose role in BT is not well established. An in-depth expression analysis of transcriptional and protein levels in patients carrying a range of mutations associated with BT phenotypes indicated that APEX1 was increased in both BTI and BTM. Furthermore, higher amounts of HMGB1 was found in the plasma of BTI patients. Our findings suggest that these proteins have important roles in BT and could represent new targets for further studies aiming to improve the management of the disease.

Complex taxonomy of the 'brush tail' peregrine earthworm Pontoscolex corethrurus.

Pontoscolex corethrurus is the most widespread earthworm species in tropical and sub-tropical zones and one of the most studied in soil science. Although, ecological interactions of P. corethrurus with its environment are well documented, the taxonomic status of the species remains unclear. In this study, we investigated phylogenetic relationships within the genus Pontoscolex, in particular focusing on morphologically indistinguishable (i.e., cryptic) lineages. A total of 792 specimens collected from 25 different countries and islands all over the world were analyzed using two mitochondrial (COI and 16S rDNA) and two nuclear (internal transcribed spacers 2 and 28S rDNA) markers, and a total of 11 morphological characters both internal and external were investigated in all genetically characterized lineages. A large-scale multilocus sequence data matrix was also obtained for Pontoscolex spp. specimens using the Anchored Hybrid Enrichment (AHE) method. Multilocus phylogenetic and phylogenomic analyses, combined with species delimitation methods; including single locus (mPTP, ABGD) and multilocus (BPP) approaches, revealed congruent results. Four cryptic species were supported within the P. corethrurus species complex, and four potentially new species within the genus Pontoscolex. One widespread lineage (L1), within P. corethrurus complex was observed in the current population of Fritz Müller's garden where P. corethrurus was first described in 1856. Cryptic lineages were observed in sympatry at several localities. This, in combination with observed heteroplasmy in COI gene in one population raises an important question of reproductive isolation between these species.