Search

We here present CLUES2, a full-likelihood method to infer natural selection from sequence data that is an extension of the method CLUES. We make several substantial improvements to the CLUES method that greatly increases both its applicability and its speed. We add the ability to use ancestral recombination graphs on ancient data as emissions to the underlying hidden Markov model, which enables CLUES2 to use both temporal and linkage information to make estimates of selection coefficients. We also fully implement the ability to estimate distinct selection coefficients in different epochs, which allows for the analysis of changes in selective pressures through time, as well as selection with dominance. In addition, we greatly increase the computational efficiency of CLUES2 over CLUES using several approximations to the forward-backward algorithms and develop a new way to reconstruct historic allele frequencies by integrating over the uncertainty in the estimation of the selection coefficients. We illustrate the accuracy of CLUES2 through extensive simulations and validate the importance sampling framework for integrating over the uncertainty in the inference of gene trees. We also show that CLUES2 is well-calibrated by showing that under the null hypothesis, the distribution of log-likelihood ratios follows a χ2 distribution with the appropriate degrees of freedom. We run CLUES2 on a set of recently published ancient human data from Western Eurasia and test for evidence of changing selection coefficients through time. We find significant evidence of changing selective pressures in several genes correlated with the introduction of agriculture to Europe and the ensuing dietary and demographic shifts of that time. In particular, our analysis supports previous hypotheses of strong selection on lactase persistence during periods of ancient famines and attenuated selection in more modern periods.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

In various organisms, sequencing of selectively bred lines at apparent selection limits has demonstrated that genetic variation can remain at many loci, implying that evolution at the genetic level may continue even if the population mean phenotype remains constant. We compared selection signatures at generations 22 and 61 of the High Runner mouse experiment, which includes 4 replicate lines bred for voluntary wheel-running behavior (HR) and 4 non-selected control (C) lines. Previously, we reported multiple regions of differentiation between the HR and C lines, based on whole-genome sequence data for 10 mice from each line at generation 61, which was >31 generations after selection limits had been reached in all HR lines. Here, we analyzed pooled sequencing data from ~20 mice for each of the 8 lines at generation 22, around when HR lines were reaching limits. Differentiation analyses of allele frequencies at ~4.4 million SNP loci used the regularized T-test and detected 258 differentiated regions with FDR = 0.01. Comparable analyses involving pooling generation 61 individual mouse genotypes into allele frequencies by line produced only 11 such regions, with almost no overlap among the largest and most statistically significant peaks between the two generations. These results implicate a sort of genetic churn that continues at loci relevant for running. Simulations indicate that loss of statistical power due to random genetic drift and sampling error are insufficient to explain the differences in selection signatures. The 13 differentiated regions at generation 22 with strict culling measures include 79 genes related to a wide variety of functions. Gene ontology identified pathways related to olfaction and vomeronasal pathways as being overrepresented, consistent with generation 61 analyses, despite those specific regions differing between generations. Genes Dspp and Rbm24 are also identified as potentially explaining known bone and skeletal muscle differences, respectively, between the linetypes.

Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

Rampant C-to-U RNA editing drives the mutation and evolution of SARS-CoV-2. While much attention has been paid to missense mutations, the C-to-U events leading to AUG and thus creating novel ORFs were uninvestigated. By utilizing the public time-course mutation data from the worldwide SARS-CoV-2 population, we systematically identified the "AUG-gain mutations" caused by C-to-U RNA editing. Synonymous mutations were of special focus. A total of 58 synonymous C-to-U sites are able to create out-of-frame AUG in coding sequence (CDS). These 58 synonymous sites showed significantly higher allele frequency (AF) and increasing rate (dAF/dt) than other C-to-U synonymous sites in the SARS-CoV-2 population, suggesting that these 58 AUG-gain events conferred additional benefits to the virus and are subjected to positive selection. The 58 predicted new ORFs created by AUG-gain events showed the following advantages compared to random expectation: they have longer lengths, higher codon adaptation index (CAI), higher Kozak scores, and higher tRNA adaptation index (tAI). The 58 putatively novel ORFs have high expressibility and are very likely to be functional, providing an explanation for the positive selection on the 58 AUG-gain mutations. Our study proposed a possible mechanism of the emergence of de novo genes in SARS-CoV-2. This idea should be helpful in studying the mutation and evolution of SARS-CoV-2. [ABSTRACT FROM AUTHOR]

Background: PLA2G6-associated neurodegeneration (PLAN) comprises three diseases with overlapping features: infantile neuroaxonal dystrophy (INAD), atypical neuroaxonal dystrophy (atypical NAD), and PLA2G6-related dystonia-parkinsonism. INAD is an early onset disease characterized by progressive loss of vision, muscular control, and mental skills. The prevalence of PLA2G6-associated diseases has not been previously calculated. Methods: To provide the most accurate prevalence estimate, we utilized two independent approaches: database-based approach which included collecting variants from ClinVar, Human Gene Mutation Database (HGMD) and high confidence predicted loss-of-function (pLoF) from gnomAD (Rare Genomes Project Genetic Prevalence Estimator; GeniE), and literature-based approach which gathered variants through Mastermind Genomic Search Engine (Genomenon, Inc). Genetic prevalence of PLAN was calculated based on allele frequencies from gnomAD, assuming Hardy–Weinberg equilibrium. Results: In the PLA2G6 gene, our analysis found 122 pathogenic, 82 VUS, and 15 variants with conflicting interpretations (pathogenic vs VUS) between two approaches. Allele frequency was available for 58 pathogenic, 42 VUS, and 15 conflicting variants in gnomAD database. If pathogenic and/or conflicting variants are included, the overall genetic prevalence was estimated to be between 1 in 987,267 to 1 in 1,570,079 pregnancies, with the highest genetic prevalence in African/African-American (1 in 421,960 to 1 in 365,197) and East-Asian (1 in 683,978 to 1 in 190,771) populations. Conclusion: Our estimates highlight the significant underdiagnosis of PLA2G6-associated neurodegeneration and underscores the need for increased awareness and diagnostic efforts. Furthermore, our study revealed a higher carrier frequency of PLA2G6 variants in African and Asian populations, stressing the importance of expanded genetic sequencing in non-European populations to ensure accurate and comprehensive diagnosis. Future research should focus on confirming our findings and implementing expanded sequencing strategies to facilitate maximal and accurate diagnosis, particularly in non-European populations. [ABSTRACT FROM AUTHOR]

High heterogeneity in clinical benefit characterizes the use of immune checkpoint inhibitors (ICIs) in non-small cell lung cancer (NSCLC). We prospectively enrolled 113 advanced NSCLC patients treated with ICIs and performed liquid biopsy at the time of ICI start (T1), after 3 weeks (T2) and at the time of radiological evaluation (T3). Molecular variables were associated with outcome endpoints: cfDNA quantification, its dynamic change (∆T2–T1), variant allele frequency (VAF) of the gene with the highest frequency detected at baseline with NGS (maxVAF) and its dynamic change (∆T2–T1). At multivariate analysis, PD-L1 negativity, T1 cfDNA, cfDNA increase (∆T2–T1), and T2 VAF were significantly associated with shorter progression-free survival (PFS); PD-L1 negativity, squamous histology, T1 cfDNA, cfDNA (∆T2–T1) increase, and T2 maxVAF affected overall survival (OS). Among high PD-L1 expressing patients treated in first-line, elevated T2 maxVAF and cfDNA increase (∆T2–T1) correlated with worse PFS; higher T2 maxVAF and cfDNA increase (∆T2–T1) with worse OS. Derived integrated models were used to build nomograms and categorize different risk groups. [ABSTRACT FROM AUTHOR]

Introduction: In collaboration with the Minnesota Hmong community, we have previously discovered significant differences in allele frequencies for key Single Nucleotide Variations (SNVs) within Very Important Pharmacogenes (VIPs) between Hmong and East Asians. Recognizing the potential clinical implications of these observed differences, we sought to validate these observations in a Hmong cohort residing in California, the state with the largest Hmong population in the US. Robust validation of these differences would affect motivation for clinicians treating individuals who identify as Hmong to consider pharmacogenomic (PGx) testing as a means to improve clinical decision making when using therapeutic agents in this unique population. Method: Guided by California Hmong community leaders and utilizing the basic approach of community-based participatory research, demographic, clinical information and a buccal swab was obtained from Hmong adults residing in California. A commercial PGx testing panel was performed on these samples and specific allele frequencies of interest were compared between California and Minnesota Hmong. Allele frequency differences between California Hmong, East Asians and Europeans, were also compared. Return-of-PGx-results and presentations of group data were made to members of the Hmong along with PGx educational sessions to help interpret the observations. Results: In 118 California Hmong who completed the study, the allele frequencies for SNV's were similar to previous Minnesota Hmong results. Furthermore, out of the 18 SNVs that were not previously reported in Hmong, allele frequencies were statistically different in 38% (7/18) of SNVs comparing California Hmong to East Asians, and in 77.8% (14/18) SNVs comparing California Hmong to Europeans. Conclusion: These results validate the original study's findings that Hmong people living in different US locations have similar allele frequencies for key PGx genes. Further, for many of these PGx genes, their allele frequencies are significantly different compared to either East Asians or Europeans. Clinicians should consider these important differences when prescribing medications for people who identify as Hmong. [ABSTRACT FROM AUTHOR]

Identifying genetic drivers of chronic diseases is necessary for drug discovery. Here, we develop a machine learning-assisted genetic priority score, which we call ML-GPS, that incorporates genetic associations with predicted disease phenotypes to enhance target discovery. First, we construct gradient boosting models to predict 112 chronic disease phecodes in the UK Biobank and analyze associations of predicted and observed phenotypes with common, rare, and ultra-rare variants to model the allelic series. We integrate these associations with existing evidence using gradient boosting with continuous feature encoding to construct ML-GPS, training it to predict drug indications in Open Targets and externally testing it in SIDER. We then generate ML-GPS predictions for 2,362,636 gene-phecode pairs. We find that the use of predicted phenotypes, which identify substantially more genetic associations than observed phenotypes across the allele frequency spectrum, significantly improves the performance of ML-GPS. ML-GPS increases coverage of drug targets, with the top 1% of all scores providing support for 15,077 gene-phecode pairs that previously had no support. ML-GPS can also identify well-known target-disease relationships, promising targets without indicated drugs, and targets for several drugs in clinical trials, including LRRK2 inhibitors for Parkinson's disease and olpasiran for cardiovascular disease. Here, the authors introduce ML-GPS, a machine learning framework that prioritizes drug targets for 112 chronic diseases and integrates genetic associations with predicted phenotypes. [ABSTRACT FROM AUTHOR]

Background: The Csangos are an East-Central European ethnographic group living mainly in east of Transylvania in Romania. Traditionally, ethnography distinguishes three Csango subpopulations, the Moldavian, Gyimes and Burzenland Csangos. In our previous study we found that the Moldavian Csangos have East Asian/Siberian Turkic ancestry components that might be unique in the East-Central European region and might help to better understand the history of Hungarian speaking ethnic groups of the area. Since then, we obtained further Csango samples from Moldavia and from a distinct region of Gyimes, which two Csango subgroups are traditionally different since they live in a degree of isolation not only from other people but also from each other. Here we present the first genomic analysis of Gyimes Csangos, which intended to compare the genetic makeup of these two Csango subgroups using both allele-frequency and haplotype-based methods. The main goal of the study was to investigate the genetic isolation of the Csangos on a genome-wide SNP basis and to assess the isolation of Gyimes Csangos, which in contrast to the Moldavians was not yet studied. Results: Our results show that these two Csango groups show slight differences from each other. We confirmed the genetic isolation of Moldavian Csangos and revealed that Gyimes Csangos have a similar, but detectably weaker isolation. In the case of Gyimes Csangos we detected also a stronger East European or presumably Asian derived ancestry. Conclusion: The Gyimes Csangos show a degree of genetic isolation comparable to that of the Moldavians. The Asian ancestry that differentiates the Moldavian Csango people from the other East-Central European populations may be present in the Gyimes Csangos in an even higher degree, since Gyimes Csango individuals show a more significant share from that ancestry component. [ABSTRACT FROM AUTHOR]

Background: Thioester-containing protein 1 (TEP1) is a crucial component of mosquitoes' natural resistance to parasites. To effectively combat malaria, there is a need to better understand how TEP1 polymorphism affects phenotypic traits during infections. Therefore, the purpose of this study was to determine the Tep1 genotype frequency in malaria vector populations from south-western Ethiopia and investigate its effect on Plasmodium oocyst development in Anopheles arabiensis populations. Methods: Using standard dippers, Anopheles mosquito larvae were collected from aquatic habitats in Asendabo, Arjo Dedessa, and Gambella in 2019 and 2020. Collected larvae were reared to adults and identified morphologically. Female An. gambiae s.l. were allowed to feed on infected blood containing the same number of gametocytes obtained from P. falciparum and P. vivax gametocyte-positive individuals using indirect membrane feeding methods. Polymerase Chain Reaction (PCR) was used to identify An. gambiae s.l. sibling species. Three hundred thirty An. gambiae s.l. were genotyped using Restricted Fragment Length Polymorphism (RFLP) PCR and sub samples were sequenced to validate the TEP1 genotyping. Results: Among the 330 samples genotyped, two TEP1 alleles, TEP1*S1 (82% frequency) and TEP1*R1 (18% frequency), were identified. Three equivalent genotypes, TEP1*S1/S1, TEP1*R1/R1, and TEP1*S1/R1, had mean frequencies of 65.15%, 2.12%, and 32.73%, respectively. The nucleotide diversity was ranging from 0.36554 to 0. 46751 while haplotype diversity ranged from 0.48871 to 0.63161, across all loci. All sample sites had positive Tajima's D and Fu's Fs values. There was a significant difference in the TEP1 allele frequency and genotype frequency among mosquito populations (p < 0.05), except populations of Anopheles arabiensis from Asendabo and Gambella (p > 0.05). In addition, mosquitoes with the TEP1 *RR genotype were susceptible and produced fewer Plasmodium oocysts than mosquitoes with the TEP1 *SR and TEP1 *SS genotypes. Conclusion: The alleles identified in populations of An. arabiensis were TEP1*R1 and TEP1*S1. There was no significant variation in TEP1*R1 allele frequency between the high and low transmission areas. Furthermore, An. arabiensis carrying the TEP1*R1 allele was susceptible to Plasmodium infection. Further studies on vector-parasite interactions, particularly on the TEP1 gene, are required for vector control techniques. [ABSTRACT FROM AUTHOR]

Simple Summary: This study focused on 50 Hu sheep from Gansu to analyze the genetic diversity and kinship relationships within the population and to construct family pedigrees. Candidate genes related to growth development and high reproductive performance in Hu sheep were identified. These findings provide a reference for future Hu sheep breeding, genetic evaluation, and population utilization. This research is designed to examine the genetic diversity and kinship among Hu sheep, as well as to discover genes associated with crucial economic traits. A selection of 50 unrelated adult male Hu sheep underwent genotyping with the SNP50K BeadChip. Seven indicators of genetic diversity were assessed based on high-quality SNP data: effective population size (Ne), polymorphic information content (PIC), polymorphic marker ratio (PN), expected heterozygosity (He), observed heterozygosity (Ho), effective number of alleles, and minor allele frequency (MAF). Plink software was employed to compute the IBS genetic distance matrix and detect runs of homozygosity (ROHs), while the G matrix and principal component analysis were performed using GCTA software. Selective sweep analysis was carried out using ROH, Pi, and Tajima's D methodologies. This study identified a total of 64,734 SNPs, of which 56,522 SNPs remained for downstream analysis after quality control. The population displayed relatively high genetic diversity. The 50 Hu sheep were ultimately grouped into 12 distinct families, with families 6, 8, and 10 having the highest numbers of individuals, each consisting of 6 sheep. Furthermore, a total of 294 ROHs were detected, with the majority having lengths between 1 and 5 Mb, and the inbreeding coefficient FROH was 0.01. In addition, 41, 440, and 994 candidate genes were identified by ROH, Pi, and Tajima's D methods, respectively, with 3 genes overlapping (BMPR1B, KCNIP4, and FAM13A). These results offer valuable insights for future Hu sheep breeding, genetic assessment, and population management. [ABSTRACT FROM AUTHOR]

Caenorhabditis elegans is one of the most popular model organisms used to genetically dissect complex biological phenomena. One common technique used routinely in the C. elegans laboratory is the generation of strains carrying combinations of genetic mutations via classical genetic crosses. Here, we have developed a simple and convenient application to quickly identify useful genetic markers (phenotypical and fluorescent) and their chromosomal positions to aid in the development of genetic cross strategies. The user-friendly software identifies and prioritizes markers with the least genetic distance to a gene of interest, as well as displays the strain name, ease of scoring, nature of the marker (fluorescent transgene or phenotypic information), mating efficiency, and number of available alleles. In addition, recombination frequencies between the gene of interest and each genetic marker are calculated automatically. The application, called "SoMarker," is designed for both MacOS and Windows environments and is available to freely download and modify through open-source software. [ABSTRACT FROM AUTHOR]

Single nucleotide polymorphisms (SNPs) of the IL-16 gene have been reported to influence the risk of several cancers, but their role in ovarian cancer (OC) has not been studied. Using the restriction fragment length polymorphism (PCR-RFLP) method, we examined four IL-16 SNPs: rs11556218 (T > G), rs4778889 (T > C), rs4072111 (C > T), and rs1131445 (T > C) in blood samples from 413 women of Central European descent, including 200 OC patients and 213 healthy controls. Among the patients, 62% were postmenopausal, 84.5% were diagnosed in late stages (FIGO IIb-IV), and 73.5% had high-grade serous OC (HGSOC). Minor allele frequencies in controls were 9.2% for rs11556218 (G allele), 13.7% for rs4778889 (C allele), 10.4% for rs4072111 (T allele), and 32.3% for rs1131445 (C allele). We found significant associations of rs11556218 (G vs. T allele: OR 2.76, 95% CI 1.84–4.14, p < 0.0001) with elevated OC risk in the whole cohort (p < 0.001) and in both premenopausal (p < 0.001) and postmenopausal (p = 0.001) subgroups. These associations remained significant across heterozygote (p < 0.001), dominant (p < 0.001), and overdominant (p < 0.001) models. IL-16 rs4778889 was associated with OC risk predominantly in premenopausal women (p < 0.0001 in almost all models). In the whole cohort, the C allele was associated with OC risk (OR 1.54, CI 95% 1.06–2.23, p = 0.024), and the association of rs4778889 was significant in dominant (p = 0.019), overdominant (p = 0.033), and heterozygote (p = 0.027) models. Furthermore, rs4778889 was linked with HGSOC (p = 0.036) and endometriosis-related OC subtypes (p = 0.002). No significant associations were found for rs4072111 or rs1131445 (p = 0.81 or 0.47, respectively). In conclusion, rs11556218 and rs4778889 SNPs are associated with OC risk, especially in premenopausal women. [ABSTRACT FROM AUTHOR]

Background: The clustered regulatory interspaced short palindromic repeats (CRISPR)-Cas13a system has strong potential for highly sensitive detection of exogenous sequences. The detection of KRASG12 point mutations with low allele frequencies may prove powerful for the formal diagnosis of pancreatic ductal adenocarcinoma (PDAC). Results: We implemented preamplification of KRAS alleles (wild-type and mutant) to reveal the presence of mutant KRAS with CRISPR-Cas13a. The discrimination of KRASG12D from KRASWT was poor for the generic KRAS preamplification templates and depended on the crRNA design, the secondary structure of the target templates, and the nature of the mismatches between the guide and the templates. To improve the specificity, we used an allele-specific PCR preamplification method called CASPER (Cas13a Allele-Specific PCR Enzyme Recognition). CASPER enabled specific and sensitive detection of KRASG12D with low DNA input. CASPER detected KRAS mutations in DNA extracted from patients' pancreatic ultrasound-guided fine-needle aspiration fluid. Conclusion: CASPER is easy to implement and is a versatile and reliable method that is virtually adaptable to any point mutation. [ABSTRACT FROM AUTHOR]

Objective: To develop a high‐throughput sequencing panel for the diagnosis of developmental and epileptic encephalopathy in Tunisia and to clarify the frequency of disease‐causing genes in this region. Methods: We developed a custom panel for next‐generation sequencing of the coding sequences of 116 genes in individuals with developmental and epileptic encephalopathy from the Tunisian population. Segregation analyses and in silico studies have been conducted to assess the identified variants' pathogenicity. Results: We report 12 pathogenic variants in SCN1A, CHD2, CDKL5, SZT2, KCNT1, GNAO1, PCDH19, MECP2, GRIN2A, and SYNGAP1 in patients with developmental and epileptic encephalopathy. Five of these variants are novel: "c.149delA, p.(Asn50MetfsTer26)" in CDKL5; "c.3616C > T, p.(Arg1206Ter)" in SZT2; "c.111_113del, p.(Leu39del)" in GNAO1; "c.1435G>C, p.(Asp479His)" in PCDH19; and "c.2143delC, p.(Arg716GlyfsTer10)" in SYNGAP1. Additionally, for four of our patients, the genetic result facilitated the choice of the appropriate treatment. Significance: This is the first report of a custom gene panel to identify genetic variants implicated in developmental and epileptic encephalopathy in the Tunisian population as well as the North African region (Tunisia, Egypt, Libya, Algeria, Morocco) with a diagnostic rate of 30%. This high‐throughput sequencing panel has considerably improved the rate of positive diagnosis of developmental and epileptic encephalopathy in the Tunisian population, which was less than 15% using Sanger sequencing. The benefit of genetic testing in these patients was approved by both physicians and parents. [ABSTRACT FROM AUTHOR]

Stenotrophomonas maltophilia causes challenging infections in immunocompromised patients, exhibiting increasing resistance to multiple antimicrobials and possessing various virulence genes, including emerging resistance to trimethoprim‐sulfamethoxazole. A total of 80 clinical isolates of S. maltophilia were collected from multiple hospitals in Tehran, Iran. This study conducted an analysis of antibiotic susceptibility by disc diffusion method and E‐test assay, resistance and virulence gene frequencies were examined by PCR‐sequencing, and multilocus sequencing typing (MLST) was performed for strain typing. Across the tested isolates, we observed notably high resistance rates for imipenem 80 (100%), meropenem 78(97.5%), and ceftazidime 72 (90%), while trimethoprim‐sulfamethoxazole (SXT) showed a lower resistance rate of 2 (2.5%). Minocycline and levofloxacin demonstrated the highest susceptibility rates, with 70 (87.5%) and 80 (100%), respectively. The prevalence of antibiotic resistance genes blaL1, and blaL2 was 71 (88.75%) and 76 (95%), respectively. Additionally, the PCR analysis revealed that the frequency of virulence genes (fliC, virB, papD, pilU, hlyIII, stmPr1, and stmPr2) was 78 (97.5%), 77 (96.25%), 58 (72.5%), 77 (96.2%), 76 (95%), 31 (38.75%), and 80 (100%), respectively. Resistance to SXT isolate belong to the sequence type (ST15) and exhibits allelic profiles of (10, 29, 21, 21, 32, 32, and 10). The data obtained from our investigation have indicated that SXT remains an efficacious antibiotic and also highlighted the importance of effective management, identification of resistant isolates, and typing methods to address the global prevalence of antibiotic resistance in S. maltophilia. [ABSTRACT FROM AUTHOR]

DNA methyltransferase 3 A mutations (DNMT3AMT) are frequent in myeloid neoplasia (MN) and mostly heterozygous. However, cases with multiple DNMT3AMT can be also encountered but their clinical and genetic landscape remains unexplored. We retrospectively analyzed 533 cases with DNMT3AMT identified out of 5,603 consecutive MNs, of whom 8.4% had multiple DNMT3AMT hits. They were most frequent in acute myeloid leukemia (AML) with R882 variant accounting for 13.3% of the multi-hits. Multiple DNMT3AMT more likely coincided with IDH2 (P = 0.005) and ETV6 (P = 0.044) mutations compared to patients with single DNMT3AMT. When the sum of variant allele frequencies (VAFs) for multiple DNMT3AMT exceeded 60%, we found a significant positive clonal burden correlation of the two DNMT3A variants (P < 0.0001) suggesting that they occurred in biallelic configuration. AML patients with biallelic DNMT3A inactivation (n = 52) presented with older age (P = 0.029), higher leukocytes (P < 0.0001) and peripheral blast counts (P = 0.0001) and significantly poorer survival rate (5.6% vs. 47.6% at 2 years; P = 0.002) than monoallelic DNMT3AMT. Multivariate analysis identified biallelic DNMT3AMT (HR 2.65; P = 0.001), male gender (HR 2.05; P = 0.014) and adverse genetic alteration according to the European LeukemiaNet 2022 classification (HR 1.84; P = 0.028) as independent adverse factors for survival, whereas intensive chemotherapy (HR 0.47; P = 0.011) favorably influenced outcomes. Longitudinal molecular analysis of 12 cases with biallelic DNMT3AMT demonstrated that such clones persisted or expanded in 9 relapsed or transformed cases (75%) suggesting the early origin of biallelic hits with strong leukemogenic potential. Our study describes the likelihood that biallelic DNMT3AMT, while rare, are indeed compatible with clonal expansion and thus questions the applicability of synthetic lethality strategies. [ABSTRACT FROM AUTHOR]

Hybridization is a common technique in the breeding of new varieties. In this study, we obtained F1 and F2 hybrids through interspecific hybridization between Siniperca chuatsi and Siniperca scherzeri and backcross offspring by backcrossing between F1 and Siniperca chuatsi. We analyzed the morphological traits and genetic diversity of the parents and their progeny. The results showed that although most morphological traits of the progeny were not significantly different from those of their parents, the F1 hybrids exhibited the most pronounced morphological variations. The F2 and backcross hybrids appeared morphological differentiation within their populations. The CPL/CPW ratio emerged as a reliable marker for identifying hybrids morphologically. The expected heterozygosity for F1 (0.6264), F2 (0.6677), and backcross offspring (0.5260) were intermediate between Siniperca chuatsi (0.7388) and Siniperca scherzeri (0.3910), indicating a relatively high level of genetic diversity in the offspring. Additionally, the private allele frequency of F1 and F2 was higher than backcross offspring. One SSR marker (G5220) was identified as a specific molecular marker to distinguish the hybrid F1 population from the parents, hybrid F2, and backcross F1 populations by analyzing specific alleles. The result of Ne, He, I, and PIC confirmed that F2 maintained the genetic diversity of F1. We speculated that both F1 and F2 have higher breeding potential. This study provides essential baseline data for genetic breeding of mandarin fish that will be instrumental in the future. [ABSTRACT FROM AUTHOR]

Streptococcus suis is a major swine pathogen and zoonotic agent, causing important economic losses to the porcine industry. Here, we used genomics approaches to characterize 251 S. suis isolates recovered from diseased pigs across Belgium, France, Germany, Hungary, the Netherlands, Spain, and the United Kingdom. We identified 13 serotypes, being serotypes 9 and 2 the most prevalent, and 34 sequence types (STs), including 16 novel STs, although ST16 and ST1 dominated the strain population. Phylogenetic analysis revealed complex genetic relationships, notable geographic clustering, and potential differential capacity for capsular switching among serotype 9 isolates. We found antimicrobial resistance (AMR) genes in 85.3% of the isolates, with high frequencies of genes conferring resistance to tetracyclines and macrolides. Specifically, 49.4% of the isolates harbored the tetO gene, and 64.9% possessed the ermB gene. Additionally, we observed a diverse array of virulence-associated genes (VAGs), including the classical VAGs mrp, epf, and sly, with variable presence across different genotypes. The high genetic diversity among European S. suis isolates highlights the importance of targeted antimicrobial use and flexible vaccine strategies. Rapid strain characterization is crucial for optimizing swine health management, enabling tailored interventions like the development of autovaccines to mitigate S. suis infections. [ABSTRACT FROM AUTHOR]

G‐protein‐coupled receptors (GPCRs) are the largest superfamily in the human genome and the major targets for the market drugs. Recent massive genomics studies revealed numerous natural variations in the general population. 54KJPN is the most extensive Japanese population genomics study, curating the whole genome sequences from about 54,000 individuals. Here, by analyzing 390 non‐olfactory GPCR genes in the 54KJPN dataset, we annotated 25,443 missense single‐nucleotide variations. Among them, we found 120 major variations that appear with an allele frequency greater than 0.5, including variations that occurred on posttranslational modification sites. Structural alignment of GPCRs using the generic numbering system in the GPCRdb reveals enrichment of alterations in the conserved arginine residue within the DRY motif, which contributes to downstream G‐protein signaling. A comparison with the worldwide 1000 Genomes Project (1KGP) dataset found 23 variations that were present exclusively in the 54KJPN dataset. This study will be the basis for future pharmacogenomics studies for the Japanese population. [ABSTRACT FROM AUTHOR]

Simple Summary: Canine degenerative myelopathy is a fatal neurodegenerative disorder that affects the spinal cord. It is a late-onset disease, with the first clinical signs becoming evident later in life at approximately 8 years of age. The aim of this study was to retrospectively evaluate allelic and genotypic frequencies of the c.118G > A and c.52A > T mutations in an Italian canine population of 1667 dogs belonging to 84 breeds. The results showed that 65.47% (n. 893) of the subjects were clear, 25.59% (n. 349) were heterozygous carriers, and 8.94% (n. 122) were homozygous for the risk allele. The highest frequency was observed in Pembroke Welsh Corgis (55.49%). Canine degenerative myelopathy is a fatal neurodegenerative disorder that affects the spinal cord. It is a late-onset disease, with symptoms becoming evident later in life at approximately 8 years of age. The principal aim of this study was to retrospectively evaluate allelic and genotypic frequencies of the c.118G > A and c.52A > T mutations located on the SOD1 gene in an Italian canine population to provide detailed information on the prevalence of the mutations in the country. The genetic data of different breeds were collected through DNA tests over a nine-year period in the Italian canine population. For each dog, the breed, sex, age, and DNA test results were recorded. Allelic and genotypic frequencies were calculated. A total of 1667 DNA tests for the c.118G > A and c.52A > T mutations were carried out on 84 breeds. For the analysis of prevalence, only breeds counting more than 20 subjects have been considered, for a total of 1410 DNA tests obtained from 13 different breeds. In the population tested for the c.118G > A mutation, 65.47% (n. 893) of the subjects were clear, 25.59% (n. 349) were heterozygous carriers, and 8.94% (n. 122) were homozygous for the mutated allele. The mutation showed the highest frequency in Pembroke Welsh Corgis (55.49%) and the lowest frequencies in Poodles (6.32%) and Australian Shepherds (7.14%). The allelic frequency of the c.52A > T mutation was 7.61% in the Bernese Mountain dog. Neither variant differed between females and males in genotypic frequencies. The present study provides insights into the allelic and genotypic frequencies of canine degenerative myelopathy in different dog breeds in Italy. [ABSTRACT FROM AUTHOR]

Summary: Background: Azithromycin (AZ) is a widely used antibiotic. The aim of this study was to characterise the clinical features, outcomes, and HLA association in patients with drug‐induced liver injury (DILI) due to AZ. Methods: The clinical characteristics of individuals with definite, highly likely, or probable AZ‐DILI enrolled in the US Drug‐Induced Liver Injury Network (DILIN) were reviewed. HLA typing was performed using an Illumina MiSeq platform. The allele frequency (AF) of AZ‐DILI cases was compared to population controls, other DILI cases, and other antibiotic‐associated DILI cases. Results: Thirty cases (4 definite, 14 highly likely, 12 probable) of AZ‐DILI were enrolled between 2004 and 2022 with a median age of 46 years, 83% white, and 60% female. Median duration of AZ treatment was 5 days. Latency was 18.5 days. 73% were jaundiced at presentation. The injury pattern was hepatocellular in 60%, cholestatic in 27%, and mixed in 3%. Ten cases (33%) were severe or fatal; 90% of these were hepatocellular. Two patients required liver transplantation. One patient with chronic liver disease died of hepatic failure. Chronic liver injury developed in 17%, of which 80% had hepatocellular injury at onset. HLA‐DQA1*03:01 was significantly more common in AZ‐DILI versus population controls and amoxicillin‐clavulanate DILI cases (AF: 0.29 vs. 0.11, p = 0.001 and 0.002, respectively). Conclusion: Azithromycin therapy can lead to rapid onset of severe hepatic morbidity and mortality in adult and paediatric populations. Hepatocellular injury and younger age were associated with worse outcomes. HLA‐DQA1*03:01 was significantly more common in AZ cases compared to controls. [ABSTRACT FROM AUTHOR]

The genetic factors of stroke in South Asians are largely unexplored. Exome-wide sequencing and association analysis (ExWAS) in 75 K Pakistanis identified NM_000435.3(NOTCH3):c.3691 C > T, encoding the missense amino acid substitution p.Arg1231Cys, enriched in South Asians (alternate allele frequency = 0.58% compared to 0.019% in Western Europeans), and associated with subcortical hemorrhagic stroke [odds ratio (OR) = 3.39, 95% confidence interval (CI) = [2.26, 5.10], p = 3.87 × 10−9), and all strokes (OR [CI] = 2.30 [1.77, 3.01], p = 7.79 × 10−10). NOTCH3 p.Arg231Cys was strongly associated with white matter hyperintensity on MRI in United Kingdom Biobank (UKB) participants (effect [95% CI] in SD units = 1.1 [0.61, 1.5], p = 3.0 × 10−6). The variant is attributable for approximately 2.0% of hemorrhagic strokes and 1.1% of all strokes in South Asians. These findings highlight the value of diversity in genetic studies and have major implications for genomic medicine and therapeutic development in South Asian populations. Stroke is a multifactorial disease influenced by genetic and environmental factors. Here, the authors apply exome-wide association analysis to find rare coding variants associated with stroke in a Pakistani cohort, finding a significant association of a variant in NOTCH3 that is highly enriched in South Asians. [ABSTRACT FROM AUTHOR]

Background: Wilson disease (WD) is an autosomal recessive inherited disease caused by ATP7B variants and characterized by copper metabolism defects. However, children with WD are often asymptomatic, making the clinical diagnosis difficult. Therefore, more accurate methods are required for clinical diagnosis. The objective of this study was to highlight the phenotypic and genetic characteristics of children with WD in northeast China. Methods: We retrospectively analyzed the clinical data and gene sequencing results of 65 children with WD from January 1, 2014, to December 31, 2022, at the Shengjing Hospital of China Medical University. All data refer to the time of diagnosis before treatment. Results: The median age at diagnosis was 5 years (range 1.2–15 years). In 50 cases (50/65, 76.9%) patients, routine physical examinations revealed only abnormal liver function. However, they had a significantly negative (p < 0.05) Kayser–Fleischer ring (KF). Children with acute liver failure had significantly increased 24 h urinary copper excretion (p < 0.05). We detected 46 genetic variants of ATP7B, including seven novel variants. The most frequent variant was p.R778L with an allele frequency of 38.7%. Phenotype-genotype correlation analysis suggested that p.R778L was significantly associated with lower serum ceruloplasmin levels and higher zinc levels (p < 0.05). The loss-of-function (LOF) variant was associated with significantly lower albumin levels (p < 0.05). Conclusion: Most children with WD are asymptomatic, which makes early diagnosis of WD difficult. Therefore, clinical and laboratory characteristics as well as genetic testing are essential. p.R778L is the most frequent variant of ATP7B in China and may play an important role in lowering serum ceruloplasmin levels. [ABSTRACT FROM AUTHOR]

Introduction: Owing to advances in high-throughput genome sequencing, QTLSeq mapping of salt tolerance traits is a major platform for identifying soil-salinity tolerance QTLs to accelerate marker-assisted selection for salt-tolerant rice varieties. We performed QTL-BSA-Seq in the seedling stage of rice from a genetic cross of the extreme salt-sensitive variety, IR29, and "Jao Khao" (JK), a Thai salt-tolerant variety. Methods: A total of 462 F2 progeny grown in soil and treated with 160 mM NaCl were used as the QTL mapping population. Two high- and low-bulk sets, based on cell membrane stability (CMS) and tiller number at the recovery stage (TN), were equally sampled. The genomes of each pool were sequenced, and statistical significance of QTL was calculated using QTLseq and G prime (G') analysis, which is based on calculating the allele frequency differences or D(SNP index). Results: Both methods detected the overlapping interval region, wherein CMSbulk was mapped at two loci in the 38.41-38.85 Mb region with 336 SNPs on chromosome 1 (qCMS1) and the 26.13-26.80 Mb region with 1,011 SNPs on chromosome 3 (qCMS3); the D(SNP index) peaks were -0.2709 and 0.3127, respectively. TN-bulk was mapped at only one locus in the overlapping 38.26-38.95 Mb region on chromosome 1 with 575 SNPs (qTN1) and a D(SNP index) peak of -0.3544. These identified QTLs in two different genetic backgrounds of segregating populations derived from JK were validated. The results confirmed the colocalization of the qCMS1 and qTN1 traits on chromosome 1. Based on the CMS trait, qCMS1/qTN1 stably expressed 6%-18% of the phenotypic variance in the two validation populations, while qCMS1/qTN1 accounted for 16%-20% of the phenotypic variance in one validation population based on the TN trait. Conclusion: The findings confirm that the CMS and TN traits are tightly linked to the long arm of chromosome 1 rather than to chromosome 3. The validated qCMS-TN1 QTL can be used for gene/QTL pyramiding in marker-assisted selection to expedite breeding for salt resistance in rice at the seedling stage. [ABSTRACT FROM AUTHOR]

The regulatory serine protease, complement factor I (FI), in conjunction with one of its cofactors (FH, C4BP, MCP, or CR1), plays an essential role in controlling complement activity through inactivation of C3b and C4b. The functional impact by missense variants in the CFI gene, particularly those with minor allele frequencies of 0.01% to 0.1%, is infrequently studied. As such, these variants are typically classified as variants of uncertain significance (VUS) when they are identified by clinical testing. Herein, we utilized a minigene splicing assay to assess the functional impact of 36 ultra-rare variants of CFI. These variants were selected based on their minor allele frequencies (MAF) and their association with low-normal FI levels. Four variants lead to aberrant splicing-one 5' consensus splice site (NM_000204.5: c.1429G>C, p. Asp477His) and three exonic changes (c.355G>A, p.Gly119Arg; c.472G>A, p.Gly158Arg; and c.950G>A, p.Arg317Gln)- enabling their reclassification to likely pathogenic (LP) or pathogenic (P) based on ACMG guidelines. These findings underscore the value of functional assays, such as the minigene assay, in assessing the clinical relevance of rare variants in CFI. [ABSTRACT FROM AUTHOR]

Motivation Selective sweeps can successfully be distinguished from neutral genetic data using summary statistics and likelihood-based methods that analyze single nucleotide polymorphisms (SNPs). However, these methods are sensitive to confounding factors, such as severe population bottlenecks and old migration. By virtue of machine learning, and specifically convolutional neural networks (CNNs), new accurate classification models that are robust to confounding factors have been recently proposed. However, such methods are more computationally expensive than summary-statistic-based ones, yielding them impractical for processing large-scale genomic data. Moreover, SNP data are frequently preprocessed to improve classification accuracy, further exacerbating the long analysis times. Results To this end, we propose a 1D CNN-based model, dubbed FAST-NN, that does not require any preprocessing while using only derived allele frequencies instead of summary statistics or raw SNP data, thereby yielding a sample-size-invariant, scalable solution. We evaluated several data fusion approaches to account for the variance of the density of genetic diversity across genomic regions (a selective sweep signature), and performed an extensive neural architecture search based on a state-of-the-art reference network architecture (SweepNet). The resulting model, FAST-NN, outperforms the reference architecture by up to 12% inference accuracy over all challenging evolutionary scenarios with confounding factors that were evaluated. Moreover, FAST-NN is between 30× and 259× faster on a single CPU core, and between 2.0× and 6.2× faster on a GPU, when processing sample sizes between 128 and 1000 samples. Our work paves the way for the practical use of CNNs in large-scale selective sweep detection. Availability and implementation https://github.com/SjoerdvandenBelt/FAST-NN [ABSTRACT FROM AUTHOR]

Mitochondria play a crucial role in numerous biological processes; however, limited methods and research have focused on revealing mitochondrial heterogeneity at the single‐cell level. In this study, we optimized the DNBelab C4 single‐cell ATAC (assay for transposase‐accessible chromatin) sequencing workflow for single‐cell mitochondrial sequencing (C4_mtscATAC‐seq). We validated the effectiveness of our C4_mtscATAC‐seq protocol by sequencing the HEK‐293T cell line with two biological replicates, successfully capturing both mitochondrial content (~68% of total sequencing data) and open chromatin status simultaneously. Subsequently, we applied C4_mtscATAC‐seq to investigate two mouse tissues, spleen and bone marrow, obtained from two mice aged 2 months and two mice aged 23 months. Our findings revealed higher mitochondrial DNA (mtDNA) content in young tissues compared to more variable mitochondrial content in aged tissues, consistent with higher activity scores of nuclear genes associated with mitochondrial replication and transcription in young tissues. We detected a total of 22, 15, and 21 mtDNA mutations in the young spleen, aged spleen, and bone marrow, respectively, with most variant allele frequencies (VAF) below 1%. Moreover, we observed a higher number of mtDNA mutations with higher VAF in aged tissues compared to young tissues. Importantly, we identified three mtDNA variations (m.9821A>T, m.15219T>C, and m.15984C>T) with the highest VAF in both aged spleen and aged bone marrow. By comparing cells with and without these mtDNA variations, we analyzed differential open chromatin status to identify potential genes associated with these mtDNA variations, including transcription factors such as KLF15 and NRF1. Our study presents an alternative single‐cell mitochondrial sequencing method and provides crude insights into age‐related single‐cell mitochondrial variations. [ABSTRACT FROM AUTHOR]

Motivation Metagenome-Assembled Genomes (MAGs) or Single-cell Amplified Genomes (SAGs) are often incomplete, with sequences missing due to errors in assembly or low coverage. This presents a particular challenge for the identification of true gene frequencies within a microbial population, as core genes missing in only a few assemblies will be mischaracterized by current pangenome approaches. Results Here, we present CELEBRIMBOR, a Snakemake pangenome analysis pipeline which uses a measure of genome completeness to automatically adjust the frequency threshold at which core genes are identified, enabling accurate core gene identification in MAGs and SAGs. Availability and implementation CELEBRIMBOR is published under open source Apache 2.0 licence at https://github.com/bacpop/CELEBRIMBOR and is available as a Docker container from this repository. Supplementary material is available in the online version of the article. [ABSTRACT FROM AUTHOR]

There is a collective push to diversify human genetic studies by including underrepresented populations. However, analyzing DNA sequence reads involves the initial step of aligning the reads to the GRCh38/hg38 reference genome which is inadequate for non-European ancestries. In this study, using long-read sequencing technology, we constructed de novo genome assemblies from two indigenous Americans from Arizona (IAZ). Each assembly included ∼17 Mb of DNA sequence not present [nonreference sequence (NRS)] in hg38, which consists mostly of repeat elements. Forty NRSs totaling 240 kb were uniquely anchored to the hg38 primary assembly generating a modified hg38-NRS reference genome. DNA sequence alignment and variant calling were then conducted with whole-genome sequencing (WGS) sequencing data from 387 IAZ using both the hg38 and modified hg38-NRS reference maps. Variant calling with the hg38-NRS map identified ∼50,000 single-nucleotide variants present in at least 5% of the WGS samples which were not detected with the hg38 reference map. We also directly assessed the NRSs positioned within genes. Seventeen NRSs anchored to regions including an identical 187 bp NRS found in both de novo assemblies. The NRS is located in HCN2 79 bp downstream of Exon 3 and contains several putative transcriptional regulatory elements. Genotyping of the HCN2 -NRS revealed that the insertion is enriched in IAZ (minor allele frequency = 0.45) compared to other reference populations tested. This study shows that inclusion of population-specific NRSs can dramatically change the variant profile in an underrepresented ethnic groups and thereby lead to the discovery of previously missed common variations. [ABSTRACT FROM AUTHOR]

Frequencies of pharmacogenetic (PGx) variants are known to differ substantially across populations but much of the available PGx literature focuses on one or a few population groups, often defined in nonstandardized ways, or on a specific gene or variant. Guidelines produced by the Clinical Pharmacogenetic Implementation Consortium (CPIC) provide consistent methods of literature extraction, curation, and reporting, including comprehensive curation of allele frequency data across nine defined "biogeographic groups" from the PGx literature. We extracted data from 23 CPIC guidelines encompassing 19 genes to compare the sizes of the populations from each group and allele frequencies of altered function alleles across groups. The European group was the largest in the curated literature for 16 of the 19 genes, while the American and Oceanian groups were the smallest. Nearly 200 alleles were detected in nonreference groups that were not reported in the largest (reference) group. The genes CYP2B6 and CYP2C9 were more likely to have higher frequencies of altered function alleles in nonreference groups compared to the reference group, while the genes CYP4F2, DPYD, SLCO1B1, and UGT1A1 were less likely to have higher frequencies in nonreference groups. PGx allele frequencies and function differ substantially across nine biogeographic groups, all but two of which are underrepresented in available PGx data. Awareness of these differences and increased efforts to characterize the breadth of global PGx variation are needed to ensure that implementation of PGx‐guided drug selection does not further widen existing health disparities among populations currently underrepresented in PGx data. [ABSTRACT FROM AUTHOR]

In this study, genome-wide variability and selection signatures in Kangal Akkaraman sheep were assessed by 238.103 bi-allelic single nucleotide polymorphisms (SNPs) recovered from genotyping-by-sequencing (GBS) libraries processed in Illumina HiSeq X Ten instrument. Summary statistics of genetic diversity such as minor allele frequency (MAF), observed (HO) and expected (HE) heterozygosity were estimated at 0.32, 0.29, and 0.30, respectively. A declining trend in effective population size was observed through generations in which the current population was estimated at 978 individuals 150 generations ago. 608 of 647 runs of homozygosity (ROH) islands were classified into =2 Mb. Strong selection signals were identified in thirteen genomic intervals overlapped with 17 protein-coding genes. The sheep quantitative trait locus (QTL) database confirmed that these genomic regions were associated with economically important traits such as milk content (KCNH5, KCNH7, LRP1B, SNAPC1, and SYT16) and fleece yield (CCDC85A, EFEMP1, and PPP4R3B), parasite resistance (MMS22L and KLHL32), fat deposition in the tail (JAZF1, TAXIBP1, EVX1, and HOXA13), and water-holding capacity (KLHL1 and DACH1). This study implies that the Kangal Akkaraman sheep will play a vital role in developing some genotypes tolerant to environmental challenges, parasite infections, fat deposition, and water-holding capacity in the future. Still, the other native sheep should be screened to identify genomic regions under selection practices using high-density genetic data obtained from next-generation sequencing (NGS) platforms. [ABSTRACT FROM AUTHOR]

Despite advancements in radiologic, laboratory, and pathological evaluations, differentiating between benign and malignant bile duct strictures remains a diagnostic challenge. Recent developments in massive parallel sequencing (MPS) have introduced new opportunities for early cancer detection and management, but these techniques have not yet been rigorously applied to biliary samples. We prospectively evaluated the Oncomine Comprehensive Assay (OCA) and the Oncomine Pan-Cancer Cell-Free Assay (OPCCFA) using biliary brush cytology and bile fluid obtained via endoscopic retrograde cholangiopancreatography from patients with bile duct strictures. The diagnostic performance of MPS testing was assessed and compared to the pathological findings of biliary brush cytology and primary tissue. Mutations in TP53, BRAF, CTNNB1, SMAD4, and K-/N-RAS identified in biliary brush cytology samples were also detected in the corresponding bile fluid samples from patients with extrahepatic cholangiocarcinoma. These mutations were also identified in the bile fluid samples, but with variant allele frequencies lower than those in the corresponding biliary brush cytology samples. In control patients diagnosed with gallstones, neither the biliary brush cytology samples nor the bile fluid samples showed any pathogenic mutations classified as tier 1 or 2. Our study represents a prospective investigation into the role of MPS-based molecular testing in evaluating bile duct strictures. MPS-based molecular testing shows promise in identifying actionable genomic alterations, potentially enabling the stratification of patients for targeted chemotherapeutic treatments. Future research should focus on integrating OCA and OPCCFA testing, as well as similar MPS-based assays, into existing surveillance and management protocols for patients with bile duct strictures. [ABSTRACT FROM AUTHOR]

Background & Objective: Hypovitaminosis D has shown to be linked with T2DM development and control in numerous studies. The association of SNPs in genes related to VitD metabolism with T2DM has not been sufficiently studied. Consequently, our aim in the present study was to explore the association between genetic variants in genes connected with VitD, mainly a SNP in GC (rs1155563), a SNP in DHCR7 (rs12785878) and a SNP in CYP2R1 (rs10500804) with glycaemic parameters in females with T2DM in Saudi Arabia. Methods: The cross-sectional study included 149 females (age 38-52 years) with T2DM from Jeddah, Saudi Arabia (September 2022-March 2023). Blood was extracted from the participants for biochemical tests including measuring VitD [25(OH)D] concentration, parameters of glycaemia (HbA1c, insulin, fasting glucose and insulin sensitivity indices including HOMA2-IR and HOMA2-%β), and for genomic DNA isolation. Sanger DNA sequencing was used to screen for VitD genetic polymorphisms (rs1155563, rs12785878 and rs10500804). Results: Minor allele frequency for rs1155563C, rs12785878T and rs10500804G was 0.21, 0.23 and 0.37, respectively. Levels of 25(OH)D and glycaemic parameters as well did not show any significant difference between the genotypes of each SNP. Conclusion: This study showed lack of association of rs1155563 in GC, rs12785878 in DHCR7 and rs10500804 in CYP2R1 with VitD level primarily and with glycaemic parameters secondarily. Additional research is required to explore further other VitD genetic polymorphisms influencing T2DM which might lead consequently to genetically-based personalized management for T2DM. [ABSTRACT FROM AUTHOR]

The list of standard abbreviations for JDS is available at adsa.org/jds-abbreviations-24. Nonstandard abbreviations are available in the Notes. This study aimed to evaluate the impact of copy number variants (CNV) on 13 reproduction and 12 disease traits in Holstein cattle. Intensity signal files containing log R ratio and B allele frequency information from 13,730 Holstein animals genotyped with a 95K SNP panel, and 8,467 Holstein animals genotyped with a 50K SNP panel were used to identify the CNVs. Subsequently, the identified CNVs were validated using whole-genome sequence data from 126 animals, resulting in 870 high-confidence copy number variant regions (CNVR) on 12,131 animals. Out of these, 54 CNVR had frequencies higher than or equal to 1% in the population and were used in the genome-wide association analysis (one CNVR at a time, including the G matrix). Results revealed that 4 CNVR were significantly associated with at least one of the traits analyzed in this study. Specifically, 2 CNVR were associated with 3 reproduction traits (i.e., calf survival, first service to conception, and nonreturn rate), and 2 CNVR were associated with 2 disease traits (i.e., metritis and retained placenta). These CNVR harbored genes implicated in immune response, cellular signaling, and neuronal development, supporting their potential involvement in these traits. Further investigations to unravel the mechanistic and functional implications of these CNVR on the mentioned traits are warranted. [ABSTRACT FROM AUTHOR]

Purpose: Aging is characterized by the gradual physiological changes and alterations that accumulate over time in the human body. The combination of obesity and ageing can lead to an increased risk of serious health issues or death. Single nucleotide variants (SNVs) in the Klotho gene were commonly studied, including that in type 2 diabetes mellitus (T2DM). Aim: The aim of this study is to examine the possible effect of SNVs in Klotho on the obese population in Saudi Arabia using middle-aged participants with and without T2DM. Methods: This study consists of 100 controls and 100 obesity patients, in which 50 had T2DM and the remaining 50 were obese without T2DM. Genotyping was performed with PCR, and Sanger sequencing analysis was used to validate the molecular association. Results: In this study, rs1207568 (p = 0.001– 0.003) and rs9527025 (p = 0.001– 0.00004) SNVs were associated with obesity cases. However, none of the genotypes or allele frequencies showed a positive association with the rs564481 SNV (p = 0.344– 0.881). The multiple linear regression model showed that waist and hip were associated (p = 0.01– 0.02). ANOVA analysis showed age (p = 0.04), hip (p = 0.002), SBP, and TC (p = 0.02) were associated. Finally, SNV (rs1207568 and rs95270250) and obesity (p < 0.001) associations were confirmed through gene multifactor dimensionality reduction analysis with gene–gene interaction, dendrogram, and graphical depletion method. Conclusion: This study concludes that rs1207568 and rs9527025 SNVs are associated with obesity in the Saudi population. Additional genetical statistics showed significant association between dependent and independent variables. SNVs in Klotho play a role in the Saudi population's susceptibility to obesity. [ABSTRACT FROM AUTHOR]

The black-striped pipefish (Syngnathus abaster) is an Atlanto-Mediterranean marine and estuarine fish species which recently expanded its range into freshwaters. iPBS genetic markers analysis revealed the population divergence of this fish species from different localities in Ukraine. In the course of molecular genetic analysis using four inter-primer binding site (iPBS) primers, 82 amplification fragments were detected. The difference between populations depends on the frequency of a particular allele allows us to separate the population from the Tylihul Estuary and all freshwater ones. The samples from the Dnipro River basin (Stuhna River and Dnipro Reservoir) related to different subclusters, which may indicate that the introduction of the pipefish into the Dnipro basin could have occurred several times and in different ways: 1. Natural range expansion; 2. Artificial acclimatisation; 3. Introduction with river transport. The differences in allele frequencies suggest the existence of populations with different origins in the Dnipro River, which can be explained by the introduction of this species from the delta areas of the Danube and Dniester (population from the Stuhna River) and self-settlement upstream of the Dnipro River (population from the Dnipro Reservoir). The black-striped pipefish (Syngnathus abaster) is an Atlanto-Mediterranean marine and estuarine fish species which recently expanded its range into freshwaters. iPBS genetic markers analysis revealed the population divergence of this fish species from different localities in Ukraine. In the course of molecular genetic analysis using four inter-primer binding site (iPBS) primers, 82 amplification fragments were detected. The difference between populations depends on the frequency of a particular allele allows us to separate the population from the Tylihul Estuary and all freshwater ones. The samples from the Dnipro River basin (Stuhna River and Dnipro Reservoir) related to different subclusters, which may indicate that the introduction of the pipefish into the Dnipro basin could have occurred several times and in different ways: 1. Natural range expansion; 2. Artificial acclimatisation; 3. Introduction with river transport. The differences in allele frequencies suggest the existence of populations with different origins in the Dnipro River, which can be explained by the introduction of this species from the delta areas of the Danube and Dniester (population from the Stuhna River) and self-settlement upstream of the Dnipro River (population from the Dnipro Reservoir). [ABSTRACT FROM AUTHOR]

Aging is a nearly inescapable trait among organisms yet lifespan varies tremendously across different species and spans several orders of magnitude in vertebrates alone. This vast phenotypic diversity is driven by distinct evolutionary trajectories and tradeoffs that are reflected in patterns of diversification and constraint in organismal genomes. Age-specific impacts of selection also shape allele frequencies in populations, thus impacting disease susceptibility and environment-specific mortality risk. Further, the mutational processes that spawn this genetic diversity in both germline and somatic cells are strongly influenced by age and life history. We discuss recent advances in our understanding of the evolution of aging and lifespan at organismal, population, and cellular scales, and highlight outstanding questions that remain unanswered.

Transcriptome prediction models built with data from European-descent individuals are less accurate when applied to different populations because of differences in linkage disequilibrium patterns and allele frequencies. We hypothesized that methods that leverage shared regulatory effects across different conditions, in this case, across different populations, may improve cross-population transcriptome prediction. To test this hypothesis, we made transcriptome prediction models for use in transcriptome-wide association studies (TWASs) using different methods (elastic net, joint-tissue imputation [JTI], matrix expression quantitative trait loci [Matrix eQTL], multivariate adaptive shrinkage in R [MASHR], and transcriptome-integrated genetic association resource [TIGAR]) and tested their out-of-sample transcriptome prediction accuracy in population-matched and cross-population scenarios. Additionally, to evaluate model applicability in TWASs, we integrated publicly available multiethnic genome-wide association study (GWAS) summary statistics from the Population Architecture using Genomics and Epidemiology (PAGE) study and Pan-ancestry genetic analysis of the UK Biobank (PanUKBB) with our developed transcriptome prediction models. In regard to transcriptome prediction accuracy, MASHR models performed better or the same as other methods in both population-matched and cross-population transcriptome predictions. Furthermore, in multiethnic TWASs, MASHR models yielded more discoveries that replicate in both PAGE and PanUKBB across all methods analyzed, including loci previously mapped in GWASs and loci previously not found in GWASs. Overall, our study demonstrates the importance of using methods that benefit from different populations' effect size estimates in order to improve TWASs for multiethnic or underrepresented populations.

The accurate determination of allele frequencies is crucially important across a wide range of problems in genetics, such as developing population genetic models, making inferences from genome-wide association studies, determining genetic risk for diseases, as well as other scientific and medical applications. Furthermore, understanding how allele frequencies change over time in populations is central to ascertaining their evolutionary dynamics. We present a precise, efficient, and economical method (FREQ-Seq2) for quantifying the relative frequencies of different alleles at loci of interest in mixed population samples. Through the creative use of paired barcode sequences, we exponentially increased the throughput of the original FREQ-Seq method from 48 to 2,304 samples. FREQ-Seq2 can be targeted to specific genomic regions of interest, which are amplified using universal barcoded adapters to generate Illumina sequencing libraries. Our enhanced method, available as a kit along with open-source software for analyzing sequenced libraries, enables the detection and removal of errors that are undetectable in the original FREQ-Seq method as well as other conventional methods for allele frequency quantification. Finally, we validated the performance of our sequencing-based approach with a highly multiplexed set of control samples as well as a competitive evolution experiment in Escherichia coli and compare the latter to estimates derived from manual colony counting. Our analyses demonstrate that FREQ-Seq2 is flexible, inexpensive, and produces large amounts of data with low error, low noise, and desirable statistical properties. In summary, FREQ-Seq2 is a powerful method for quantifying allele frequency that provides a versatile approach for profiling mixed populations.

Over 80% of genetic studies in the Parkinson's disease (PD) field have been conducted on individuals of European descent. There is a social and scientific imperative to understand the genetic basis of PD across global populations for therapeutic development and deployment. PD etiology is impacted by genetic and environmental factors that are variable by ancestry and region, emphasising the need for worldwide programs to gather large numbers of patients to identify novel candidate genes and risk loci involved in disease. Only a handful of documented genetic assessments have investigated families with PD in AfrAbia, which comprises the member nations of the Arab League and the African Union, with very limited cohort and case-control studies reported. This review article summarises prior research on PD genetics in AfrAbia, highlighting gaps and challenges. We discuss the etiological risk spectrum in the context of historical interactions, highlighting allele frequencies, penetrance, and the clinical manifestations of known genetic variants in the AfrAbian PD patient community. [ABSTRACT FROM AUTHOR]

Objective To study the RHD genotyping of primary screening RhD negative population in Beijing, and analyze the serological characteristics of different genotypes and the distribution of Rh phenotypes. Methods From August 2022 to January 2024, 204 RhD negative samples in our hospital were initially screened by microplate or microcolumn agglutination method. RhD antigen phenotype was identified by saline tube method, and RhD negative confirmation test was performed by microcolumn gel-indirect antiglobulin technology.RHD genotyping was performed by polymerasing chain reaction-sequence specific primmer (PCR-SSP) method and sequencing technology, and the phenotypic frequency and gene frequency were calculated. Results The distribution of Rh serologic phenotypes in 204 RhD negative samples were 112 cases of ccee (54.90%) > 73 cases of Ccee (35.78%) > 9 cases of CCee (4.41%) > 6 cases of ccEe (2.94%) > 4 cases of CcEe (1.96%). A total of 7 RHD genotypes were detected in this study, including 129 cases of RHD* 01N.01 (63.24%), 44 cases of RHD* 01EL.01 (21.57%), 26 cases of RHD* 01N.03 (12.75%), 2 cases of RHD* 01N.16 (0.98%), 1 case of RHD* 06.03.01 (0.49%), 1 case of weak RHD* 15 (0.49%) and 1 case of RHD* 13.01/RHD* 01N.01 (0.49%). The frequency and phenotype distribution of Rh blood group antigen gene were consistent with Hardy-Weinberg equilibrium (P>0.05). Conclusion The RHD gene is polymorphic in the primary screening RhD negative population in Beijing, and the most common RhD negative genotype is RHD* 01N.01, and the RHD variant is RHD* 01EL.01.This study may provide data to support the clinical implementation of genotype-matched blood transfusions.

Abstract Background Factor XI (FXI) deficiency is an autosomal hemorrhagic disorder characterized by reduced plasma FXI levels. Multiple ancestral variants in the F11 gene have been identified in Ashkenazi Jews and other selected European populations. However, there are few reports of predominant variants in Chinese and/or East Asian populations. The aim of this study is to characterize the genotypes and phenotypes of FXI deficiency and identify the predominant variants. Results Of the 41 FXI-deficient patients, 39 exhibited severe FXI defects, considerably more than those with partial defects. The APTT levels showed a negative correlation with FXI activity levels (coefficient=-0.584, P

Introduction: The correct labeling of a genetic variant as pathogenic is important as breeding decisions based on incorrect DNA tests can lead to the unwarranted exclusion of animals, potentially compromising the long-term health of a population. In human medicine, the American college of Medical Genetics (ACMG) guidelines provide a framework for variant classification. This study aims to apply these guidelines to six genetic variants associated with hypertrophic cardiomyopathy (HCM) in certain cat breeds and to propose a modified criterion for variant classification. Methods: Genetic samples were sourced from five cat breeds: Maine Coon, Sphynx, Ragdoll, Devon Rex, and British Short- and Longhair. Allele frequencies were determined, and in the subset with phenotypes available, odds ratios to determine the association with HCM were calculated. In silico evaluation followed with joint evidence and data from other publications assisting in the classification of each variant. Results: Two variants, MYBPC3:c.91G  >  C [A31P] and MYBPC3:c.2453C  >  T [R818W], were designated as pathogenic. One variant, MYH7:c.5647G  >  A [E1883K], was found likely pathogenic, while the remaining three were labeled as variants of unknown significance. Discussion: Routine genetic testing is advised solely for the MYBPC3:c.91G  >  C [A31P] in the Maine Coon and MYBPC3:c.2453C  >  T [R818W] in the Ragdoll breed. The human ACMG guidelines serve as a suitable foundational tool to ascertain which variants to include; however, refining them for application in veterinary medicine might be beneficial. [ABSTRACT FROM AUTHOR]

Digital PCR (dPCR) is a powerful method for highly sensitive and precise quantification of nucleic acids. However, designing and optimizing new multiplex dPCR assays using target sequence specific probes remains cumbersome, since fluorescent signals must be optimized for every new target panel. As a solution, we established a generic fluorogenic 6-plex reporter set, based on mediator probe technology, that decouples target detection from signal generation. This generic reporter set is compatible with different target panels and thus provides already optimized fluorescence signals from the start of new assay development. Generic reporters showed high population separability in a colorimetric 6-plex mediator probe dPCR, due to their tailored fluorophore and quencher selection. These reporters were further tested using different KRAS, NRAS and BRAF single-nucleotide polymorphisms (SNP), which are frequent point mutation targets in liquid biopsy. We specifically quantified SNP targets in our multiplex approach down to 0.4 copies per microliter (cp/µL) reaction mix, equaling 10 copies per reaction, on a wild-type background of 400 cp/µL for each, equaling 0.1% variant allele frequencies. We also demonstrated the design of an alternative generic reporter set from scratch in order to give detailed step-by-step guidance on how to systematically establish and optimize novel generic reporter sets. Those generic reporter sets can be customized for various digital PCR platforms or target panels with different degrees of multiplexing. [ABSTRACT FROM AUTHOR]

The HLA genes are associated with various autoimmune pathologies, with the control of the immune response also being significant in organs and cells transplantation. The aim of the study is to identify the HLA-A, HLA-B, and HLA-C alleles frequencies in the analyzed Romanian cohort. We performed HLA typing using next-generation sequencing (NGS) in a Romanian cohort to estimate class I HLA allele frequencies up to a six-digit resolution. A total of 420 voluntary donors from the National Registry of Voluntary Hematopoietic Stem Cell Donors (RNDVCSH) were included in the study for HLA genotyping. Peripheral blood samples were taken and brought to the Fundeni Clinical Institute during 2020–2021. HLA genotyping was performed using the Immucor Mia Fora NGS MFlex kit. A total of 109 different alleles were detected in 420 analyzed samples, out of which 31 were for HLA-A, 49 for HLA-B, and 29 for HLA-C. The most frequent HLA-A alleles were HLA-A*02:01:01 (26.11%), HLA-A*01:01:01 (12.5%), HLA-A*24:02:01 (11.67%), HLA-A*03:01:01 (9.72%), HLA-A*11:01:01, and HLA-A*32:01:01 (each with 8.6%). For the HLA-B locus, the most frequent allele was HLA-B*18:01:01 (11.25%), followed by HLA-B*51:01:01 (10.83%) and HLA-B*08:01:01 (7.78%). The most common HLA-C alleles were HLA-C*07:01:01 (17.36%), HLA-C*04:01:01 (13.47%), and HLA-C*12:03:01 (10.69%). Follow-up studies are ongoing for confirming the detected results. [ABSTRACT FROM AUTHOR]

Introduction: Spillover events of Mycoplasma ovipneumoniae have devastating effects on the wild sheep populations. Multilocus sequence typing (MLST) is used to monitor spillover events and the spread of M. ovipneumoniae between the sheep populations. Most studies involving the typing of M. ovipneumoniae have used Sanger sequencing. However, this technology is time-consuming, expensive, and is not well suited to efficient batch sample processing. Methods: Our study aimed to develop and validate an MLST workflow for typing of M. ovipneumoniae using Nanopore Rapid Barcoding sequencing and multiplex polymerase chain reaction (PCR). We compare the workflow with Nanopore Native Barcoding library preparation and Illumina MiSeq amplicon protocols to determine the most accurate and cost-effective method for sequencing multiplex amplicons. A multiplex PCR was optimized for four housekeeping genes of M. ovipneumoniae using archived DNA samples (N = 68) from nasal swabs. Results: Sequences recovered from Nanopore Rapid Barcoding correctly identified all MLST types with the shortest total workflow time and lowest cost per sample when compared with Nanopore Native Barcoding and Illumina MiSeq methods. Discussion: Our proposed workflow is a convenient and effective method for strain typing of M. ovipneumoniae and can be applied to other bacterial MLST schemes. The workflow is suitable for diagnostic settings, where reduced hands-on time, cost, and multiplexing capabilities are important. [ABSTRACT FROM AUTHOR]