Evaluation of respondent-driven sampling in seven studies of people who use drugs from rural populations: findings from the Rural Opioid Initiative.

BACKGROUND: Accurate prevalence estimates of drug use and its harms are important to characterize burden and develop interventions to reduce negative health outcomes and disparities. Lack of a sampling frame for marginalized/stigmatized populations, including persons who use drugs (PWUD) in rural settings, makes this challenging. Respondent-driven sampling (RDS) is frequently used to recruit PWUD. However, the validity of RDS-generated population-level prevalence estimates relies on assumptions that should be evaluated. METHODS: RDS was used to recruit PWUD across seven Rural Opioid Initiative studies between 2018-2020. To evaluate RDS assumptions, we computed recruitment homophily and design effects, generated convergence and bottleneck plots, and tested for recruitment and degree differences. We compared sample proportions with three RDS-adjusted estimators (two variations of RDS-I and RDS-II) for five variables of interest (past 30-day use of heroin, fentanyl, and methamphetamine; past 6-month homelessness; and being positive for hepatitis C virus (HCV) antibody) using linear regression with robust confidence intervals. We compared regression estimates for the associations between HCV positive antibody status and (a) heroin use, (b) fentanyl use, and (c) age using RDS-1 and RDS-II probability weights and no weights using logistic and modified Poisson regression and random-effects meta-analyses. RESULTS: Among 2,842 PWUD, median age was 34 years and 43% were female. Most participants (54%) reported opioids as their drug of choice, however regional differences were present (e.g., methamphetamine range: 4-52%). Many recruitment chains were not long enough to achieve sample equilibrium. Recruitment homophily was present for some variables. Differences with respect to recruitment and degree varied across studies. Prevalence estimates varied only slightly with different RDS weighting approaches, most confidence intervals overlapped. Variations in measures of association varied little based on weighting approach. CONCLUSIONS: RDS was a useful recruitment tool for PWUD in rural settings. However, several violations of key RDS assumptions were observed which slightly impacts estimation of proportion although not associations.

A computationally designed antigen eliciting broad humoral responses against SARS-CoV-2 and related sarbecoviruses.

The threat of spillovers of coronaviruses associated with the severe acute respiratory syndrome (SARS) from animals to humans necessitates vaccines that offer broader protection from sarbecoviruses. By leveraging a viral-genome-informed computational method for selecting immune-optimized and structurally engineered antigens, here we show that a single antigen based on the receptor binding domain of the spike protein of sarbecoviruses elicits broad humoral responses against SARS-CoV-1, SARS-CoV-2, WIV16 and RaTG13 in mice, rabbits and guinea pigs. When administered as a DNA immunogen or by a vector based on a modified vaccinia virus Ankara, the optimized antigen induced vaccine protection from the Delta variant of SARS-CoV-2 in mice genetically engineered to express angiotensin-converting enzyme 2 and primed by a viral-vector vaccine (AZD1222) against SARS-CoV-2. A vaccine formulation incorporating mRNA coding for the optimized antigen further validated its broad immunogenicity. Vaccines that elicit broad immune responses across subgroups of coronaviruses may counteract the threat of zoonotic spillovers of betacoronaviruses.

The IRDiRC Chrysalis Task Force: making rare disease research attractive to companies.

Background: The International Rare Diseases Research Consortium (IRDiRC) is an international initiative that aims to use research to facilitate rapid diagnosis and treatment of rare diseases. Objective: IRDiRC launched the Chrysalis Task Force to identify key financial and nonfinancial factors that make rare disease research and development attractive to companies. Methods: The Chrysalis Task Force was comprised of thought leaders from companies, patient advocacy groups, regulatory agencies, and research funders. The Task Force created a survey that was distributed to companies of different sizes with varied investment portfolios and interests in rare disease research. Based on the survey results, the Task Force then conducted targeted interviews. Results: The survey and interview respondents identified several factors that make rare disease research and development attractive (e.g. a good understanding of the underlying biology) as well as barriers (e.g. absence of an advocacy organization representing the affected community's needs). The concept of Return On Investment allowed the exploration of factors that were weighed differently by survey and interview respondents, depending on a number of intrinsic and extrinsic issues. Conclusions: The Chrysalis Task Force identified factors attributable to rare disease research and development that may be of interest to and actionable by funders, academic researchers, patients and their families, companies, regulators, and payers in the medium term to short term. By addressing the identified challenges, involved parties may seek solutions to significantly advance the research and development of treatments for rare diseases.

Making rare disease research attractive to companies The International Rare Diseases Research Consortium (IRDiRC) is an international initiative that aims to speed the diagnosis and treatment of rare diseases through research. The IRDiRC Chrysalis Task Force, comprised of thought leaders from companies, patient advocacy groups, regulatory agencies, and research funders, identified key factors that make rare disease research and development attractive to companies. The Task Force distributed a survey to companies with varied investment portfolios and interests in rare disease research, followed by in-depth interviews based on the survey results. The survey and interview respondents identified both attractive factors and barriers to rare disease research and development. The concept of Return On Investment was used to frame discussion of factors that companies weighed differently, depending on a number of issues that were a function of both the company itself and outside factors. The identified challenges can be addressed by funders, academic researchers, patients and their families, companies, regulators, and payers, which hopefully will lead to significant advances in the research and development of treatments for rare diseases.

Metagenomic surveillance uncovers diverse and novel viral taxa in febrile patients from Nigeria.

Effective infectious disease surveillance in high-risk regions is critical for clinical care and pandemic preemption; however, few clinical diagnostics are available for the wide range of potential human pathogens. Here, we conduct unbiased metagenomic sequencing of 593 samples from febrile Nigerian patients collected in three settings: i) population-level surveillance of individuals presenting with symptoms consistent with Lassa Fever (LF); ii) real-time investigations of outbreaks with suspected infectious etiologies; and iii) undiagnosed clinically challenging cases. We identify 13 distinct viruses, including the second and third documented cases of human blood-associated dicistrovirus, and a highly divergent, unclassified dicistrovirus that we name human blood-associated dicistrovirus 2. We show that pegivirus C is a common co-infection in individuals with LF and is associated with lower Lassa viral loads and favorable outcomes. We help uncover the causes of three outbreaks as yellow fever virus, monkeypox virus, and a noninfectious cause, the latter ultimately determined to be pesticide poisoning. We demonstrate that a local, Nigerian-driven metagenomics response to complex public health scenarios generates accurate, real-time differential diagnoses, yielding insights that inform policy.

Use of Next-Generation Sequencing in a State-Wide Strategy of HIV-1 Surveillance: Impact of the SARS-COV-2 Pandemic on HIV-1 Diagnosis and Transmission.

BACKGROUND: The ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic posed an unpreceded threat to the management of other pandemics such as human immunodeficiency virus-1 (HIV-1) in the United States. The full impact of the SARS-CoV-2 pandemic on the HIV-1 pandemic needs to be evaluated. METHODS: All individuals with newly reported HIV-1 diagnoses from NC State Laboratory of Public Health were enrolled in this prospective observational study, 2018-2021. We used a sequencing-based recency assay to identify recent HIV-1 infections and to determine the days postinfection (DPI) for each person at the time of diagnosis. RESULTS: Sequencing used diagnostic serum samples from 814 individuals with new HIV-1 diagnoses spanning this 4-year period. Characteristics of individuals diagnosed in 2020 differed from those in other years. People of color diagnosed in 2021 were on average 6 months delayed in their diagnosis compared to those diagnosed in 2020. There was a trend that genetic networks were more known for individuals diagnosed in 2021. We observed no major integrase resistance mutations over the course of the study. CONCLUSIONS: SARS-CoV-2 pandemic may contribute to the spread of HIV-1. Public health resources need to focus on restoring HIV-1 testing and interrupting active, ongoing, transmission.

Model design for nonparametric phylodynamic inference and applications to pathogen surveillance.

Inference of effective population size from genomic data can provide unique information about demographic history and, when applied to pathogen genetic data, can also provide insights into epidemiological dynamics. The combination of nonparametric models for population dynamics with molecular clock models which relate genetic data to time has enabled phylodynamic inference based on large sets of time-stamped genetic sequence data. The methodology for nonparametric inference of effective population size is well-developed in the Bayesian setting, but here we develop a frequentist approach based on nonparametric latent process models of population size dynamics. We appeal to statistical principles based on out-of-sample prediction accuracy in order to optimize parameters that control shape and smoothness of the population size over time. Our methodology is implemented in a new R package entitled mlesky. We demonstrate the flexibility and speed of this approach in a series of simulation experiments and apply the methodology to a dataset of HIV-1 in the USA. We also estimate the impact of non-pharmaceutical interventions for COVID-19 in England using thousands of SARS-CoV-2 sequences. By incorporating a measure of the strength of these interventions over time within the phylodynamic model, we estimate the impact of the first national lockdown in the UK on the epidemic reproduction number.

Phylogenomic analysis uncovers a 9-year variation of Uganda influenza type-A strains from the WHO-recommended vaccines and other Africa strains.

Genetic characterisation of circulating influenza viruses directs annual vaccine strain selection and mitigation of infection spread. We used next-generation sequencing to locally generate whole genomes from 116 A(H1N1)pdm09 and 118 A(H3N2) positive patient swabs collected across Uganda between 2010 and 2018. We recovered sequences from 92% (215/234) of the swabs, 90% (193/215) of which were whole genomes. The newly-generated sequences were genetically and phylogenetically compared to the WHO-recommended vaccines and other Africa strains sampled since 1994. Uganda strain hemagglutinin (n = 206), neuraminidase (n = 207), and matrix protein (MP, n = 213) sequences had 95.23-99.65%, 95.31-99.79%, and 95.46-100% amino acid similarity to the 2010-2020 season vaccines, respectively, with several mutated hemagglutinin antigenic, receptor binding, and N-linked glycosylation sites. Uganda influenza type-A virus strains sequenced before 2016 clustered uniquely while later strains mixed with other Africa and global strains. We are the first to report novel A(H1N1)pdm09 subclades 6B.1A.3, 6B.1A.5(a,b), and 6B.1A.6 (± T120A) that circulated in Eastern, Western, and Southern Africa in 2017-2019. Africa forms part of the global influenza ecology with high viral genetic diversity, progressive antigenic drift, and local transmissions. For a continent with inadequate health resources and where social distancing is unsustainable, vaccination is the best option. Hence, African stakeholders should prioritise routine genome sequencing and analysis to direct vaccine selection and virus control.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses.

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Humoral and cellular immune responses to Lassa fever virus in Lassa fever survivors and their exposed contacts in Southern Nigeria.

Elucidating the adaptive immune characteristics of natural protection to Lassa fever (LF) is vital in designing and selecting optimal vaccine candidates. With rejuvenated interest in LF and a call for accelerated research on the Lassa virus (LASV) vaccine, there is a need to define the correlates of natural protective immune responses to LF. Here, we describe cellular and antibody immune responses present in survivors of LF (N = 370) and their exposed contacts (N = 170) in a LASV endemic region in Nigeria. Interestingly, our data showed comparable T cell and binding antibody responses from both survivors and their contacts, while neutralizing antibody responses were primarily seen in the LF survivors and not their contacts. Neutralizing antibody responses were found to be cross-reactive against all five lineages of LASV with a strong bias to Lineage II, the prevalent strain in southern Nigeria. We demonstrated that both T cell and antibody responses were not detectable in peripheral blood after a decade in LF survivors. Notably LF survivors maintained high levels of detectable binding antibody response for six months while their contacts did not. Lastly, as potential vaccine targets, we identified the regions of the LASV Glycoprotein (GP) and Nucleoprotein (NP) that induced the broadest peptide-specific T cell responses. Taken together this data informs immunological readouts and potential benchmarks for clinical trials evaluating LASV vaccine candidates.

Genome-wide mutational analysis of Chikungunya strains from 2016 to 2017 outbreak of central India: An attempt to elucidate the immunological basis for outbreak.

Chikungunya re-emerged in India in 2016-2017, as the first major outbreak since 2006. In our previous study, we undertook partial E1 gene sequencing and phylogenetic/mutational analysis of strains from the 2016-2017 outbreak of Chikungunya in central India and reported important mutations associated with the outbreak. This study was performed to validate the previous findings and to identify key mutations that had emerged throughout the entire genome of Chikungunya virus that could be driving the enormity of this outbreak. The phylogenetic analysis revealed the closeness of our isolates with ECSA genotype, specifically with the Singapore 2015 strain. We found 2 mutations in C and E2 genes, which were present in our isolates but were non-existent during the period of 2010-2016. Furthermore, re-emergence of Arg amino acid in place of stop codon in nsP3 gene and Thr at E2:312 positions was observed after 2011. We also used computational tools to assess the effect of the identified mutations on the T cell and B cell epitopes that could influence the protective immune response against this infection.

Pneumococcal within-host diversity during colonization, transmission and treatment.

Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.

Lessons from the pandemic: Responding to emerging zoonotic viral diseases-a Keystone Symposia report.

The COVID-19 pandemic caught the world largely unprepared, including scientific and policy communities. On April 10-13, 2022, researchers across academia, industry, government, and nonprofit organizations met at the Keystone symposium "Lessons from the Pandemic: Responding to Emerging Zoonotic Viral Diseases" to discuss the successes and challenges of the COVID-19 pandemic and what lessons can be applied moving forward. Speakers focused on experiences not only from the COVID-19 pandemic but also from outbreaks of other pathogens, including the Ebola virus, Lassa virus, and Nipah virus. A general consensus was that investments made during the COVID-19 pandemic in infrastructure, collaborations, laboratory and manufacturing capacity, diagnostics, clinical trial networks, and regulatory enhancements-notably, in low-to-middle income countries-must be maintained and strengthened to enable quick, concerted responses to future threats, especially to zoonotic pathogens.

Differences in epidemiology of enteropathogens in children pre- and post-rotavirus vaccine introduction in Kilifi, coastal Kenya.

BACKGROUND: Kenya introduced Rotarix® (GlaxoSmithKline Biologicals, Rixensart, Belgium) vaccination into its national immunization programme beginning July 2014. The impact of this vaccination program on the local epidemiology of various known enteropathogens is not fully understood. METHODS: We used a custom TaqMan Array Card (TAC) to screen for 28 different enteropathogens in 718 stools from children aged less than 13 years admitted to Kilifi County Hospital, coastal Kenya, following presentation with diarrhea in 2013 (before vaccine introduction) and in 2016-2018 (after vaccine introduction). Pathogen positivity rate differences between pre- and post-Rotarix® vaccination introduction were examined using both univariate and multivariable logistic regression models. RESULTS: In 665 specimens (92.6%), one or more enteropathogen was detected, while in 323 specimens (48.6%) three or more enteropathogens were detected. The top six detected enteropathogens were: enteroaggregative Escherichia coli (EAggEC; 42.1%), enteropathogenic Escherichia coli (EPEC; 30.2%), enterovirus (26.9%), rotavirus group A (RVA; 24.8%), parechovirus (16.6%) and norovirus GI/GII (14.4%). Post-rotavirus vaccine introduction, there was a significant increase in the proportion of samples testing positive for EAggEC (35.7% vs. 45.3%, p = 0.014), cytomegalovirus (4.2% vs. 9.9%, p = 0.008), Vibrio cholerae (0.0% vs. 2.3%, p = 0.019), Strongyloides species (0.8% vs. 3.6%, p = 0.048) and Dientamoeba fragilis (2.1% vs. 7.8%, p = 0.004). Although not reaching statistical significance, the positivity rate of adenovirus 40/41 (5.8% vs. 7.3%, p = 0.444), norovirus GI/GII (11.2% vs. 15.9%, p = 0.089), Shigella species (8.7% vs. 13.0%, p = 0.092) and Cryptosporidium spp. (11.6% vs. 14.7%, p = 0.261) appeared to increase post-vaccine introduction. Conversely, the positivity rate of sapovirus decreased significantly post-vaccine introduction (7.8% vs. 4.0%, p = 0.030) while that of RVA appeared not to change (27.4% vs. 23.5%, p = 0.253). More enteropathogen coinfections were detected per child post-vaccine introduction compared to before (mean: 2.7 vs. 2.3; p = 0.0025). CONCLUSIONS: In this rural Coastal Kenya setting, childhood enteropathogen infection burden was high both pre- and post-rotavirus vaccination introduction. Children who had diarrheal admissions post-vaccination showed an increase in coinfections and changes in specific enteropathogen positivity rates. This study highlights the utility of multipathogen detection platforms such as TAC in understanding etiology of childhood acute gastroenteritis in resource-limited regions.

Molecular Epidemiology and Evolutionary Dynamics of Human Influenza Type-A Viruses in Africa: A Systematic Review.

Genomic characterization of circulating influenza type-A viruses (IAVs) directs the selection of appropriate vaccine formulations and early detection of potentially pandemic virus strains. However, longitudinal data on the genomic evolution and transmission of IAVs in Africa are scarce, limiting Africa's benefits from potential influenza control strategies. We searched seven databases: African Journals Online, Embase, Global Health, Google Scholar, PubMed, Scopus, and Web of Science according to the PRISMA guidelines for studies that sequenced and/or genomically characterized Africa IAVs. Our review highlights the emergence and diversification of IAVs in Africa since 1993. Circulating strains continuously acquired new amino acid substitutions at the major antigenic and potential N-linked glycosylation sites in their hemagglutinin proteins, which dramatically affected vaccine protectiveness. Africa IAVs phylogenetically mixed with global strains forming strong temporal and geographical evolution structures. Phylogeographic analyses confirmed that viral migration into Africa from abroad, especially South Asia, Europe, and North America, and extensive local viral mixing sustained the genomic diversity, antigenic drift, and persistence of IAVs in Africa. However, the role of reassortment and zoonosis remains unknown. Interestingly, we observed substitutions and clades and persistent viral lineages unique to Africa. Therefore, Africa's contribution to the global influenza ecology may be understated. Our results were geographically biased, with data from 63% (34/54) of African countries. Thus, there is a need to expand influenza surveillance across Africa and prioritize routine whole-genome sequencing and genomic analysis to detect new strains early for effective viral control.

Whole-genome analysis to determine the rate and patterns of intra-subtype reassortment among influenza type-A viruses in Africa.

Influenza type-A viruses (IAVs) present a global burden of human respiratory infections and mortality. Genome reassortment is an important mechanism through which epidemiologically novel influenza viruses emerge and a core step in the safe reassortment-incompetent live-attenuated influenza vaccine development. Currently, there are no data on the rate, spatial and temporal distribution, and role of reassortment in the evolution and diversification of IAVs circulating in Africa. We aimed to detect intra-subtype reassortment among Africa pandemic H1N1pdm09 (2009-10), seasonal H1N1pdm09 (2011-20), and seasonal H3N2 viruses and characterize the genomic architecture and temporal and spatial distribution patterns of the resulting reassortants. Our study was nested within the Uganda National Influenza Surveillance Programme. Next-generation sequencing was used to generate whole genomes (WGs) from 234 H1N1pdm09 (n = 116) and H3N2 (n = 118) viruses sampled between 2010 and 2018 from seven districts in Uganda. We combined our newly generated WGs with 658 H1N1pdm09 and 1131 H3N2 WGs sampled between 1994 and 2020 across Africa and identified reassortants using an automated Graph Incompatibility Based Reassortment Finder software. Viral reassortment rates were estimated using a coalescent reassortant constant population model. Phylogenetic analysis was used to assess the effect of reassortment on viral genetic evolution. We observed a high frequency of intra-subtype reassortment events, 12 · 4 per cent (94/758) and 20 · 9 per cent (256/1,224), and reassortants, 13 · 3 per cent (101/758) and 38 · 6 per cent (472/1,224), among Africa H1N1pdm09 and H3N2 viruses, respectively. H1N1pdm09 reassorted at higher rates (0.1237-0.4255) than H3N2 viruses (0 · 00912-0.0355 events/lineage/year), a case unique to Uganda. Viral reassortants were sampled in 2009 through 2020, except in 2012. 78 · 2 per cent (79/101) of H1N1pdm09 reassortants acquired new non-structural, while 57 · 8 per cent (273/472) of the H3N2 reassortants had new hemagglutinin (H3) genes. Africa H3N2 viruses underwent more reassortment events involving larger reassortant sets than H1N1pdm09 viruses. Viruses with a specific reassortment architecture circulated for up to five consecutive years in specific countries and regions. The Eastern (Uganda and Kenya) and Western Africa harboured 84 · 2 per cent (85/101) and 55 · 9 per cent (264/472) of the continent's H1N1pdm09 and H3N2 reassortants, respectively. The frequent reassortment involving multi-genes observed among Africa IAVs showed the intracontinental viral evolution and diversification possibly sustained by viral importation from outside Africa and/or local viral genomic mixing and transmission. Novel reassortant viruses emerged every year, and some persisted in different countries and regions, thereby presenting a risk of influenza outbreaks in Africa. Our findings highlight Africa as part of the global influenza ecology and the advantage of implementing routine whole-over partial genome sequencing and analyses to monitor circulating and detect emerging viruses. Furthermore, this study provides evidence and heightens our knowledge on IAV evolution, which is integral in directing vaccine strain selection and the update of master donor viruses used in recombinant vaccine development.

Intersection of Syphilis and Human Immunodeficiency Virus (HIV) Networks to Identify Opportunities to Enhance HIV Prevention.

BACKGROUND: Human immunodeficiency virus (HIV) and syphilis infection continue at disproportionate rates among minority men who have sex with men (MSM) in the United States. The integration of HIV genetic clustering with partner services can provide important insight into local epidemic trends to guide interventions and control efforts. METHODS: We evaluated contact networks of index persons defined as minority men and transgender women diagnosed with early syphilis and/or HIV infection between 2018 and 2020 in 2 North Carolina regions. HIV clusters were constructed from pol sequences collected through statewide surveillance. A combined "HIV-risk" network, which included persons with any links (genetic or sexual contact) to HIV-positive persons, was evaluated by component size, demographic factors, and HIV viral suppression. RESULTS: In total, 1289 index persons were identified and 55% named 1153 contacts. Most index persons were Black (88%) and young (median age 30 years); 70% had early syphilis and 43% had prevalent HIV infection. Most people with HIV (65%) appeared in an HIV cluster. The combined HIV-risk network (1590 contact network and 1500 cluster members) included 287 distinct components; however, 1586 (51%) were in a single component. Fifty-five percent of network members with HIV had no evidence of viral suppression. Overall, fewer index persons needed to be interviewed to identify 1 HIV-positive member without viral suppression (1.3 vs 4.0 for contact tracing). CONCLUSIONS: Integration of HIV clusters and viral loads illuminate networks with high HIV prevalence, indicating recent and ongoing transmission. Interventions intensified toward these networks may efficiently reach persons for HIV prevention and care re-engagement.

VGEA: an RNA viral assembly toolkit.

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user's choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

Molecular epidemiology of dengue in a setting of low reported endemicity: Kupang, East Nusa Tenggara province, Indonesia.

BACKGROUND: Most regions in Indonesia experience annual dengue epidemics. However, the province of East Nusa Tenggara has consistently reported low incidence. We conducted a dengue molecular epidemiology study in Kupang, the capital of the province. METHODS: Dengue patients were recruited from May 2016 to September 2017. Dengue virus (DENV) screening was performed using NS1 and immunoglobulin G (IgG)/IgM detection. Serotype was determined using reverse transcription polymerase chain reaction and the envelope genes were sequenced to infer the genetic identity and phylogeny. RESULTS: From 119 patients, dengue was confirmed in 62 (52%). Compared with official data, underreporting of dengue incidence was observed. The majority (36%) of patients were children <10 y of age. Most patients (80%) experienced mild fever. All serotypes were detected, with DENV-3 as the predominant (57%). Kupang DENV-1 isolate was classified as genotype IV, an old and endemic strain, DENV-2 as cosmopolitan, DENV-3 as genotype I and DENV-4 as genotype II. Most isolates showed relatively low evolutionary rates and are closely related with strains from Bali and Timor Leste. CONCLUSIONS: The low dengue incidence was most likely caused by sustained local circulation of endemic viruses. This study provides information on the epidemiology of dengue in a low-endemicity setting that should help future mitigation and disease management.

HIV-1 Transmission linkages among persons with incident infection to inform public health surveillance.

BACKGROUND: We evaluated features of HIV transmission networks involving persons diagnosed during incident HIV infection (IHI) to assess network-based opportunities to curtail onward transmission. METHODS: Transmission networks were constructed using partial pol sequences reported to North Carolina surveillance among persons with recent (2014-2018) and past (<2014) HIV diagnoses. IHI were defined as documented acute infections or seroconversion. Demographic and virologic features of HIV genetic clusters (<1.5% pairwise genetic distance) involving ≥ 1 IHI were assessed. Persons with viral genetic links and who had diagnoses >90 days prior to an IHI were further characterized. We assessed named partner outcomes among IHI index persons using contact tracing data. FINDINGS: Of 4,405 HIV diagnoses 2014-2018 with sequences, there were 323 (7%) IHI index persons; most were male (88%), Black (65%), young (68% <30 years), and reported sex with men (MSM) risk (79%). Index persons were more likely to be cluster members compared to non-index persons diagnosed during the same period (72% vs. 49%). In total, 162 clusters were identified involving 233 IHI, 577 recent diagnoses, and 163 past diagnoses. Most IHI cases (53%) had viral linkages to ≥1 previously diagnosed person without evidence of HIV viral suppression in the year prior to the diagnosis of the IHI index. In contact tracing, only 53% IHI cases named an HIV-positive contact, resulting in 0.5 previously diagnosed persons detected per IHI investigated. When combined with viral analyses, the detection rate of viremic previously diagnosed persons increased to 1.3. INTERPRETATION: Integrating public health with molecular epidemiology, revealed that more than half of IHI have viral links to persons with previously diagnosed unsuppressed HIV infection which was largely unrecognized by traditional contact tracing. Enhanced partner services to support engagement and retention in HIV care and improved case finding supported by rapid phylogenetic analysis are tools to substantially reduce onward HIV transmission.

Phylogenetic Characterization of Crimean-Congo Hemorrhagic Fever Virus Detected in African Blue Ticks Feeding on Cattle in a Ugandan Abattoir.

Crimean-Congo hemorrhagic fever virus (CCHFV) is the most geographically widespread of the tick-borne viruses. However, African strains of CCHFV are poorly represented in sequence databases. In addition, almost all sequence data collected to date have been obtained from cases of human disease, while information regarding the circulation of the virus in tick and animal reservoirs is severely lacking. Here, we characterize the complete coding region of a novel CCHFV strain, detected in African blue ticks (Rhipicephalus (Boophilus) decoloratus) feeding on cattle in an abattoir in Kampala, Uganda. These cattle originated from a farm in Mbarara, a major cattle-trading hub for much of Uganda. Phylogenetic analysis indicates that the newly sequenced strain belongs to the African genotype II clade, which predominantly contains the sequences of strains isolated from West Africa in the 1950s, and South Africa in the 1980s. Whilst the viral S (nucleoprotein) and L (RNA polymerase) genome segments shared >90% nucleotide similarity with previously reported genotype II strains, the glycoprotein-coding M segment shared only 80% nucleotide similarity with the next most closely related strains, which were derived from ticks in Western India and Northern China. This genome segment also displayed a large number of non-synonymous mutations previously unreported in the genotype II strains. Characterization of this novel strain adds to our limited understanding of the natural diversity of CCHFV circulating in both ticks and in Africa. Such data can be used to inform the design of vaccines and diagnostics, as well as studies exploring the epidemiology and evolution of the virus for the establishment of future CCHFV control strategies.