DNA Barcoding Analysis of Trinidad Haemagogus Mosquitoes Reveals Evidence for Putative New Species

Background: Haemagogus janthinomys is a primary sylvan vector of yellow fever virus and the emerging Mayaro virus. However, despite its medical importance, there is a dearth of data on the molecular taxonomy of this mosquito species. Methods: In this study, DNA barcoding analysis was performed on 64 adult female mosquitoes from Trinidad morphologically identified as Hg. janthinomys. The mitochondrial cytochrome c oxidase I (COI) gene and ribosomal DNA internal transcribed spacer 2 (ITS2) region of the mosquitoes were PCR amplified and sequenced, and molecular phylogenies inferred. Results: The BLASTN analysis showed that only 20% (n = 13/66) of COI sequences had high similarity (>99% identity) to Hg. janthinomys and the remaining sequences had low similarity (<90% identity) to reference GenBank sequences. Phylogenetic analysis of COI sequences revealed the presence of four strongly supported groups, with one distinct clade that did not align with any reference sequences. Corresponding ITS2 sequences for samples in this distinct COI group clustered into three clades. Conclusions: These molecular findings suggest the existence of a putative new Haemagogus mosquito species and underscore the need for further, more in-depth investigations into the taxonomy and classification of the Haemagogus genus.

Correction: Implementation of Genomic Surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings

[This corrects the article DOI: 10.1371/journal.pgph.0001455.].

Toward a global virus genomic surveillance network

The COVID-19 pandemic galvanized the field of virus genomic surveillance, demonstrating its utility for public health. Now, we must harness the momentum that led to increased infrastructure, training, and political will to build a sustainable global genomic surveillance network for other epidemic and endemic viruses. We suggest a generalizable modular sequencing framework wherein users can easily switch between virus targets to maximize cost-effectiveness and maintain readiness for new threats. We also highlight challenges associated with genomic surveillance and when global inequalities persist. We propose solutions to mitigate some of these issues, including training and multilateral partnerships. Exploring alternatives to clinical sequencing can also reduce the cost of surveillance programs. Finally, we discuss how establishing genomic surveillance would aid control programs and potentially provide a warning system for outbreaks, using a global respiratory virus (RSV), an arbovirus (dengue virus), and a regional zoonotic virus (Lassa virus) as examples.

Implementation of genomic surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings

The COVID-19 pandemic highlighted the importance of global genomic surveillance to monitor the emergence and spread of SARS-CoV-2 variants and inform public health decision-making. Until December 2020 there was minimal capacity for viral genomic surveillance in most Caribbean countries. To overcome this constraint, the COVID-19: Infectious disease Molecular epidemiology for PAthogen Control & Tracking (COVID-19 IMPACT) project was implemented to establish rapid SARS-CoV-2 whole genome nanopore sequencing at The University of the West Indies (UWI) in Trinidad and Tobago (T&T) and provide needed SARS-CoV-2 sequencing services for T&T and other Caribbean Public Health Agency Member States (CMS). Using the Oxford Nanopore Technologies MinION sequencing platform and ARTIC network sequencing protocols and bioinformatics pipeline, a total of 3610 SARS-CoV-2 positive RNA samples, received from 17 CMS, were sequenced in-situ during the period December 5th 2020 to December 31st 2021. Ninety-one Pango lineages, including those of five variants of concern (VOC), were identified. Genetic analysis revealed at least 260 introductions to the CMS from other global regions. For each of the 17 CMS, the percentage of reported COVID-19 cases sequenced by the COVID-19 IMPACT laboratory ranged from 0·02% to 3·80% (median = 1·12%). Sequences submitted to GISAID by our study represented 73·3% of all SARS-CoV-2 sequences from the 17 CMS available on the database up to December 31st 2021. Increased staffing, process and infrastructural improvement over the course of the project helped reduce turnaround times for reporting to originating institutions and sequence uploads to GISAID. Insights from our genomic surveillance network in the Caribbean region directly influenced non-pharmaceutical countermeasures in the CMS countries. However, limited availability of associated surveillance and clinical data made it challenging to contextualise the observed SARS-CoV-2 diversity and evolution, highlighting the need for development of infrastructure for collecting and integrating genomic sequencing data and sample-associated metadata.

Phylogenetic characterization of Orthobunyaviruses isolated from Trinidad shows evidence of natural reassortment

The genus Orthobunyavirus is a diverse group of viruses in the family Peribunyaviridae, recently classified into 20 serogroups, and 103 virus species. Although most viruses within these serogroups are phylogenetically distinct, the absence of complete genome sequences has left several viruses incompletely characterized. Here we report the complete genome sequences for 11 orthobunyaviruses isolated from Trinidad, French Guiana, Guatemala, and Panama that were serologically classified into six serogroups and 10 species. Phylogenetic analyses of these 11 newly derived sequences indicate that viruses belonging to the Patois, Capim, Guama, and Group C serocomplexes all have a close genetic origin. We show that three of the 11 orthobunyaviruses characterized (belonging to the Group C and Bunyamwera serogroups) have evidence of histories of natural reassortment through the M genome segment. Our data also suggests that two distinct lineages of Group C viruses concurrently circulate in Trinidad and are transmitted by the same mosquito vectors. This study also highlights the importance of complementing serological identification with nucleotide sequencing when characterizing orthobunyaviruses.

Genomics-informed outbreak investigations of SARS-CoV-2 using civet

The scale of data produced during the SARS-CoV-2 pandemic has been unprecedented, with more than 13 million sequences shared publicly at the time of writing. This wealth of sequence data provides important context for interpreting local outbreaks. However, placing sequences of interest into national and international context is difficult given the size of the global dataset. Often outbreak investigations and genomic surveillance efforts require running similar analyses again and again on the latest dataset and producing reports. We developed civet (cluster investigation and virus epidemiology tool) to aid these routine analyses and facilitate virus outbreak investigation and surveillance. Civet can place sequences of interest in the local context of background diversity, resolving the query into different 'catchments' and presenting the phylogenetic results alongside metadata in an interactive, distributable report. Civet can be used on a fine scale for clinical outbreak investigation, for local surveillance and cluster discovery, and to routinely summarise the virus diversity circulating on a national level. Civet reports have helped researchers and public health bodies feedback genomic information in the appropriate context within a timeframe that is useful for public health.

Saliva-based methods for SARS-CoV-2 testing in low- and middle-income countries

As the coronavirus disease 2019 (COVID-19) continues to disproportionately affect low- and middle-income countries, the need for simple, accessible and frequent diagnostic testing grows. In lower-resource settings, case detection is often limited by a lack of available testing for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To address global inequities in testing, alternative sample types could be used to increase access to testing by reducing the associated costs. Saliva is a sensitive, minimally invasive and inexpensive diagnostic sample for SARS-CoV-2 detection that is appropriate for asymptomatic surveillance, symptomatic testing and at-home collection. Saliva testing can lessen two major challenges faced by lower- and middle-income countries: constrained resources and overburdened health workers. Saliva sampling enables convenient self-collection and requires fewer resources than swab-based methods. However, saliva testing for SARS-CoV-2 diagnostics has not been implemented on a large scale in low- and middle-income countries. While numerous studies based in these settings have demonstrated the usefulness of saliva sampling, there has been insufficient attention on optimizing its implementation in practice. We argue that implementation science research is needed to bridge this gap between evidence and practice. Low- and middle-income countries face many barriers as they continue their efforts to provide mass COVID-19 testing in the face of substantial inequities in global access to vaccines. Laboratories should look to replicate successful approaches for sensitive detection of SARS-CoV-2 in saliva, while governments should act to facilitate mass testing by lifting restrictions that limit implementation of saliva-based methods.

Comparison of surveillance trapping methods to monitor Culicoides biting midge activity in Trinidad, West Indies

Culicoides biting midges (Diptera: Ceratopogonidae) are biting nuisances and arbovirus vectors of both public health and veterinary significance in Trinidad. We compared sampling methods to define the behaviour and bionomics of adult Culicoides populations at a commercial dairy goat farm. Three static trap designs were compared: (a) Centre for Disease Control (CDC) downdraft UV trap; (b) CDC trap with an incandescent bulb and (c) CDC trap with semiochemical lure consisting of R-(-)-1-octen-3-ol and CO₂ (no bulb). Sweep netting was used to define diel periodicity. A total of 30,701 biting midges were collected using static traps, dominated by female Culicoides furens (>70% of trap collections across all three designs). There was no significant difference in the Margalef's index between the three traps; however, trap designs A and C collected a significantly greater number of individuals than trap B, and trap C gained highest species richness. The greatest species richness and abundance of Culicoides collected by sweep net was observed between 6:00 and 6:15 pm and notable differences in the crepuscular activity pattern of several species were identified. Comparative data on Culicoides species richness, abundance, sex and reproductive status is discussed and can be used to improve surveillance strategies, research designs and risk management.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

A high-throughput screening assay to identify inhibitory antibodies targeting alphavirus release

BACKGROUND

Several studies have demonstrated neutralizing antibodies to be highly effective against alphavirus infection in animal models, both prophylactically and remedially. In most studies, neutralizing antibodies have been evaluated for their ability to block viral entry in vitro but recent evidence suggests that antibody inhibition through other mechanisms, including viral budding/release, significantly contributes to viral control in vivo for a number of alphaviruses.

RESULTS

We describe a BSL-2, cell-based, high-throughput screening system that specifically screens for inhibitors of alphavirus egress using chikungunya virus (CHIKV) and Mayaro virus (MAYV) novel replication competent nano-luciferase (nLuc) reporter viruses. Screening of both polyclonal sera and memory B-cell clones from CHIKV immune individuals using the optimized assay detected several antibodies that display potent anti-budding activity.

CONCLUSIONS

We describe an "anti-budding assay" to specifically screen for inhibitors of viral egress using novel CHIKV and MAYV nLuc reporter viruses. This BSL-2 safe, high-throughput system can be utilized to explore neutralizing "anti-budding" antibodies to yield potent candidates for CHIKV and MAYV therapeutics and prophylaxis.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing provides critical information to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments and vaccines, and guide public health responses. To investigate the spatiotemporal heterogeneity in the global SARS-CoV-2 genomic surveillance, we estimated the impact of sequencing intensity and turnaround times (TAT) on variant detection in 167 countries. Most countries submit genomes >21 days after sample collection, and 77% of low and middle income countries sequenced <0.5% of their cases. We found that sequencing at least 0.5% of the cases, with a TAT <21 days, could be a benchmark for SARS-CoV-2 genomic surveillance efforts. Socioeconomic inequalities substantially impact our ability to quickly detect SARS-CoV-2 variants, and undermine the global pandemic preparedness.

Characterization of the virome associated with Haemagogus mosquitoes in Trinidad, West Indies

Currently, there are increasing concerns about the possibility of a new epidemic due to emerging reports of Mayaro virus (MAYV) fever outbreaks in areas of South and Central America. Haemagogus mosquitoes, the primary sylvan vectors of MAYV are poorly characterized and a better understanding of the mosquito's viral transmission dynamics and interactions with MAYV and other microorganisms would be important in devising effective control strategies. In this study, a metatranscriptomic based approach was utilized to determine the prevalence of RNA viruses in field-caught mosquitoes morphologically identified as Haemagogus janthinomys from twelve (12) forest locations in Trinidad, West Indies. Known insect specific viruses including the Phasi Charoen-like and Humaiata-Tubiacanga virus dominated the virome of the mosquitoes throughout sampling locations while other viruses such as the avian leukosis virus, MAYV and several unclassified viruses had a narrower distribution. Additionally, assembled contigs from the Ecclesville location suggests the presence of a unique uncharacterized picorna-like virus. Mapping of RNA sequencing reads to reference mitochondrial sequences of potential feeding host animals showed hits against avian and rodent sequences, which putatively adds to the growing body of evidence of a potentially wide feeding host-range for the Haemagogus mosquito vector.

Isolation of a novel insect-specific flavivirus with immunomodulatory effects in vertebrate systems

We describe the isolation and characterization of a novel insect-specific flavivirus (ISFV), tentatively named Aripo virus (ARPV), that was isolated from Psorophora albipes mosquitoes collected in Trinidad. The ARPV genome was determined and phylogenetic analyses showed that it is a dual host associated ISFV, and clusters with the main mosquito-borne flaviviruses. ARPV antigen was significantly cross-reactive with Japanese encephalitis virus serogroup antisera, with significant cross-reactivity to Ilheus and West Nile virus (WNV). Results suggest that ARPV replication is limited to mosquitoes, as it did not replicate in the sandfly, culicoides or vertebrate cell lines tested. We also demonstrated that ARPV is endocytosed into vertebrate cells and is highly immunomodulatory, producing a robust innate immune response despite its inability to replicate in vertebrate systems. We show that prior infection or coinfection with ARPV limits WNV-induced disease in mouse models, likely the result of a robust ARPV-induced type I interferon response.

Novel quaranjavirus and other viral sequences identified from ticks parasitizing hunted wildlife in Trinidad and Tobago

Hunters are at a higher risk for exposure to zoonotic pathogens due to their close interactions with wildlife and arthropod vectors. In this study, high throughput sequencing was used to explore the viromes of two tick species, Amblyomma dissimile and Haemaphysalis juxtakochi, removed from hunted wildlife in Trinidad and Tobago. We identified sequences from 3 new viral species, from the viral families Orthomyxoviridae, Chuviridae and Tetraviridae in A. dissimile.

Generation of recombinant hyperimmune globulins from diverse B-cell repertoires

Plasma-derived polyclonal antibody therapeutics, such as intravenous immunoglobulin, have multiple drawbacks, including low potency, impurities, insufficient supply, and batch-to-batch variation. Here we describe a microfluidics and molecular genomics strategy for capturing diverse mammalian antibody repertoires to create recombinant multivalent hyperimmune globulins. Our method generates thousands-diverse mixtures of recombinant antibodies, enriched for specificity and activity against therapeutic targets. Each hyperimmune globulin product comprised thousands to tens of thousands of antibodies derived from convalescent or vaccinated human donors, or immunized mice. Using this approach, we generated hyperimmune globulins with potent neutralizing activity against Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) in under three months, Fc-engineered hyperimmune globulins specific for Zika virus that lacked antibody-dependent enhancement of disease, and hyperimmune globulins specific for lung pathogens present in patients with primary immune deficiency. To address the limitations of rabbit-derived anti-thymocyte globulin (ATG), we generated a recombinant human version and demonstrated its efficacy in mice against graft-versus-host disease.

Serological Evidence for Henipa-like and Filo-like Viruses in Trinidad Bats

Bat-borne zoonotic pathogens belonging to the family Paramxyoviridae, including Nipah and Hendra viruses, and the family Filoviridae, including Ebola and Marburg viruses, can cause severe disease and high mortality rates on spillover into human populations. Surveillance efforts for henipaviruses and filoviruses have been largely restricted to the Old World; however, recent studies suggest a potentially broader distribution for henipaviruses and filoviruses than previously recognized. In the current study, we screened for henipaviruses and filoviruses in New World bats collected across 4 locations in Trinidad near the coast of Venezuela. Bat tissue samples were screened using previously established reverse-transcription polymerase chain reaction assays. Serum were screened using a multiplex immunoassay to detect antibodies reactive with the envelope glycoprotein of viruses in the genus Henipavirus and the family Filoviridae. Serum samples were also screened by means of enzyme-linked immunosorbent assay for antibodies reactive with Nipah G and F glycoproteins. Of 84 serum samples, 28 were reactive with ≥1 henipavirus glycoprotein by ≥1 serological method, and 6 serum samples were reactive against ≥1 filovirus glycoproteins. These data provide evidence of potential circulation of viruses related to the henipaviruses and filoviruses in New World bats.

Determinants of dengue virus dispersal in the Americas

Dengue viruses (DENVs) are classified into four serotypes, each of which contains multiple genotypes. DENV genotypes introduced into the Americas over the past five decades have exhibited different rates and patterns of spatial dispersal. In order to understand factors underlying these patterns, we utilized a statistical framework that allows for the integration of ecological, socioeconomic, and air transport mobility data as predictors of viral diffusion while inferring the phylogeographic history. Predictors describing spatial diffusion based on several covariates were compared using a generalized linear model approach, where the support for each scenario and its contribution is estimated simultaneously from the data set. Although different predictors were identified for different serotypes, our analysis suggests that overall diffusion of DENV-1, -2, and -3 in the Americas was associated with airline traffic. The other significant predictors included human population size, the geographical distance between countries and between urban centers and the density of people living in urban environments.

Characterization of novel, pathogenic field strains of infectious bronchitis virus (IBV) in poultry in Trinidad and Tobago

Avian coronaviruses, including infectious bronchitis virus (IBV) and turkey coronavirus (TCoV), are economically important viruses affecting poultry worldwide. IBV is responsible for causing severe losses to the commercial poultry sector globally. The objectives of this study were to identify the viruses that were causing outbreaks of severe respiratory disease in chickens in Trinidad and Tobago (T&T) and to characterize the strains. Swab samples were collected from birds showing severe respiratory signs in five farms on the island of Trinidad. Samples were tested for the presence of IBV, as well as avian influenza virus (AIV), Newcastle disease virus (NDV) and avian metapneumovirus (aMPV) by real-time reverse transcription polymerase chain reaction (qRT-PCR). All samples from the five farms tested negative for AIV, NDV and aMPV; however, samples from clinically affected birds in all five of the farms tested positive for IBV. Genetic data revealed the presence of TCoV in chickens on two of the farms. Interestingly, these two farms had never reared turkeys. Phylogenetic analysis showed that IBV S1 sequences formed two distinct clusters. Two sequences grouped with vaccine strains within the GI-1 lineage, whereas three sequences grouped together, but separately from other defined lineages, forming a likely new lineage of IBV. Pairwise comparison revealed that the three unique variant strains within the distinct lineage of IBV were significantly different in their S1 nucleotide coding regions from viruses in the closest lineage (16% difference) and locally used vaccine strains (>20% difference). Results also suggested that one of the samples was a recombinant virus, generated from a recombination event between a Trinidad virus of the GI-1 lineage and a Trinidad virus of the newly defined lineage. Many amino acid differences were also observed between the S1 coding regions of the circulating field and vaccine strains, indicating that the IBV vaccines may not be protective. Vaccine-challenge studies are however needed to prove this.

The Serological Prevalence of Rabies Virus-Neutralizing Antibodies in the Bat Population on the Caribbean Island of Trinidad

: Rabies virus (RABV) is the only lyssavirus known to be present within the Caribbean. The island of Trinidad, is richly diverse in chiropteran fauna and endemic for bat-transmitted rabies with low RABV isolation rates observed in this population. We aimed to determine the seroprevalence of rabies virus neutralizing antibodies (RVNA) in light of spatio-temporal and bat demographic factors to infer the extent of natural exposure to RABV in the Trinidadian bat population. RVNA titers were determined by the RABV micro-neutralization test on 383 bat samples representing 21 species, comprising 30.9% of local bat diversity, from 31 locations across the island over 5 years. RVNA was positively detected in 33 samples (8.6%) representing 6 bat species (mainly frugivorous) with titers ranging from 0.1 to 19 IU/mL (mean 1.66 IU/mL). The analyses based on a multivariable binomial generalised linear mixed-effects model showed that bat age and year of capture were significant predictors of seropositivity. Thus, juvenile bats were more likely to be seropositive when compared to adults (estimate 1.13; p = 0.04) which may suggest early exposure to the RABV with possible implications for viral amplification in this population. Temporal variation in rabies seropositivity, 2012-2014 versus 2015-2017 (estimate 1.07; p = 0.03) may have been related to the prevailing rabies epizootic situation. Regarding other factors investigated, RVNA was found in bats from both rural and non-rural areas, as well as in both hematophagous and non-hematophagous bat species. The most common seropositive species, Artibeusjamaicensisplanirostris is ubiquitous throughout the island which may potentially facilitate human exposure. The findings of this study should be factored into public health assessments on the potential for rabies transmission by non-hematophagous bats in Trinidad.

Xylitol and sorbitol effects on the microbiome of saliva and plaque

Chewing gum containing xylitol may help prevent caries by reducing levels of mutans streptococci (MS) and lactobacilli in saliva and plaque. Very little is known about other species which are possibly beneficial to oral health. In this study, we employed high-throughput sequencing of the 16S rRNA gene to profile microbial communities of saliva and plaque following short-term consumption of xylitol and sorbitol containing chewing gum. Participants (n = 30) underwent a washout period and were randomly assigned to one of two groups. Each group chewed either xylitol or sorbitol gum for three weeks, before undergoing a second four-week washout period after which they switched to the alternate gum for three weeks. Analysis of samples collected before and after each intervention identified distinct plaque and saliva microbial communities that altered dependent on the order in which gum treatments were given. Neither the xylitol nor sorbitol treatments significantly affected the bacterial composition of plaque. Lactobacilli were undetected and the number of Streptococcus mutans sequence reads was very low and unaffected by either xylitol or sorbitol. However, sorbitol affected several other streptococcal species in saliva including increasing the abundance of S. cristatus, an oral commensal shown to inhibit bacteria associated with chronic periodontitis.

Rabies in the Caribbean: A Situational Analysis and Historic Review

Rabies virus is the only Lyssavirus species found in the Americas. In discussions about rabies, Latin America and the Caribbean are often grouped together. Our study aimed to independently analyse the rabies situation in the Caribbean and examine changes in rabies spatiotemporal epidemiology. A questionnaire was administered to the 33 member countries and territories of the Caribbean Animal Health Network (CaribVET) to collect current data, which was collated with a literature review. Rabies was endemic in ten Caribbean localities, with the dog, mongoose, and vampire bat identified as enzootic reservoirs. The majority of animal cases occurred in Puerto Rico, the Dominican Republic, and Haiti, while human cases only consistently occurred in the latter two areas. Rabies vaccination was conducted for high-risk animal populations with variable coverage, and rabies diagnostic capacities varied widely throughout the region. Illegal importation and natural migration of animals may facilitate the introduction of rabies virus variants into virus-naïve areas. Passive surveillance, together with enhanced methods and serological screening techniques, can therefore be of value. The insularity of the Caribbean makes it ideal for conducting pilot studies on reservoir host population management. Best practice guidelines developed for these reservoir hosts can be individually modified to the epidemiological status and available resources within each locality.

Serological evidence for eight globally important poultry viruses in Trinidad & Tobago

Viruses affecting poultry cause significant levels of disease leading to severe economic losses among poultry farmers worldwide. The Americas region continues to be vulnerable to the spread of poultry viruses across the continents and Caribbean island chains. In Trinidad and Tobago (T&T) there is limited information on the viruses circulating in poultry. Many flock are vulnerable to infection and there are occasional outbreaks of disease resulting in high levels of morbidity and mortality. This study aims to identify important viruses of poultry circulating in T&T through a broad-based surveillance approach. Serum samples from 29 layer farms in Trinidad and 14 layer farms in Tobago were collected from the eldest laying hens. Samples were tested from unvaccinated birds for antibodies by enzyme-linked immunosorbent assay (ELISA) against Avian influenza virus (AIV), Infectious bronchitis virus (IBV), Newcastle disease virus (NDV), Infectious laryngotracheitis virus (ILTV), Avian pneumovirus (APV), Infectious bursal disease virus (IBDV), Fowl adenovirus Gp1 (FADV) and Egg drop syndrome virus (EDSV). In Trinidad, the estimated true seroprevalence levels of antibodies were 0% (CI 95%: 0-0%) for AIV, 100% (CI 95%: 97-100%) for IBV, 79.8% (CI 95%: 70.6-86.9%) for NDV, 1% (CI 95%: 0-2.6%) for ILTV, 67.55% (CI 95%: 62.3-72.4%) for APV, 94.93% (CI 95%: 88.0-98.6%) for IBDV, 100% (CI 95%: 99.7-100%) for FADV and 67.8% (CI 95%: 62.4-72.8%) for EDSV. In Tobago, seroprevalence levels were 0% (CI 95%: 0-0%) for AIV, 100% (CI 95%: 95.6-100%) for IBV, 80.5% (CI 95%: 70.1-88.5%) for NDV, 29.9% (CI 95%: 20.8-40.6%) for ILTV, 100% (CI 95%: 97.7-100%) for APV, 97.1% (CI95%: 89.9-100%) for IBDV, 100% (CI 95%: 97.5-100%) for FADV and 100% (CI 95%: 99-100%) for EDSV. The results reveal strong evidence for the circulation of IBV, NDV, APV, IBDV, FADV and EDSV in layer poultry on both islands, as well as ILTV in Tobago.

Seroprevalence of economically important viral pathogens in swine populations of Trinidad and Tobago, West Indies

The objective of this study was to evaluate the seroprevalence and identify the strains of swine influenza virus (SwIV), as well as the seroprevalence of porcine parvovirus (PPV), transmissible gastroenteritis virus (TGEV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine respiratory coronavirus (PRCV), porcine circovirus type 2 (PCV-2), and classical swine fever virus (CSFV) in pigs in Trinidad and Tobago (T&T). Blood samples (309) were randomly collected from pigs at farms throughout T&T. Serum samples were tested for the presence of antibodies to the aforementioned viruses using commercial ELISA kits, and the circulating strains of SwIV were identified by the hemagglutination inhibition test (HIT). Antibodies against SwIV were detected in 114 out of the 309 samples (37%). Out of a total of 26 farms, 14 tested positive for SwIV antibodies. HI testing revealed high titers against the A/sw/Minnesota/593/99 H3N2 strain and the pH1N1 2009 pandemic strain. Antibodies against PPV were detected in 87 out of the 309 samples (28%), with 11 out of 26 farms testing positive for PPV antibodies. Antibodies against PCV-2 were detected in 205 out of the 309 samples tested (66%), with 25 out of the 26 farms testing positive for PCV-2 antibodies. No antibodies were detected in any of the tested pigs to PRRSV, TGEV, PRCV, or CSFV.

The History of Rabies in Trinidad: Epidemiology and Control Measures

Vampire bat-transmitted rabies was first recognized in Trinidad during a major outbreak reported in 1925. Trinidad is the only Caribbean island with vampire bat-transmitted rabies. We conducted a literature review to describe the changing epidemiology of rabies in Trinidad and give a historical perspective to rabies prevention and control measures on the island. The last human case of rabies occurred in 1937 and although no case of canine-transmitted rabies was reported since 1914, sporadic outbreaks of bat-transmitted rabies still occur in livestock to date. Over the last century, seven notable epidemics were recorded in Trinidad with the loss of over 3000 animals. During the 1950s, several measures were effectively adopted for the prevention and control of the disease which led to a significant reduction in the number of cases. These measures include: vampire bat population control, livestock vaccination, and animal surveillance. However, due to lapses in these measures over the years (e.g., periods of limited vampire control and incomplete herd vaccination), epidemics have occurred. In light of the significant negative impact of rabies on animal production and human health, rabies surveillance in Trinidad should be enhanced and cases evaluated towards the design and implementation of more evidence-based prevention and control programs.

Understanding the evolution and spread of chikungunya virus in the Americas using complete genome sequences

Local transmission of chikungunya virus (CHIKV) was first detected in the Americas in December 2013, after which it spread rapidly throughout the Caribbean islands and American mainland, causing a major chikungunya fever epidemic. Previous phylogenetic analysis of CHIKV from a limited number of countries in the Americas suggests that an Asian genotype strain was responsible, except in Brazil where both Asian and East/Central/South African (ECSA) lineage strains were detected. In this study, we sequenced thirty-three complete CHIKV genomes from viruses isolated in 2014 from fourteen Caribbean islands, the Bahamas and two mainland countries in the Americas. Phylogenetic analyses confirmed that they all belonged to the Asian genotype and clustered together with other Caribbean and mainland sequences isolated during the American outbreak, forming an 'Asian/American' lineage defined by two amino acid substitutions, E2 V368A and 6K L20M, and divided into two well-supported clades. This lineage is estimated to be evolving at a mean rate of 5 × 10^-4 substitutions per site per year (95% higher probability density, 2.9-7.9 × 10^-4) and to have arisen from an ancestor introduced to the Caribbean (most likely from Oceania) in about March 2013, 9 months prior to the first report of CHIKV in the Americas. Estimation of evolutionary rates for individual gene regions and selection analyses indicate that (in contrast to the Indian Ocean Lineage that emerged from the ECSA genotype followed by adaptive evolution and with a significantly higher substitution rate) the evolutionary dynamics of the Asian/American lineage are very similar to the rest of the Asian genotype and natural selection does not appear to have played a major role in its emergence. However, several codon sites with evidence of positive selection were identified within the non-structural regions of Asian genotype sequences outside of the Asian/American lineage.

Molecular Characterisation of Chikungunya Virus Infections in Trinidad and Comparison of Clinical and Laboratory Features with Dengue and Other Acute Febrile Cases

Local transmission of Chikungunya virus (CHIKV) was first documented in Trinidad and Tobago (T&T) in July 2014 preceding a large epidemic. At initial presentation, it is difficult to distinguish chikungunya fever (CHIKF) from other acute undifferentiated febrile illnesses (AUFIs), including life-threatening dengue disease. We characterised and compared dengue virus (DENV) and CHIKV infections in 158 patients presenting with suspected dengue fever (DF) and CHIKF at a major hospital in T&T, and performed phylogenetic analyses on CHIKV genomic sequences recovered from 8 individuals. The characteristics of patients with and without PCR-confirmed CHIKV were compared using Pearson’s χ² and student’s t-tests, and adjusted odds ratios (aORs) and 95% confidence intervals (CIs) were determined using logistic regression. We then compared signs and symptoms of people with RT-qPCR-confirmed CHIKV and DENV infections using the Mann-Whitney U, Pearson’s χ² and Fisher’s exact tests. Among the 158 persons there were 8 (6%) RT-qPCR-confirmed DENV and 30 (22%) RT-qPCR-confirmed CHIKV infections. Phylogenetic analyses showed that the CHIKV strains belonged to the Asian genotype and were most closely related to a British Virgin Islands strain isolated at the beginning of the 2013/14 outbreak in the Americas. Compared to persons who were RT-qPCR-negative for CHIKV, RT-qPCR-positive individuals were significantly more likely to have joint pain (aOR: 4.52 [95% CI: 1.28–16.00]), less likely to be interviewed at a later stage of illness (days post onset of fever—aOR: 0.56 [0.40–0.78]) and had a lower white blood cell count (aOR: 0.83 [0.71–0.96]). Among the 38 patients with RT-qPCR-confirmed CHIKV or DENV, there were no significant differences in symptomatic presentation. However when individuals with serological evidence of recent DENV or CHIKV infection were included in the analyses, there were key differences in clinical presentation between CHIKF and other AUFIs including DF, which can be used to triage patients for appropriate care in the clinical setting.

Chikungunya virus (CHIKV) recently emerged in the Americas and caused a major epidemic of chikungunya fever (CHIKF). While not usually life threatening, CHIKF is a debilitating and often chronic illness resulting in major morbidity and economic losses. It is difficult to distinguish CHIKF from other viral illnesses that cause acute fevers including dengue fever (DF), an important consideration since DF can be life-threatening and early identification and treatment of cases is key to reducing mortality. In this study we investigated individuals presenting to a major hospital in Trinidad and Tobago (T&T) with suspected DF or CHIKF, and identified signs and symptoms that distinguish these two illnesses. We also recovered complete genome sequences for CHIKV and show that the etiologic strain in T&T is a closely related descendent of the strain first isolated from the British Virgin Islands at the beginning of the outbreak in the Americas.

Improved serotype-specific dengue virus detection in Trinidad and Tobago using a multiplex, real-time RT-PCR

Dengue virus (DENV) transmission occurs throughout the Caribbean, though laboratory confirmation and epidemiologic surveillance are limited by the availability of serotype-specific molecular diagnostics. In this study, we show that a serotype-specific DENV multiplex, real-time reverse transcriptase-PCR (RT-PCR) detected DENV RNA in significantly more samples (82/182) than a reference hemi-nested RT-PCR (57/182; P=0.01).

Seroepidemiology of selected alphaviruses and flaviviruses in bats in Trinidad

A serosurvey of antibodies against selected flaviviruses and alphaviruses in 384 bats (representing 10 genera and 14 species) was conducted in the Caribbean island of Trinidad. Sera were analysed using epitope‐blocking enzyme‐linked immunosorbent assays (ELISAs) specific for antibodies against West Nile virus (WNV), Venezuelan equine encephalitis virus (VEEV) and eastern equine encephalitis virus (EEEV), all of which are zoonotic viruses of public health significance in the region. Overall, the ELISAs resulted in the detection of VEEV‐specific antibodies in 11 (2.9%) of 384 bats. Antibodies to WNV and EEEV were not detected in any sera. Of the 384 sera, 308 were also screened using hemagglutination inhibition assay (HIA) for antibodies to the aforementioned viruses as well as St. Louis encephalitis virus (SLEV; which also causes epidemic disease in humans), Rio Bravo virus (RBV), Tamana bat virus (TABV) and western equine encephalitis virus (WEEV). Using this approach, antibodies to TABV and RBV were detected in 47 (15.3%) and 3 (1.0%) bats, respectively. HIA results also suggest the presence of antibodies to an undetermined flavivirus(es) in 8 (2.6%) bats. Seropositivity for TABV was significantly (P < 0.05; χ²) associated with bat species, location and feeding preference, and for VEEV with roost type and location. Differences in prevalence rates between urban and rural locations were statistically significant (P < 0.05; χ²) for TABV only. None of the aforementioned factors was significantly associated with RBV seropositivity rates.

Characterization of a novel Negevirus and a novel Bunyavirus isolated from Culex (Culex) declarator mosquitoes in Trinidad

Pools of mosquitoes were tested for insect-specific viruses using cytopathic effect (CPE) assays on Aedes albopictus (C6/36) cells. Illumina sequencing of RNA from pool TR7094, which produced extensive CPE 2 days post-infection, yielded the complete genome sequences of a previously unknown Bunyavirus, designated Cumuto virus (CUMV), and a second virus designated Wallerfield virus (WALV). WALV shared highest amino acid identity (60.1 %) with Dezidougou virus from Côte d'Ivoire, a positive-sense, single-strand RNA, insect-specific virus belonging to the newly proposed genus Negevirus associated with mosquitoes and phlebotomine sandflies. The S, M and L segments of CUMV were most closely related to those of Gouleako virus, also from Côte d'Ivoire (amino acid identities of 36 %, 38% and 54 % respectively). Neither virus produced CPE on vertebrate cells, or illness in newborn mice. Isolation and characterization of these viruses increase our knowledge of the geographical distribution, diversity and host range of mosquito-specific bunyaviruses and negeviruses.

Evolutionary history and phylogeography of rabies viruses associated with outbreaks in Trinidad

Bat rabies is an emerging disease of public health significance in the Americas. The Caribbean island of Trinidad experiences periodic outbreaks within the livestock population. We performed molecular characterisation of Trinidad rabies virus (RABV) and used a Bayesian phylogeographic approach to investigate the extent to which outbreaks are a result of in situ evolution versus importation of virus from the nearby South American mainland. Trinidadian RABV sequences were confirmed as bat variant and clustered with Desmodus rotundus (vampire bat) related sequences. They fell into two largely temporally defined lineages designated Trinidad I and II. The Trinidad I lineage which included sequences from 1997-2000 (all but two of which were from the northeast of the island) was most closely related to RABV from Ecuador (2005, 2007), French Guiana (1990) and Venezuela (1993, 1994). Trinidad II comprised sequences from the southwest of the island, which clustered into two groups: Trinidad IIa, which included one sequence each from 2000 and 2007, and Trinidad IIb including all 2010 sequences. The Trinidad II sequences were most closely related to sequences from Brazil (1999, 2004) and Uruguay (2007, 2008). Phylogeographic analyses support three separate RABV introductions from the mainland from which each of the three Trinidadian lineages arose. The estimated dates for the introductions and subsequent lineage expansions suggest periods of in situ evolution within Trinidad following each introduction. These data also indicate co-circulation of Trinidad lineage I and IIa during 2000. In light of these findings and the likely vampire bat origin of Trinidadian RABV, further studies should be conducted to investigate the relationship between RABV spatiotemporal dynamics and vampire bat population ecology, in particular any movement between the mainland and Trinidad.

Serological evidence of flaviviruses and alphaviruses in livestock and wildlife in Trinidad

Seroprevalence rates of selected arboviruses in animal populations in Trinidad were determined using serum samples collected between 2006 and 2009 from horses (n=506), cattle (n=163), sheep (n=198), goats (n=82), pigs (n=184), birds (n=140), rodents (n=116), and other vertebrates (n=23). The sera were screened for antibodies to West Nile virus (WNV), St. Louis encephalitis virus (SLEV), Ilheus virus (ILHV), Bussuquara virus (BSQV), Venezuelan equine encephalitis virus (VEEV), eastern equine encephalitis virus (EEEV), and western equine encephalitis virus (WEEV), using hemagglutination inhibition assay (HIA) and epitope-blocking enzyme-linked immunosorbent assays (ELISA). Antibodies to SLEV were detected in a total of 49 (9.7%) horses, 8 (4.9%) cattle, 1 (1.2%) goat, 2 (1.4%) wild birds, and 3 (2.2%) wild rodents by both methods. In contrast, antibodies to EEEV, VEEV, and WNV were detected only in horses, at rates of 4.3%, 0.8%, and 17.2%, respectively, by ELISA, and IgM capture ELISA was WNV-positive in 3 (0.6%) of these sera. Among locally bred unvaccinated horses that had never left Trinidad, seroprevalence rates against WNV were 12.1% and 17.2% by ELISA and HIA, respectively. The presence of WNV- and SLEV-specific antibodies in a representative sample of horse sera that were both ELISA- and HIA-seropositive was confirmed by plaque reduction neutralization testing (PRNT). Antibodies to ILHV and BSQV were not detected in any of the serum samples tested (i.e., sera from horses, other livestock, and wild birds in the case of ILHV, and wild mammals in the case of BSQV). The data indicate the presence of WNV in Trinidad, and continuing low-level circulation of SLEV, EEEV, and VEEV.

Phylogeography and population dynamics of dengue viruses in the Americas

Changes in Dengue virus (DENV) disease patterns in the Americas over recent decades have been attributed, at least in part, to repeated introduction of DENV strains from other regions, resulting in a shift from hypoendemicity to hyperendemicity. Using newly sequenced DENV-1 and DENV-3 envelope (E) gene isolates from 11 Caribbean countries, along with sequences available on GenBank, we sought to document the population genetic and spatiotemporal transmission histories of the four main invading DENV genotypes within the Americas and investigate factors that influence the rate and intensity of DENV transmission. For all genotypes, there was an initial invasion phase characterized by rapid increases in genetic diversity, which coincided with the first confirmed cases of each genotype in the region. Rapid geographic dispersal occurred upon each genotype's introduction, after which individual lineages were locally maintained, and gene flow was primarily observed among neighboring and nearby countries. There were, however, centers of viral diversity (Barbados, Puerto Rico, Colombia, Suriname, Venezuela, and Brazil) that were repeatedly involved in gene flow with more distant locations. For DENV-1 and DENV-2, we found that a "distance-informed" model, which posits that the intensity of virus movement between locations is inversely proportional to the distance between them, provided a better fit than a model assuming equal rates of movement between all pairs of countries. However, for DENV-3 and DENV-4, the more stochastic "equal rates" model was preferred.

Human papillomavirus genotypes and their prevalence in a cohort of women in Trinidad

OBJECTIVE

Human papillomavirus (HPV) genotypes and their relative prevalences were determined in a cohort of 310 sexually active women in Trinidad, West Indies.

METHODS

Cervical samples were collected with Ayre's spatulas and endocervical brushes. Samples were used for the conventional Papanicolaou test and for determining HPV genotypes by amplification of a section of the viral L1 gene, followed by DNA sequencing and probe hybridization.

RESULTS

HPV infections were identified in 126 of 310 (40.6%) women. Of them, 83 (65.8%) were infected with high-risk HPV, 16 (12.7%) with low-risk HPV, and 27 (21.4%) with HPV types of unknown risk. HPV 52 (12.7%) was the most frequently occurring high-risk type, followed by HPV 66 (10.3%), HPV 16 (9.5%), and HPV 18 (8.6%). High-risk types HPV 16 and HPV 66 were each found in 3 (20.0%) and HPV 18 was found in 1 (6.6%) of the 15 women with abnormal cytology.

CONCLUSIONS

Cervical HPV prevalence and heterogeneity of HPV genotypes are high in this Trinidad cohort. The relative importance of HPV genotypes in the development of cervical lesions needs further investigation in Trinidad in order to better understand the epidemiology of HPV infections as well as to determine the role of HPV testing in the screening, prevention, and control of cervical cancer. This pilot study provided important information on the prevalence of HPV genotypes, which will be used in future nationwide studies.

Isolation and characterization of sylvatic mosquito-borne viruses in Trinidad: enzootic transmission and a new potential vector of Mucambo virus

Mosquito surveillance was carried out in three forested regions of Trinidad during July 2007-March 2009. A total of 185,397 mosquitoes representing at least 46 species was collected, divided into pools of 1-50 mosquitoes according to species and sex, and screened for arboviruses using cytopathic effect assays on Vero cell monolayers. Eighty-five viruses were isolated, including members of the genera Alphavirus (Mucambo virus; MUCV) and Orthobunyavirus (Caraparu, Oriboca, Bimiti, and Wyeomyia viruses). Species of the Culex subgenus Melanoconion accounted for 56% of the total number of mosquitoes collected and 97% of the viruses isolated; Cx. (Mel.) portesi accounted for 92% of virus isolations. Our results also implicate for the first time Aedes (Ochlerotatus) hortator as a potential vector of MUCV. Phylogenetic analyses of 43 MUCV strains suggest population subdivision within Trinidad, consistent with the hypothesis of enzootic maintenance in localized rodent populations.

Yellow fever virus maintenance in Trinidad and its dispersal throughout the Americas

Trinidad, like many other American regions, experiences repeated epizootics of yellow fever virus (YFV). However, it is unclear whether these result from in situ evolution (enzootic maintenance) or regular reintroduction of YFV from the South American mainland. To discriminate between these hypotheses, we carried out a Bayesian phylogeographic analysis of over 100 prM/E gene sequences sampled from 8 South American countries. These included newly sequenced isolates from the recent 2008-2009 Trinidad epizootic and isolates derived from mainland countries within the last decade. The results indicate that the most recent common ancestor of the 2008-2009 epizootic existed in Trinidad 4.2 years prior to 2009 (95% highest probability density [HPD], 0.5 to 9.0 years). Our data also suggest a Trinidad origin for the progenitor of the 1995 Trinidad epizootic and support in situ evolution of YFV between the 1979 and 1988-1989 Trinidad epizootics. Using the same phylogeographic approach, we also inferred the historical spread of YFV in the Americas. The results suggest a Brazilian origin for YFV in the Americas and an overall dispersal rate of 182 km/year (95% HPD, 52 to 462 km/year), with Brazil as the major source population for surrounding countries. There is also strong statistical support for epidemiological links between four Brazilian regions and other countries. In contrast, while there were well-supported epidemiological links within Peru, the only statistically supported external link was a relatively weak link with neighboring Bolivia. Lastly, we performed a complete analysis of the genome of a newly sequenced Trinidad 2009 isolate, the first complete genome for a genotype I YFV isolate.

Meiotic recombination generates rich diversity in NK cell receptor genes, alleles, and haplotypes

Natural killer (NK) cells contribute to the essential functions of innate immunity and reproduction. Various genes encode NK cell receptors that recognize the major histocompatibility complex (MHC) Class I molecules expressed by other cells. For primate NK cells, the killer-cell immunoglobulin-like receptors (KIR) are a variable and rapidly evolving family of MHC Class I receptors. Studied here is KIR3DL1/S1, which encodes receptors for highly polymorphic human HLA-A and -B and comprises three ancient allelic lineages that have been preserved by balancing selection throughout human evolution. While the 3DS1 lineage of activating receptors has been conserved, the two 3DL1 lineages of inhibitory receptors were diversified through inter-lineage recombination with each other and with 3DS1. Prominent targets for recombination were D0-domain polymorphisms, which modulate enhancer function, and dimorphism at position 283 in the D2 domain, which influences inhibitory function. In African populations, unequal crossing over between the 3DL1 and 3DL2 genes produced a deleted KIR haplotype in which the telomeric "half" was reduced to a single fusion gene with functional properties distinct from its 3DL1 and 3DL2 parents. Conversely, in Eurasian populations, duplication of the KIR3DL1/S1 locus by unequal crossing over has enabled individuals to carry and express alleles of all three KIR3DL1/S1 lineages. These results demonstrate how meiotic recombination combines with an ancient, preserved diversity to create new KIR phenotypes upon which natural selection acts. A consequence of such recombination is to blur the distinction between alleles and loci in the rapidly evolving human KIR gene family.

Isolation and phylogenetic analysis of Mucambo virus (Venezuelan equine encephalitis complex subtype IIIA) in Trinidad

In the 1950s and 1960s, alphaviruses in the Venezuelan equine encephalitis (VEE) antigenic complex were the most frequently isolated arboviruses in Trinidad. Since then, there has been very little research performed with these viruses. Herein, we report on the isolation, sequencing, and phylogenetic analyses of Mucambo virus (MUCV; VEE complex subtype IIIA), including 6 recently isolated from Culex (Melanoconion) portesi mosquitoes and 11 previously isolated in Trinidad and Brazil. Results show that nucleotide and amino acid identities across the complete structural polyprotein for the MUCV isolates were 96.6-100% and 98.7-100%, respectively, and the phylogenetic tree inferred for MUCV was highly geographically- and temporally-structured. Bayesian analyses suggest that the sampled MUCV lineages have a recent common ancestry of approximately 198 years (with a 95% highest posterior density (HPD) interval of 63-448 years) prior to 2007, and an overall rate of evolution of 1.28 x 10(-4) substitutions/site/yr.

Detection and phylogenetic analysis of group 1 coronaviruses in South American bats

Bat coronaviruses (Bt-CoVs) are thought to be the precursors of severe acute respiratory syndrome coronavirus. We detected Bt-CoVs in 2 bat species from Trinidad. Phylogenetic analysis of the RNA-dependent RNA polymerase gene and helicase confirmed them as group 1 coronaviruses.

Characterization of Culex Flavivirus (Flaviviridae) strains isolated from mosquitoes in the United States and Trinidad

Recent reports indicate that flaviviruses similar to the cell fusing agent virus (CFAV) naturally infect a wide variety of mosquito species. These newly recognized insect-specific viruses comprise a distinct CFAV complex within the genus Flavivirus. Here, we describe the isolation and characterization of nine strains of Culex flavivirus (Cx FV), a member of the CFAV complex, from mosquitoes collected in the United States (East Texas) and Trinidad. Phylogenetic analyses of the envelope protein gene sequences of these nine mosquito isolates with those of other CFAV complex flaviviruses in GenBank indicate that the U.S. isolates group with CxFV isolates from Asia (Japan and Indonesia), while the Trinidad isolates are more similar to CxFV isolates from Central America. A discussion follows on the possible biological significance of the CFAV complex flaviviruses.

Evolution and dispersal of St. Louis encephalitis virus in the Americas

Using a Bayesian coalescent approach on a dataset of 73 envelope gene sequences we estimated substitution rates and dates of divergence for St. Louis encephalitis virus (SLEV) in the Americas. We found significant rate heterogeneity among lineages, such that "relaxed" molecular clock models were much better supported than a strict molecular clock. The mean substitution rate estimated for all SLEV was 4.1x10(-4)substitutions/site/year (95% HPD 2.5-5.7)-higher than previous estimates that relied on the less well-suited strict clock. Mean substitution rates for individual lineages varied from 3.7x10(-4) to 7.2x10(-4)substitutions/site/year. For the first time we also assessed the magnitude and direction of viral gene flow within the Americas. The overall direction of gene flow during the period represented by the phylogeny is from South to North, and the region between 15 degrees N and 30 degrees N latitude appears to be the major source of virus for the rest of North America, which is consistent with migratory birds returning to their northern breeding grounds having acquired infection while wintering in the region of the Gulf of Mexico.

Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans

Interactions of killer cell immunoglobulin-like receptors (KIRs) with major histocompatibility complex (MHC) class I ligands diversify natural killer cell responses to infection. By analyzing sequence variation in diverse human populations, we show that the KIR3DL1/S1 locus encodes two lineages of polymorphic inhibitory KIR3DL1 allotypes that recognize Bw4 epitopes of protein">HLA-A and HLA-B and one lineage of conserved activating KIR3DS1 allotypes, also implicated in Bw4 recognition. Balancing selection has maintained these three lineages for over 3 million years. Variation was selected at D1 and D2 domain residues that contact HLA class I and at two sites on D0, the domain that enhances the binding of KIR3D to HLA class I. HLA-B variants that gained Bw4 through interallelic microconversion are also products of selection. A worldwide comparison uncovers unusual KIR3DL1/S1 evolution in modern sub-Saharan Africans. Balancing selection is weak and confined to D0, KIR3DS1 is rare and KIR3DL1 allotypes with similar binding sites predominate. Natural killer cells express the dominant KIR3DL1 at a high frequency and with high surface density, providing strong responses to cells perturbed in Bw4 expression.

Invasion and maintenance of dengue virus type 2 and type 4 in the Americas

Dengue virus type 4 (DENV-4) was first reported in the Americas in 1981, where it caused epidemics of dengue fever throughout the region. In the same year, the region's first epidemic of dengue hemorrhagic fever was reported, caused by an Asian strain of dengue virus type 2 (DENV-2) that was distinct from the American subtype circulating previously. Despite the importance of these epidemics, little is known about the rates or determinants of viral spread among island and mainland populations or their directions of movement. We employed a Bayesian coalescent approach to investigate the transmission histories of DENV-2 and DENV-4 since their introduction in 1981 and a parsimony method to assess patterns of strain migration. For both viruses there was an initial invasion phase characterized by an exponential increase in the number of DENV lineages, after which levels of genetic diversity remained constant despite reported fluctuations in DENV-2 and DENV-4 activity. Strikingly, viral lineage numbers increased far more rapidly for DENV-4 than DENV-2, indicative of a more rapid rate of exponential population growth in DENV-4 or a higher rate of geographic dispersal, allowing this virus to move more effectively among localities. We propose that these contrasting dynamics may reflect underlying differences in patterns of host immunity. Despite continued gene flow along particular transmission routes, the overall extent of viral traffic was less than expected under panmixis. Hence, DENV in the Americas has a clear geographic structure that maintains viral diversity between outbreaks.

Phylogeography and molecular evolution of dengue 2 in the Caribbean basin, 1981-2000

We sequenced the envelope (E) genes of 59 DEN-2 isolates collected from ten Caribbean islands, six South American countries, and two Central American countries between 1981 and 2000, a period characterized by hyperendemicity and increased incidence of severe dengue. Fifty-two isolates belonged to "American/Asian" subtype IIIb, possessing a characteristic polar residue at envelope aa position 390 (N [n = 48] or S [n = 4]) common to that group. Six isolates from Trinidad (1981), Honduras (1991 [4]), and El Salvador (1987) fell into the "Native American" subtype V (D at aa 390), and one from Honduras (1986) belonged to "Asian" subtype I. The data suggest that after its first isolation in the Caribbean in 1981, genotype IIIb spread throughout the Americas and effectively replaced subtype V throughout the Caribbean basin. The strain also evolved into several distinct lineages, based on substitutions in the E glycoprotein (amino acids 91 and 131), two of which were still in circulation in 2000. Interestingly, a molecular clock did not fit the data well, suggesting that other sources of rate variation, such as differential selection or differences in effective population sizes, may exist among lineages. Our results indicate the importance of large temporal- and geographical-scale phylogenetic studies in understanding disease dynamics, particularly where replacements between regions can occur.

SNP haplotypes and allele frequencies show evidence for disruptive and balancing selection in the human leukocyte receptor complex

The human leukocyte receptor complex (LRC) of Chromosome 19q13.4 encodes polymorphic and highly homologous genes that are expressed by cells of the immune system and regulate their function. There is an enormous diversity at the LRC, most particularly the variable number of killer cell immunoglobulin-like receptor (KIR) genes. KIR have been associated with several disease processes due to their interaction with polymorphic human leukocyte antigen class I molecules. We have assessed haplotype compositions, linkage disequilibrium patterns and allele frequencies in two Caucasoid population samples (n=54, n=100), using a composite of single-nucleotide polymorphism (SNP) markers and high-resolution, allele-specific molecular genotyping. Particular KIR loci segregated with SNP and other markers, forming two blocks that were separated by a region with a greater history of recombination. The KIR haplotype composition and allele frequency distributions were consistent with KIR having been subject to balancing selection (Watterson's F: P=0.001). In contrast, there was a high inter-population heterogeneity measure for the LRC-encoded leukocyte immunoglobulin-like receptor A3 (LILRA3), indicating pathogen-driven disruptive selection (Wright's FST=0.32). An assessment of seven populations representative of African, Asian and Caucasoid ethnic groups (total n=593) provided little evidence for long-range LRC haplotypes. The different natural selection pressures acting on each locus may have contributed to a lack of linkage disequilibrium between them.

Pfcrt and pfmdr1 alleles associated with chloroquine resistance in Plasmodium falciparum from Guyana, South America

Using DNA extracted from 112 parasitised blood blots, we screened for the population marker of chloroquine resistance (CQR) pfcrt K76T in Plasmodium falciparum infections from Guyana. Pfmdr1 mutations S1034C, N1042D, and D1246Y also associated with CQR were surveyed as well in 15 isolates for which the in vitro responses to CQ were known. Results indicate that the pfcrt K76T is ubiquitous in this environment, and confirmatory sequencing of codons 72 and 76 revealed two novel allelic sequences SVMIT and RVMNT in addition to the previously identified CVMNT and SVMNT haplotypes. The frequency of the pfcrt K76T despite its presence in both CQR and CQS (chloroquine sensitive) infections measured in vivo and in vitro, suggests that it is a useful population marker in this low-transmission setting of sweeping CQR.

Analysis of Fc gamma receptor II (CD32) polymorphism in populations of African and South Asian ancestry reveals east-west geographic gradients of allele frequencies

Analysis of FcgammaRIIA alleles in Pakistanis and in Trinidadians of South Asian, African and mixed ancestry revealed no significant differences between Trinidadian South Asians and Pakistanis. H131 homozygotes were more common among Trinidadian South Asians than among Africans and those of mixed ancestry. Comparison with other populations revealed east-west geographic gradients of allele frequencies.

Molecular evolution and phylogeny of dengue type 4 virus in the Caribbean

We sequenced the E gene and adjacent prM/M and NS1 junctions (1940 bp) of 48 Dengue-4 (DEN-4) isolates collected between 1981 and 1999 from 8 Caribbean islands and from 7 South and Central American countries. Phylogenetic analysis confirms a single introduction in the early 1980s and a high degree of gene flow resulting in a pattern of evolution defined more by time period than geographic origin, especially within the Caribbean basin. A modern Caribbean clade consisting of four distinct lineages has arisen, comprised of isolates from Caribbean islands and nearby regions of South America. This clade is defined by three amino acid substitutions in the E (aa 163 and 351) and NS1 (aa 52) proteins. These findings highlight the importance of migration and gene flow in dengue viral change and suggest that efforts to understand disease dynamics in the Caribbean basin need to focus at regional, rather than local scales.

A comparison of HLA-DR and -DQ allele and haplotype frequencies in Trinidadian populations of African, South Asian, and mixed ancestry

Using polymerase chain reaction-sequence-specific primer (PCR-SSP) typing, this study determined the frequency of human leukocyte antigen (HLA) DR- and -DQ alleles and haplotypes in individuals of African (n = 75), South Asian (n = 98), and mixed (n = 102) ancestry from the Caribbean island of Trinidad. We detected 19 different haplotypes containing DRB3, 8 containing DRB4, 6 containing DRB5, and 6 different haplotypes without DRB3/4/5 genes. Twenty-nine haplotypes were identified in Africans, 24 in the South Asians, and 32 in the mixed group. We detected significant differences between the populations, principally at the DQA1 and DQB1 loci, although the allele frequency for DRB1*0901 was highest in the Africans (p(c) < 0.05). Trinidad African and mixed groups were generally more diverse than the South Asians and displayed a wider range of DRB1-DQB1 associations; DQB1*02 and DQB1*0301 each associated with five to six different DRB1 alleles in the Africans and mixed group but only two in South Asians. In the Africans and the mixed group, DQB1*04 was found with DRB1*0302 and DRB1*04, but only with DRB1*08 in the South Asians. Trinidad Africans revealed consistencies with populations in Western, Central, and Northern Africa, but differed considerably from individual populations on the African continent. Trinidad South Asians displayed similar allele frequencies and associations to other populations from Northern India.

Natural killer cell immunoglobulin-like receptor (KIR) locus profiles in African and South Asian populations

Natural killer (NK) and some T cells express killer cell immunoglobulin-like receptors (KIRs), which interact with HLA class I expressed by target cells and consequently regulate cytolytic activity. The number of KIR loci can vary and so a range of genetic profiles is observed. We have determined the KIR genetic profiles from one African (n = 62) and two South Asian (n = 108, n = 78) populations. Several of the KIRs are present at significantly different frequencies between the two major ethnic groups (eg KIR2DS4 gene frequency 0.82 African, 0.47 S Asian. Pc < 1 x 10(-6)) and this is due to uneven distribution of two KIR haplotype families 'A' and 'B'. All three populations described here displayed a greater degree of diversity of KIR genetic profiles than other populations investigated, which indicates further complexity of underlying haplotypes; in this respect we describe two individuals who appear homozygous for a large deletion including the previously ubiquitous 2DL4. We have also reanalysed three populations that we studied previously, for the presence of a KIR which is now known to be an indicator of the 'B' haplotype. South Asians had the highest overall frequencies of all KIR loci characteristic of 'B' haplotypes (Pc < 0.0001 to < 0.004). Furthermore, gene frequency independent deviances in the linkage disequilibrium were apparent between populations.

A CD45 polymorphism associated with abnormal splicing is absent in African populations

The CD45 antigen is essential for normal antigen receptor-mediated signalling in lymphocytes, and different patterns of splicing of CD45 are associated with distinct functions in lymphocytes. Abnormal CD45 splicing has been recognized in humans, caused by a C77G transversion in the gene encoding CD45 (PTPRC). Recently the C77G polymorphism has been associated with multiple sclerosis and increased susceptibility to HIV-1 infection. These studies suggest that the regulation of CD45 splicing may be critical for the proper function of the immune system. Because of these data we examined the frequency of the C77G allele in African and Asian populations from countries with high or low prevalence of HIV infection. Here we report that the variant CD45 C77G allele is absent in African populations. We further show that populations living in the Pamir mountains of Central Asia have a very high prevalence of the C77G variant.