Lining the small intestine with mycobacteriophages protects from Mycobacterium avium subsp. paratuberculosis and eliminates fecal shedding.

Johne's disease (JD), a chronic, infectious enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), affects wild and domestic ruminants. There is no cure or effective prevention, and current vaccines have substantial limitations, leaving this disease widespread in all substantial dairy industries causing economic, and animal welfare implications. Mycobacteriophages (MPs) have been gaining interest in recent years and are proposed as a promising solution to curtailing MAP infection. Using a well-validated infection model, we have demonstrated the preventative potential of MPs to protect dairy calves against MAP infection. Calves were supplemented daily with a phage cocktail from birth till weaning at 2 m of age and inoculated with MAP at 2 wk of age. Infection status was measured for 4.5 mo through blood, fecal, and postmortem tissue samples. Our findings highlight the remarkable efficacy of orally administered MPs. Notably, fecal shedding of MAP was entirely eliminated within 10 wk, in contrast to the infected control group where shedding continued for the entirety of the trial period. Postmortem tissue culture analysis further supported the effectiveness of MPs, with only 1 out of 6 animals in the phage-treated group testing positive for MAP colonized tissues compared to 6 out of 6 animals in the infected control group. Additionally, plaque assay results demonstrated the ability of phages to persist within the intestinal tract. Collectively, these results underscore the potential of orally administered MP cocktails as a highly effective intervention strategy to combat JD in dairy calves and by extension in the dairy industry.

Phylogenomics and pervasive genome-wide phylogenetic discordance among fin whales (Balaenoptera physalus).

Phylogenomics has the power to uncover complex phylogenetic scenarios across the genome. In most cases, no single topology is reflected across the entire genome as the phylogenetic signal differs among genomic regions due to processes, such as introgression and incomplete lineage sorting. Baleen whales are among the largest vertebrates on Earth with a high dispersal potential in a relatively unrestricted habitat, the oceans. The fin whale (Balaenoptera physalus) is one of the most enigmatic baleen whale species, currently divided into four subspecies. It has been a matter of debate whether phylogeographic patterns explain taxonomic variation in fin whales. Here we present a chromosome-level whole genome analysis of the phylogenetic relationships among fin whales from multiple ocean basins. First, we estimated concatenated and consensus phylogenies for both the mitochondrial and nuclear genomes. The consensus phylogenies based upon the autosomal genome uncovered monophyletic clades associated with each ocean basin, aligning with the current understanding of subspecies division. Nevertheless, discordances were detected in the phylogenies based on the Y chromosome, mitochondrial genome, autosomal genome and X chromosome. Furthermore, we detected signs of introgression and pervasive phylogenetic discordance across the autosomal genome. This complex phylogenetic scenario could be explained by a puzzle of introgressive events, not yet documented in fin whales. Similarly, incomplete lineage sorting and low phylogenetic signal could lead to such phylogenetic discordances. Our study reinforces the pitfalls of relying on concatenated or single locus phylogenies to determine taxonomic relationships below the species level by illustrating the underlying nuances which some phylogenetic approaches may fail to capture. We emphasize the significance of accurate taxonomic delineation in fin whales by exploring crucial information revealed through genome-wide assessments.

Lacticaseibacillus parahuelsenbergensis sp. nov., Lacticaseibacillus styriensis sp. nov. and Lacticaseibacillus zeae subsp. silagei subsp. nov., isolated from different grass and corn silage.

Four rod-shaped, non-motile, non-spore-forming, facultative anaerobic, Gram-stain-positive lactic acid bacteria, designated as EB0058T, SCR0080, LD0937T and SCR0063T, were isolated from different corn and grass silage samples. The isolated strains were characterized using a polyphasic approach and EB0058T and SCR0080 were identified as Lacticaseibacillus zeae by 16S rRNA gene sequence analysis. Based on whole-genome sequence-based characterization, EB0058T and SCR0080 were separated into a distinct clade from Lacticaseibacillus zeae DSM 20178T, together with CECT9104 and UD2202, whose genomic sequences are available from NCBI GenBank. The average nucleotide identity (ANI) values within the new subgroup are 99.9 % and the digital DNA-DNA hybridization (dDDH) values are 99.3-99.9 %, respectively. In contrast, comparison of the new subgroup with publicly available genomic sequences of L. zeae strains, including the type strain DSM 20178T, revealed dDDH values of 70.2-72.5 % and ANI values of 96.2-96.6 %. Based on their chemotaxonomic, phenotypic and phylogenetic characteristics, EB0058T and SCR0080 represent a new subspecies of L. zeae. The name Lacticaseibacillus zeae subsp. silagei subsp. nov. is proposed with the type strain EB0058T (=DSM 116376T=NCIMB 15474T). According to the results of 16S rRNA gene sequencing, LD0937T and SCR0063T are members of the Lacticaseibacillus group. The dDDH value between the isolates LD0937T and SCR0063T was 67.6 %, which is below the species threshold of 70 %, clearly showing that these two isolates belong to different species. For both strains, whole genome-sequencing revealed that the closest relatives within the Lacticaseibacillus group were Lacticaseibacillus huelsenbergensis DSM 115425 (dDDH 66.5 and 65.9 %) and Lacticaseibacillus casei DSM 20011T (dDDH 64.1 and 64.9 %). Based on the genomic, chemotaxonomic and morphological data obtained in this study, two novel species, Lacticaseibacillus parahuelsenbergensis sp. nov. and Lacticaseibacillus styriensis sp. nov. are proposed and the type strains are LD0937T (=DSM 116105T=NCIMB 15471T) and SCR0063T (=DSM 116297T=NCIMB 15473T), respectively.

A lineage perspective on hominin taxonomy and evolution.

An uncritical reliance on the phylogenetic species concept has led paleoanthropologists to become increasingly typological in their delimitation of new species in the hominin fossil record. As a practical matter, this approach identifies species as diagnosably distinct groups of fossils that share a unique suite of morphological characters but, ontologically, a species is a metapopulation lineage segment that extends from initial divergence to eventual extinction or subsequent speciation. Working from first principles of species concept theory, it is clear that a reliance on morphological diagnosabilty will systematically overestimate species diversity in the fossil record; because morphology can evolve within a lineage segment, it follows that early and late populations of the same species can be diagnosably distinct from each other. We suggest that a combination of morphology and chronology provides a more robust test of the single-species null hypothesis than morphology alone.

Streptococcus equi subspecies equi from strangles suspected equines: molecular detection, antibiogram profiles and risk factors.

Strangles, caused by Streptococcus equi subspecies equi, is a highly infectious disease of equines causing major health issues and financial losses. The aim of the study was to detect the presence of the SeM gene in Streptococcus equi isolated from equine suspected of having strangles. A cross-sectional study design was conducted from July to December 2022 in five districts of the central Gondar zone, Ethiopia. One-hundred sixty swab samples were taken from animals that had been clinically suspected. The SeM gene was detected using polymerase chain reaction, and the antimicrobial susceptibility test was performed using the Kirby-Bauer disc diffusion method. The binary logistic regression model was employed to test for statistical significance. In 31.87% (51/160) of the samples, Streptococcus equi species were isolated, and 31.37% (16/51) of these species carried the SeM gene. There was a significant amount of tetracycline (81.5%), erythromycin (81.5%), and vancomycin (75.5%) resistance among the 16 isolates. Strangles were more likely to be present in animals who shared feed containers (AOR = 7.59; 95% CI = 1.44-39.93), drank from the same water troughs (AOR = 7.74; 95% CI = 1.44-41.01), and spent the night together (AOR = 5.97; 95% CI 1.41-25.37). The findings of this study showed that the research areas harboured Streptococcus equi subspecies equi. Sharing feed containers and water troughs were potential sources of strangles infection; thus, these containers need to be cleaned regularly.

Pyelonephritis with bacteremia caused by Salmonella Choleraesuis in a Japanese patient with carcinoma of unknown primary origin: A case report.

Salmonella enterica subspecies enterica serovar Choleraesuis (S. Choleraesuis) is a nontyphoidal Salmonella pathogen that causes swine paratyphoids. S. Choleraesuis is a zoonotic pathogen transmitted to humans via contaminated food and causes sepsis. Here, we report a rare case of pyelonephritis caused by S. Choleraesuis in a Japanese patient with a carcinoma of unknown primary origin. On the day of admission, the patient was diagnosed with pyelonephritis associated with ureteral stent obstruction. He had no history of raw pork consumption or gastrointestinal symptoms. Gram-negative rods were isolated from urine and blood cultures, identified as Salmonella enterica subsp. enterica using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. The serological typing results were O7: -: 1 and 5; however, the serotypes could not be determined. The isolate was identified as S. Choleraesuis using multilocus sequence typing, nucleotide sequence analysis of the fliC gene, and biochemical examination. Four days after a 14-day course of intravenous piperacillin-tazobactam (9 g/day), the patient showed relapse of the condition. Subsequently, the patient was treated with intravenous ceftriaxone (2 g/day) and oral amoxicillin (1000 mg/day) for 14 days each; recurrence was not observed. This novel case of pyelonephritis with bacteremia was caused by S. Choleraesuis in Japan. Conventional testing methods could not identify the serotypes; however, the case highlights the importance of adopting advanced diagnostic techniques based on molecular biology to ensure accurate pathogen identification.

Phylogeography, genetic diversity, and intraspecific genetic structure of the black-horned capuchin (Sapajus nigritus).

Recent molecular studies have clarified the overarching taxonomy of capuchin monkeys, but intraspecific genetic diversity remains unexplored for most capuchin species. One example is Sapajus nigritus, the southernmost capuchin monkey, found in Brazil and Argentina; its phenotypic diversity has been recognized as two geographic subspecies, but the intraspecific genetic structure of this taxon is poorly known. Here, we sampled across most of this species' geographic distribution, producing a newly sequenced data set for genetic analyses that included 78 individuals from 14 populations. We investigated the intraspecific diversity, genetic structure, and evolutionary history using three mitochondrial markers. Our results indicated that S. nigritus populations exhibited high levels of genetic structure. We found strong support for two monophyletic clades within this species with a deep phylogenetic split, and clear separation from other related taxa. Vicariance events seem to have played a prevalent role in shaping S. nigritus genetic differentiation. The Paraíba do Sul River may have driven the deep divergence between southern and northern clades, whereas the Tietê River may have had a weaker, more recent effect on the divergence of populations within the southern clade.

Phylogenomics of the Olea europaea complex using 15 whole genomes supports recurrent genetic admixture together with differentiation into seven subspecies.

BACKGROUND: The last taxonomic account of Olea recognises six subspecies within Olea europaea L., including the Mediterranean olive tree (subsp. europaea) and five other subspecies (laperrinei, guanchica, maroccana, cerasiformis, and cuspidata) distributed across the Old World, including Macaronesian islands. The evolutionary history of this monophyletic group (O. europaea complex) has revealed a reticulated scenario involving hybridization and polyploidization events, leading to the presence of a polyploid series associated with the subspecies. However, how the polyploids originated, and how the different subspecies contributed to the domestication of the cultivated olive are questions still debated. Tracing the recent evolution and genetic diversification of the species is key for the management and preservation of its genetic resources. To study the recent history of the O. europaea complex, we compared newly sequenced and available genomes for 27 individuals representing the six subspecies. RESULTS: Our results show discordance between current subspecies distributions and phylogenomic patterns, which support intricate biogeographic patterns. The subspecies guanchica, restricted to the Canary Islands, is closely related to subsp. europaea, and shows a high genetic diversity. The subsp. laperrinei, restricted now to high mountains of the Sahara desert, and the Canarian subsp. guanchica contributed to the formation of the allotetraploid subsp. cerasiformis (Madeira islands) and the allohexaploid subsp. maroccana (western Sahara region). Our phylogenomic data support the recognition of one more taxon (subsp. ferruginea) for the Asian populations, which is clearly segregated from the African subsp. cuspidata. CONCLUSIONS: In sum, the O. europaea complex underwent several processes of hybridization, polyploidy, and geographical isolation resulting in seven independent lineages with certain morphological traits recognised into subspecies.

Alcohol Consumption, High-Density Lipoprotein Particles and Subspecies, and Risk of Cardiovascular Disease: Findings from the PREVEND Prospective Study.

The associations of HDL particle (HDL-P) and subspecies concentrations with alcohol consumption are unclear. We aimed to evaluate the interplay between alcohol consumption, HDL parameters and cardiovascular disease (CVD) risk. In the PREVEND study of 5151 participants (mean age, 53 years; 47.5% males), self-reported alcohol consumption and HDL-P and subspecies (small, medium, and large) by nuclear magnetic resonance spectroscopy were assessed. Hazard ratios (HRs) with 95% CIs for first CVD events were estimated. In multivariable linear regression analyses, increasing alcohol consumption increased HDL-C, HDL-P, large and medium HDL, HDL size, and HDL subspecies (H3P, H4P, H6 and H7) in a dose-dependent manner. During a median follow-up of 8.3 years, 323 first CVD events were recorded. Compared with abstainers, the multivariable adjusted HRs (95% CIs) of CVD for occasional to light, moderate, and heavy alcohol consumers were 0.72 (0.55-0.94), 0.74 (0.54-1.02), and 0.65 (0.38-1.09), respectively. These associations remained consistent on additional adjustment for each HDL parameter. For CVD, only HDL-C was associated with a statistically significant decreased risk of CVD in a fully adjusted analysis (HR 0.84, 95% CI 0.72-0.97 per 1 SD increment). For coronary heart disease, HDL-C, HDL-P, medium HDL, HDL size, and H4P showed inverse associations, whereas HDL-C and HDL size modestly increased stroke risk. Except for H6P, alcohol consumption did not modify the associations between HDL parameters and CVD risk. The addition of HDL-C, HDL size, or H4P to a CVD risk prediction model containing established risk factors improved risk discrimination. Increasing alcohol consumption is associated with increased HDL-C, HDL-P, large and medium HDL, HDL size, and some HDL subspecies. Associations of alcohol consumption with CVD are largely independent of HDL parameters. The associations of HDL parameters with incident CVD are generally not attenuated or modified by alcohol consumption.

A Rough Colony Morphology of Mycobacterium abscessus Is Associated With Cavitary Pulmonary Disease and Poor Clinical Outcome.

BACKGROUND: The Mycobacterium abscessus complex (MABC) is a difficult to treat mycobacterium with two distinct morphologies: smooth and rough. As the clinical implications are unclear, we explored the morphology of MABC in relation to disease and outcome. METHODS: We performed a retrospective multicenter cohort study including patients with confirmed MABC in Sweden, 2009-2020, with treatment outcome as the primary outcome. MABC colony morphology was determined by light microscopy on Middlebrook 7H10 agar plates. RESULTS: Of the 71 MABC isolates, a defined morphology could be determined for 63 isolates, of which 40 were smooth (56%) and 23 were rough (32%). Immunosuppression, pulmonary disease, and cavitary lesion on chest radiographs were significantly associated with a rough isolate morphology. Participants with smooth isolates had more favorable treatment outcomes (12/14, 86%) compared to those with rough isolates (3/10, 30%). In an age-adjusted logistic regression, rough morphology of MABC was associated to lower odds of clinical cure compared to smooth morphology (adjusted odds ratio, 0.12; P = .049). CONCLUSIONS: Study participants with rough MABC colony morphology of isolates had a worse clinical outcome compared to those with smooth isolates. The biological mechanisms should be further characterized and colony morphology of MABC taken into account during clinical management.

Cold-adaptive evolution at the reproductive stage in Geng/japonica subspecies reveals the role of OsMAPK3 and OsLEA9.

Cold stress at the reproductive stage severely affects the production and geographic distribution of rice. The Geng/japonica subpopulation gradually developed stronger cold adaptation than the Xian/indica subpopulation during the long-term domestication of cultivated rice. However, the evolutionary path and natural alleles underlying the cold adaptability of intra-Geng subspecies remain largely unknown. Here, we identified MITOGEN-ACTIVATED PROTEIN KINASE 3 (OsMAPK3) and LATE EMBRYOGENESIS ABUNDANT PROTEIN 9 (OsLEA9) as two important regulators for the cold adaptation of Geng subspecies from a combination of transcriptome analysis and genome-wide association study. Transgenic validation showed that OsMAPK3 and OsLEA9 confer cold tolerance at the reproductive stage. Selection and evolution analysis suggested that the Geng version of OsMAPK3 (OsMAPK3Geng ) directly evolved from Chinese Oryza rufipogon III and was largely retained in high-latitude and high-altitude regions with low temperatures during domestication. Later, the functional nucleotide polymorphism (FNP-776) in the Kunmingxiaobaigu and Lijiangxiaoheigu version of the OsLEA9 (OsLEA9KL ) promoter originated from novel variation of intra-Geng was selected and predominantly retained in temperate Geng to improve the adaptation of Geng together with OsMAPK3Geng to colder climatic conditions in high-latitude areas. Breeding potential analysis suggested that pyramiding of OsMAPK3Geng and OsLEA9KL enhanced the cold tolerance of Geng and promotes the expansion of cultivated rice to colder regions. This study not only highlights the evolutionary path taken by the cold-adaptive differentiation of intra-Geng, but also provides new genetic resources for rice molecular breeding in low-temperature areas.

Long read isoform sequencing reveals hidden transcriptional complexity between cattle subspecies.

The Iso-Seq method of full-length cDNA sequencing is suitable to quantify differentially expressed genes (DEGs), transcripts (DETs) and transcript usage (DTU). However, the higher cost of Iso-Seq relative to RNA-seq has limited the comparison of both methods. Transcript abundance estimated by RNA-seq and deep Iso-Seq data for fetal liver from two cattle subspecies were compared to evaluate concordance. Inter-sample correlation of gene- and transcript-level abundance was higher within technology than between technologies. Identification of DEGs between the cattle subspecies depended on sequencing method with only 44 genes identified by both that included 6 novel genes annotated by Iso-Seq. There was a pronounced difference between Iso-Seq and RNA-seq results at transcript-level wherein Iso-Seq revealed several magnitudes more transcript abundance and usage differences between subspecies. Factors influencing DEG identification included size selection during Iso-Seq library preparation, average transcript abundance, multi-mapping of RNA-seq reads to the reference genome, and overlapping coordinates of genes. Some DEGs called by RNA-seq alone appear to be sequence duplication artifacts. Among the 44 DEGs identified by both technologies some play a role in immune system, thyroid function and cell growth. Iso-Seq revealed hidden transcriptional complexity in DEGs, DETs and DTU genes between cattle subspecies previously missed by RNA-seq.

Human Tularemia Epididymo-Orchitis Caused by Francisella tularensis Subspecies holartica, Austria.

A previously healthy man in Austria had tularemia epididymo-orchitis develop, leading to unilateral orchiectomy. Francisella tularensis subspecies holartica was detected by 16S rRNA gene sequencing analysis of inflamed granulomatous testicular tissue. Clinicians should suspect F. tularensis as a rare etiologic microorganism in epididymo-orchitis patients with relevant risk factors.

Clonal Expansion of Multidrug-Resistant Streptococcus dysgalactiae Subspecies equisimilis Causing Bacteremia, Japan, 2005-2021.

Incidence of Streptococcus dysgalactiae subspecies equisimilis (SDSE) bacteremia is increasing in the Kyoto-Shiga region of Japan. We retrospectively analyzed clinical features of SDSE bacteremia and conducted comparative genomic analyses of isolates collected from 146 bacteremia episodes among 133 patients during 2005-2021. Of those patients, 7.7% required vasopressor support, and 7.0% died while in the hospital. The prevalence of isolates resistant to erythromycin, minocycline, and clindamycin increased from 8.6% during 2005-2017 to 21.6% during 2018-2021. Our genomic analysis demonstrated that sequence type 525 and clonal complex 25 were predominant in SDSE isolates collected during 2018-2021. In addition, those isolates had acquired 2 antimicrobial-resistance genes, ermB and tetM, via Tn916-like integrative and conjugative elements (ICEs). Phylogenetic analysis revealed clonal distribution of Tn916-like ICEs in SDSE isolates. Our findings suggest that Tn916-like ICEs contributed to the emergence and recent increase of multidrug-resistant SDSE bacteremia in this region of Japan.

Photobacterium damselae subspecies damselae Pneumonia in Dead, Stranded Bottlenose Dolphin, Eastern Mediterranean Sea.

Photobacterium damselae subspecies damselae, an abundant, generalist marine pathogen, has been reported in various cetaceans worldwide. We report a bottlenose dolphin in the eastern Mediterranean Sea that was found stranded and dead. The dolphin had a severe case of chronic suppurative pneumonia and splenic lymphoid depletion caused by this pathogen.

Severe Streptococcus equi Subspecies zooepidemicus Outbreak from Unpasteurized Dairy Product Consumption, Italy.

During November 2021-May 2022, we identified 37 clinical cases of Streptococcus equi subspecies zooepidemicus infections in central Italy. Epidemiologic investigations and whole-genome sequencing showed unpasteurized fresh dairy products were the outbreak source. Early diagnosis by using sequencing technology prevented the spread of life-threatening S. equi subsp. zooepidemicus infections.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica).

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here, we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (<20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (YD) (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Identification of Mycobacterium abscessus Subspecies by MALDI-TOF Mass Spectrometry and Machine Learning.

Mycobacterium abscessus is one of the most common and pathogenic nontuberculous mycobacteria (NTM) isolated in clinical laboratories. It consists of three subspecies: M. abscessus subsp. abscessus, M. abscessus subsp. bolletii, and M. abscessus subsp. massiliense. Due to their different antibiotic susceptibility pattern, a rapid and accurate identification method is necessary for their differentiation. Although matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) has proven useful for NTM identification, the differentiation of M. abscessus subspecies is challenging. In this study, a collection of 325 clinical isolates of M. abscessus was used for MALDI-TOF MS analysis and for the development of machine learning predictive models based on MALDI-TOF MS protein spectra. Overall, using a random forest model with several confidence criteria (samples by triplicate and similarity values >60%), a total of 96.5% of isolates were correctly identified at the subspecies level. Moreover, an improved model with Spanish isolates was able to identify 88.9% of strains collected in other countries. In addition, differences in culture media, colony morphology, and geographic origin of the strains were evaluated, showing that the latter had an impact on the protein spectra. Finally, after studying all protein peaks previously reported for this species, two novel peaks with potential for subspecies differentiation were found. Therefore, machine learning methodology has proven to be a promising approach for rapid and accurate identification of subspecies of M. abscessus using MALDI-TOF MS.

Differences in the suitability of published honey bee (Apis mellifera) reference genes between the African subspecies Apis mellifera scutellata and European derived Apis mellifera.

Quantitative real-time polymerase chain reaction (qPCR) is a method widely used to determine changes and differences in gene expression. As target gene expression is most often quantified relative to the expression of reference genes, the validation of suitable reference genes is of critical importance. In practice, however, such validation might not be thoroughly conducted if the same species and the same tissue or body parts are used for qPCR experiments. Here we show, that qPCR reference genes published for workers of European honey bee (Apis mellifera) subspecies fail to be stably expressed in workers of the African subspecies Apis mellifera scutellata. This is the case even when the sampled workers are in the same life stage, the same organ was dissected and the same reagents were used. Thus, reference genes need to be thoroughly re-tested before they can be used as suitable references even when the only thing that changes is the subspecies used.

Evolutionary history of an island endemic, the Azorean common quail.

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.