Pairwise Sequence Identity Research Articles

Genomic Sequence of the Threespine Stickleback Iridovirus (TSIV) from Wild Gasterosteus aculeatus in Stormy Lake, Alaska.

The threespine stickleback iridovirus (TSIV), a double-stranded DNA virus, was the first megalocytivirus detected in wild North American fishes. We report a second occurrence of TSIV in threespine stickleback (Gasterosteus aculeatus) from Stormy Lake, Alaska, and assemble a nearly complete genome of TSIV. The 115-kilobase TSIV genome contains 94 open reading frames (ORFs), with 91 that share homology with other known iridoviruses. We identify three ORFs that likely originate from recent lateral gene transfers from a eukaryotic host and one ORF with homology to B22 poxvirus proteins that likely originated from a lateral gene transfer between viruses. Phylogenetic analysis of 24 iridovirus core genes and pairwise sequence identity analysis support TSIV as a divergent sister taxon to other megalocytiviruses and a candidate for a novel species designation. Screening of stickleback collected from Stormy Lake before and after a 2012 rotenone treatment to eliminate invasive fish shows 100% positivity for TSIV in the two years before treatment (95% confidence interval: 89-100% prevalence) and 0% positivity for TSIV in 2024 after treatment (95% confidence interval: 0 to 3.7% prevalence), suggesting that the rotenone treatment and subsequent crash and reestablishment of the stickleback population is associated with loss of TSIV.

Characterization of garlic virus A, garlic virus D, and onion yellow dwarf virus infecting onion

Background: Allium cepa is one of the major vegetable crops cultivated in Saudi Arabia. This economically important crop faces several biotic stresses which adversely affect the quality and quantity of its production. Several virus groups (potyviruses, carlaviruses, tospoviruses, and allexiviruses) have been reported infecting Allium crops.Methods: During the growing season of 2021-2022, a total of 81 onion samples exhibiting virus-like symptoms were collected from two different geographical regions in Saudi Arabia. The serological technique (ELISA) was used to detect the important allexiviruses and potyviruses. RT-PCR amplification of partial genome sequence was done using degenerate primers for allexiviruses and potyviruses and the phylogenetic trees were constructed using different bioinformatic tools.Results: The results obtained from ELISA tests showed that 26% and 32% of onion samples were positive with both Garlic Virus A (GarV-A) and Onion yellow dwarf virus (OYDV) respectively. RT-PCR amplification and sequencing results showed that two allexiviruses, GarV-A, garlic virus D (GarV-D), and one Potyvirus (OYDV) were detected in both regions. Sequence data were deposited in the GenBank database with accession numbers, OQ397545, OQ397546 for GarV-A, OQ397547 for GarV-D, and OQ397548, OQ397549 for OYDV, sequentially. Phylogenetic tree analysis of these virus isolates showed making clades with closely related isolates of their respective viruses. Pairwise nucleotide sequence identity showed their similarity with GarV-A, GarV-D, and OYDV isolates reported earlier in the GenBank.Conclusion: To the best of our knowledge, these two distinct allexiviruses (GarV-A, GarV-D) and one Potyvirus (OYDV) were isolated for the first time from an onion crop in Saudi Arabia.Keywords: Allexiviruses; Potyviruses; Serological Detection; Allium cepa; RT-PCR; Sequence

Mitochondrial DNA sequence-based identification of two subterranean termite species, from Riyadh Province, Kingdom of Saudi Arabia.

Termites are economically important wood-destroying and agricultural pests. The termite fauna almost consists of 2900 described species in 286 genera worldwide. In the present study, hundreds of termite samples from 42 different locations in the Riyadh province were collected. These samples were previously used for morphometric identification and reported two subterranean termite species, Coptotermes heimi and Psammotermes hypostoma, in the family Rhinotermitidae. In the present study, these samples were analysed using DNA barcoding with the mitochondrial cytochrome c oxidase subunit 1 gene to confirm the conventional taxonomical identification on a molecular basis. The obtained COI gene sequences of all 42 termite specimens were submitted to GenBank (accession numbers: ON529959-ON529969, OP825131-OP825132, and OP890882-OP890910). Eleven of the 42 samples were thus identified as C. heimi and the remaining 31 samples as P. hypostoma, which were phylogenetically analysed. All the 11 C. heimi sequences were grouped in a single clade, indicating close relatedness. While 31 sequences of P. hypostoma constituted two clades in the phylogenetic tree. Pairwise nucleotide sequence identity and divergence analysis showed that C. heimi sequences showed high nucleotide identities of 87.6-99.5% and less divergence ranging from 0.5% to 13.6%. Similarly, sequences of P. hypostoma also showed high nucleotide identity of 78.6-100% and low divergence among them ranging from 0-10.7%. A further application, significance, and shortcomings of COI-based DNA barcoding have been discussed. DNA barcoding using the COI gene is a reliable tool to distinguish C. heimi and P. hypostoma genotypes.

A study on the development of O2O operation service platform using augmented reality

An attempt of sequence selection is made for multiple alignments of transmembrane proteins in order to detect the structural similarity in the protein families. Treatment for the sequence selection, which is based on the pairwise sequence identities, is applied to ten sequence data sets of functional group of transmembrane proteins extracted from SWISS-PROT. The treatment excludes the sequences with low sequence identity from the divergent data sets effectively. Obtained multiple alignments are evaluated using two newly developed indices. The selected sequences, which are well aligned with few inserted gaps, seem to contain enough information for extracting the structural features of transmembrane functional groups.

2024 taxonomy updatefor the family Circoviridae.

Circovirids have a circular single-stranded DNA genome packed into a small icosahedral capsid. They are classified within two genera, Circovirus and Cyclovirus, in the family Circoviridae (phylum Cressdnaviricota, class Arfiviricetes, order Cirlivirales). Over the last five years, a number of new circovirids have been identified, and, as a result, 54 new species have been created for their classification based on the previously established species demarcation criterion, namely, that viruses classified into different species share less than 80% genome-wide pairwise sequence identity. Of note, one of the newly created species includes a circovirus that was identified in human hepatocytes and suspected of causing liver damage. Furthermore, to comply with binomial species nomenclature, all new and previously recognized species have been (re)named in binomial format with a freeform epithet. Here, we provide a summary of the properties of circovirid genomes and their classification as of June 2024 (65 species in the genus Circovirus and 90 species in the genus Cyclovirus). Finally, we provide reference datasets of the nucleotide and amino acid sequences representing each of the officially recognized circovirid species to facilitate further classification of newly discovered members of the Circoviridae.

Graphical models for identifying pore-forming proteins.

Pore-forming toxins (PFTs) are proteins that form lesions in biological membranes. Better understanding of the structure and function of these proteins will be beneficial in a number of biotechnological applications, including the development of new pest control methods in agriculture. When searching for new pore formers, existing sequence homology-based methods fail to discover truly novel proteins with low sequence identity to known proteins. Search methodologies based on protein structures would help us move beyond this limitation. As the number of known structures for PFTs is very limited, it's quite challenging to identify new proteins having similar structures using computational approaches like deep learning. In this article, we therefore propose a sample-efficient graphical model, where a protein structure graph is first constructed according to consensus secondary structures. A semi-Markov conditional random fields model is then developed to perform protein sequence segmentation. We demonstrate that our method is able to distinguish structurally similar proteins even in the absence of sequence similarity (pairwise sequence identity < 0.4)-a feat not achievable by traditional approaches like HMMs. To extract proteins of interest from a genome-wide protein database for further study, we also develop an efficient framework for UniRef50 with 43 million proteins.

Computational analysis of human medium-chain dehydrogenases/reductases revealing substrate- and coenzyme-binding characteristics

The medium-chain dehydrogenase/reductase (MDR) superfamily has more than 600,000 members in UniProt as of March 2023. As the family has been growing, the proportion of functionally characterized proteins has been falling behind. The aim of this project was to investigate the binding pockets of nine different MDR protein families based on sequence conservation patterns and three-dimensional structures of members within the respective families.A search and analysis methodology was developed. Using this, a total of 2000 eukaryotic MDR sequences belonging to nine different families were identified. The pairwise sequence identities within each of the families were 80–90 % for the mammalian sequences, like the levels observed for alcohol dehydrogenase, another MDR family. Twenty conserved residues were identified in the coenzyme part of the binding site by matching structural and conservation data of all nine protein families. The conserved residues in the substrate part of the binding pocket varied between the nine MDR families, implying divergent functions for the different families. Studying each family separately made it possible to identify multiple conserved residues that are expected to be important for substrate binding or catalysis of the enzymatic reaction.By combining structural data with the conservation of the amino acid residues in each protein, important residues in the binding pocket were identified for each of the nine MDRs. The obtained results add new positions of interest for the binding and activity of the enzyme family as well as fit well to earlier published data. Three distinct types of binding pockets were identified, containing no, one, or two tyrosine residues.

Remote homolog detection places insect chemoreceptors in a cryptic protein superfamily spanning the tree of life

Many proteins exist in the so-called "twilight zone" of sequence alignment, where low pairwise sequence identity makes it difficult to determine homology and phylogeny.1,2 As protein tertiary structure is often more conserved,3 recent advances in ab initio protein folding have made structure-based identification of putative homologs feasible.4,5,6 We present a pipeline for the identification and characterization of distant homologs and apply it to 7-transmembrane-domain ion channels (7TMICs), a protein group founded by insect odorant and gustatory receptors. Previous sequence and limited structure-based searches identified putatively related proteins, mainly in other animals and plants.7,8,9,10 However, very few 7TMICs have been identified in non-animal, non-plant taxa. Moreover, these proteins' remarkable sequence dissimilarity made it uncertain whether disparate 7TMIC types (Gr/Or, Grl, GRL, DUF3537, PHTF, and GrlHz) are homologous or convergent, leaving their evolutionary history unresolved. Our pipeline identified thousands of new 7TMICs in archaea, bacteria, and unicellular eukaryotes. Using graph-based analyses and protein language models to extract family-wide signatures, we demonstrate that 7TMICs have structure and sequence similarity, supporting homology. Through sequence- and structure-based phylogenetics, we classify eukaryotic 7TMICs into two families (Class-A and Class-B), which are the result of a gene duplication predating the split(s) leading to Amorphea (animals, fungi, and allies) and Diaphoretickes (plants and allies). Our work reveals 7TMICs as a cryptic superfamily, with origins close to the evolution of cellular life. More generally, this study serves as a methodological proof of principle for the identification of extremely distant protein homologs.

A distinct strain of tomato leaf curl New Delhi virus that causes mosaic disease in ash gourd and other cucurbitaceous crops.

Ash gourd (Benincasa hispida) is a cucurbitaceous crop cultivated as an edible vegetable rich in vitamins, minerals, dietary fibers and antioxidants. In a field survey conducted in the Udumalpet region of Tamil Nadu during 2019, the incidence of mosaic disease on ash gourd crop was observed to be 75%. The DNA-A and DNA-B components of begomovirus genome have been identified as associated with this disease. Both the cloned DNA-A and DNA-B genomic components shared highest pairwise sequence identities with the isolates of tomato leaf curl New Delhi virus (ToLCNDV), a bipartite begomovirus. Recombinant analysis showed that both the components are possibly evolved through intra-species recombination between ToLCNDV isolates. Tomato leaf curl Bangladesh betasatellite (ToLCBB) is not naturally associated with this sample. The results of infectivity studies on ash gourd and other cucurbitaceous crops demonstrates the Koch's postulates, when co-inoculation of DNA-A and DNA-B of ToLCNDV was undertaken. However, the inoculation of non-cognate ToLCBB along with DNA-A and DNA-B enhances the symptom expression and reduces the time taken for symptom development. Thus, Koch's postulates were proved for these virus complexes on cucurbitaceous crops. Furthermore, an enhanced accumulation of DNA-A component was detected in the cucurbits co-inoculated with ToLCNDV and ToLCBB. This report highlights the importance of investigating the spread of these disease complexes with other cucurbitaceous crops in India.

Monitoring the Spread of Grapevine Viruses in Vineyards of Contrasting Agronomic Practices: A Metagenomic Investigation.

This study investigated the transmission of grapevine viruses, specifically grapevine red blotch virus (GRBV) and grapevine Pinot gris virus (GPGV), in vineyards in Niagara Region, Ontario, Canada. Forty sentinel vines that were confirmed free of GRBV and GPGV by both high-throughput sequencing (HTS) and endpoint polymerase chain reaction (PCR) were introduced to two vineyards (one organic and one conventional) that were heavily infected with both GRBV and GPGV. Four months post-introduction, the sentinel vines were relocated to a phytotron. The HTS results from 15 months post-introduction revealed a widespread infection of GPGV among the sentinel vines but did not detect any GRBV. The GPGV infection rate of sentinel vines in the organic vineyard (13/18) was higher than in the conventional vineyard (1/19). The possibility of an alternative viral reservoir was assessed by testing the most abundant plants in between rows (Medicago sativa, Trifolium repens, Cirsium arvense and Taraxacum officinale), perennial plants in border areas (Fraxinus americana, Ulmus americana, Rhamnus cathartica) and wild grape (unknown Vitis sp.). The HTS result showed that cover crops and perennial plants did not harbor any grapevine viruses, while 4/5 wild grapes tested positive for GPGV but not GRBV. A pairwise sequence identity analysis revealed high similarities between the GPGV isolates found in the established vines on the vineyard and the newly contracted GPGV isolates in the sentinel vines, implicating a recent transmission event. This work provides novel insights into the spread of grapevine viruses in Niagara Region and is also the first direct proof of the spread of GPGV in natural vineyard conditions in North America.

Interaction of S100A6 Protein with the Four-Helical Cytokines.

S100 is a family of over 20 structurally homologous, but functionally diverse regulatory (calcium/zinc)-binding proteins of vertebrates. The involvement of S100 proteins in numerous vital (patho)physiological processes is mediated by their interaction with various (intra/extra)cellular protein partners, including cell surface receptors. Furthermore, recent studies have revealed the ability of specific S100 proteins to modulate cell signaling via direct interaction with cytokines. Previously, we revealed the binding of ca. 71% of the four-helical cytokines via the S100P protein, due to the presence in its molecule of a cytokine-binding site overlapping with the binding site for the S100P receptor. Here, we show that another S100 protein, S100A6 (that has a pairwise sequence identity with S100P of 35%), specifically binds numerous four-helical cytokines. We have studied the affinity of the recombinant forms of 35 human four-helical cytokines from all structural families of this fold to Ca2+-loaded recombinant human S100A6, using surface plasmon resonance spectroscopy. S100A6 recognizes 26 of the cytokines from all families of this fold, with equilibrium dissociation constants from 0.3 nM to 12 µM. Overall, S100A6 interacts with ca. 73% of the four-helical cytokines studied to date, with a selectivity equivalent to that for the S100P protein, with the differences limited to the binding of interleukin-2 and oncostatin M. The molecular docking study evidences the presence in the S100A6 molecule of a cytokine-binding site, analogous to that found in S100P. The findings argue the presence in some of the promiscuous members of the S100 family of a site specific to a wide range of four-helical cytokines. This unique feature of the S100 proteins potentially allows them to modulate the activity of the numerous four-helical cytokines in the disorders accompanied by an excessive release of the cytokines.

Dynamic hybridization between two spleenworts, Asplenium incisum and Asplenium ruprechtii in Korea.

Natural hybridization between Asplenium incisum and A. ruprechtii has been observed in Northeast Asia and its allotetraploid species, A. castaneoviride, was reported. However, the hybridization process between the parental species and the origin of the allotetraploid taxon remains obscure. Additionally, the systematic affinities of the recently described hybrid A. bimixtum, considered to have originated from the hybridization of A. ruprechtii, A. trichomanes, and A. incisum, is unresolved owing to its similarity to A. castaneoviride. The goals of this study were to (1) investigate the hybridization between A. ruprechtii and A. incisum; (2) verify the origin of A. castaneoviride occurring in Korea, whether it independently arose from 2x sterile hybrids; and (3) elucidate the reliability of identifying A. bimixtum. Three genotypes, A. incisum, A. ruprechtii, and their hybrid, were identified based on the nuclear gene pgiC sequence and finally divided them into six types by ploidy levels: diploid A. incisum, A. ruprechtii, and four hybrid types (diploid A. × castaneoviride, triploid A. × castaneoviride, allotetraploid A. castaneoviride, and A. bimixtum). In the analyses of plastid DNA, all hybrids had an A. ruprechtii-type rbcL gene. In addition, the four plastomes of A. ruprechtii and the hybrids had high pairwise sequence identities greater than 98.48%. They increased up to 99.88% when a large deletion of A. x castaneoriviride (2x) collected from Buramsan populations was ignored. Notably, this large deletion was also found in triploid A. × castaneoviride and allotetraploid A. castaneoviride in the same populations. Sequence data of the nuclear and plastid genes showed that hybridization is unidirectional, and A. ruprechtii is the maternal parent. The large deletion of rpoC2-rps2 commonly found in the different ploidy hybrids of the Buramsan population suggests that the allotetraploid A. castaneoviride can be created independently from sterile hybrids. We assume that both polyploidization driving allopolyploidy and minority cytotype exclusion took place independently in the population, since A castaenoviride co-occurs with A. ruprechtii in small populations. Furthermore, it was also observed that an enlarged noncoding region in fern organelle (ENRIFO) of the plastome was found in the genus Asplenium.

Characterization, Amplification, and Phylogenetic Analysis of Gossypium herbaceum Using rbcL Molecular Marker

People have used medicinal plants for centuries to produce traditional remedies that greatly interest modern health care. One of these plants, Gossypium herbaceum or commonly called Arabian cotton, has been used in various medicinal applications. Scientists are turning to DNA barcoding, a molecular technique that identifies species using standardized DNA regions. Objective: To evaluate samples of Gossypium herbaceum and their physical properties. Methods: DNA was extracted from the plant material, and its quality and quantity were checked. Using PCR and gel electrophoresis, amplification of the RBCL gene was done. Purification of the PCR products was done for DNA sequencing. After that, all the DNA sequences were compared with the available DNA sequences in public databases. The relationship between Gossypium herbaceum and other related species was evaluated using the neighbour-joining method for phylogenetic analysis. Results: The results showed a high percentage of pairwise nucleotide sequence identity with Gossypium richmond and Gossypium hirsutum. The study demonstrated the potential of DNA barcoding using the rbcL gene as a reliable method for identifying and confirming Gossypium herbaceum. Also, this study provides valuable insight into the phylogenetic relationships of this medicinal plant species. Conclusions: The findings support the conservation and appropriate use of medicinal plants and highlight the importance of ensuring the authenticity and quality of herbal products.

Identification and Molecular Characterization of a Novel Carlavirus Infecting Chrysanthemum morifolium in China.

Chrysanthemum (Chrysanthemum morifolium) is an important ornamental and medicinal plant suffering from many viruses and viroids worldwide. In this study, a new carlavirus, tentatively named Chinese isolate of Carya illinoinensis carlavirus 1 (CiCV1-CN), was identified from chrysanthemum plants in Zhejiang Province, China. The genome sequence of CiCV1-CN was 8795 nucleotides (nt) in length, with a 68-nt 5'-untranslated region (UTR) and a 76-nt 3'-UTR, which contained six predicted open reading frames (ORFs) that encode six corresponding proteins of various sizes. Phylogenetic analyses based on full-length genome and coat protein sequences revealed that CiCV1-CN is in an evolutionary branch with chrysanthemum virus R (CVR) in the Carlavirus genus. Pairwise sequence identity analysis showed that, except for CiCV1, CiCV1-CN has the highest whole-genome sequence identity of 71.3% to CVR-X6. At the amino acid level, the highest identities of predicted proteins encoded by the ORF1, ORF2, ORF3, ORF4, ORF5, and ORF6 of CiCV1-CN were 77.1% in the CVR-X21 ORF1, 80.3% in the CVR-X13 ORF2, 74.8% in the CVR-X21 ORF3, 60.9% in the CVR-BJ ORF4, 90.2% in the CVR-X6 and CVR-TX ORF5s, and 79.4% in the CVR-X21 ORF6. Furthermore, we also found a transient expression of the cysteine-rich protein (CRP) encoded by the ORF6 of CiCV1-CN in Nicotiana benthamiana plants using a potato virus X-based vector, which can result in a downward leaf curl and hypersensitive cell death over the time course. These results demonstrated that CiCV1-CN is a pathogenic virus and C. morifolium is a natural host of CiCV1.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences.

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Nearest neighbor search on embeddings rapidly identifies distant protein relations.

Since 1992, all state-of-the-art methods for fast and sensitive identification of evolutionary, structural, and functional relations between proteins (also referred to as "homology detection") use sequences and sequence-profiles (PSSMs). Protein Language Models (pLMs) generalize sequences, possibly capturing the same constraints as PSSMs, e.g., through embeddings. Here, we explored how to use such embeddings for nearest neighbor searches to identify relations between protein pairs with diverged sequences (remote homology detection for levels of <20% pairwise sequence identity, PIDE). While this approach excelled for proteins with single domains, we demonstrated the current challenges applying this to multi-domain proteins and presented some ideas how to overcome existing limitations, in principle. We observed that sufficiently challenging data set separations were crucial to provide deeply relevant insights into the behavior of nearest neighbor search when applied to the protein embedding space, and made all our methods readily available for others.

Surveillance of avian influenza viruses in live bird markets of Shandong province from 2013 to 2019.

Avian influenza viruses (AIVs) seriously affect the poultry industry and pose a great threat to humans. Timely surveillance of AIVs is the basis for preparedness of the virus. This study reported the long-term surveillance of AIVs in the live bird market (LBM) of 16 cities in Shandong province from 2013 to 2019. A total of 29,895 samples were obtained and the overall positive rate of AIVs was 9.7%. The H9 was found to be the most predominant subtype in most of the time and contributed most to the monthly positve rate of AIVs as supported by the univariate and multivariate analysis, while H5 and H7 only circulated in some short periods. Then, the whole-genome sequences of 62 representative H9N2 viruses including one human isolate from a 7-year-old boy in were determined and they were genetically similar to each other with the median pairwise sequence identities ranging from 0.96 to 0.98 for all segments. The newly sequenced viruses were most similar to viruses isolated in chickens in mainland China, especially the provinces in Eastern China. Phylogenetic analysis showed that these newly sequenced H9N2 viruses belonged to the same clade for all segments except PB1. Nearly all of these viruses belonged to the G57 genotype which has dominated in China since 2010. Finally, several molecular markers associated with human adaptation, mammalian virulence, and drug resistance were identified in the newly sequenced H9N2 viruses. Overall, the study deepens our understanding of the epidemic and evolution of AIVs and provides a basis for effective control of AIVs in China.

Complete mitochondrial genome of the branching octocoral Paramuricea grayi (Johnson, 1861), phylogenetic relationships and divergence analysis

The Gray’s sea fan, Paramuricea grayi (Johnson, 1861), typically inhabits deep littoral and circalittoral habitats of the eastern temperate and tropical Atlantic Ocean. Along the Iberian Peninsula, where P. grayi is a dominant constituent of circalittoral coral gardens, two segregating lineages (yellow and purple morphotypes) were recently identified using single-copy nuclear orthologues. The mitochondrial genomes of 9 P. grayi individuals covering both color morphotypes were assembled from RNA-seq data, using samples collected at three sites in southern (Sagres and Tavira) and western (Cape Espichel) Portugal. The complete circular mitogenome is 18,668 bp in length, has an A + T-rich base composition (62.5%) and contains the 17 genes typically found in Octocorallia: 14 protein-coding genes (atp6, atp8, cob, cox1-3, mt-mutS, nad1-6, and nad4L), the small and large subunit rRNAs (rns and rnl), and one transfer RNA (trnM). The mitogenomes were nearly identical for all specimens, though we identified a noteworthy polymorphism (two SNPs 9 bp apart) in the mt-mutS of one purple individual that is shared with the sister species P. clavata. The mitogenomes of the two species have a pairwise sequence identity of 99.0%, with nad6 and mt-mutS having the highest rates of non-synonymous substitutions.

Revisiting Papillomavirus Taxonomy: A Proposal for Updating the Current Classification in Line with Evolutionary Evidence.

Papillomaviruses infect a wide array of animal hosts and are responsible for roughly 5% of all human cancers. Comparative genomics between different virus types belonging to specific taxonomic groupings (e.g., species, and genera) has the potential to illuminate physiological differences between viruses with different biological outcomes. Likewise, extrapolation of features between related viruses can be very powerful but requires a solid foundation supporting the evolutionary relationships between viruses. The current papillomavirus classification system is based on pairwise sequence identity. However, with the advent of metagenomics as facilitated by high-throughput sequencing and molecular tools of enriching circular DNA molecules using rolling circle amplification, there has been a dramatic increase in the described diversity of this viral family. Not surprisingly, this resulted in a dramatic increase in absolute number of viral types (i.e., sequences sharing <90% L1 gene pairwise identity). Many of these novel viruses are the sole member of a novel species within a novel genus (i.e., singletons), highlighting that we have only scratched the surface of papillomavirus diversity. I will discuss how this increase in observed sequence diversity complicates papillomavirus classification. I will propose a potential solution to these issues by explicitly basing the species and genera classification on the evolutionary history of these viruses based on the core viral proteins (E1, E2, and L1) of papillomaviruses. This strategy means that it is possible that a virus identified as the closest neighbor based on the E1, E2, L1 phylogenetic tree, is not the closest neighbor based on L1 nucleotide identity. In this case, I propose that a virus would be considered a novel type if it shares less than 90% identity with its closest neighbors in the E1, E2, L1 phylogenetic tree.

Genetic diversity of imported PRRSV-2 strains, 2005-2020, Hungary.

Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) remains sporadic in Europe. In this study, we investigated the molecular epidemiology of PRRSV-2 infections encompassing 15 years in Hungary. Partial (423 bp long) ORF5 sequences (n = 44) from 20 Hungarian pig herds were analyzed. The study strains fell into two genetic lineages, L1 and L5, being L5 strains more prevalent (88.6 vs. 11.4%). Pairwise sequence identities within Hungarian representative PRRSV-2 strains ranged between 84.7 to 100% (nucleotide, nt) and 85 to 100% (amino acid, aa). When compared with reference strains, identity values fell between 87 and 100% (L1, nt 87–91%, aa 87–93%, reference strain IAF-exp91; L5, nt 87–100%, aa 88–100%, reference strain Ingelvac MLV). Epidemiologic examination implied that the majority of L5 strains were imported repeatedly from other European countries where Ingelvac MLV was approved for routine use. The emergence of L1 strains was thought to be associated with a single introduction and subsequent dissemination between pig farms of a large integrator. Results presented here contribute to a better understanding of the epizootiology of PRRSV-2 infections and shed light on the genetic diversity of viral strains in non-endemic countries.