Genome Alignment Research Articles

An updated reference genome sequence and annotation reveals gene losses and gains underlying naked mole-rat biology.

The naked mole-rat (NMR; Heterocephalus glaber ) is a eusocial subterranean rodent with a highly unusual set of physiological traits that has attracted great interest amongst the scientific community. However, the genetic basis of most of these traits has not been elucidated. To facilitate our understanding of the molecular mechanisms underlying NMR physiology and behaviour, we generated a long-read chromosomal-level genome assembly of the NMR. This genome was subsequently annotated and incorporated into multiple whole genome alignments in the Ensembl database. Our long-read assembly identified thousands of repeats and genes that were previously unassembled in the NMR and improved the results of routinely used short-read sequencing-based experiments such as RNA-seq, snRNA-seq, and ATAC-seq. We identified several spermatozoa related gene losses that may underlie the unique degenerative sperm phenotype in NMRs ( IRGC , FSCB , AKAP3 , MROH2B , CATSPER1 , DCDC2C , ATP1A4 , TEKT5 , and ZAN ), and an additional gene loss related to the established NK-cell absence in NMRs (PILRB). We resolved several tandem duplications in genes related to pathways underlying unique NMR adaptations including hypoxia tolerance, oxidative stress, and nervous system protection ( TINF2 , TCP1 , KYAT1 ). Lastly, we describe our ongoing efforts to generate a reference telomere-to-telomere assembly in the NMR which includes the resolution of complex gene families. This new reference genome should accelerate the discovery of the genetic underpinnings of NMR physiology and adaptation.

An efficient deep learning method for amino acid substitution model selection.

Amino acid substitution models play an important role in studying the evolutionary relationships among species from protein sequences. The amino acid substitution model consists of a large number of parameters; therefore, it is estimated from hundreds or thousands of alignments. Both general models and clade-specific models have been estimated and widely used in phylogenetic analyses. The maximum likelihood method is normally used to select the best fit model for a specific protein alignment under the study. A number of studies have discussed theoretical concerns as well as computational burden of the maximum likelihood methods in model selection. Recently, machine learning methods have been proposed for selecting nucleotide models. In this paper, we propose a method to measure substitution rates among amino acids (called summary statistics) from protein alignments to efficiently train a deep learning network of so-called ModelDetector for detecting amino acid substitution models. The ModelDetector network was trained from 2,246,400 alignments on a computer with 8 cores (without GPU) in about 3.3 hours. Experiments on simulation data showed that the accuracy of the ModelDetector was comparable with that of the maximum likelihood method ModelFinder. It was orders of magnitude faster than the maximum likelihood method in inferring amino acid substitution models and able to analyze genome alignments with millions of sites in minutes. The results indicate that the deep learning network can play as a promising tool for amino acid substitution model selection.

Genetic Variation Study of Several Romanian Pepper (Capsicum annuum L.) Varieties Revealed by Molecular Markers and Whole Genome Resequencing.

Numerous varieties of Capsicum annuum L. with multiple valuable traits, such as adaptation to biotic and abiotic stress factors, can be found in south-east Romania, well known for vegetable cultivation and an important area of biodiversity conservation. To obtain useful information about sustainable agriculture, management, and conservation of local pepper varieties, we analyzed the genetic diversity and conducted deep molecular characterization using whole genome resequencing (WGS) for variant/mutation detection. The pepper varieties used in the present study were registered by VRDS in the ISTIS catalog between 1974 and 2019 and maintained in conservative selection; however, no studies have been published yet using WGS analysis in order to characterize this specific germplasm. The genome sequences, annotation, and alignments provided in this study offer essential resources for genomic research as well as for future breeding efforts using the C. annuum local varieties.

Synoptic Variation Drives Genetic Diversity and Transmission Mode of Airborne DNA Viruses in Urban Space.

Airborne viruses are ubiquitous and play critical roles in maintaining ecosystem balance, however, they remain unexplored. Here, it is aimed to demonstrate that highly diverse airborne viromes carry out specific metabolic functions and use different transmission modes under different air quality conditions. A total of 263.5-Gb data are collected from 13 air samples for viral metagenomic analysis. After assembly and curation, a total of 12 484 viral contigs (1.5-184.2 kb) are assigned to 221 genus-level clades belonging to 47 families, 19 orders, and 15 classes. The composition of viral communities is influenced by weather conditions, with the main biomarker being Caudoviricetes. The most dominant viruses in these air samples belong to the dsDNA Caudoviricetes (54.0%) and ssDNA Repensiviricetes (31.2%) classes. Twelve novel candidate viruses are identified at the order/family/genus levels by alignment of complete genomes and core genes. Notably, Caudoviricetes are highly prevalent in cloudy and smoggy air, whereas Repensiviricetesare highly dominant in sunny and rainy air. Diverse auxiliary metabolic genes of airborne viruses are mainly involved in deoxynucleotide synthesis, implying their unique roles in atmosphere ecosystem. These findings deepen the understanding of the meteorological impacts on viral composition, transmission mode, and ecological roles in the air that we breathe.

Tn4661-mediated transfer of bla CTX-M-15 from Klebsiella michiganensis to an outbreak clone of Pseudomonas aeruginosa.

Carriage of CTX-M-type extended-spectrum β-lactamase (ESBL) is rare in Pseudomonas aeruginosa. During routine surveillance of an endemic ST-621 P. aeruginosa at a large hospital, isolate MRSN 100690 carrying bla CTX-M-15 was cultured from a patient (P2). This was the first detection of this ESBL in the endemic ST-621 lineage. All 1 488 bacterial isolates collected from the same facility in the 12 months prior to the incidence of 100 690 were screened for the presence of bla CTX-M-15. A set of 183 isolates was identified, in which corresponding patient metadata was evaluated for spatiotemporal overlaps with P2. The resulting three isolates, along with 100 690, were long-read sequenced using the Oxford Nanopore MinION platform to determine a potential donor of bla CTX-M-15. The screen revealed a single Klebsiella michiganensis isolate, MRSN 895358, which carried an IncA/C2 plasmid harbouring bla CTX-M-15. Notably, the patient harbouring 895358, P1, occupied the same hospital room as P2 9 months prior. Genomic alignment revealed that both isolates shared an identical 80.8 kb region containing the IncA/C2 plasmid replicon and bla CTX-M-15. This region was plasmid bound in 895 358, but chromosomally bound in 100 690 due to Tn4661-mediated transposition. ESBL bla CTX-M-15 was acquired and subsequently integrated into the chromosome of a ST-621 P. aeruginosa, likely initiated by plasmid transfer from a K. michiganensis strain.

Identification of two effective sex-specific DNA markers in silver arowana (Osteoglossum bicirrhosum)

The silver arowana (Osteoglossum bicirrhosum), often referred to as a living fossil, presents significant challenges in sex identification based on morphological characteristics, with its sex determination mechanism remaining poorly understood. This issue poses considerable obstacles to its captive breeding efforts. Given the species' substantial value for both consumption and ornamental purposes, it is a critical candidate for aquaculture. Through genome sequencing, read mapping, and alignment analysis of five male and five female silver arowanas, along with one mixed pool of males, we successfully identified two sex-specific SNPs and confirmed a ZW model of sex determination in this ancient species. These two SNPs demonstrated remarkable accuracy, achieving 100 % success in genetic sex identification, as validated by Sanger sequencing of PCR products from 135 individuals. Both sex-specific SNPs are located approximately 26 kb upstream of the foxl2 gene, indicating a potentially pivotal role for the foxl2 gene in the sex-determination process of the silver arowana. Consequently, we have developed precise molecular markers for sex identification in the silver arowana, significantly enhancing reproductive efficiency and promoting further industrial development. This achievement lays a foundational basis for investigating the sex-determination mechanism of the silver arowana and other osteoglossid fish, thereby advancing aquaculture practices within this species.

Molecular Characterisation of US6 (gD) and US8 (gE) Glycoprotein Genes Located on a Unique Short Region of Bovine Alphaherpesvirus-1 (BoAHV-1) Genome

Background: Bovine alphaherpesvirus type 1 (BoAHV-1) is a significant pathogen causing respiratory diseases and reproductive issues in cattle and buffaloes, resulting in substantial economic losses globally. This study provides a comprehensive analysis of the genetic characteristics of BoAHV-1 isolates from India, focusing on glycoprotein genes US6 (gD) and US8(gE). Methods: Ten archival isolates that were stored at the ICAR-NIVEDI BSL-2 facility were chosen, revived on MDBK cell lines and subjected for molecular characterization through PCR amplification of partial US6 (gD) and US8 (gE) gene. Genome alignment and phylogenetic analysis were performed by GeneTool and Mega 11. Result: Phylogenetic analysis of both genes revealed two cluster formation, BoAHV-1.1, BoAHV-1.2 in one cluster and BoAHV-5, Bubaline alphaherpesvirus (BuAHV-1) in another cluster. All ten isolates displayed close relation to BoAH 1.1. Seven isolates from each phylogenetic analysis showed 100% genome identity to vaccine strains from USA- Bovishield Gold MLV, Pyramid IBR MLV, Switzerland Cooper strain, ATCC Los Angeles strain, previously reported BoAHV 1.1 strains of India and Egypt. Genome identity of BoAHV-1.1 against BoAHV-1.2 was 97.70 to 98.70% for gD gene and 99.0 to 99.70% for gE gene respectively. The analysis highlights the higher genome identity among BoAHV-1 subtypes in the gE gene, suggesting its potential for developing more sensitive immunoassays for seroepidemiology. The developed assay could also be exploited for detecting BoAHV-5 infections. The findings contribute to the knowledge base essential for development of companion ELISA for implementing Differentiating Infected from Vaccinated Animals (DIVA) vaccination strategies.

Prokrustean Graph: A substring index for rapid k-mer size analysis.

Despite the widespread adoption of -mer-based methods in bioinformatics, understanding the influence of -mer sizes remains a persistent challenge. Selecting an optimal -mer size or employing multiple -mer sizes is often arbitrary, application-specific, and fraught with computational complexities. Typically, the influence of -mer size is obscured by the outputs of complex bioinformatics tasks, such as genome analysis, comparison, assembly, alignment, and error correction. However, it is frequently overlooked that every method is built above a well-defined -mer-based object like Jaccard Similarity, de Bruijn graphs, -mer spectra, and Bray-Curtis Dissimilarity. Despite these objects offering a clearer perspective on the role of -mer sizes, the dynamics of -mer-based objects with respect to -mer sizes remain surprisingly elusive. This paper introduces a computational framework that generalizes the transition of -mer-based objects across -mer sizes, utilizing a novel substring index, the Prokrustean graph. The primary contribution of this framework is to compute quantities associated with -mer-based objects for all -mer sizes, where the computational complexity depends solely on the number of maximal repeats and is independent of the range of -mer sizes. For example, counting vertices of compacted de Bruijn graphs for can be accomplished in mere seconds with our substring index constructed on a gigabase-sized read set. Additionally, we derive a space-efficient algorithm to extract the Prokrustean graph from the Burrows-Wheeler Transform. It becomes evident that modern substring indices, mostly based on longest common prefixes of suffix arrays, inherently face difficulties at exploring varying -mer sizes due to their limitations at grouping co-occurring substrings. We have implemented four applications that utilize quantities critical in modern pangenomics and metagenomics. The code for these applications and the construction algorithm is available at https://github.com/KoslickiLab/prokrustean.

KPC variants conferring resistance to ceftazidime-avibactam in Pseudomonas aeruginosa strains

BackgroundThis study aimed to characterize three KPC variants (KPC-33, KPC-100, and KPC-201) obtained from a clinical isolate of Pseudomonas aeruginosa (#700), along with two induced strains C109 and C108. MethodsGenomic DNAs of #700 (ST235), C109 (ST463), and C108 (ST1076) were sequenced using Illumina and Oxford Nanopore technologies. The transferability and stability of the plasmid was assessed through conjugation experiments and plasmid stability experiments, respectively. Minimum inhibitory concentrations of bacterial strains were determined using broth microdilution methods. In vitro induction was performed using ceftazidime-avibactam (CZA) at concentrations of 6/4 µg/ml. Linear genomic alignments were visualized using Easyfig, and protein structure modeling of the novel KPC variant (KPC-201) was conducted using PyMol. ResultsThe plasmids carrying the KPC variants in the three CZA-resistant strains (C109, C108, and #700) had sizes of 39,251 bp (KPC-100), 394,978 bp (KPC-201), and 48,994 bp (KPC-33). All three plasmids belonged to the IncP-like incompatibility (Inc) groups, and the plasmid exhibited relatively high plasmid stability, KPC-33 and KPC-201-harboring plasmids were successfully transferred to the recipient strain P. aeruginosa PAO1rifR. The genetic environments of the three blaKPC genes differed from each other. The mobile elements of the three blaKPC genes were as follows, TnAS1-IS26-ΔISKpn27-blaKPC-33-ISKpn6-IS26, IS6-ΔISKpn27-blaKPC-100-ISKpn6-IS26-Tn3-IS26, and IS6100-ISKpn27-blaKPC-201-ISKpn6-TnAS1. Notably, the length of ΔISKpn27 upstream of the blaKPC-33 and blaKPC-100 genes were remarkably short, measuring 114 bp and 56 bp, respectively, deviating significantly from typical lengths associated with ISKpn27 elements. Moreover, the novel KPC variant, KPC-201, featured a deletion of amino acids LDR at positions 161–163 in KPC-3, resulting in a looser pocket structure contributing to its avibactam resistance. ConclusionsKPC-201, identified as a novel KPC variant, exhibits resistance to CZA. The presence of multiple mobile elements surrounding the blaKPC-variant genes on stable plasmids is concerning. Urgent preventive measures are crucial to curb its dissemination in clinical settings.

Long-term persistence of diverse clones shapes the transmission landscape of invasive Listeria monocytogenes

BackgroundThe foodborne bacterium Listeria monocytogenes (Lm) causes a range of diseases, from mild gastroenteritis to invasive infections that have high fatality rate in vulnerable individuals. Understanding the population genomic structure of invasive Lm is critical to informing public health interventions and infection control policies that will be most effective especially in local and regional communities.MethodsWe sequenced the whole draft genomes of 936 Lm isolates from human clinical samples obtained in a two-decade active surveillance program across 58 counties in New York State, USA. Samples came mostly from blood and cerebrospinal fluid. We characterized the phylogenetic relationships, population structure, antimicrobial resistance genes, virulence genes, and mobile genetic elements.ResultsThe population is genetically heterogenous, consisting of lineages I–IV, 89 clonal complexes, 200 sequence types, and six known serogroups. In addition to intrinsic antimicrobial resistance genes (fosX, lin, norB, and sul), other resistance genes tetM, tetS, ermG, msrD, and mefA were sparsely distributed in the population. Within each lineage, we identified clusters of isolates with ≤ 20 single nucleotide polymorphisms in the core genome alignment. These clusters may represent isolates that share a most recent common ancestor, e.g., they are derived from the same contamination source or demonstrate evidence of transmission or outbreak. We identified 38 epidemiologically linked clusters of isolates, confirming eight previously reported disease outbreaks and the discovery of cryptic outbreaks and undetected chains of transmission, even in the rarely reported Lm lineage III (ST3171). The presence of animal-associated lineages III and IV may suggest a possible spillover of animal-restricted strains to humans. Many transmissible clones persisted over several years and traversed distant sites across the state.ConclusionsOur findings revealed the bacterial determinants of invasive listeriosis, driven mainly by the diversity of locally circulating lineages, intrinsic and mobile antimicrobial resistance and virulence genes, and persistence across geographical and temporal scales. Our findings will inform public health efforts to reduce the burden of invasive listeriosis, including the design of food safety measures, source traceback, and outbreak detection.

Comparative plastome analyses and evolutionary relationships of 25 East Asian species within the medicinal plant genus Scrophularia (Scrophulariaceae)

BackgroudScrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear.MethodsIn this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes.Results and discussionThe conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia.

KegAlign: Optimizing pairwise alignments with diagonal partitioning.

Our ability to generate sequencing data and assemble it into high quality complete genomes has rapidly advanced in recent years. These data promise to advance our understanding of organismal biology and answer longstanding evolutionary questions. Multiple genome alignment is a key tool in this quest. It is also the area which is lagging: today we can generate genomes faster than we can construct and update multiple alignments containing them. The bottleneck is in considerable computational time required to generate accurate pairwise alignments between divergent genomes, an unavoidable precursor to multiple alignments. This step is typically performed with lastZ, a very sensitive and yet equally slow tool. Here we describe an optimized GPU-enabled pairwise aligner KegAlign. It incorporates a new parallelization strategy, diagonal partitioning, with the latest features of modern GPUs. With KegAlign a typical human/mouse alignment can be computed in under 6 hours on a machine containing a single NVidia A100 GPU and 80 CPU cores without the need for any pre-partitioning of input sequences: a ~150× improvement over lastZ. While other pairwise aligners can complete this task in a fraction of that time, none achieves the sensitivity of KegAlign's main alignment engine, lastZ, and thus may not be suitable for comparing divergent genomes. In addition to providing the source code and a Conda package for KegAlign we also provide a Galaxy workflow that can be readily used by anyone.

Study of Endogenous Viruses in the Strawberry Plants.

Endogenous viral elements (EVEs) have been reported to exist widely in the genomes of eukaryotic organisms, and they are closely associated with the growth, development, genetics, adaptation, and evolution of their hosts. In this study, two methods-homologous sequence search and genome alignment-were used to explore the endogenous viral sequences in the genomes of Fragaria species. Results revealed abundant endogenous pararetroviruses (EPRVs) in the genomes of Fragaria species, including 786 sequences belonging to five known taxa such as Caulimovirus and other unclassified taxa. Differences were observed in the detected EPRVs between the two methods, with the homologous sequence search having a greater number of EPRVs. On the contrary, genome alignment identified various types and sources of virus-like sequences. Furthermore, through genome alignment, a 267-bp sequence with 95% similarity to the gene encoding the aphid-transmitted protein of Strawberry vein banding virus (Caulimovirus venafragariae) was discovered in the F. chiloensis genome, which was likely a recent insertion. In addition, the statistical analysis of the genome alignment results indicated a remarkably higher abundance of virus-like sequences in the genomes of polyploid strawberries compared with diploid ones. Moreover, the differences in virus-like sequences were observed between the genomes of Fragaria species and those of their close relatives. This study enriched the diversity of viruses that infect strawberries, and laid a theoretical foundation for further research on the origin of endogenous viruses in the strawberry genome, host-virus interactions, adaptation, evolution, and their functions.

Microbial Forensics: Comparison of MLVA Results According to NGS Methods, and Forensic DNA Analysis Using MLVA

Microbial forensics is a scientific discipline for analyzing evidence related to biological crimes by identifying the origin of microorganisms. Multiple locus variable number tandem repeat analysis(MLVA) is one of the microbiological analysis methods used to specify subtypes within a species based on the number of tandem repeat in the genome, and advances in next generation sequencing(NGS) technology have enabled <i>in silico</i> anlysis of full-length whole genome sequences. In this paper, we analyzed unknown samples provided by Robert Koch Institute(RKI) through The United Nations Secretary-General’s Mechanism(UNSGM)’s external quality assessment exercise(EQAE) project, which we officially participated in 2023. We confirmed that the 3 unknown samples were <i>B. anthracis</i> through nucleic acid isolation and genetic sequence analysis studies. MLVA results on 32 loci of <i>B. anthracis</i> were analysed by using genome sequences obtained from NGS(NextSeq and MinION) and Sanger sequencing. The MLVA typing using short-reads based NGS platform(NextSeq) showed a high probability of causing assembly error when a size of the tandem repeats was grater than 200 bp, while long-reads based NGS platform(MinION) showed higher accuracy than NextSeq, although insertion and deletion was observed. We also showed hybrid assembly can correct most indel error caused by MinION. Based on the MLVA results, genetic identification was performed compared to the 2,975 published MLVA databases of <i>B. anthracis</i>, and MLVA results of 10 strains were identical with 3 unkonwn samples. As a result of whole genome alignment of the 10 strains and 3 unknown samples, all samples were identified as <i>B. anthracis</i> strain A4564 which is associated with injectional anthrax isolates in heroin users.

Comparative plastome and mitogenome analyses indicate that the marine prasinophyte green algae Pycnococcus provasolii and Pseudoscourfieldia marina (Pseudoscourfieldiophyceae class nov., Chlorophyta) represent morphotypes of the same species.

The marine prasinophyte green algae Pycnococcus provasolii and Pseudoscourfieldia marina represent the only extant genera and known species of the Pycnococcaceae. However, their taxonomic status needs to be reassessed, owing to the very close relationship inferred from previous sequence comparisons of individual genes. Although Py. provasolii and Ps. marina are morphologically different, their plastid rbcL and nuclear small subunit rRNA genes were observed to be nearly or entirely identical in sequence, thus leading to the hypothesis that they represent distinct growth forms or alternate life-cycle stages of the same organism. To evaluate this hypothesis, we used organelle genomes as molecular markers. The plastome and mitogenome of Ps. marina UIO 007 were sequenced and compared with those available for two isolates of Py. provasolii (CCMP 1203 and CCAP 190/2). The Ps. marina organelle genomes proved to be almost identical in size and had the same gene content and gene order as their Py. provasolii counterparts. Single nucleotide substitutions and insertions/deletions were localized using genome-scale sequence alignments. Over 99.70% sequence identities were observed in all pairwise comparisons of plastomes and mitogenomes. Alignments of both organelle genomes revealed that Ps. marina UIO 007 is closer to Py. provasolii CCAP 190/2 than are the two Py. provasolii strains to one another. Therefore, our results are not consistent with the placement of Ps. marina and Py. provasolii strains into distinct genera. We propose a taxonomic revision of the Pycnococcaceae and the erection of a new class of Chlorophyta, the Pseudoscourfieldiophyceae.

Geographic Distribution of Rabies Virus and Genomic Sequence Alignment of Wild and Vaccine Strains, Kenya.

Rabies, a viral disease that causes lethal encephalitis, kills ≈59,000 persons worldwide annually, despite availability of effective countermeasures. Rabies is endemic in Kenya and is mainly transmitted to humans through bites from rabid domestic dogs. We analyzed 164 brain stems collected from rabid animals in western and eastern Kenya and evaluated the phylogenetic relationships of rabies virus (RABV) from the 2 regions. We also analyzed RABV genomes for potential amino acid changes in the vaccine antigenic sites of nucleoprotein and glycoprotein compared with RABV vaccine strains commonly used in Kenya. We found that RABV genomes from eastern Kenya overwhelmingly clustered with the Africa-1b subclade and RABV from western Kenya clustered with Africa-1a. We noted minimal amino acid variances between the wild and vaccine virus strains. These data confirm minimal viral migration between the 2 regions and that rabies endemicity is the result of limited vaccine coverage rather than limited efficacy.

Analysis of the recombination and evolution of the new type mutant pseudorabies virus XJ5 in China

Pseudorabies have caused enormous economic losses in China’s pig industry and have recurred on many large pig farms since late 2011. The disease is caused by highly pathogenic, antigenic variant pseudorabies virus (vPRV) strains. Our laboratory isolated a pseudorabies virus in 2015 and named it XJ5. The pathogenic ability of this mutant strain was much stronger than that of the original isolate. After we sequenced its whole genome (GenBank accession number: OP512542), we found that its overall structure was not greatly changed compared with that of the previous strain Ea (KX423960.1). The whole genome alignment showed that XJ5 had a strong genetic relationship with the strains isolated in China after 2012 reported in GenBank. Based on the isolation time of XJ5 and the mutation and recombination analysis of programs, we found that the whole genome homology of XJ5 and other strains with Chinese isolates was greater than 95%, while the homology with strains outside Asia was less than 94%, which indicated that there may be some recombination and mutation patterns. We found that virulent PRV isolates emerged successively in China in 2011 and formed two different evolutionary clades from foreign isolates. At the same time, this may be due to improper immunization and the presence of wild strains in the field, and recent reports have confirmed that Bartha vaccine strains recombine with wild strains to obtain new pathogenic strains. We performed genetic evolution analysis of XJ5 isolated and sequenced in our laboratory to trace its possible mutations and recombination. We found that XJ5 may be the result of natural mutation of a virus in a branch of mutant strains widely existing in China.

Deciphering the complete chloroplast genome sequence of Meconopsis torquata Prain: Insights into genome structure, comparative analysis and phylogenetic relationship

In the present study, we have characterized the complete chloroplast (Cp) genome of Meconopsis torquata Prain (family Papaveraceae), revealing the plastome size of 153,290 bp, and a GC content of 38.72 %. The cp genome features the typical circular quadripartite structure found in flowering plants, including a pair of inverted repeat regions (25,816 bp), isolated by a small single-copy region (17,740 bp) and a large single-copy (83,918 bp). Genome annotation revealed 132 genes: 87 protein-coding genes, 37 tRNAs and eight rRNAs. This comparative study demonstrated that the genome structure, gene number and GC ratio are consistent with several other cp genomes of Meconopsis and Papaver genera. A total of 120 SSRs were detected in the plastome, the majority (111) of which were mononucleotide repeats. Among the longer repeats, palindromic sequences were most common, followed by forward, reverse, and complement repeats. The whole genome alignment revealed the conserved nature of the inverted repeat region over single-copy zones. Nucleotide diversity unveiled hypervariable sites (ycf1, rps16, accD, atpB and psbD) in both the small and large single-copy regions, which could be useful for designing molecular markers for taxonomic identification. Phylogenetic analysis revealed a close alliance of M. torquata with other Meconopsis species, such as M. pinnatifolia and M. paniculata, with strong bootstrap support. Molecular dating suggests that M. torquata originated during the Tortonian age of the Miocene epoch of the Cenozoic era. These findings provide valuable insights for biological research, especially in understanding the genetic and evolutionary divergence within the Papaveraceae family.

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

Prediction of cancer drug sensitivity based on genomic feature distribution alignment and drug structure information

In recent years, precision medicine has demonstrated wide applications in cancer therapy, and the focus of precision medicine lies in accurately predicting the responses of different patients to drug treatment. We propose a model for predicting cancer drug sensitivity based on genomic feature distribution alignment and drug structure information. This model initially aligns the genomic features from cell lines with those from patients and removes noise from gene expression data. Subsequently, it integrates drug structure features and employs multi-task learning to predict the drug sensitivity of patients. The experimental results on the genomics of drug sensitivity in cancer (GDSC) dataset indicates that this method achieved a reduced mean square error of 0.905 2, an increased correlation coefficient of 0.875 4, and an enhanced accuracy rate of 0.836 0 which significantly outperformed the recently published methods. On the cancer genome atlas (TCGA) dataset, this method demonstrates an improved average recall rate of 0.571 4 and an increased F1-score of 0.658 0 in predicting drug sensitivity, exhibiting excellent generalization performance. The result demonstrates the potential of this method to assist in the selection of clinical treatment plans in the future.