Species-specific Sequences Research Articles

Major Histocompatibility Complex Class II Genes Allele Diversity in Landlocked Seals.

The allelic diversity of exon 2 (DQB gene) and exon 3 (DRB gene) of major histocompatibility complex class II was studied for the first time in two species of the landlocked pinnipeds, Baikal (N = 79) and Caspian (N = 32) seals, and these were in compared with the widespread Arctic species, the ringed seal (N = 13). The analysis of the second exon comprising the antigen-binding region revealed high allelic diversity in all three species but the pattern of the diversity was the most specific for the Baikal seal. This species differs from the other two by the smallest number of alleles in the population, yet they have the largest number of alleles per individual and by the maximum similarity of individual genotypes. Presumably, this specificity is a consequence of the spatial and temporal homogeneity of the Lake Baikal environment. Analysis of the third exon encoding the conserved β2-domain showed that the Baikal seal differs by the greatest number of amino acid sequences per individual, while the Caspian seal has the lowest number of variants. A single variant of the β2-domain, the same as in the ringed seal, predominates in the Caspian seal, whereas in the Baikal seal the two other variants predominate. At the same time, three species-specific amino acid sequences were observed among minor variants in the Caspian seal, while only one was found in the Baikal seal. This fact may suggest a longer period of independent evolution in the Caspian seal compared to the Baikal seal.

Circadian plasticity evolves through regulatory changes in a neuropeptide gene.

Many organisms, including cosmopolitan drosophilids, show circadian plasticity, varying their activity with changing dawn-dusk intervals1. How this behaviour evolves is unclear. Here we compare Drosophila melanogaster with Drosophila sechellia, an equatorial, ecological specialist that experiences minimal photoperiod variation, to investigate the mechanistic basis of circadian plasticity evolution2. D. sechellia has lost the ability to delay its evening activity peak time under long photoperiods. Screening of circadian mutants in D. melanogaster/D. sechellia hybrids identifies a contribution of the neuropeptide pigment-dispersing factor (Pdf) to this loss. Pdf exhibits species-specific temporal expression, due in part to cis-regulatory divergence. RNA interference and rescue experiments in D. melanogaster using species-specific Pdf regulatory sequences demonstrate that modulation of this neuropeptide's expression affects the degree of behavioural plasticity. The Pdf regulatory region exhibits signals of selection in D. sechellia and across populations of D. melanogaster from different latitudes. We provide evidence that plasticity confers a selective advantage for D. melanogaster at elevated latitude, whereas D. sechellia probably suffers fitness costs through reduced copulation success outside its range. Our findings highlight this neuropeptide gene as a hotspot locus for circadian plasticity evolution that might have contributed to both D. melanogaster's global distribution and D. sechellia's specialization.

Species-specific satellite DNA composition dictates de novo formation of PRC1-mediated pericentric heterochromatin.

Pericentromeres are heterochromatic regions adjacent to centromeres that ensure accurate chromosome segregation. Despite their conserved function, they are composed of rapidly evolving A/T-rich satellite DNA. This paradoxical observation is partially resolved by epigenetic mechanisms that maintain heterochromatin, independent of specific DNA sequences. However, it is unclear how satellite DNA sequence variation impacts de novo formation of pericentric heterochromatin, which is initially absent from paternal chromosomes in the zygote. Here we show that satellite variation has functional consequences for zygotic heterochromatin formation, recruitment of the Chromosome Passenger Complex (CPC), and interactions with spindle microtubules. In M. musculus zygotes, Polycomb Repressive Complex 1 (PRC1) is recruited to pericentric satellites by its AT-hook domain, which binds runs of A/T nucleotides, to generate H2AK119ub1 (H2Aub) heterochromatin. By fertilizing M. musculus eggs with sperm from other mouse species, we show that species-specific satellite sequences differ in their ability to recruit PRC1 and form H2Aub. This satellite-DNA mediated increase in PRC1 heterochromatin leads to reduced CPC recruitment and increased microtubule forces on kinetochores. Our results provide a direct link between satellite DNA composition and pericentromere function in the zygote, when epigenetic pathways maintaining pericentromere heterochromatin are absent.

Differences in seasonal dynamics and pyrethroid resistance development among Anopheles Hyrcanus group species

BackgroundThe Anopheles Hyrcanus group, which transmits Plasmodium vivax, consists of six confirmed species in South Korea. An epidemiological study revealed differences in the seasonal occurrence patterns of each species. Pyrethroid resistance in An. sinensis dates back to the early 2000s, whereas information on pyrethroid resistance in other species is lacking despite their greater significance in malaria epidemiology.MethodsAnopheles mosquitoes were collected from two malaria-endemic regions in South Korea for 2 years and their knockdown resistance (kdr) mutations were genotyped. The larval susceptibility to λ-cyhalothrin was compared in six Anopheles species and its seasonal changes in three species were investigated. The full-length sequences of the voltage-sensitive sodium channel (VSSC) were compared across six species to evaluate potential target-site insensitivity. The contribution of the kdr mutation to phenotypic resistance was confirmed by comparing median lethal time (LT50) to λ-cyhalothrin between populations of Anopheles belenrae with distinct genotypes.ResultsThe composition and seasonal occurrence of rare species (Anopheles kleini, Anopheles lestri, and Anopheles sineroides) varied considerably, whereas An. sinensis occurs continuously throughout the season. A kdr mutation in the form of heterozygous allele was newly identified in An. belenrae, An. lesteri, An. pullus, and An. sineroides. The baseline susceptibility to λ-cyhalothrin was the highest in An. belenrae, followed by An. lesteri, An. sineroides, An. kleini, An. pullus, and An. sinensis, with median lethal concentration (LC50) values ranging from 6.0- to 73.5-fold higher than that of An. belenrae. The susceptibility of An. sinensis and An. pullus varied by season, whereas that of An. belenrae remained stable. The kdr-heterozygous An. belenare population exhibited 5.1 times higher LT50 than that of the susceptible population. Species-specific VSSC sequence differences were observed among the six species.ConclusionsOur findings suggest that the status and extent of pyrethroid resistance vary among Anopheles Hyrcanus group species. While An. sinensis, the predominant species, developed a considerable level of pyrethroid resistance through kdr mutation, the resistance levels of other species appeared to be less pronounced. Large-scale monitoring is crucial to fully understand species-specific seasonal occurrence and resistance status for effective management strategies, considering the ongoing impact of climate change on their vectorial capacity.Graphical

Construction of core genome multi-locus sequence typing schemes for population structure analyses of Nocardia species

Nocardia, a member of the Actinobacteria phylum, populates diverse habitats globally, with certain species being the cause of various clinical infections in humans. There is paucity of data regarding the population structure of this genus and of established genomic-based phylogenetic methods. We examined the whole genome sequences of 193 isolates spanning five major pathogenic Nocardia species sourced from public databases, encompassing diverse geographic regions. Using the chewBBACA pipeline, a species-specific core genome multilocus sequence typing (cgMLST) schema was created for N. cyriacigeorgica, N. farcinica, N. brasiliensis, N. wallacei, and N. abscessus. Additional genomic features that were examined included virulence factor (VF) profile, total length and open-reading frame count, the core genome length and core gene count, and GC content. Our findings indicated that: (i) N. brasiliensis diverges significantly from the other four species, underscoring its distinct evolutionary trajectory; (ii) the population structures of all species were polyclonal, with phylogenetic clustering occurring in the minority of isolates; (iii) clonal complexes were largely restricted to specific geographical locations, rather than demonstrating a global distribution, and (iv) initial evidence suggests no direct common-source transmission amongst the studied strains. Our study establishes a comprehensive genome-based phylogenetic methodology for population structure of Nocardia species.

Detection of Podosphaera macularis in Air Samples by Quantitative PCR.

Detection and quantification of pathogen propagules in the air or other environmental samples is facilitated by culture-independent assays. We developed a quantitative PCR assay for the hop powdery mildew fungus, Podosphaera macularis, for detection of the organism from air samples. The assay uses primers and a TaqMan probe designed to target species-specific sequences in the 28S large subunit of the nuclear ribosomal DNA. Analytical sensitivity was not affected by the presence of an exogenous internal control or potential PCR inhibitors associated with DNA extracted from soil. The level of quantification of the assay was between 200 and 350 conidia when DNA was extracted from a fixed number of conidia. The assay amplified all isolates of P. macularis tested and had minimal cross-reactivity with other Podosphaera species when assayed with biologically relevant quantities of DNA. Standard curves generated independently in two other laboratories indicated that assay sensitivity was qualitatively similar and reproducible. All laboratories successfully detected eight unknown isolates of P. macularis and correctly discriminated Pseudoperonospora humuli and a water control. The usefulness of the assay for air sampling for late-season inoculum of P. macularis was demonstrated in field studies in 2019 and 2020. In both years, airborne populations of P. macularis in hop yards were detected consistently and increased during bloom and cone development.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo.

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Modeling 0.6 million genes for the rational design of functional cis-regulatory variants and de novo design of cis-regulatory sequences

Rational design of plant cis-regulatory DNA sequences without expert intervention or prior domain knowledge is still a daunting task. Here, we developed PhytoExpr, a deep learning framework capable of predicting both mRNA abundance and plant species using the proximal regulatory sequence as the sole input. PhytoExpr was trained over 17 species representative of major clades of the plant kingdom to enhance its generalizability. Via input perturbation, quantitative functional annotation of the input sequence was achieved at single-nucleotide resolution, revealing an abundance of predicted high-impact nucleotides in conserved noncoding sequences and transcription factor binding sites. Evaluation of maize HapMap3 single-nucleotide polymorphisms (SNPs) by PhytoExpr demonstrates an enrichment of predicted high-impact SNPs in cis-eQTL. Additionally, we provided two algorithms that harnessed the power of PhytoExpr in designing functional cis-regulatory variants, and de novo creation of species-specific cis-regulatory sequences through in silico evolution of random DNA sequences. Our model represents a general and robust approach for functional variant discovery in population genetics and rational design of regulatory sequences for genome editing and synthetic biology.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo.

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Simultaneous determination of small molecules and proteins in wastewater-based epidemiology

Wastewater-based epidemiology (WBE) aims to understand a population’s consumption habits, exposure to chemicals, and the prevalence of specific diseases or pathogens. This is achieved by the chemical or biological/genomic determination of biomarkers (e.g., excreted metabolic products), which are in urban wastewater generated by that population. WBE has been mostly linked to the determination of small molecules of human origin using liquid-chromatography mass spectrometry (LC-MS). In this Perspective, we provide a state-of-the-art and critical evaluation of further developments in the information achieved by determining small molecules as well as the most promising analytical techniques to enlarge the information obtained. By simultaneously monitoring small and large molecules we can comprehensively trace the population’s health by their consumption of prescribed pharmaceuticals and illegal drugs, as well as by the amount of excreted macromolecule biomarkers such as peptides and proteins. Moreover, species-specific protein sequences allow us to monitor animal populations reflecting farming and slaughterhouse activities (poultry, pigs…) or pest occurrences (rats). To this end, the capability of proteomic studies using high-resolution tandem mass spectrometry is highlighted and compared in the context of other advances in the broader field of high-resolution mass spectrometry (HRMS).

Молекулярные методы аутентификации молочного сырья в сыре: обзор

The modern dairy industry needs authentic raw milk in cheese production to provide high-quality product control and consumer protection. Unfortunately, unreliable cheese producers sometimes replace goat’s or sheep’s milk with cow’s milk. This process is associated with seasonality and the high cost of milk from small ruminants. The article reviews the methods molecular genetics used to detect milk adulteration in cheese making. It highlights the problems and prospects of applying molecular genetics to complex food matrices, such as cheese. Modern PCR methods prove quite efficient in identifying species-specific sequences by using mitochondrial DNA (mtDNA) as markers. Multiplex PCR, quantitative PCR (qPCR), and digital droplet PCR (ddPCR) have their advantages and shortcomings. The multiplex PCR method presupposes a simultaneous amplification of multiple DNA targets in a single tube while the qPCR method provides a quantitative profile. These methods make it possible to maintain high industry standards by ensuring the authenticity of cheeses and cheese products. In addition, molecular methods provide supply chain stability, increase consumer confidence, and advance dairy technology.

Nonribosomal Peptide Synthetase Specific Genome Amplification Using Rolling Circle Amplification for Targeted Gene Sequencing.

Next-generation sequencing has transformed the acquisition of vast amounts of genomic information, including the rapid identification of target gene sequences in metagenomic databases. However, dominant species can sometimes hinder the detection of rare bacterial species. Therefore, a highly sensitive amplification technique that can selectively amplify bacterial genomes containing target genes of interest was developed in this study. The rolling circle amplification (RCA) method can initiate amplification from a single locus using a specific single primer to amplify a specific whole genome. A mixed cell suspension was prepared using Pseudomonas fluorescens ATCC17400 (targeting nonribosomal peptide synthetase [NRPS]) and Escherichia coli (non-target), and a specific primer designed for the NRPS was used for the RCA reaction. The resulting RCA product (RCP) amplified only the Pseudomonas genome. The NRPS was successfully amplified using RCP as a template from even five cells, indicating that the single-priming RCA technique can specifically enrich the target genome using gene-specific primers. Ultimately, this specific genome RCA technique was applied to metagenomes extracted from sponge-associated bacteria, and NRPS sequences were successfully obtained from an unknown sponge-associated bacterium. Therefore, this method could be effective for accessing species-specific sequences of NRPS in unknown bacteria, including viable but non-culturable bacteria.

Exploring mosquito abundance and Plasmodium infection through nested-PCR: implications for disease surveillance and control

The Plasmodium is responsible for malaria which poses a major health threat, globally. This study is based on the estimation of the relative abundance of mosquitoes, and finding out the correlations of meteorological parameters (temperature, humidity and rainfall) with the abundance of mosquitoes. In addition, this study also focused on the use of nested PCR (species-specific nucleotide sequences of 18S rRNA genes) to explore the Plasmodium spp. in female Anopheles. In the current study, the percentage relative abundance of Culex mosquitoes was 57.65% and Anopheles 42.34% among the study areas. In addition, the highest number of mosquitoes was found in March in district Mandi Bahauddin at 21 °C (Tmax = 27, Tmin = 15) average temperature, 69% average relative humidity and 131 mm rainfall, and these climatic factors were found to affect the abundance of the mosquitoes, directly or indirectly. Molecular analysis showed that overall, 41.3% of the female Anopheles pools were positive for genus Plasmodium. Among species, the prevalence of Plasmodium (P.) vivax (78.1%) was significantly higher than P. falciparum (21.9%). This study will be helpful in the estimation of future risk of mosquito-borne diseases along with population dynamic of mosquitoes to enhance the effectiveness of vector surveillance and control programs.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

KmerDB: A database encompassing the set of genomic and proteomic sequence information for each species

The decrease in sequencing expenses has facilitated the creation of reference genomes and proteomes for an expanding array of organisms. Nevertheless, no established repository that details organism-specific genomic and proteomic sequences of specific lengths, referred to as kmers, exists to our knowledge. In this article, we present kmerDB, a database accessible through an interactive web interface that provides kmer-based information from genomic and proteomic sequences in a systematic way. kmerDB currently contains 202,340,859,107 base pairs and 19,304,903,356 amino acids, spanning 54,039 and 21,865 reference genomes and proteomes, respectively, as well as 6,905,362 and 149,305,183 genomic and proteomic species-specific sequences, termed quasi-primes. Additionally, we provide access to 5,186,757 nucleic and 214,904,089 peptide sequences absent from every genome and proteome, termed primes. kmerDB features a user-friendly interface offering various search options and filters for easy parsing and searching. The service is available at: www.kmerdb.com.

Real-time PCR method based on single-copy nuclear DNA sequences for the quantitative detection of pork adulteration in processed beef products

Meat adulteration, which adversely impacts food safety, economy, and religious beliefs, has come to be considered as a matter of global concern. The lack of a standard method that can be employed to quantitatively detect the pork content in processed meat products has led to demands for an accurate and sensitive method capable of detecting the proportion of pork in processed beef products. In this study, we aimed to develop a real-time PCR-based method to accurately and quantitatively detect pork adulteration in processed beef. We used bovine and porcine species-specific sequences to design primers and probes that specifically amplified bovine and porcine genes. Our real-time PCR method was evaluated via quantitative detection of reference samples and found to be highly accurate. Sensitivity analyses indicated that the limit of detection for beef as well as for pork was as low as five copies, while the limit of quantification for beef and for pork was 0.1%. To compare qualitative and quantitative results, beef samples cross-contaminated with pork were analyzed using both quantitative and standard qualitative methods. No undeclared meat components were found among the 25 processed meat products that were analyzed to evaluate the applicability of the quantitative method. However, the beef content in some processed beef products was very low, even approaching zero. This real-time PCR method may be reliably applied for quantitatively detect pork adulteration in processed beef products. Thus, this method of ours may help promote healthy and sustainable development of meat products.

The Chaperone NASP Contributes to de Novo Deposition of the Centromeric Histone Variant CENH3 in Arabidopsis Early Embryogenesis.

The centromere is an essential chromosome region where the kinetochore is formed to control equal chromosome distribution during cell division. The centromere-specific histone H3 variant CENH3 (also called CENP-A) is a prerequisite for the kinetochore formation. Since CENH3 evolves rapidly, associated factors, including histone chaperones mediating the deposition of CENH3 on the centromere, are thought to act through species-specific amino acid sequences. The functions and interaction networks of CENH3 and histone chaperons have been well-characterized in animals and yeasts. However, molecular mechanisms involved in recognition and deposition of CENH3 are still unclear in plants. Here, we used a swapping strategy between domains of CENH3 of Arabidopsis thaliana and the liverwort Marchantia polymorpha to identify specific regions of CENH3 involved in targeting the centromeres and interacting with the general histone H3 chaperone, nuclear autoantigenic sperm protein (NASP). CENH3's LoopN-α1 region was necessary and sufficient for the centromere targeting in cooperation with the α2 region and was involved in interaction with NASP in cooperation with αN, suggesting a species-specific CENH3 recognition. In addition, by generating an Arabidopsis nasp knock-out mutant in the background of a fully fertile GFP-CENH3/cenh3-1 line, we found that NASP was implicated for de novo CENH3 deposition after fertilization and thus for early embryo development. Our results imply that the NASP mediates the supply of CENH3 in the context of the rapidly evolving centromere identity in land plants.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

Analysis of Whole-Genome as a Novel Strategy for Animal Species Identification.

Survival crises stalk many animals, especially endangered and rare animals. Accurate species identification plays a pivotal role in animal resource conservation. In this study, we developed an animal species identification method called Analysis of whole-GEnome (AGE), which identifies species by finding species-specific sequences through bioinformatics analysis of the whole genome and subsequently recognizing these sequences using experimental technologies. To clearly demonstrate the AGE method, Cervus nippon, a well-known endangered species, and a closely related species, Cervus elaphus, were set as model species, without and with published genomes, respectively. By analyzing the whole genomes of C. nippon and C. elaphus, which were obtained through next-generation sequencing and online databases, we built specific sequence databases containing 7,670,140 and 570,981 sequences, respectively. Then, the species specificities of the sequences were confirmed experimentally using Sanger sequencing and the CRISPR-Cas12a system. Moreover, for 11 fresh animal samples and 35 commercially available products, our results were in complete agreement with those of other authoritative identification methods, demonstrating AGE's precision and potential application. Notably, AGE found a mixture in the 35 commercially available products and successfully identified it. This study broadens the horizons of species identification using the whole genome and sheds light on the potential of AGE for conserving animal resources.

Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966).

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.