Large-Scale Comparative Analyses of Tick Genomes Elucidate Their Genetic Diversity and Vector Capacities.

Among arthropod vectors, ticks transmit the most diverse human and animal pathogens, leading to an increasing number of new challenges worldwide. Here we sequenced and assembled high-quality genomes of six ixodid tick species and further resequenced 678 tick specimens to understand three key aspects of ticks: genetic diversity, population structure, and pathogen distribution. We explored the genetic basis common to ticks, including heme and hemoglobin digestion, iron metabolism, and reactive oxygen species, and unveiled for the first time that genetic structure and pathogen composition in different tick species are mainly shaped by ecological and geographic factors. We further identified species-specific determinants associated with different host ranges, life cycles, and distributions. The findings of this study are an invaluable resource for research and control of ticks and tick-borne diseases.

Identifying SARS-CoV-2-related coronaviruses in Malayan pangolins.

The ongoing outbreak of viral pneumonia in China and across the world is associated with a new coronavirus, SARS-CoV-21. This outbreak has been tentatively associated with a seafood market in Wuhan, China, where the sale of wild animals may be the source of zoonotic infection2. Although bats are probable reservoir hosts for SARS-CoV-2, the identity of any intermediate host that may have facilitated transfer to humans is unknown. Here we report the identification of SARS-CoV-2-related coronaviruses in Malayan pangolins (Manis javanica) seized in anti-smuggling operations in southern China. Metagenomic sequencing identified pangolin-associated coronaviruses that belong to two sub-lineages of SARS-CoV-2-related coronaviruses, including one that exhibits strong similarity in the receptor-binding domain to SARS-CoV-2. The discovery of multiple lineages of pangolin coronavirus and their similarity to SARS-CoV-2 suggests that pangolins should be considered as possible hosts in the emergence of new coronaviruses and should be removed from wet markets to prevent zoonotic transmission.

Robust expansion of phylogeny for fast-growing genome sequence data.

Massive sequencing of SARS-CoV-2 genomes has urged novel methods that employ existing phylogenies to add new samples efficiently instead of de novo inference. 'TIPars' was developed for such challenge integrating parsimony analysis with pre-computed ancestral sequences. It took about 21 seconds to insert 100 SARS-CoV-2 genomes into a 100k-taxa reference tree using 1.4 gigabytes. Benchmarking on four datasets, TIPars achieved the highest accuracy for phylogenies of moderately similar sequences. For highly similar and divergent scenarios, fully parsimony-based and likelihood-based phylogenetic placement methods performed the best respectively while TIPars was the second best. TIPars accomplished efficient and accurate expansion of phylogenies of both similar and divergent sequences, which would have broad biological applications beyond SARS-CoV-2. TIPars is accessible from https://tipars.hku.hk/ and source codes are available at https://github.com/id-bioinfo/TIPars.

Amino acid substitution L232F in non-structural protein 6 identified as a possible human-adaptive mutation in clade B MERS coronaviruses.

IMPORTANCE: Viral host adaptation plays an important role in inter-species transmission of coronaviruses and influenza viruses. Multiple human-adaptive mutations have been identified in influenza viruses but not so far in MERS-CoV that circulates widely in dromedary camels in the Arabian Peninsula leading to zoonotic transmission. Here, we analyzed clade B MERS-CoV sequences and identified an amino acid substitution L232F in nsp6 that repeatedly occurs in human MERS-CoV. Using a loss-of-function reverse genetics approach, we found the nsp6 L232F conferred increased viral replication competence in vitro, in cultures of the upper human respiratory tract ex vivo, and in lungs of mice infected in vivo. Our results showed that nsp6 L232F may be an adaptive mutation associated with zoonotic transmission of MERS-CoV. This study highlighted the capacity of MERS-CoV to adapt to transmission to humans and also the need for continued surveillance of MERS-CoV in camels.

Characterization of an Emergent Chicken H3N8 Influenza Virus in Southern China: a Potential Threat to Public Health.

Although influenza A viruses of several subtypes have occasionally infected humans, to date only those of the H1, H2, and H3 subtypes have led to pandemics and become established in humans. The detection of two human infections by avian H3N8 viruses in April and May of 2022 raised pandemic concerns. Recent studies have shown the H3N8 viruses were introduced into humans from poultry, although their genesis, prevalence, and transmissibility in mammals have not been fully elucidated. Findings generated from our systematic influenza surveillance showed that this H3N8 influenza virus was first detected in chickens in July 2021 and then disseminated and became established in chickens over wider regions of China. Phylogenetic analyses revealed that the H3 HA and N8 NA were derived from avian viruses prevalent in domestic ducks in the Guangxi-Guangdong region, while all internal genes were from enzootic poultry H9N2 viruses. The novel H3N8 viruses form independent lineages in the glycoprotein gene trees, but their internal genes are mixed with those of H9N2 viruses, indicating continuous gene exchange among these viruses. Experimental infection of ferrets with three chicken H3N8 viruses showed transmission through direct contact and inefficient transmission by airborne exposure. Examination of contemporary human sera detected only very limited antibody cross-reaction to these viruses. The continuing evolution of these viruses in poultry could pose an ongoing pandemic threat. IMPORTANCE A novel H3N8 virus with demonstrated zoonotic potential has emerged and disseminated in chickens in China. It was generated by reassortment between avian H3 and N8 virus(es) and long-term enzootic H9N2 viruses present in southern China. This H3N8 virus has maintained independent H3 and N8 gene lineages but continues to exchange internal genes with other H9N2 viruses to form novel variants. Our experimental studies showed that these H3N8 viruses were transmissible in ferrets, and serological data suggest that the human population lacks effective immunological protection against it. With its wide geographical distribution and continuing evolution in chickens, other spillovers to humans can be expected and might lead to more efficient transmission in humans.

ggmsa: a visual exploration tool for multiple sequence alignment and associated data.

The identification of the conserved and variable regions in the multiple sequence alignment (MSA) is critical to accelerating the process of understanding the function of genes. MSA visualizations allow us to transform sequence features into understandable visual representations. As the sequence-structure-function relationship gains increasing attention in molecular biology studies, the simple display of nucleotide or protein sequence alignment is not satisfied. A more scalable visualization is required to broaden the scope of sequence investigation. Here we present ggmsa, an R package for mining comprehensive sequence features and integrating the associated data of MSA by a variety of display methods. To uncover sequence conservation patterns, variations and recombination at the site level, sequence bundles, sequence logos, stacked sequence alignment and comparative plots are implemented. ggmsa supports integrating the correlation of MSA sequences and their phenotypes, as well as other traits such as ancestral sequences, molecular structures, molecular functions and expression levels. We also design a new visualization method for genome alignments in multiple alignment format to explore the pattern of within and between species variation. Combining these visual representations with prime knowledge, ggmsa assists researchers in discovering MSA and making decisions. The ggmsa package is open-source software released under the Artistic-2.0 license, and it is freely available on Bioconductor (https://bioconductor.org/packages/ggmsa) and Github (https://github.com/YuLab-SMU/ggmsa).

Australia as a global sink for the genetic diversity of avian influenza A virus.

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, on the available data we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.

The Predictive Value of Gut Microbiota Composition for Sustained Immunogenicity following Two Doses of CoronaVac.

CoronaVac immunogenicity decreases with time, and we aimed to investigate whether gut microbiota associate with longer-term immunogenicity of CoronaVac. This was a prospective cohort study recruiting two-dose CoronaVac recipients from three centres in Hong Kong. We collected blood samples at baseline and day 180 after the first dose and used chemiluminescence immunoassay to test for neutralizing antibodies (NAbs) against the receptor-binding domain (RBD) of wild-type SARS-CoV-2 virus. We performed shotgun metagenomic sequencing performed on baseline stool samples. The primary outcome was the NAb seroconversion rate (seropositivity defined as NAb ≥ 15AU/mL) at day 180. Linear discriminant analysis [LDA] effect size analysis was used to identify putative bacterial species and metabolic pathways. A univariate logistic regression model was used to derive the odds ratio (OR) of seropositivity with bacterial species. Of 119 CoronaVac recipients (median age: 53.4 years [IQR: 47.8-61.3]; male: 39 [32.8%]), only 8 (6.7%) remained seropositive at 6 months after vaccination. Bacteroides uniformis (log10LDA score = 4.39) and Bacteroides eggerthii (log10LDA score = 3.89) were significantly enriched in seropositive than seronegative participants. Seropositivity was associated with B. eggerthii (OR: 5.73; 95% CI: 1.32-29.55; p = 0.022) and B. uniformis with borderline significance (OR: 3.27; 95% CI: 0.73-14.72; p = 0.110). Additionally, B. uniformis was positively correlated with most enriched metabolic pathways in seropositive vaccinees, including the superpathway of adenosine nucleotide de novo biosynthesis I (log10LDA score = 2.88) and II (log10LDA score = 2.91), as well as pathways related to vitamin B biosynthesis, all of which are known to promote immune functions. In conclusion, certain gut bacterial species (B. eggerthii and B. uniformis) and metabolic pathways were associated with longer-term CoronaVac immunogenicity.

Tracking the Genomic Footprints of SARS-CoV-2 Transmission.

There is considerable public and scientific interest in the origin, spread, and evolution of SARS-CoV-2. Lu et al. recently conducted genomic sequencing and analysis of SARS-CoV-2 in Guangdong, revealing its early transmission out of Hubei and shedding light on the effectiveness of controlling local transmission chains.

Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia.

BACKGROUND: The initial cases of novel coronavirus (2019-nCoV)-infected pneumonia (NCIP) occurred in Wuhan, Hubei Province, China, in December 2019 and January 2020. We analyzed data on the first 425 confirmed cases in Wuhan to determine the epidemiologic characteristics of NCIP. METHODS: We collected information on demographic characteristics, exposure history, and illness timelines of laboratory-confirmed cases of NCIP that had been reported by January 22, 2020. We described characteristics of the cases and estimated the key epidemiologic time-delay distributions. In the early period of exponential growth, we estimated the epidemic doubling time and the basic reproductive number. RESULTS: Among the first 425 patients with confirmed NCIP, the median age was 59 years and 56% were male. The majority of cases (55%) with onset before January 1, 2020, were linked to the Huanan Seafood Wholesale Market, as compared with 8.6% of the subsequent cases. The mean incubation period was 5.2 days (95% confidence interval [CI], 4.1 to 7.0), with the 95th percentile of the distribution at 12.5 days. In its early stages, the epidemic doubled in size every 7.4 days. With a mean serial interval of 7.5 days (95% CI, 5.3 to 19), the basic reproductive number was estimated to be 2.2 (95% CI, 1.4 to 3.9). CONCLUSIONS: On the basis of this information, there is evidence that human-to-human transmission has occurred among close contacts since the middle of December 2019. Considerable efforts to reduce transmission will be required to control outbreaks if similar dynamics apply elsewhere. Measures to prevent or reduce transmission should be implemented in populations at risk. (Funded by the Ministry of Science and Technology of China and others.).

Toward a Universal Unit for Quantification of Antibiotic Resistance Genes in Environmental Samples.

Surveillance of antibiotic resistance genes (ARGs) has been increasingly conducted in environmental sectors to complement the surveys in human and animal sectors under the "One-Health" framework. However, there are substantial challenges in comparing and synthesizing the results of multiple studies that employ different test methods and approaches in bioinformatic analysis. In this article, we consider the commonly used quantification units (ARG copy per cell, ARG copy per genome, ARG density, ARG copy per 16S rRNA gene, RPKM, coverage, PPM, etc.) for profiling ARGs and suggest a universal unit (ARG copy per cell) for reporting such biological measurements of samples and improving the comparability of different surveillance efforts.

COVID-19 Vaccine Allergy Safety Track (VAS-Track) pathway: real-world outcomes on vaccination rates and antibody protection.

BACKGROUND: Misdiagnosed vaccine-related "allergies" lead to unnecessary vaccine deferrals and incomplete vaccinations, leaving patients unprotected against COVID-19. To overcome limitations and queues for Allergist assessment, the "VAS-Track" pathway was developed to evaluate patients via a multi-disciplinary triage model including nurses, non-specialists, and Allergists. OBJECTIVE: We assessed the effectiveness and safety of VAS-Track and evaluate its real-world impact in terms of vaccination rates and COVID-19 protection. METHODS: Patients referred to VAS-Track between September 2021 and March 2022 were recruited. Subgroup analysis was performed with prospective pre- and post-clinic antibody levels. RESULTS: Nurse-assisted screening identified 10,412 (76%) referrals as inappropriate. 369 patients were assessed by VAS-Track. Overall, 100% of patients were recommended to complete vaccination and 332 (90%) completed their primary series. No patients reported any significant allergic reactions following subsequent vaccination. Vaccination completion rates between patients seen by non-specialists and additional Allergist review were similar (90% vs. 89%, p = 0.617). Vaccination rates were higher among patients with prior history of immediate-type reactions (odds ratio: 2.43, p = 0.025). Subgroup analysis revealed that only 20% (56/284) of patients had seropositive COVID-19 neutralizing antibody levels (≥ 15 AU/mL) prior to VAS-Track, which increased to 92% after vaccine completion (pre-clinic antibody level 6.0 ± 13.5 AU/mL vs. post-clinic antibody level 778.8 ± 337.4 AU/mL, p > 0.001). CONCLUSIONS: A multi-disciplinary allergy team was able to streamline our COVID-19 VAS services, enabling almost all patients to complete their primary series, significantly boosting antibody levels and real-world COVID-19 protection. We propose similar multidisciplinary models to be further utilized, especially in the settings with limited allergy services.

Recombinant BA.1/BA.2 SARS-CoV-2 Virus in Arriving Travelers, Hong Kong, February 2022.

We studied SARS-CoV-2 genomes from travelers arriving in Hong Kong during November 2021-February 2022. In addition to Omicron and Delta variants, we detected a BA.1/BA.2 recombinant with a breakpoint near the 5' end of the spike gene in 2 epidemiologically linked case-patients. Continued surveillance for SARS-CoV-2 recombinants is needed.

Surveillance of Rodent Pests for SARS-CoV-2 and Other Coronaviruses, Hong Kong.

We report surveillance conducted in 217 pestiferous rodents in Hong Kong for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We did not detect SARS-CoV-2 RNA but identified 1 seropositive rodent, suggesting exposure to a virus antigenically similar to SARS-CoV-2. Potential exposure of urban rodents to SARS-CoV-2 cannot be ruled out.

Predictors of rapidly progressive interstitial lung disease and mortality in patients with autoantibodies against melanoma differentiation-associated protein 5 dermatomyositis.

OBJECTIVE: Anti-melanoma differentiation-associated protein 5 (MDA5)-positive DM is associated with rapidly progressive interstitial lung disease (RP-ILD) and high mortality. This multicentre retrospective study aimed to identify predictors of mortality and RP-ILD. METHODS: Anti-MDA5-positive DM patients were identified from the Hong Kong Myositis Registry and the Clinical Data Analysis and Reporting System. Clinical characteristics were reviewed. Risk factors for mortality and RP-ILD were identified. RESULTS: Among the 116 recruited patients, 100 (86.2%) had ILD, 47 (40.5%) had RP-ILD and 44 (37.9%) patients died. Cox regression analysis revealed RP-ILD [hazard ratio (HR) 9.735 (95% CI 3.905, 24.272)], age >52 years [HR 4.750 (95% CI 1.692, 13.333)], ferritin level >2800 pmol/l [HR 3.042 (95% CI 1.323, 6.997)] and lactate dehydrogenase (LDH) >400 IU/l [HR 2.290 (95% CI 1.009, 5.198)] were independent predictors of mortality. With regard to RP-ILD, analyses showed that potential predictors at baseline included age >50 years [HR 2.640 (95% CI 1.277, 5.455)], LDH >300 IU/l [HR 3.189 (95% CI 1.469, 6.918)], fever [HR 1.903 (95% CI 0.956, 3.790)] and neutrophil:lymphocyte ratio >7.0 [HR 1.967 (95% CI 0.942, 4.107)]. We proposed a prediction model based on fever, LDH, age and white cell count (FLAW) to stratify the risk of development of RP-ILD. The probability of RP-ILD in a patient with a score of 4 was 100%. A small internal validation cohort showed the odds of RP-ILD with FLAW scores of 0, 1, 2 and 3 were 0%, 0%, 42.9% and 75%, respectively. CONCLUSIONS: Anti-MDA5-associated RP-ILD is significantly associated with poor survival rates. The FLAW model maybe useful to predict the development of RP-ILD.

Genotype and organ effect on the occupancy of phyllosphere prokaryotes in different rice landraces.

Bacteria community provides essential ecological services to rice plants. The bacterial diversity of rice varies across host plant genotype and organs. This study employed 16S rDNA amplicon sequencing to characterise the bacterial community associated with three rice landraces using leaf blade and stem samples. The prokaryotic community found in these rice landraces comprised of two kingdoms, 12 phyla, 25 classes, 40 orders, 80 families, and 118 genera. Proteobacteria (53.9%) was the most abundant phylum. The most abundant genus was an undefined genus under Cyanobacteria (33.0%). Homogeneity of prokaryotic community was observed across the three rice landraces, which may suggest a high similarity in biological and genetical properties of the rice landraces. The difference in prokaryotic composition between leaf blade and stem was depicted based on principal coordinate analysis. This study observed that the prokaryotic inhabitants in rice plants is predominantly determined by rice plant organs.

Association between inflammatory cytokines and anti-SARS-CoV-2 antibodies in hospitalized patients with COVID-19.

BACKGROUND: COVID-19 patients may experience "cytokine storm" when human immune system produces excessive cytokines/chemokines. However, it remains unclear whether early responses of inflammatory cytokines would lead to high or low titers of anti-SARS-CoV-2 antibodies. METHODS: This retrospective study enrolled a cohort of 272 hospitalized patients with laboratory-confirmed SARS-CoV-2. Laboratory assessments of serum cytokines (IL-2R, IL-6, IL-8, IL-10, TNF-α), anti-SARS-CoV-2 IgG/IgM antibodies, and peripheral blood biomarkers were conducted during hospitalization. RESULTS: At hospital admission, 36.4% patients were severely ill, 51.5% patients were ≥ 65 years, and 60.3% patients had comorbidities. Higher levels of IL-2R and IL-6 were observed in older patients (≥65 years). Significant differences of IL-2R (week 2 to week ≥5 from symptom onset), IL-6 (week 1 to week ≥5), IL-8 (week 2 to week ≥5), and IL-10 (week 1 to week 3) were observed between moderately-ill and severely ill patients. Anti-SARS-CoV-2 IgG titers were significantly higher in severely ill patients than in moderately ill patients, but such difference was not observed for IgM. High titers of early-stage IL-6, IL-8, and TNF-α (≤2 weeks after symptom onset) were positively correlated with high titers of late-stage IgG (≥5 weeks after symptom onset). Deaths were mostly observed in severely ill older patients (45.9%). Survival analyses revealed risk factors of patient age, baseline COVID-19 severity, and baseline IL-6 that affected survival time, especially in severely ill older patients. CONCLUSION: Early responses of elevated cytokines such as IL-6 reflect the active immune responses, leading to high titers of IgG antibodies against COVID-19.

Highly Pathogenic Avian Influenza A(H5N6) Virus Clade 2.3.4.4h in Wild Birds and Live Poultry Markets, Bangladesh.

Migratory birds play a major role in spreading influenza viruses over long distances. We report highly pathogenic avian influenza A(H5N6) viruses in migratory and resident ducks in Bangladesh. The viruses were genetically similar to viruses detected in wild birds in China and Mongolia, suggesting migration-associated dissemination of these zoonotic pathogens.

Treeio: An R Package for Phylogenetic Tree Input and Output with Richly Annotated and Associated Data.

Phylogenetic trees and data are often stored in incompatible and inconsistent formats. The outputs of software tools that contain trees with analysis findings are often not compatible with each other, making it hard to integrate the results of different analyses in a comparative study. The treeio package is designed to connect phylogenetic tree input and output. It supports extracting phylogenetic trees as well as the outputs of commonly used analytical software. It can link external data to phylogenies and merge tree data obtained from different sources, enabling analyses of phylogeny-associated data from different disciplines in an evolutionary context. Treeio also supports export of a phylogenetic tree with heterogeneous-associated data to a single tree file, including BEAST compatible NEXUS and jtree formats; these facilitate data sharing as well as file format conversion for downstream analysis. The treeio package is designed to work with the tidytree and ggtree packages. Tree data can be processed using the tidy interface with tidytree and visualized by ggtree. The treeio package is released within the Bioconductor and rOpenSci projects. It is available at https://www.bioconductor.org/packages/treeio/.

Incorporating heterogeneous sampling probabilities in continuous phylogeographic inference - Application to H5N1 spread in the Mekong region.

MOTIVATION: The potentially low precision associated with the geographic origin of sampled sequences represents an important limitation for spatially explicit (i.e. continuous) phylogeographic inference of fast-evolving pathogens such as RNA viruses. A substantial proportion of publicly available sequences is geo-referenced at broad spatial scale such as the administrative unit of origin, rather than more precise locations (e.g. geographic coordinates). Most frequently, such sequences are either discarded prior to continuous phylogeographic inference or arbitrarily assigned to the geographic coordinates of the centroid of their administrative area of origin for lack of a better alternative. RESULTS: We here implement and describe a new approach that allows to incorporate heterogeneous prior sampling probabilities over a geographic area. External data, such as outbreak locations, are used to specify these prior sampling probabilities over a collection of sub-polygons. We apply this new method to the analysis of highly pathogenic avian influenza H5N1 clade data in the Mekong region. Our method allows to properly include, in continuous phylogeographic analyses, H5N1 sequences that are only associated with large administrative areas of origin and assign them with more accurate locations. Finally, we use continuous phylogeographic reconstructions to analyse the dispersal dynamics of different H5N1 clades and investigate the impact of environmental factors on lineage dispersal velocities. AVAILABILITY AND IMPLEMENTATION: Our new method allowing heterogeneous sampling priors for continuous phylogeographic inference is implemented in the open-source multi-platform software package BEAST 1.10. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.