Antigenic cartography using hamster sera identifies SARS-CoV-2 JN.1 evasion seen in human XBB.1.5 booster sera.

Antigenic assessments of SARS-CoV-2 variants inform decisions to update COVID-19 vaccines. Primary infection sera are often used for assessments, but such sera are rare due to population immunity from SARS-CoV-2 infections and COVID-19 vaccinations. Here, we show that neutralization titers and breadth of matched human and hamster pre-Omicron variant primary infection sera correlate well and generate similar antigenic maps. The hamster antigenic map shows modest antigenic drift among XBB sub-lineage variants, with JN.1 and BA.4/BA.5 variants within the XBB cluster, but with five to six-fold antigenic differences between these variants and XBB.1.5. Compared to sera following only ancestral or bivalent COVID-19 vaccinations, or with post-vaccination infections, XBB.1.5 booster sera had the broadest neutralization against XBB sub-lineage variants, although a five-fold titer difference was still observed between JN.1 and XBB.1.5 variants. These findings suggest that antibody coverage of antigenically divergent JN.1 could be improved with a matched vaccine antigen.

Antigenic distance between primary and secondary dengue infections correlates with disease risk.

Many pathogens continuously change their protein structure in response to immune-driven selection, resulting in weakened protection even in previously exposed individuals. In addition, for some pathogens, such as dengue virus, poorly targeted immunity is associated with increased risk of severe disease through a mechanism known as antibody-dependent enhancement. However, it remains unclear whether the antigenic distances between an individual's first infection and subsequent exposures dictate disease risk, explaining the observed large-scale differences in dengue hospitalizations across years. Here, we develop a framework that combines detailed antigenic and genetic characterization of viruses with details on hospitalized cases from 21 years of dengue surveillance in Bangkok, Thailand, to identify the role of the antigenic profile of circulating viruses in determining disease risk. We found that the risk of hospitalization depended on both the specific order of infecting serotypes and the antigenic distance between an individual's primary and secondary infections, with risk maximized at intermediate antigenic distances. These findings suggest that immune imprinting helps determine dengue disease risk and provide a pathway to monitor the changing risk profile of populations and to quantifying risk profiles of candidate vaccines.

Metagenomics analysis reveals presence of the Merida-like virus in Georgia.

Arbovirus surveillance is fundamental for the discovery of novel viruses and prevention of febrile vector-borne illnesses. Vector-borne pathogens can rapidly expand and adapt in new geographic and environmental conditions. In this study, metagenomic surveillance was conducted to identify novel viruses in the Country of Georgia. A total of 521 mosquitoes were captured near a military training facility and pooled from species Culex pipiens (Linnaeus) (87%) and Aedes albopictus (Skuse) (13%). We decided to further analyze the Culex pipiens mosquitoes, due to the more extensive number of samples collected. Our approach was to utilize an unbiased total RNA-seq for pathogen discovery in order to explore the mosquito virome. The viral reads from this analysis were mostly aligned to Insect-specific viruses from two main families, the Iflaviridae; a positive-stranded RNA virus and the Rhabdoviridae; a negative- and single-stranded RNA virus. Our pathogen discovery analysis revealed viral reads aligning to the Merida-like virus Turkey (MERDLVT) strain among the Rhabdoviridae. To further validate this result, we conducted a BLAST sequence comparison analysis of our samples with the MERDLVT strain. Our positive samples aligned to the MERDLVT strain with 96-100% sequence identity and 99.7-100% sequence coverage. A bootstrapped maximum-likelihood phylogenetic tree was used to evaluate the evolutionary relationships among these positive pooled specimens with the (MERDLVT) strain. The Georgia samples clustered most closely with two strains from Turkey, the Merida-like virus KE-2017a isolate 139-1-21 and the Merida-like virus Turkey isolate P431. Collectively, these results show the presence of the MERDLVT strain in Georgia.

Antigenic diversity and dengue disease risk.

Many pathogens continuously change their protein structure in response to immune-driven selection, resulting in weakened protection. In addition, for some pathogens such as dengue virus, poorly targeted immunity is associated with increased risk of severe disease, through a mechanism known as antibody-dependent enhancement. However, it remains a mystery whether the antigenic distance between an individual's first infection and subsequent exposures dictate disease risk, explaining the observed large-scale differences in dengue hospitalisations across years. Here we develop an inferential framework that combines detailed antigenic and genetic characterisation of viruses, and hospitalised cases from 21 years of surveillance in Bangkok, Thailand to identify the role of the antigenic profile of circulating viruses in determining disease risk. We find that the risk of hospitalisation depends on both the specific order of infecting serotypes and the antigenic distance between an individual's primary and secondary infections, with risk maximised at intermediate antigenic distances. These findings suggest immune imprinting helps determine dengue disease risk, and provides a pathway to monitor the changing risk profile of populations and to quantifying risk profiles of candidate vaccines.

Transmission of SARS-CoV-2 among recruits in a US Army training environment: a brief report.

BACKGROUND: In 2020, preventive measures were implemented to mitigate the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) among 600-700 recruits arriving weekly at a basic combat training (BCT) facility in the southern United States. Trainees were sorted into companies and platoons (cocoons) at arrival, tested, quarantined for 14 days with daily temperature and respiratory-symptom monitoring and retested before release into larger groups for training where symptomatic testing was conducted. Nonpharmaceutical measures, such as masking, and social distancing, were maintained throughout quarantine and BCT. We assessed for SARS-CoV-2 transmission in the quarantine milieu. METHODS: Nasopharyngeal (NP) swabs were collected at arrival and at the end of quarantine and blood specimens at both timepoints and at the end of BCT. Epidemiological characteristics were analyzed for transmission clusters identified from whole-genome sequencing of NP samples. RESULTS: Among 1403 trainees enrolled from 25 August to 7 October 2020, epidemiological analysis identified three transmission clusters (n = 20 SARS-CoV-2 genomes) during quarantine, which spanned five different cocoons. However, SARS-CoV-2 incidence decreased from 2.7% during quarantine to 1.5% at the end of BCT; prevalence at arrival was 3.3%. CONCLUSIONS: These findings suggest layered SARS-CoV-2 mitigation measures implemented during quarantine minimized the risk of further transmission in BCT.

Clustered cases of human adenovirus types 4, 7, and 14 infections in US Department of Defense Beneficiaries during the 2018-2019 season.

Human adenoviruses (HAdV) are genetically diverse and can infect a number of tissues with severities varied from mild to fatal. HAdV types 3, 4, 7, 11, 14, 21, and 55 were associated with acute respiratory illnesses outbreaks in the United States and in other countries. The risk of outbreaks can be effectively controlled by HAdV vaccination or mitigated by screening and preventive measures. During the influenza season 2018-2019, the DoD Global Respiratory Pathogen Surveillance Program (DoDGRS) received 24 300 respiratory specimens. HAdV samples that produced positive cytopathic effects in viral cultivation were subjected to next-generation sequencing for genome sequence assembly, genome typing, whole genome phylogeny, and sequence comparative analyses. A variety of HAdV types were identified in this study, including HAdV types 1-7, 14, 55, and 56. HAdV types 4, 7, and 14 were found in clustered cases in Colorado, Florida, New York, and South Carolina. Comparative sequence analyses of these isolates revealed the emergence of novel genetic mutations despite the stability of adenovirus genomes. Genomic surveillance of HAdV suggested possible undetected outbreaks and shed light on prevalence, genetic divergence, and viral evolution of HAdV. Continued surveillance will inform risk assessment and countermeasures.

Intrahost SARS-CoV-2 k-mer Identification Method (iSKIM) for Rapid Detection of Mutations of Concern Reveals Emergence of Global Mutation Patterns.

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM ("intrahost SARS-CoV-2 k-mer identification method") to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

Discovery of Rickettsia spp. in mosquitoes collected in Georgia by metagenomics analysis and molecular characterization.

Arthropods have a broad and expanding worldwide presence and can transmit a variety of viral, bacterial, and parasite pathogens. A number of Rickettsia and Orientia species associated with ticks, fleas, lice, and mites have been detected in, or isolated from, patients with febrile illness and/or animal reservoirs throughout the world. Mosquitoes are not currently considered vectors for Rickettsia spp. pathogens to humans or to animals. In this study, we conducted a random metagenome next-generation sequencing (NGS) of 475 pools of Aedes, Culex, and Culiseta species of mosquitoes collected in Georgia from 2018 to 2019, identifying rickettsial gene sequences in 33 pools of mosquitoes. We further confirmed the findings of the Rickettsia by genus-specific quantitative PCR (qPCR) and multi-locus sequence typing (MLST). The NGS and MLST results indicate that Rickettsia spp. are closely related to Rickettsia bellii, which is not known to be pathogenic in humans. The results, together with other reports of Rickettsia spp. in mosquitoes and the susceptibility and transmissibility experiments, suggest that mosquitoes may play a role in the transmission cycle of Rickettsia spp.

Intrahost SARS-CoV-2 k-mer identification method (iSKIM) for rapid detection of mutations of concern reveals emergence of global mutation patterns.

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM ("intrahost SARS-CoV-2 k-mer identification method") to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

Genomic surveillance of SARS-CoV-2 in US military compounds in Afghanistan reveals multiple introductions and outbreaks of Alpha and Delta variants.

BACKGROUND: With the emergence and spread of SARS-CoV-2 variants, genomic epidemiology and surveillance have proven invaluable tools for variant tracking. Here, we analyzed SARS-CoV-2 samples collected from personnel located at the US/NATO bases across Afghanistan. RESULTS: Sequencing and phylogenetic analyses revealed at least 16 independent introductions of SARS-CoV-2 into four of these relatively isolated compounds during April and May 2021, including multiple introductions of Alpha and Delta variants. Four of the introductions resulted in sustained spread of the virus within, and in two cases between, the compounds. Three of these outbreaks, one Delta and two Alpha, occurred simultaneously. CONCLUSIONS: Even in rigorously controlled and segregated environments, SARS-CoV-2 introduction and spread may occur frequently.

New Sex Chromosomes in Lake Victoria Cichlid Fishes (Cichlidae: Haplochromini).

African cichlid fishes harbor an extraordinary diversity of sex-chromosome systems. Within just one lineage, the tribe Haplochromini, at least 6 unique sex-chromosome systems have been identified. Here we focus on characterizing sex chromosomes in cichlids from the Lake Victoria basin. In Haplochromis chilotes, we identified a new ZW system associated with the white blotch color pattern, which shows substantial sequence differentiation over most of LG16, and is likely to be present in related species. In Haplochromis sauvagei, we found a coding polymorphism in amh that may be responsible for an XY system on LG23. In Pundamilia nyererei, we identified a feminizing effect of B chromosomes together with XY- and ZW-patterned differentiation on LG23. In Haplochromis latifasciatus, we identified a duplication of amh that may be present in other species of the Lake Victoria superflock. We further characterized the LG5-14 XY system in Astatotilapia burtoni and identified the oldest stratum on LG14. This species also showed ZW differentiation on LG2. Finally, we characterized an XY system on LG7 in Astatoreochromis alluaudi. This report brings the number of distinct sex-chromosome systems in haplochromine cichlids to at least 13, and highlights the dynamic evolution of sex determination and sex chromosomes in this young lineage.

Investigation of an outbreak of COVID-19 among U.S. military personnel and beneficiaries stationed in the Republic of Korea, June-July 2021.

On 28 May 2021, leisure travel restrictions in place to control coronavirus disease 2019 (COVID-19) were eased among vaccinated U.S. military personnel and beneficiaries stationed in South Korea (USFK) allowing access to bars and clubs which were off limits. We describe results from an investigation of the largest severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as of November 2021 among USFK personnel following this change in policy. Data such as SARS-CoV-2 real-time polymerase chain reaction (RT-PCR) test results, demographic characteristics, symptom and vaccination histories, and genome sequences were analyzed. Of a total 207 new cases of COVID-19 diagnosed among USFK members from 15 June to 27 July 2021, 113 (57%) eligible cases were fully vaccinated, of whom 86 (76%) were symptomatic. RT-PCR cycling threshold values were similar among vaccinated and unvaccinated members. Whole genomic sequencing of 54 outbreak samples indicated all infections were due to the Delta variant. Phylogenetic analysis revealed two sources of SARS-CoV-2 accounted for 41% of infections among vaccinated and unvaccinated members. Vaccinated personnel were not at risk of severe illness; however, 86% experienced symptoms following infection. There were no hospitalizations among COVID-19 cases, most of whom were young military service members. Rescinded restrictions were reinstated to control the outbreak. Masking was mandated among all personnel predating U.S. national recommendations for indoor masking in high COVID-19 transmission areas. Increased vaccination with continued vigilance and extension of COVID-19 mitigation measures are warranted to contain the spread of SARS-CoV-2 variants of concern.

A comparative recombination analysis of human coronaviruses and implications for the SARS-CoV-2 pandemic.

The SARS-CoV-2 pandemic prompts evaluation of recombination in human coronavirus (hCoV) evolution. We undertook recombination analyses of 158,118 public seasonal hCoV, SARS-CoV-1, SARS-CoV-2 and MERS-CoV genome sequences using the RDP4 software. We found moderate evidence for 8 SARS-CoV-2 recombination events, two of which involved the spike gene, and low evidence for one SARS-CoV-1 recombination event. Within MERS-CoV, 229E, OC43, NL63 and HKU1 datasets, we noted 7, 1, 9, 14, and 1 high-confidence recombination events, respectively. There was propensity for recombination breakpoints in the non-ORF1 region of the genome containing structural genes, and recombination severely skewed the temporal structure of these data, especially for NL63 and OC43. Bayesian time-scaled analyses on recombinant-free data indicated the sampled diversity of seasonal CoVs emerged in the last 70 years, with 229E displaying continuous lineage replacements. These findings emphasize the importance of genomic based surveillance to detect recombination in SARS-CoV-2, particularly if recombination may lead to immune evasion.

Origin of a Giant Sex Chromosome.

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

Movement of transposable elements contributes to cichlid diversity.

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia.

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y.

Tbx2a Modulates Switching of RH2 and LWS Opsin Gene Expression.

Sensory systems are tuned by selection to maximize organismal fitness in particular environments. This tuning has implications for intraspecies communication, the maintenance of species boundaries, and speciation. Tuning of color vision largely depends on the sequence of the expressed opsin proteins. To improve tuning of visual sensitivities to shifts in habitat or foraging ecology over the course of development, many organisms change which opsins are expressed. Changes in this developmental sequence (heterochronic shifts) can create differences in visual sensitivity among closely related species. The genetic mechanisms by which these developmental shifts occur are poorly understood. Here, we use quantitative trait locus analyses, genome sequencing, and gene expression studies in African cichlid fishes to identify a role for the transcription factor Tbx2a in driving a switch between long wavelength sensitive (LWS) and Rhodopsin-like (RH2) opsin expression. We identify binding sites for Tbx2a in the LWS promoter and the highly conserved locus control region of RH2 which concurrently promote LWS expression while repressing RH2 expression. We also present evidence that a single change in Tbx2a regulatory sequence has led to a species difference in visual tuning, providing the first mechanistic model for the evolution of rapid switches in sensory tuning. This difference in visual tuning likely has important roles in evolution as it corresponds to differences in diet, microhabitat choice, and male nuptial coloration.

Diurnal variation in opsin expression and common housekeeping genes necessitates comprehensive normalization methods for quantitative real-time PCR analyses.

To determine the visual sensitivities of an organism of interest, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is often used to quantify expression of the light-sensitive opsins in the retina. While qRT-PCR is an affordable, high-throughput method for measuring expression, it comes with inherent normalization issues that affect the interpretation of results, especially as opsin expression can vary greatly based on developmental stage, light environment or diurnal cycles. We tested for diurnal cycles of opsin expression over a period of 24 hr at 1-hr increments and examined how normalization affects a data set with fluctuating expression levels using qRT-PCR and transcriptome data from the retinae of the cichlid Pelmatolapia mariae. We compared five methods of normalizing opsin expression relative to (a) the average of three stably expressed housekeeping genes (Ube2z, EF1-α and ß-actin), (b) total RNA concentration, (c) GNAT2, (the cone-specific subunit of transducin), (d) total opsin expression and (e) only opsins expressed in the same cone type. Normalizing by proportion of cone type produced the least variation and would be best for removing time-of-day variation. In contrast, normalizing by housekeeping genes produced the highest daily variation in expression and demonstrated that the peak of cone opsin expression was in the late afternoon. A weighted correlation network analysis showed that the expression of different cone opsins follows a very similar daily cycle. With the knowledge of how these normalization methods affect opsin expression data, we make recommendations for designing sampling approaches and quantification methods based upon the scientific question being examined.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes.

BACKGROUND: African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages. RESULTS: We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage. CONCLUSION: This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.