Nucleotide Diversity Research Articles

Dog10K: an integrated Dog10K database summarizing canine multi-omics.

The diversity observed in canine breed phenotypes, together with their risk for heritabily disorders of relevance to dogs and humans, makes the species an ideal subject for studies aimed at understanding the genetic basis of complex traits and human biomedical models. Dog10K is an ongoing international collaboration that aims to uncover the genetic basis of phenotypic diversity, disease, behavior, and domestication history of dogs. To best present and make the extensive data accessible and user friendly, we have established the Dog10K (http://dog10k.kiz.ac.cn/) database, a comprehensive-omics resource summarizing multiple types of data. This database integrates single nucleotide variants (SNVs) from 1987 canine genomes, de-novo mutations (DNMs) from 43 dog breeds with>40×sequence, RNA-seq data of 105057 single nuclei from hippocampus, 74067 single cells from leukocytesand 30 blood samples from published canid studies. We provide clear visualization, statistics, browse, searching, and downloading functions for all data. We have integrated three analysis tools, Selscan, LiftOver and AgeConversion, to aid researchers in custom exploration of the comprehensive-omics data. The Dog10K database will serve as a foundational platform for analyzing, presenting and utilizing canine multi-omics data.

Clinical impact of large genomic explorations at diagnosis in 198 pediatric solid tumors: a monocentric study aiming practical feasibility of precision oncology

IntroductionFaced to the growing development of collecting systematic molecular analyses in relapsed pediatric cancers to transform their targeted matched therapies, this study aimed to assess the clinical and therapeutic indications of systematic diagnostic genomic explorations performed in pediatric solid cancers to determine which type of screening and if it afford at relapse time an accurate targeted strategy.MethodsA total of 280 patients less than 22 years, referred at the University Hospitals of Strasbourg for a newly diagnosed solid tumor from January 2015 to December 2021, were prospectively genomically investigated since diagnosis. Using 7 different molecular tests going from single-gene methods (IHC, FISH, RT-PCR, Sanger sequencing, droplet digital PCR) to largescale analyses (Next-Generation sequencing, RNAsequencing and FoundationOne®CDx), we explored retrospectively the molecular findings in those pediatric solid tumors (except hematolymphoid cancers) to improve diagnosis, prognosis assessment and relapse therapeutics.ResultsOne hundred and ninety-eight patients (71%) underwent molecular biology (MB) at diagnosis. Thirty-eight different histologies were grouped into cerebral tumors (30%), sarcomas (26%, bone and soft tissues), various blastomas (27%), and other entities (17%). Over a median 40-month follow-up, the overall survival rate of patients was 85% and the relapse rate 28%. Of the 326 analyses carried out, 245 abnormalities (single nucleotide variations: 50%, fusions: 25%, copy number alteration: 20%) concerning 70 oncogenes were highlighted. The overall clinical impact rate was 84%. Broad-spectrum analyses had a higher therapeutic impact (57%) than the targeted analyses (28%). 75% of broad-spectrum tests found an actionable variant conducting 23% of patients to receive rapidly a matched targeted therapy since first relapse.ConclusionOur experience highlighted the clinical utility of molecular profiling of solid tumors as soon as at diagnosis in children to expect improving access to innovative agents at relapse.

A Comprehensive Assessment of Nutritional Value, Antioxidant Potential, and Genetic Diversity in Metapenaeus ensis from Three Different Populations.

Due to its high tolerance to salinity and temperature, as well as its strong adaptability, Metapenaeus ensis holds an important position in the Chinese aquaculture industry. However, studies on the evaluation of its germplasm resources remain insufficient. This research conducted an in-depth comparative evaluation of M. ensis from three representative farming regions in China: Sanya, Zhuhai, and Raoping. The nutritional analysis of muscle tissue showed no statistically significant differences in crude ash, moisture, and crude protein content among the populations (p > 0.05). However, significant differences were observed in crude fat and total sugar content (p < 0.05). The MeSY and MeRP populations had higher crude fat content than the MeZH population (p < 0.05), while the MeZH population exhibited the highest total sugar content. In terms of amino acid composition, the MeSY population had relatively higher total essential amino acid content and proportion, as well as higher total amino acid content, both of which were statistically significant (p < 0.05). A fatty acid composition analysis further highlighted the advantages of the MeRP population in several key fatty acids (p < 0.05). Physiological and biochemical analyses showed no significant differences among the three populations in total antioxidant capacity, superoxide dismutase activity, or catalase activity (p > 0.05). A genetic diversity analysis indicated that M. ensis has relatively low diversity, with the MeSY population showing higher SNP density and nucleotide diversity. A genetic differentiation analysis revealed significant genetic differentiation between the MeSY and MeZH populations, while differentiation between the MeZH and MeRP populations was relatively smaller. This comprehensive assessment of nutritional components, amino acids, fatty acids, antioxidant capacity, and genetic diversity highlights the advantages of germplasm resources from different regions. These findings provide valuable insights for future research on the genetic characteristics and breeding potential of M. ensis.

Optical genome mapping of structural variants in Parkinson’s disease-related induced pluripotent stem cells

BackgroundCertain structural variants (SVs) including large-scale genetic copy number variants, as well as copy number-neutral inversions and translocations may not all be resolved by chromosome karyotype studies. The identification of genetic risk factors for Parkinson’s disease (PD) has been primarily focused on the gene-disruptive single nucleotide variants. In contrast, larger SVs, which may significantly influence human phenotypes, have been largely underexplored. Optical genomic mapping (OGM) represents a novel approach that offers greater sensitivity and resolution for detecting SVs. In this study, we used induced pluripotent stem cell (iPSC) lines of patients with PD-linked SNCA and PRKN variants as a proof of concept to (i) show the detection of pathogenic SVs in PD with OGM and (ii) provide a comprehensive screening of genetic abnormalities in iPSCs.ResultsOGM detected SNCA gene triplication and duplication in patient-derived iPSC lines, which were not identified by long-read sequencing. Additionally, various exon deletions were confirmed by OGM in the PRKN gene of iPSCs, of which exon 3–5 and exon 2 deletions were unable to phase with conventional multiplex-ligation-dependent probe amplification. In terms of chromosomal abnormalities in iPSCs, no gene fusions, no aneuploidy but two balanced inter-chromosomal translocations were detected in one line that were absent in the parental fibroblasts and not identified by routine single nucleotide variant karyotyping.ConclusionsIn summary, OGM can detect pathogenic SVs in PD-linked genes as well as reveal genomic abnormalities for iPSCs that were not identified by other techniques, which is supportive for OGM’s future use in gene discovery and iPSC line screening.

Slowly progressive autosomal dominant Alport Syndrome due to COL4A3 splicing variant.

Alport syndrome is a rare genetic kidney disease caused by variants in the COL4A3/A4/A5 genes. It's characterised by progressive kidney failure, though therapies targeting Renin-Angiotensin System can delay its progression. Additionally, extrarenal manifestations may sometimes coexist. Recent advances in genetic analysis and the necessity to better clarify genotype-phenotype correlations in affected patients raises the importance of detecting even cryptic splicing variants, lying in both canonical and non-canonical splice sites variants such as last exonic nucleotide variants. These variants, often, do not cause an amino acid change but alter the snRNP proteins binding. We studied a big Italian family with Alport syndrome showing a clear dominant pattern of transmission with younger family members having only haematuria and older individuals presenting with End-Stage Kidney Failure (ESKF). Kidney biopsy showed the typical disease hallmarks. We deeply mined the data for SNV and CNV through exome sequencing on DNA from both peripheral blood samples and patients' podocytes-lineage cells. We identified an already reported synonymous variant, c.765G>A (p.(Thr255Thr)), in the last exonic nucleotide of exon 13 of the COL4A3 gene. Employing the patient's podocytes we demonstrated that this variant results in exon skipping leading to an in-frame deletion of 28 amino acids without leaky effect. According to the pattern of transmission, to the kidney biopsy and to the exome data analysis we provided further evidence that autosomal dominant Alport syndrome is a well-defined clinical entity. We also confirmed the pathogenicity of the synonymous COL4A3 variant for the first time demonstrating its role in a dominant pattern of transmission.

Phylodynamic Structure in the Botswana HIV Epidemic.

Background Studying viral sequences can provide insights into the structure of host contact networks through which the virus is transmitted. Uncovering the population structure of the HIV-1 epidemic in Botswana will help optimise public health interventions and may identify hidden sub-epidemics. We sought to determine the phylodynamic structure of the Botswana HIV-1 epidemic from viral sequence genetic data. Methods The Botswana Combination Prevention Project (BCPP) randomly sampled 20% of households in 30 villages in Botswana between 2013-2018 and tested for HIV-1. Extensive demographic data were collected from all participants and next-generation full-genome HIV-1 sequences were generated from HIV-1 positive participants (n = 4,164), 78% of whom were on antiretroviral treatment (ART). We inferred the stage of infection (< or > 1 year) among HIV-1 cases based on nucleotide diversity and clinical data using a previously trained machine learning model. We then reconstructed time-resolved gag and pol phylogenies from sequences, other Botswana cohorts and publicly available sequences that were genetically close to those from Botswana. We statistically explored phylogenies for partitions with diverging patterns of coalescence, indicating sub-epidemics, and estimated viral effective population size through time, a measure of viral incidence, for each partition. Finally, we compared the demographic makeup, clinical and geographic characteristics across partitions using χ2, ANOVA tests and Tukey analysis. Results We identified three partitions of time-resolved gag and pol phylogenies, revealing divergent patterns of coalescence and HIV-1 transmission. In both gag and pol phylogenies, partitions with persistent growth and transmission were characterised by lower treatment coverage and more recent infections when compared to other partitions. The Southern and South East regions of Botswana were over-represented in the fast-growing partitions. Conclusion Our findings suggest that transmission is slowing in segments of the population that have high ART coverage. However, recent infections are over-represented in ongoing sub-epidemics. The phylodynamic structure suggests that there are districts with higher growth and prioritising these in the deployment of public health interventions might curb new infections. Nonetheless the high mobility of Botswana residents should be taken into consideration in implementing effective interventions to combat HIV-1.

Comparative genomics of Fusarium species causing Fusarium ear rot of maize.

Fusarium ear rot (FER), caused by the fungal pathogen Fusarium verticillioides, stands as one of the most economically burdensome and pervasive diseases affecting maize worldwide. Its impact on food security is particularly pronounced due to the production of fumonisins, toxic secondary metabolites that pose serious health risks, especially for livestock. FER disease severity is complex and polygenic, with few resistance (R) genes being identified for use in breeding resistant varieties. While FER is the subject of several breeding programs, only a few studies have investigated entire populations of F. verticillioides with corresponding virulence data to better understand and characterize the pathogenomics. Here, we sequenced and compared the genomes of 50 Fusarium isolates (43 F. verticillioides and 7 other Fusarium spp.) that were used to inoculate a diverse maize population. Our objectives were to elucidate the genome size and composition of F. verticillioides, explore the variable relationship between fumonisin production and visual disease severity, and shed light on the phylogenetic relationships among the isolates. Additionally, we conducted a comparative analysis of the nucleotide variants (SNPs) and the isolates' effectoromes to uncover potential genetic determinants of pathogenicity. Our findings revealed several promising leads, notably the association of certain gene groups, such as pectate lyase, with disease severity. These genes should be investigated further as putative alleles for breeding resistant maize varieties. We suggest that, beyond validation of the alleles identified in this study, researchers validate each phenotypic dataset on an individual basis, particularly if considering fumonisin concentrations and when using diverse populations. Our study underscores the importance of genomic analysis in tackling FER and offers insights that could inform the development of resilient maize cultivars. By leveraging advances in genomics and incorporating pathogen populations into breeding programs, resistance to FER can be advanced.

Behavioral screening reveals a conserved residue in Y-Box RNA-binding protein required for associative learning and memory in C. elegans.

RNA-binding proteins (RBPs) regulate translation and plasticity which are required for memory. RBP dysfunction has been linked to a range of neurological disorders where cognitive impairments are a key symptom. However, of the 2,000 RBPs in the human genome, many are uncharacterized with regards to neurological phenotypes. To address this, we used the model organism C. elegans to assess the role of 20 conserved RBPs in memory. We identified eight previously uncharacterized memory regulators, three of which are in the C. elegans Y-Box (CEY) RBP family. Of these, we determined that cey-1 is the closest ortholog to the mammalian Y-Box (YBX) RBPs. We found that CEY-1 is both necessary in the nervous system for memory ability and sufficient to promote memory. Leveraging human datasets, we found both copy number variation losses and single nucleotide variants in YBX1 and YBX3 in individuals with neurological symptoms. We identified one predicted deleterious YBX3 variant of unknown significance, p.Asn127Tyr, in two individuals with neurological symptoms. Introducing this variant into endogenous cey-1 locus caused memory deficits in the worm. We further generated two humanized worm lines expressing human YBX3 or YBX1 at the cey-1 locus to test evolutionary conservation of YBXs in memory and the potential functional significance of the p.Asn127Tyr variant. Both YBX1/3 can functionally replace cey-1, and introduction of p.Asn127Tyr into the humanized YBX3 locus caused memory deficits. Our study highlights the worm as a model to reveal memory regulators and identifies YBX dysfunction as a potential new source of rare neurological disease.

Stabilizing selection and adaptation shape cis and trans gene expression variation in C. elegans.

An outstanding question in the evolution of gene expression is the relative influence of neutral processes versus natural selection, including adaptive change driven by directional selection as well as stabilizing selection, which may include compensatory dynamics. These forces shape patterns of gene expression variation within and between species, including the regulatory mechanisms governing expression in cis and trans . In this study, we interrogate intraspecific gene expression variation among seven wild C. elegans strains, with varying degrees of genomic divergence from the reference strain N2, leveraging this system's unique advantages to comprehensively evaluate gene expression evolution. By capturing allele-specific and between-strain changes in expression, we characterize the regulatory architecture and inheritance mode of gene expression variation within C. elegans and assess their relationship to nucleotide diversity, genome evolutionary history, gene essentiality, and other biological factors. We conclude that stabilizing selection is a dominant influence in maintaining expression phenotypes within the species, and the discovery that genes with higher overall expression tend to exhibit fewer expression differences supports this conclusion, as do widespread instances of cis differences compensated in trans . Moreover, analyses of human expression data replicate our finding that higher expression genes have less variable expression. We also observe evidence for directional selection driving expression divergence, and that expression divergence accelerates with increasing genomic divergence. To provide community access to the data from this first analysis of allele-specific expression in C. elegans , we introduce an interactive web application, where users can submit gene-specific queries to view expression, regulatory pattern, inheritance mode, and other information: https://wildworm.biosci.gatech.edu/ase/ .

ScTML: a pan-cancer single-cell landscape of multiple mutation types.

Investigating mutations, including single nucleotide variations (SNVs), gene fusions, alternative splicing and copy number variations (CNVs), is fundamental to cancer study. Recent computational methods and biological research have demonstrated the reliability and biological significance of detecting mutations from single-cell transcriptomic data. However, there is a lack of a single-cell-level database containing comprehensive mutation information in all types of cancer. Establishing a single-cell mutation landscape from the huge emerging single-cell transcriptomic data can provide a critical resource for elucidating the mechanisms of tumorigenesis and evolution. Here, we developed scTML (http://sctml.xglab.tech/), the first database offering a pan-cancer single-cell landscape of multiple mutation types. It includes SNVs, insertions/deletions, gene fusions, alternative splicing and CNVs, along with gene expression, cell states and other phenotype information. The data are from 74 datasets with 2 582 633 cells, including 35 full-length (Smart-seq2) transcriptomic single-cell datasets (all publicly available data with raw sequencing files), 23 datasets from 10X technology and 16 spatial transcriptomic datasets. scTML enables users to interactively explore multiple mutation landscapes across tumors or cell types, analyze single-cell-level mutation-phenotype associations and detect cell subclusters of interest. scTML is an important resource that will significantly advance deciphering intra-tumor and inter-tumor heterogeneity, and how mutations shape cell phenotypes.

DNAJC13 influences responses of the extended reward system to conditioned stimuli: a genome-wide association study.

Reward system dysfunction is implicated in the pathogenesis of major psychiatric disorders. We conducted a genome-wide association study (GWAS) to identify genes that influence activation strength of brain regions within the extended reward system in humans. A homogeneous sample of 214 participants was genotyped and underwent functional magnetic resonance imaging (fMRI). All subjects performed the 'desire-reason dilemma' (DRD) paradigm allowing systematic investigation of systems-level mechanisms of reward processing in humans. As a main finding, we identified the single nucleotide variant rs113408797 in the DnaJ Heat Shock Protein Family Member C13 gene [DNAJC13], alias Receptor-Mediated Endocytosis 8 [RME-8], that was associated with the activation strength of the ventral tegmental area (VTA; p = 2.50E-07) and the nucleus accumbens (NAcc; p = 5.31E-05) in response to conditioned reward stimuli. Moreover, haplotype analysis assessing the information across the entire DNAJC13 locus demonstrated an impact of a five-marker haplotype on VTA activation (p = 3.21E-07), which further corroborates a link between this gene and reward processing. The present findings provide first direct empirical evidence that genetic variation of DNAJC13 influences neural responses within the extended reward system to conditioned stimuli. Further studies are required to investigate the role of this gene in the pathogenesis and pathophysiology of neuropsychiatric disorders.

Conserved Plastid Genomes of Pourthiaea Trees: Comparative Analyses and Phylogenetic Relationship

The genus Pourthiaea Decne., a deciduous woody group with high ornamental value, belongs to the family Rosaceae. Here, we reported newly sequenced plastid genome sequences of Pourthiaea beauverdiana (C. K. Schneid.) Hatus., Pourthiaea parvifolia E. Pritz., Pourthiaea villosa (Thunb.) Decne., and Photinia glomerata Rehder &amp; E. H. Wilson. The plastomes of these three Pourthiaea species shared the typical quadripartite structures, ranging in size from 159,903 bp (P. parvifolia) to 160,090 bp (P. beauverdiana). The three Pourthiaea plastomes contained a pair of inverted repeat regions (26,394–26,399 bp), separated by a small single-copy region (19,304–19,322 bp) and a large single-copy region (87,811–87,973 bp). A total of 113 unique genes were predicted for the three Pourthiaea plastomes, including four ribosomal RNA genes, 30 transfer RNA genes, and 79 protein-coding genes. Analyses of inverted repeat/single-copy boundary, mVISTA, nucleotide diversity, and genetic distance showed that the plastomes of 13 Pourthiaea species (including 10 published plastomes) are highly conserved. The number of simple sequence repeats and long repeat sequences is similar among 13 Pourthiaea species. The three non-coding regions (trnT-GGU-psbD, trnR-UCU-atpA, and trnH-GUG-psbA) were the most divergent. Only one plastid protein-coding gene, rbcL, was under positive selection. Phylogenetic analyses based on 78 shared plastid protein-coding sequences and 29 nrDNA sequences strongly supported the monophyly of Pourthiaea. As for the relationship with other genera in our phylogenies, Pourthiaea was sister to Malus in plastome phylogenies, while it was sister to the remaining genera in nrDNA phylogenies. Furthermore, significant cytonuclear discordance likely stems from hybridization events within Pourthiaea, reflecting complex evolutionary dynamics within the genus. Our study provides valuable genetic insights for further phylogenetic, taxonomic, and species delimitation studies in Pourthiaea, as well as essential support for horticultural improvement and conservation of the germplasm resources.

Comprehensive pan-cancer analysis of FUTs family as prognostic and immunity markers based on multi-omics data

BackgroundThe dysregulation of fucosyltransferases (FUTs) contributes to alterations in fucosylated epitope expression, which serve as distinctive features of cancer cells. Nonetheless, a comprehensive elucidation of the prognostic biological marker and therapeutic target of the FUTs family in pan-cancer remains elusive.MethodsOver 10,000 individuals' profiling information was examined, including information on 750 small molecule drugs, 33 types of cancer, and 24 types of immune cells. We focused on POFUT2's function and applied GSVA (Gene Set Variation Analysis) to calculate the FUT score. Survival and cancer pathways were found to be correlated with this score. After deriving a signature via univariate Cox and LASSO regression, we generated and analyzed the ROC curve and developed a nomogram.ResultsOur comprehensive analysis revealed epigenetic, genomic, and immunogenomic changes in FUTs, particularly POFUT2, resulting in aberrant expression. Elevated frequencies of CNV (Copy number variation), SNV (Single Nucleotide Variant), and hypermethylation were observed in FUTs. Additionally, the survival of patients with various types of cancers may be predicted by FUT expression. Immune response and prognosis in numerous types of cancer were found to be strongly linked to aberrant POFUT2 expression. Pathway analysis unveiled the role of FUTs in apoptosis, epithelial-to-mesenchymal transition (EMT), cell cycle, DNA damage response, RAS/MAPK, TSC/mTOR, PI3K/AKT, AR, ER, and RTK. A prognostic index for patients diagnosed with adrenocortical carcinoma (ACC) was established by applying a risk model incorporating nine FUTs and based on the findings of the GSVA.ConclusionsFUTs, particularly POFUT2, emerge as candidate targets for improving the outcomes of immune therapy. The significance of aberrant MUC12 expression, cancer immune therapy, and patient survival in the context of diverse malignancies is enhanced by the strong correlation observed among these factors. Our five-gene risk signature provides patients with ACC with an independent prognostic indicator, emphasizing the critical function of these genes in inhibiting the immune system's response in ACC.

Molecular characterization and genetic variability of Toxocara vitulorum from naturally infected buffalo calves for the first time in Bangladesh.

Toxocara vitulorum is one of the deadliest parasite of buffalo calves in Bangladesh. This study was conducted to explore genetic variability within and among the T. vitulorum populations in buffalo calves of Bangladesh. Genomic DNA was extracted, ITS2, COX1 and NAD1 gene were amplified and sequenced. Distinct 29 ITS2, 21 unique NAD1 and 24 COX1 genotypes were detected among the T. vitulorum of different geographic regions. These three gene genotypes similarities ranged from 97 to 99%, when these were compared to best hit scoring T. vitulorum sequences retrieved from GenBank. A total of 12 and 6 unique haplotypes were detected for COX1 and NAD1 gene sequences. The average nucleotide and haplotype diversity for COX1 and NAD1 were 0.0931 & 0.89493 and 0.00658 & 0.77895 respectively and the recorded values were more dispersed than previously published values. The pairwise Nst values ranged from −0.050 to 0.602 and Fst from −0.050 to 0.600 between all the T. vitulorum genotypes indicated huge genetic differentiation which were reportedly higher than other published reports Fst values. This is the first report of T. vitulorum on the basis of COX1 gene in Bangladesh. The study findings will be helpful for further extensive epidemiological studies regarding anthelmintic resistance, control and prevention of T. vitulorum infection in buffalo calves.

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries.

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

KBG Syndrome in 16 Indian Individuals.

We aimed to describe the clinical and genetic characteristics of 16 individuals with KBG syndrome (KBGS) from 13 Indian families. We retrospectively analyzed the clinical details of individuals with KBGS harboring a likely pathogenic/pathogenic variant in ANKRD11. We also analyzed their facial gestalt using Face2Gene and recorded the top three differential disorders suggested by the application. The most frequent clinical features observed in our cohort were as follows: learning and intellectual disability-14/15 (93%), skeletal abnormalities-14/15 (93%), postnatal short stature-13/15 (87%), brachydactyly-11/15 (73%), and characteristic facial appearance-13/15 (87%). We identified 12 single nucleotide variants (SNVs), including six recurrent and six novel variants, and a copy number variant in the 16q24.3 region encompassing ANKRD11 gene. The novel variants were as follows: p.(Gln1236Ter), p.(Asp884ThrfsTer93), p.(Arg1466GlyfsTer87), p.(Tyr2056Ter), p.(Leu955TrpfsTer22), and p.(Lys766ArgfsTer10). The identified SNVs in ANKRD11 clustered around exon 9. We observed a high concordance of Face2Gene in predicting KBGS.

Exploring Distribution and Evolution of Pi-ta Haplotypes in Rice Landraces across Different Rice Cultivation Regions in Yunnan.

Background: Rice blast, caused by Magnaporthe oryzae, seriously damages the yield and quality of rice worldwide. Pi-ta is a durable resistance gene that combats M. oryzae carrying AVR-Pita1. However, the distribution of the Pi-ta gene in rice germplasms in Yunnan Province has been inadequately studied. Methods: We analyzed the potential molecular evolution pattern of Pi-ta alleles by examining the diversity in the coding sequence (CDS) among rice varieties. Results: The results revealed that 95% of 405 rice landraces collected from different ecological regions in Yunnan Province carry Pi-ta alleles. We identified 17 nucleotide variation sites in the CDS regions of the Pi-ta gene across 385 rice landraces. These variations led to the identification of 28 Pi-ta haplotypes, encoding 12 novel variants. Among these, 5 Pi-ta haplotypes (62 rice landraces) carried R alleles. The evolutionary cluster and network of the Pi-ta haplotypes suggested that the Pi-ta S alleles were the ancestral alleles, which could potentially evolve into R variants through base substitution. Conclusions: This study suggests that Pi-ta alleles are diverse in the rice landraces in Yunnan, and the Pi-ta sites resistant to blast evolved from the susceptible plants of the rice landraces. These results provide the basis for breeding resistant varieties.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands.

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Machine learning classification of archaea and bacteria identifies novel predictive genomic features

BackgroundArchaea and Bacteria are distinct domains of life that are adapted to a variety of ecological niches. Several genome-based methods have been developed for their accurate classification, yet many aspects of the specific genomic features that determine these differences are not fully understood. In this study, we used publicly available whole-genome sequences from bacteria (N=2546\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=2546$$\\end{document}) and archaea (N=109\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=109$$\\end{document}). From these, a set of genomic features (nucleotide frequencies and proportions, coding sequences (CDS), non-coding, ribosomal and transfer RNA genes (ncRNA, rRNA, tRNA), Chargaff’s, topological entropy and Shannon’s entropy scores) was extracted and used as input data to develop machine learning models for the classification of archaea and bacteria.ResultsThe classification accuracy ranged from 0.993 (Random Forest) to 0.998 (Neural Networks). Over the four models, only 11 examples were misclassified, especially those belonging to the minority class (Archaea). From variable importance, tRNA topological and Shannon’s entropy, nucleotide frequencies in tRNA, rRNA and ncRNA, CDS, tRNA and rRNA Chargaff’s scores have emerged as the top discriminating factors. In particular, tRNA entropy (both topological and Shannon’s) was the most important genomic feature for classification, pointing at the complex interactions between the genetic code, tRNAs and the translational machinery.ConclusionstRNA, rRNA and ncRNA genes emerged as the key genomic elements that underpin the classification of archaea and bacteria. In particular, higher nucleotide diversity was found in tRNA from bacteria compared to archaea. The analysis of the few classification errors reflects the complex phylogenetic relationships between bacteria, archaea and eukaryotes.

Simultaneous Genotyping of Three SNVs, rs5471, rs5472, and rs2000999 Involved in Serum Haptoglobin Levels by Fluorescent Probe-Based Melting Curve Analysis.

Haptoglobin (Hp) is a hemoglobin-binding acute-phase serum protein. Several single nucleotide variations (SNVs) within the Hp gene (HP) or Hp-related protein gene (HPR), such as rs5471 (A>C) and rs5472 (A>G) in HP promoter region and rs2000999 (G>A) in intron 2 of HRP, are suggested to correlate with the serum Hp levels. To determine these three SNVs simultaneously, a genotyping assay based on duplex dual-labeled fluorescent probes was developed. The method was then validated by analyzing genomic DNA from 121 Ghanaian and two Japanese subjects who had been previously genotyped for rs5471, rs5472, and rs2000999. Both rs5471 and rs5472 could be determined as haplotypes with a single FAM-labeled fluorescent probe, and rs2000999 could be genotyped with a HEX-labeled fluorescent probe. The results obtained with the present method were consistent with the previous results except for those of three Ghanaian subjects. All three subjects appear to have multiple HPR copy number variants characteristic of African populations, which may have led to incorrect results during previous genotyping. This method allows us to genotype these three SNVs in a relatively large number of samples, especially in African populations where rs5471 is uniquely distributed.