Full-length Genome Sequences Research Articles

Insights into genomic sequence diversity of the SAG surface antigen superfamily in geographically diverse Eimeria tenella isolates

Eimeria tenella is among the protozoan parasites that cause the infectious disease coccidiosis in chickens, incurring huge economic losses to the global poultry industry. Surface antigens (EtSAGs) involved in host-parasite interaction are potential targets for control strategies. However, the occurrence of genetic diversity for EtSAGs in field populations is unknown, as is the risk of such diversity to the efficacy of EtSAG-based control approaches. Here, the extent of EtSAG genetic diversity and its implications on protein structure and function is assessed. Eighty-seven full-length EtSAG genomic sequences were identified from E. tenella genome assemblies of isolates sampled from continents including North America (United States), Europe (United Kingdom), Asia (Malaysia and Japan) and Africa (Nigeria). Limited diversity was observed in the EtSAG sequences. However, distinctive patterns of polymorphism were identified between EtSAG subfamilies, suggesting functional differences among these antigen families. Polymorphisms were sparsely distributed across isolates, with a small number of variants exclusive to specific geographical regions. These findings enhance our understanding of EtSAGs, particularly in elucidating functional differences among the antigens that could inform the development of more effective and long-lasting anticoccidial control strategies.

Complete genome sequence of Kozakia baliensis LTH 7215, a newly isolated strain from grape must.

We report the complete genome sequence of the novel Kozakia baliensis strain LTH 7215, isolated from grape must. This strain exhibits a mucoid phenotype on solid medium, associated with exopolysaccharide production. The high-quality and full-length genome sequence was obtained using PacBio Sequencing technology, refined with Illumina sequencing.

The Challenging Life of Mutators: How Pseudomonas aeruginosa Survives between Persistence and Evolution in Cystic Fibrosis Lung.

Cystic fibrosis (CF) is a life-threatening genetic disease characterised by chronic lung infections sustained by opportunistic pathogens such as Pseudomonas aeruginosa. During the chronic long-lasting lung infections, P. aeruginosa adapts to the host environment. Hypermutability, mainly due to defects in the DNA repair system, resulting in an increased spontaneous mutation rate, represents a way to boost the rapid adaptation frequently encountered in CF P. aeruginosa isolates. We selected 609 isolates from 51 patients with CF chronically colonised by P. aeruginosa to study, by full-length genome sequencing, the longitudinal evolution of the bacterium. We recovered at least one hypermutable (mutator) isolate in 57% of patients. By combining genomic information and phenotypic analyses, we followed the evolutionary pathways of the P. aeruginosa mutator strains, identifying their contribution to multi-drug resistance and the emergence of new sub-lineages. By implementing patient clinical data, we observed that mutators preferentially follow a specific evolutionary trajectory in patients with a negative clinical outcome and that maintenance antibiotic polytherapy, based on alternating molecules, apparently reduces the occurrence of hypermutability. Finally, we draw attention to the possibility that modulator-induced changes in the pulmonary environment may be associated with the onset of hypermutability.

Genotypic and phenotypic characterisation of respiratory syncytial virus after nirsevimab breakthrough infections: a large, multicentre, observational, real-world study

Nirsevimab, a long-acting monoclonal antibody, has been approved for the prevention of respiratory syncytial virus (RSV) infection in infants. In France, more than 210 000 single doses were administered in infants younger than 1 year during the 2023-24 season. In this context, the selection and spread of escape variants might be a concern. Here, we aimed to characterise RSV associated with breakthrough infection. We did a multicentre, national, observational study in France during the 2023-24 RSV season in RSV-infected infants (aged <1 year) who either received or did not receive a dose of nirsevimab before their first RSV season. We excluded infants with insufficient information about nirsevimab treatment or without parental consent. We used respiratory samples collected in each laboratory for full-length RSV RNA sequencing to analyse changes in the nirsevimab binding site Ø. We tested clinical RSV isolates for neutralisation by nirsevimab. We analysed F candidate substitutions by fusion-inhibition assay. Of the 695 RSV infected infants, we analysed 545 (78%) full-length RSV genome sequences: 260 (48%) from nirsevimab-treated breakthrough infections (236 [91%] RSV-A and 24 [9%] RSV-B) and 285 (52%) from untreated RSV-infected infants (236 [83%] RSV-A and 49 [17%] RSV-B). Analysis of RSV-A did not reveal any substitution in site Ø known to be associated with resistance to nirsevimab. Two (8%) of 24 RSV-B breakthrough infections had resistance-associated substitutions: F:N208D (dominant resistance-associated substitution) and a newly described F:I64M plus F:K65R combination (minority resistance-associated substitution), both of which induced high levels of resistance in the fusion-inhibition assay. This study is, to the best of our knowledge, the largest genotypic and phenotypic surveillance study of nirsevimab breakthrough infections to date. Nirsevimab breakthrough variants remain very rare despite the drug's widespread use. The detection of resistance-associated substitutions in the RSV-B F protein highlights the importance of active molecular surveillance. ANRS Maladies Infectieuses Emergentes and the French Ministry of Health and Prevention.

Genomic full-length sequence of the HLA-B*40:96 allele was identified by full length group-specific sequencing.

Genomic full-length sequence of HLA-B*40:96 was identified by a group-specific sequencing approach in a Chinese individual.

Rare BK polyomavirus subtype III detected in hematopoietic cell transplant recipient.

BK polyomavirus (BKPyV) is associated with disease in transplant patients. BKPyV genomic diversity is thought to contribute to functional differences in virulence factors. Here, we present the near full-length BKPyV genome sequence from a patient who underwent hematopoietic cell transplant and was infected with subtype III, a rare subtype.

Genotyping and phylogeographic dynamics of coxsackievirus A16

Coxsackievirus A16 (CV-A16) is one of the major pathogens of Hand, Foot and Mouth disease. Here, we analyzed 287 full-length genome sequences of CV-A16 found worldwide from 1994 to 2019 to see the genomic evolution characteristics. Full-length genome-based phylogenetic tree divided the viruses into five different genotypes, G-a to G-e. The CV-A16 strains circulating in China dominate G-a and G-c, but can also be found in other genotypes including G-b and G-e. Phylogeographic analysis showed a high diversity of CV-A16 distribution. In addition, recombination was shown to drive the genomic evolution of CV-A16 during past decades. However, the structural proteins still remain relative conserved while there is extensive genomic recombination. This study updates the phylogenetic and phylogeographic information of CV-A16 and provides insights into the genetic characteristics of CV-A16.

Molecular characterization of avian metapneumovirus subtype C associated with hydrosalpinx fluid and egg drop syndrome in Jinding ducks.

A disease called "hydrosalpinx fluid and egg drop syndrome" (HFEDS) was observed in four flocks of Jinding ducks (Anas platyrhynchos domesticus) in Northeast China during the years 2022 to 2023. Here, we investigated the possible involvement of avian metapneumovirus (AMPV) infection. Full-length genome sequencing and sequence analysis of two AMPV strains showed that they belong to Eurasian lineage of AMPV subtype C. Based on surface glycoprotein (G) sequence comparisons, the Eurasian lineage can be divided into two sublineages (E1 and E2), and sublineage E2 is circulating in Jinding duck populations in Northeast China.

New lineages of RNA viruses from clinical isolates of Rhizopus microsporus revealed by fragmented and primer-ligated dsRNA sequencing (FLDS) analysis.

Rhizopus microsporus is a species in the order Mucorales that is known to cause mucormycosis, but it is poorly understood as a host of viruses. Here, we examined 25 clinical strains of R. microsporus for viral infection with a conventional double-stranded RNA (dsRNA) assay using agarose gel electrophoresis (AGE) and the recently established fragmented and primer-ligated dsRNA sequencing (FLDS) protocol. By AGE, five virus-infected strains were detected. Then, full-length genomic sequences of 12 novel RNA viruses were revealed by FLDS, which were related to the families Mitoviridae, Narnaviridae, and Endornaviridae, ill-defined groups of single-stranded RNA (ssRNA) viruses with similarity to the established families Virgaviridae and Phasmaviridae, and the proposed family "Ambiguiviridae." All the characterized viruses, except a potential phasmavirid with a negative-sense RNA genome, had positive-sense RNA genomes. One virus belonged to a previously established species within the family Mitoviridae, whereas the other 11 viruses represented new species or even new genera. These results show that the fungal pathogen R. microsporus harbors diverse RNA viruses and extend our understanding of the diversity of RNA viruses in the fungal order Mucorales, division Mucoromycota. Identifying RNA viruses from clinical isolates of R. microsporus may expand the repertoire of natural therapeutic agents for mucormycosis in the future.IMPORTANCEThe diversity of mycoviruses in fungal hosts in the division Mucoromycota has been underestimated, mainly within the species Rhizopus microsporus. Only five positive-sense RNA genomes had previously been discovered in this species. Because current sequencing methods poorly complete the termini of genomes, we used fragmented and primer-ligated double-stranded RNA sequencing to acquire the full-length genomes. Eleven novel mycoviruses were detected in this study, including the first negative-sense RNA genome reported in R. microsporus. Our findings extend the understanding of the viral diversity in clinical strains of Mucoromycota, may provide insights into the pathogenesis and ecology of this fungus, and may offer therapeutic options.

Identification of a Novel HIV-1 Circulating Recombinant Form (CRF150_Cpx) Among Men Who Have Sex with Men in China.

Recent studies have reported increasing complexity in human immunodeficiency virus 1 (HIV-1) genotypes among men who have sex with men (MSM) in China. In an HIV-1 molecular epidemiological study conducted among MSM in Yunnan Province, China, we discovered that four samples could potentially represent a circulating recombinant form (CRF). In this study, we conducted further analysis on their nearly full-length genome (NFLG) sequences. The NFLG sequences formed a distinct monophyletic clade in the phylogenetic tree. Recombination analysis indicated that the four sequences were constructed upon the backbone of CRF149_01B, with the insertion of three CRF07_BC fragments. Consequently, they were designated as CRF150_cpx. Evolutionary analyses suggested that CRF150_cpx emerged between approximately 2014 and 2015. The identification of new CRFs not only deepens our understanding of HIV recombination but also aids in comprehending the prevalence and transmission history of HIV among specific populations.

Characterization of bean common mosaic virus isolates infecting three leguminous bean crops from South and Southeast Asia.

Bean common mosaic virus (BCMV) is causing economically important diseases in leguminous crops worldwide. In this study, BCMV isolates from country bean (CB, Lablab purpureus), yard-long bean (YLB, Vigna unguiculata) and rajma bean (RB, Phaseolus vulgaris) collected from Bangladesh, Nepal and Cambodia were characterized. Samples tested positive for BCMV in serological assays were subjected to high-throughput sequencing to generate near full-length genome sequences. In pair-wise comparisons of the polyprotein open reading frame, thirteen BCMV isolates from Bangladesh, Cambodia, and Nepal showed sequence identity of 92.1 to 98.8% at the nucleotide and 94.2 to 99% at the amino acid level among themselves and with corresponding sequences of BCMV reported previously. In phylogenetic analyses using the global BCMV sequences, they segregated into five distinct lineages, with RB isolates from Nepal clustering with US1/NL1-clade of common bean isolates from different countries, YLB isolates aligning with blackeye cowpea strain sequences reported from China, and CB isolates from Nepal and Bangladesh clustering with soybean isolates from China. One YLB isolate from Nepal was identified as a putative recombinant. None of the BCMV sequences aligned with isolates representing the RU1 or PStV clades. In grow-out tests, seed samples from local markets showed 14.3 to 38.1% transmission efficiency rate of BCMV with CB seed lots and from 9.5% to 33.3% with YLB seed lots.

Plasmid-mediated acquisition and chromosomal integration of blaCTX-M-14 in a subclade of Escherichia coli ST131-H30 clade C1.

Escherichia coli ST131 is a multidrug-resistant lineage associated with the global spread of extended-spectrum β-lactamase-producing organisms. Particularly, ST131 clade C1 is the most predominant clade in Japan, harboring blaCTX-M-14 at a high frequency. However, the process of resistance gene acquisition and spread remains unclear. Here, we performed whole-genome sequencing of 19 E. coli strains belonging to 12 STs and 12 fimH types collected between 1997 and 2016. Additionally, we analyzed the full-length genome sequences of 96 ST131-H30 clade C0 and C1 strains, including those obtained from this study and those registered in public databases, to understand how ST131 clade C1 acquired and spread blaCTX-M-14. We detected conjugative IncFII plasmids and IncB/O/K/Z plasmids carrying blaCTX-M-14 in diverse genetic lineages of E. coli strains from the 1990s to the 2010s, suggesting that these plasmids played an important role in the spread of blaCTX-M-14. Molecular phylogenetic and molecular clock analyses of the 96 ST131-H30 clade C0 and C1 strains identified 8 subclades. Strains harboring blaCTX-M-14 were clustered in subclades 4 and 5, and it was inferred that clade C1 acquired blaCTX-M-14 around 1993. All 34 strains belonging to subclade 5 possessed blaCTX-M-14 with ISEcp1 upstream at the same chromosomal position, indicating their common ancestor acquired blaCTX-M-14 in a single ISEcp1-mediated transposition event during the early formation of the subclade around 1999. Therefore, both the horizontal transfer of plasmids carrying blaCTX-M-14 to diverse genetic lineages and chromosomal integration in the predominant genetic lineage have contributed to the spread of blaCTX-M-14.

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 &lt;i&gt;in silico&lt;/i&gt; 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 &lt;i&gt;B. anthracis&lt;/i&gt; through nucleic acid isolation and genetic sequence analysis studies. MLVA results on 32 loci of &lt;i&gt;B. anthracis&lt;/i&gt; 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 &lt;i&gt;B. anthracis&lt;/i&gt;, 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 &lt;i&gt;B. anthracis&lt;/i&gt; strain A4564 which is associated with injectional anthrax isolates in heroin users.

The full-length genomic sequence of the HLA-G*01:19 allele.

Full genomic sequence shows HLA-G*01:19 differs from HLA-G*01:04:01:01 only at position 99 in exon 2.

Genomic full-length sequence of the novel HLA-B*40:540N allele.

The coding sequence of HLA-B*40:540N differs from HLA-B*40:02:01:01 by a non-synonymous nucleotide substitution in exon 3.

Highly homologous simian T-cell leukemia virus type 1 genome in Japanese macaques: a large cohort study

BackgroundSimian T-cell leukemia virus type 1 (STLV-1) is a retrovirus closely related to human T-cell leukemia virus type 1 (HTLV-1), the causative agent of adult T-cell leukemia (ATL). It has been shown that Japanese macaques (Macaca fuscata, JMs) are one of the main hosts of STLV-1 and that a high percentage of JMs (up to 60%) are infected with STLV-1; however, the molecular epidemiology of STLV-1 in JMs has not been examined.MethodsIn this study, we analyzed full-length STLV-1 genome sequences obtained from 5 independent troops including a total of 68 JMs.ResultsThe overall nucleotide heterogeneity was 4.7%, and the heterogeneity among the troops was 2.1%, irrespective of the formation of distinct subclusters in each troop. Moreover, the heterogeneity within each troop was extremely low (>99% genome homology) compared with cases of STLV-1 in African non-human primates as well as humans. It was previously reported that frequent G-to-A single-nucleotide variants (SNVs) occur in HTLV-1 proviral genomes in both ATL patients and HTLV-1 carriers, and that a G-to-A hypermutation is associated with the cellular antiviral restriction factor, Apobec3G. Surprisingly, these SNVs were scarcely observed in the STLV-1 genomes in JMs.ConclusionsTaken together, these results indicate that STLV-1 genomes in JMs are highly homologous, at least in part due to the lack of Apobec3G-dependent G-to-A hypermutation.

Identification and characterization of olfactory gene families in Macrobrachium rosenbergii based on full-length transcripts and genome sequences

The olfactory gene families include odorant binding proteins (OBPs), chemosensory proteins (CSPs), olfactory receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs). To investigate the molecular function of olfactory perception in Macrobrachium rosenbergii, we integrated the full-length transcripts and whole-genome sequences to identify the olfactory gene families. In this study, a total of 38,955 full-length transcripts with an N50 length of 3383 bp were obtained through PacBio SMRT sequencing. Through the annotation of full-length transcripts and whole-genome sequences, several olfactory gene families were identified, including 18 MrORs, 16 MrIRs, 151 MrIGluRs (ionotropic glutamate receptors), 2 MrVIGluRs (variant ionotropic glutamate receptors) and 3 MrCRs (chemosensory receptors). Notably, the CRs were first identified in prawns and shrimps. Additionally, the olfactory gene families in M. nipponense were identified, comprising 4 MnORs, 21 MnIRs, 79 MnIGluRs, 5 MnVIGluRs, 1 MnGR and 1 MnOBP, using the available whole-genome sequences. Meanwhile, the external morphology of the chemical sensory organs of M. rosenbergii was explored, and the presence of plumose setae (PS), hard thorn setae (HTS), bamboo shoot setae (BSS), soft thorn setae (STS) and aesthetascs (AE) on the antennules, HTS and BSS on the second antennae, and PS on the pereiopods were observed by scanning electron microscope. This study provides valuable insights for future functional studies into the olfactory perception of crustaceans and establishes a theoretical basis for molecular design breeding in M. rosenbergii.

Description of a new tick species, closely related to Amblyomma javanense (Supino, 1897), associated with Varanus bengalensis (Squamata: Varanidae) in Pakistan

The genus Amblyomma contains the highest percentage of reptile-associated ticks, and comprises approximately nine subgenera. One of these subgenera is Adenopleura, which also encompasses Amblyomma javanense, and its type species Amblyomma compressum. This study describes a new Amblyomma species associated with Bengal monitor lizards (Varanus bengalensis) based on morphology and its mitogenome in Khyber Pakhtunkhwa, Pakistan. Reptiles belonging to different genera were examined for Amblyomma ticks and only the monitor lizard was infested with ticks in the District Bajaur. Collected Amblyomma cf. javanense ticks were analyzed and formally described as a new species. Overall, 57 A. cf. javanense ticks were collected on monitor lizards (4/27) with a 15% prevalence of infestation, 2.1 mean abundance, and 14.3 mean intensity. Ticks comprised males (n = 23, 40%), females (n = 14, 25%) and nymphs (n = 20, 35%), while no larvae were found. BLAST analysis of A. cf. javanense sequences showed the following maximum identities; 98.25% with undetermined Amblyomma species based on 12S rRNA, 96.07% with A. javanense based on 16S rRNA, 99.56% and 90.95% with an Amblyomma sp. and A. javanense, respectively, based on ITS2. Moreover, the mitochondrial genome of A. cf. javanense showed maximum identities of 80.75%, 80.48% and 79.42% with Amblyomma testudinarium, A. javanense, and Amblyomma sp., respectively. The phylogenetic analysis of A. cf. javanense revealed that its 12S rRNA and 16S rRNA are closely related to an Amblyomma sp. and A. javanense, respectively, from Sri Lanka, its ITS2 is closely related to A. javanense from China and an Amblyomma sp. from Sri Lanka, and its mitogenome is closely related to A. javanense and Amblyomma sp. from China. The pairwise distance analysis resulted in divergence of 0–1.71% (12S rRNA), 0–17.5% (16S rRNA), 0–9.1% (ITS2) and 0–20.5% (mitochondrial genome). We also contributed the full-length mitochondrial genome sequence of A. compressum and showed that this species does not share a most recent common ancestor with A. javanense. As the subgenus Adenopleura is paraphyletic, this study could help to understand the systematics and phylogeny of this taxon.

Identification of a novel waikavirus infecting Pittosporum tobira in China.

A novel waikavirus, tentatively named "Pittosporum tobira waikavirus" (PtWV), was identified in Pittosporum tobira plants exhibiting mosaic and ringspot symptoms on foliage in Yunnan, China. The full-length genomic sequence was determined by high-throughput sequencing and rapid amplification of cDNA ends. The genome of PtWV is 12,709 nt in length and has a large open reading frame (ORF) of 11,010 nt, encoding a polyprotein, and a small ORF that encodes a 13.2-kDa bellflower vein chlorosis virus (BVCV)-like protein. Phylogenetic analysis and sequence alignment revealed that PtWV is closely related to actinidia yellowing virus 1 (AcYV1), which shares the highest amino acid (aa) sequence similarity (50.1% identity) in the Pro-RdRp region. To the best of our knowledge, this is the first report of a novel waikavirus in P. tobira.

Diverse amdoparvoviruses infection of farmed Asian badgers (Meles meles).

Amdoparvoviruses infect various carnivores, including mustelids, canids, skunks, and felids. Aleutian mink disease virus (AMDV) belongs to the prototypical species Amdoparvovirus carnivoran1. Here, we identified a novel amdoparvovirus in farmed Asian badgers (Meles meles), and we named this virus "Meles meles amdoparvovirus" (MMADV). A total of 146 clinical samples were collected from 134 individual badgers, and 30.6% (41/134) of the sampled badgers tested positive for amdoparvovirus by PCR. Viral DNA was detected in feces, blood, spleen, liver, lung, and adipose tissue from these animals. Viral sequences from eight samples were determined, five of which represented nearly full-length genome sequences (4,237-4,265 nt). Six serum samples tested positive by PCR, CIEP, and IAT, four of which had high antibody titers (> 512) against AMDV-G. Twenty-six of the 41 amdoparvovirus-positive badgers showed signs of illness, and necropsy revealed lesions in their organs. Sequence comparisons and phylogenetic analysis of the viral NS1 and VP2 genes of these badger amdoparvoviruses showed that their NS1 proteins shared 62.6%-88.8% sequence identity with known amdoparvoviruses, and they clustered phylogenetically into two related clades. The VP2 proteins shared 76.6%-97.2% identity and clustered into two clades, one of which included raccoon dog and arctic fox amdoparvovirus (RFAV), and the other of which did not include other known amdoparvoviruses. According to the NS1-protein-based criterion for parvovirus species demarcation, the MMADV isolate from farm YS should be classified as a member of a new species of the genus Amdoparvovirus. In summary, we have discovered a novel MMADV and other badger amdoparvoviruses that naturally infect Asian badgers and are possibly pathogenic in badgers.