Elucidating the Role of OXPHOS Variants in Asthenozoospermia: Insights from Whole Genome Sequencing and an In Silico Analysis.

Infertility is a global health challenge that affects an estimated 72.4 million people worldwide. Between 30 and 50% of these cases involve male factors, showcasing the complex nature of male infertility, which can be attributed to both environmental and genetic determinants. Asthenozoospermia, a condition characterized by reduced sperm motility, stands out as a significant contributor to male infertility. This study explores the involvement of the mitochondrial oxidative phosphorylation (OXPHOS) system, crucial for ATP production and sperm motility, in asthenozoospermia. Through whole-genome sequencing and in silico analysis, our aim was to identify and characterize OXPHOS gene variants specific to individuals with asthenozoospermia. Our analysis identified 680,099 unique variants, with 309 located within OXPHOS genes. Nine of these variants were prioritized due to their significant implications, such as potential associations with diseases, effects on gene expression, protein function, etc. Interestingly, none of these variants had been previously associated with male infertility, opening up new avenues for research. Thus, through our comprehensive approach, we provide valuable insights into the genetic factors that influence sperm motility, laying the foundation for future research in the field of male infertility.

Emergence of OXA-48-like producing Citrobacter species, Germany, 2011 to 2022.

BackgroundCarbapenemase-producing Enterobacterales are a public health threat worldwide and OXA-48 is the most prevalent carbapenemase in Germany and western Europe. However, the molecular epidemiology of OXA-48 in species other than Escherichia coli and Klebsiella pneumoniae remains poorly understood.AimTo analyse the molecular epidemiology of OXA-48 and OXA-48-like carbapenemases in Citrobacter species (spp.) in Germany between 2011 and 2022.MethodsData of 26,822 Enterobacterales isolates sent to the National Reference Centre (NRC) for Gram-negative bacteria were evaluated. Ninety-one Citrobacter isolates from 40 German hospitals harbouring bla OXA-48/OXA-48­like were analysed by whole genome sequencing and conjugation experiments.ResultsThe frequency of OXA-48 in Citrobacter freundii (CF) has increased steadily since 2011 and is now the most prevalent carbapenemase in this species in Germany. Among 91 in-depth analysed Citrobacter spp. isolates, CF (n = 73) and C. koseri (n = 8) were the most common species and OXA-48 was the most common variant (n = 77), followed by OXA-162 (n = 11) and OXA­181 (n = 3). Forty percent of the isolates belonged to only two sequence types (ST19 and ST22), while most other STs were singletons. The plasmids harbouring bla OXA­48 and bla OXA-162 belonged to the plasmid types IncL (n = 85) or IncF (n = 3), and plasmids harbouring bla OXA­181 to IncX3 (n = 3). Three IncL plasmid clusters (57/85 IncL plasmids) were identified, which were highly transferable in contrast to sporadic plasmids.ConclusionIn CF in Germany, OXA-48 is the predominant carbapenemase. Dissemination is likely due to distinct highly transmissible plasmids harbouring bla OXA­48 or bla OXA-48-like and the spread of the high-risk clonal lineages ST19 and ST22.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Automating the Illumina DNA library preparation kit for whole genome sequencing applications on the flowbot ONE liquid handler robot.

Whole-genome sequencing (WGS) is currently making its transition from research tool into routine (clinical) diagnostic practice. The workflow for WGS includes the highly labor-intensive library preparations (LP), one of the most critical steps in the WGS procedure. Here, we describe the automation of the LP on the flowbot ONE robot to minimize the risk of human error and reduce hands-on time (HOT). For this, the robot was equipped, programmed, and optimized to perform the Illumina DNA Prep automatically. Results obtained from 16 LP that were performed both manually and automatically showed comparable library DNA yields (median of 1.5-fold difference), similar assembly quality values, and 100% concordance on the final core genome multilocus sequence typing results. In addition, reproducibility of results was confirmed by re-processing eight of the 16 LPs using the automated workflow. With the automated workflow, the HOT was reduced to 25 min compared to the 125 min needed when performing eight LPs using the manual workflow. The turn-around time was 170 and 200 min for the automated and manual workflow, respectively. In summary, the automated workflow on the flowbot ONE generates consistent results in terms of reliability and reproducibility, while significantly reducing HOT as compared to manual LP.

Biochemical and molecular characterization of Campylobacter fetus isolates from bulls subjected to bovine genital campylobacteriosis diagnosis in Spain.

BACKGROUND: Bovine genital campylobacteriosis (BGC) is caused by Campylobacter fetus subsp. venerealis (Cfv) including its biovar intermedius (Cfvi). This sexually transmitted disease induces early reproductive failure causing considerable economic losses in the cattle industry. Using a collection of well-characterized isolates (n = 13), C. fetus field isolates (n = 64) and saprophytic isolates resembling Campylobacter (n = 75) obtained from smegma samples of breeding bulls, this study evaluated the concordance of the most used phenotypic (H2S production in cysteine medium and 1% glycine tolerance) and molecular (PCR) methods for the diagnosis of BGC and assessed possible cross-reactions in the molecular diagnostic methods. RESULTS: Characterization at the subspecies level (fetus vs. venerealis) of C. fetus isolated from bull preputial samples using phenotypic and molecular (PCR targeting nahE and ISCfe1) methods showed moderate concordance (κ = 0.462; CI: 0.256-0.669). No cross-reactions were observed with other saprophytic microaerophilic species or with other Campylobacter species that can be present in preputial samples. Whole genome sequencing (WGS) of discrepant isolates showed 100% agreement with PCR identification. For the differentiation of Cfv biovars, comparison of the H2S test (at 72 h and 5 days of incubation) and a PCR targeting the L-cysteine transporter genes showed higher concordance when H2S production was assessed after 5 days (72 h; κ = 0.553, 0.329-0.778 CI vs. 5 days; κ = 0.881, 0.631-1 CI), evidencing the efficacy of a longer incubation time. CONCLUSIONS: This study confirmed the limitations of biochemical tests to correctly identify C. fetus subspecies and biovars. However, in the case of biovars, when extended incubation times for the H2S test (5 days) were used, phenotypic identification results were significantly improved, although PCR-based methods produced more accurate results. Perfect agreement of WGS with the PCR results and absence of cross-reactions with non-C. fetus saprophytic bacteria from the smegma demonstrated the usefulness of these methods. Nevertheless, the identification of new C. fetus subspecies-specific genes would help to improve BGC diagnosis.

Whole genome sequencing analysis of Komagataeibacter nataicola reveals its potential in food waste valorisation for cellulose production.

BACKGROUND: Komagataeibacter nataicola (K. nataicola) is a gram-negative acetic acid bacterium that produces natural bacterial cellulose (BC) as a fermentation product under acidic conditions. The goal of this work was to study the complete genome of K. nataicola and gain insight into the functional genes in K. nataicola that are responsible for BC synthesis in acidic environments. METHODS AND RESULT: The pure culture of K. nataicola was obtained from yeast-glucose-calcium carbonate (YGC) agar, followed by genomic DNA extraction, and subjected to whole genome sequencing on a Nanopore flongle flow cell. The genome of K. nataicola consists of a 3,767,936 bp chromosome with six contigs and 4,557 protein coding sequences. The maximum likelihood phylogenetic tree and average nucleotide identity analysis confirmed that the bacterial isolate was K. nataicola. The gene annotation via RAST server discovered the presence of cellulose synthase, along with three genes associated with lactate utilization and eight genes involved in lactate fermentation that could potentially contribute to the increase in acid concentration during BC synthesis. CONCLUSION: A more comprehensive genome study of K. nataicola may shed light into biological pathway in BC productivity as well as benefit the analysis of metabolites generated and understanding of biological and chemical interactions in BC production later.

Impact of Whole Genome Sequencing to investigate transmission of Serratia marcescens in Neonatal Intensive Care Unit.

Newborns admitted to neonatal intensive care units (NICU) are at increased risk of health care-associated infections. Serratia marcescens represent the third most common pathogen in NICU outbreaks. Here we present an outbreak investigation performed using Whole Genome Sequencing (WGS) analyses and the control measures implemented to limit the spread of S. marcescens in the NICU of an Italian hospital. In February 2023 S. marcescens was isolated from six newborns, when in 2022 this pathogen was isolated only from two samples in the same ward. Measures for infection prevention were adopted. Routinary surveillance screening, performed with rectal swabs collected at admission and weekly thereafter, was implemented to search for S. marcescens presence. Environmental samples were collected. All the isolates, obtained from the conjunctival swab of six newborns, from rectal swab of two newborns who did not develop infections, as well as from the aerators of two faucets, were sequenced. WGS analyses showed no correlation between the isolates from newborns and environmental isolates. The implementation of the measures for infection prevention and control had enabled us to successfully control the outbreak within a short period. WGS analyses proved to be crucial in outbreak investigation to limit the spreading of the pathogens.

Understanding genetic variability: exploring large-scale copy number variants through non-invasive prenatal testing in European populations.

Large-scale copy number variants (CNVs) are structural alterations in the genome that involve the duplication or deletion of DNA segments, contributing to genetic diversity and playing a crucial role in the evolution and development of various diseases and disorders, as they can lead to the dosage imbalance of one or more genes. Massively parallel sequencing (MPS) has revolutionized the field of genetic analysis and contributed significantly to routine clinical diagnosis and screening. It offers a precise method for detecting CNVs with exceptional accuracy. In this context, a non-invasive prenatal test (NIPT) based on the sequencing of cell-free DNA (cfDNA) from pregnant women's plasma using a low-coverage whole genome MPS (WGS) approach represents a valuable source for population studies. Here, we analyzed genomic data of 12,732 pregnant women from the Slovak (9,230), Czech (1,583), and Hungarian (1,919) populations. We identified 5,062 CNVs ranging from 200 kbp and described their basic characteristics and differences between the subject populations. Our results suggest that re-analysis of sequencing data from routine WGS assays has the potential to obtain large-scale CNV population frequencies, which are not well known and may provide valuable information to support the classification and interpretation of this type of genetic variation. Furthermore, this could contribute to expanding knowledge about the central European genome without investing in additional laboratory work, as NIPTs are a relatively widely used screening method.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Whole genome sequencing of recombinant viruses obtained from co-infection and superinfection of Vero cells with modified vaccinia virus ankara vectored influenza vaccine and a naturally occurring cowpox virus.

Modified vaccinia virus Ankara (MVA) has been widely tested in clinical trials as recombinant vector vaccine against infectious diseases and cancers in humans and animals. However, one biosafety concern about the use of MVA vectored vaccine is the potential for MVA to recombine with naturally occurring orthopoxviruses in cells and hosts in which it multiplies poorly and, therefore, producing viruses with mosaic genomes with altered genetic and phenotypic properties. We previously conducted co-infection and superinfection experiments with MVA vectored influenza vaccine (MVA-HANP) and a feline Cowpox virus (CPXV-No-F1) in Vero cells (that were semi-permissive to MVA infection) and showed that recombination occurred in both co-infected and superinfected cells. In this study, we selected the putative recombinant viruses and performed genomic characterization of these viruses. Some putative recombinant viruses displayed plaque morphology distinct of that of the parental viruses. Our analysis demonstrated that they had mosaic genomes of different lengths. The recombinant viruses, with a genome more similar to MVA-HANP (>50%), rescued deleted and/or fragmented genes in MVA and gained new host ranges genes. Our analysis also revealed that some MVA-HANP contained a partially deleted transgene expression cassette and one recombinant virus contained part of the transgene expression cassette similar to that incomplete MVA-HANP. The recombination in co-infected and superinfected Vero cells resulted in recombinant viruses with unpredictable biological and genetic properties as well as recovery of delete/fragmented genes in MVA and transfer of the transgene into replication competent CPXV. These results are relevant to hazard characterization and risk assessment of MVA vectored biologicals.

Klebsiella pneumoniae exhibiting a phenotypic hyper-splitting phenomenon including the formation of small colony variants.

In this study, we characterized a Klebsiella pneumoniae strain in a patient with shrapnel hip injury, which resulted in multiple phenotypic changes, including the formation of a small colony variant (SCV) phenotype. Although already described since the 1960s, there is little knowledge about SCV phenotypes in Enterobacteriaceae. The formation of SCVs has been recognized as a bacterial strategy to evade host immune responses and compromise the efficacy of antimicrobial therapies, leading to persistent and recurrent courses of infections. In this case, 14 isolates with different resisto- and morpho-types were distinguished from the patient's urine and tissue samples. Whole genome sequencing revealed that all isolates were clonally identical belonging to the K. pneumoniae high-risk sequence type 147. Subculturing the SCV colonies consistently resulted in the reappearance of the initial SCV phenotype and three stable normal-sized phenotypes with distinct morphological characteristics. Additionally, an increase in resistance was observed over time in isolates that shared the same colony appearance. Our findings highlight the complexity of bacterial behavior by revealing a case of phenotypic "hyper-splitting" in a K. pneumoniae SCV and its potential clinical significance.

Iron-related gene mutations driving global Mycobacterium tuberculosis transmission revealed by whole-genome sequencing.

BACKGROUND: Iron plays a crucial role in the growth of Mycobacterium tuberculosis (M. tuberculosis). However, the precise regulatory mechanism governing this system requires further elucidation. Additionally, limited studies have examined the impact of gene mutations related to iron on the transmission of M. tuberculosis globally. This research aims to investigate the correlation between mutations in iron-related genes and the worldwide transmission of M. tuberculosis. RESULTS: A total of 13,532 isolates of M. tuberculosis were included in this study. Among them, 6,104 (45.11%) were identified as genomic clustered isolates, while 8,395 (62.04%) were classified as genomic clade isolates. Our results showed that a total of 12 single nucleotide polymorphisms (SNPs) showed a positive correlation with clustering, such as Rv1469 (ctpD, C758T), Rv3703c (etgB, G1122T), and Rv3743c (ctpJ, G676C). Additionally, seven SNPs, including Rv0104 (T167G, T478G), Rv0211 (pckA, A302C), Rv0283 (eccB3, C423T), Rv1436 (gap, G654T), ctpD C758T, and etgB C578A, demonstrated a positive correlation with transmission clades across different countries. Notably, our findings highlighted the positive association of Rv0104 T167G, pckA A302C, eccB3 C423T, ctpD C758T, and etgB C578A with transmission clades across diverse regions. Furthermore, our analysis identified 78 SNPs that exhibited significant associations with clade size. CONCLUSIONS: Our study reveals the link between iron-related gene SNPs and M. tuberculosis transmission, offering insights into crucial factors influencing the pathogenicity of the disease. This research holds promise for targeted strategies in prevention and treatment, advancing research and interventions in this field.

A Practical Guide to Whole Genome Sequencing in the NICU.

The neonatal period is a peak time for the presentation of genetic disorders that can be diagnosed using whole genome sequencing (WGS). While any one genetic disorder is individually rare, they collectively contribute to significant morbidity, mortality, and health-care costs. As the cost of WGS continues to decline and becomes increasingly available, the ordering of rapid WGS for NICU patients with signs or symptoms of an underlying genetic condition is now feasible. However, many neonatal clinicians are not comfortable with the testing, and unfortunately, there is a dearth of geneticists to facilitate testing for every patient that needs it. Here, we will review the science behind WGS, diagnostic capabilities, limitations of testing, time to consider testing, test initiation, interpretation of results, developing a plan of care that incorporates genomic information, and returning WGS results to families.

Whole Genome Sequencing Highlights the Pathogenic Profile in Nocardia Keratitis.

Purpose: Nocardia keratitis is a serious and sight-threatening condition. This study aims to reveal the virulence and antimicrobial resistance gene profile of Nocardia strains using whole genome sequencing. Methods: Whole-genome sequencing was performed on 23 cornea-derived Nocardia strains. Together with genomic data from the respiratory tract and the environment, 141 genomes were then utilized for phylogenetic and pan-genome analyses, followed by virulence and antibiotic resistance analysis. The correlations between virulence genes and pathogenicity were experimentally validated, including the characteristics of Nocardia colonies and clinical and histopathological evaluations of Nocardia keratitis mice models. Results: Whole-genome sequencing of 141 Nocardia strains revealed a mean of 220 virulence genes contributed to bacterial pathogenesis. The mce gene family analysis led to the categorization of strains from the cornea into groups A, B, and C. The colonies of group C had the largest diameter, height, and fastest growth rate. The size of corneal ulcers and the clinical scores showed a significant increase in mouse models induced by group C. The relative expression levels of pro-inflammatory cytokines (CD4, IFN-γ, IL-6Rα, and TNF-α) in the lesion area exhibited an increasing trend from group A to group C. Antibiotic resistance genes (ARGs) spanned nine distinct drug classes, four resistance mechanisms, and seven primary antimicrobial resistance gene families. Conclusions: Whole genome sequencing highlights the pathogenic role of mce gene family in Nocardia keratitis. Its distribution pattern may contribute to the distinct characteristics of the growth of Nocardia colonies and the clinical severity of the mice models.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Discerning the dissemination mechanisms of antibiotic resistance genes through whole genome sequencing of extended-spectrum beta-lactamase (ESBL)-producing E. coli isolated from veterinary clinics and farms in South Korea.

Extended-spectrum beta-lactamase (ESBL)-producing bacteria are resistant to most beta-lactams, including third-generation cephalosporins, limiting the treatment methods against the infections they cause. In this study, we performed whole genome sequencing of ESBL-producing E. coli to determine the mechanisms underlying the dissemination of antibiotic resistance genes. We analyzed 141 ESBL-producing isolates which had been collected from 16 veterinary clinics and 16 farms in South Korea. Long- and short-read sequencing platforms were used to obtain high-quality assemblies. The results showed that blaCTX-M is the dominant ESBL gene type found in South Korea. The spread of blaCTX-M appears to have been facilitated by both clonal spread between different host species and conjugation. Most blaCTX-M genes were found associated with diverse mobile genetic elements that may contribute to the chromosomal integration of the genes. Diverse incompatibility groups of blaCTX-M-harboring plasmids were also observed, which allows their spread among a variety of bacteria. Comprehensive whole genome sequence analysis was useful for the identification of the most prevalent types of ESBL genes and their dissemination mechanisms. The results of this study suggest that the propagation of ESBL genes can occur through clonal spread and plasmid-mediated dissemination, and that suitable action plans should be developed to prevent further propagation of these genes.

Contrasting somatic mutation patterns in aging human neurons and oligodendrocytes.

Characterizing somatic mutations in the brain is important for disentangling the complex mechanisms of aging, yet little is known about mutational patterns in different brain cell types. Here, we performed whole-genome sequencing (WGS) of 86 single oligodendrocytes, 20 mixed glia, and 56 single neurons from neurotypical individuals spanning 0.4-104 years of age and identified >92,000 somatic single-nucleotide variants (sSNVs) and small insertions/deletions (indels). Although both cell types accumulate somatic mutations linearly with age, oligodendrocytes accumulated sSNVs 81% faster than neurons and indels 28% slower than neurons. Correlation of mutations with single-nucleus RNA profiles and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed across the genome similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These stark differences suggest an assortment of active mutagenic processes in oligodendrocytes and neurons.

Opportunistic Elizabethkingia miricola Infections in Intensive Care Unit, Spain.

In 2021, we identified a cluster of Elizabethkingia miricola cases in an intensive care unit in Spain. Because E. miricola is not considered a special surveillance agent in Spain, whole-genome sequencing was not performed. The bacterial source was not identified. All Elizabethkingia species should be listed as special surveillance bacteria.

Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates.

Staphylococcus schleiferi and Staphylococcus coagulans, closely related bacterial species within the Staphylococcus genus, present a challenge in classification and diagnosis due to their close genetic proximity and overlapping phenotypic features. Moreover, our understanding of the virulence mechanisms in staphylococcal species, beyond the extensively studied Staphylococcus aureus, remains limited, underscoring the importance of using comparative data to enhance our insights into virulence within these bacterial species. This study employed a comprehensive approach, utilizing comparative genomics, to identify genomic distinctions between S. schleiferi and S. coagulans, aiming to address the challenges in the accurate classification and diagnosis of these organisms and identify unique features. Whole genome sequencing was performed on six clinical isolates, and their genomes were compared to identify variations in gene content and virulence factors. De novo assembly and annotation revealed two samples as S. coagulans and four samples as S. schleiferi. Analysis of the core genomes revealed conserved regions crucial for defining species identity, while accessory genomic elements contained unique genes, possibly impacting the pathogenicity of the species.

Population Structure and Selection Signal Analysis of Nanyang Cattle Based on Whole-Genome Sequencing Data.

With a rich breeding history, Nanyang cattle (NY cattle) have undergone extensive natural and artificial selection, resulting in distinctive traits such as high fertility, excellent meat quality, and disease resistance. This makes them an ideal model for studying the mechanisms of environmental adaptability. To assess the population structure and genetic diversity of NY cattle, we performed whole-genome resequencing on 30 individuals. These data were then compared with published whole-genome resequencing data from 432 cattle globally. The results indicate that the genetic structure of NY cattle is significantly different from European commercial breeds and is more similar to North-Central Chinese breeds. Furthermore, among all breeds, NY cattle exhibit the highest genetic diversity and the lowest population inbreeding levels. A genome-wide selection signal analysis of NY cattle and European commercial breeds using Fst, θπ-ratio, and θπ methods revealed significant selection signals in genes associated with reproductive performance and immunity. Our functional annotation analysis suggests that these genes may be responsible for reproduction (MAP2K2, PGR, and GSE1), immune response (NCOA2, HSF1, and PAX5), and olfaction (TAS1R3). We provide a comprehensive overview of sequence variations in the NY cattle genome, revealing insights into the population structure and genetic diversity of NY cattle. Additionally, we identify candidate genes associated with important economic traits, offering valuable references for future conservation and breeding efforts of NY cattle.