Several mosquito-borne diseases affecting humans are emerging or reemerging in the United States. The early detection of pathogens in mosquito populations is essential to prevent and control the spread of these diseases. In this study, we tested the potential applicability of the Lawrence Livermore Microbial Detection Array (LLMDA) to enhance biosurveillance by detecting microbes present in Aedes aegypti, Aedes albopictus, and Culex mosquitoes, which are major vector species globally, including in Texas. The sensitivity and reproducibility of the LLMDA were tested in mosquito samples spiked with different concentrations of dengue virus (DENV), revealing a detection limit of >100 but <1,000 PFU/ml. Additionally, field-collected mosquitoes from Chicago, IL, and College Station, TX, of known infection status (West Nile virus [WNV] and Culex flavivirus [CxFLAV] positive) were tested on the LLMDA to confirm its efficiency. Mosquito field samples of unknown infection status, collected in San Antonio, TX, and the Lower Rio Grande Valley (LRGV), TX, were run on the LLMDA and further confirmed by PCR or quantitative PCR (qPCR). The analysis of the field samples with the LLMDA revealed the presence of cell-fusing agent virus (CFAV) in A. aegypti populations. Wolbachia was also detected in several of the field samples (A. albopictus and Culex spp.) by the LLMDA. Our findings demonstrated that the LLMDA can be used to detect multiple arboviruses of public health importance, including viruses that belong to the Flavivirus, Alphavirus, and Orthobunyavirus genera. Additionally, insect-specific viruses and bacteria were also detected in field-collected mosquitoes. Another strength of this array is its ability to detect multiple viruses in the same mosquito pool, allowing for the detection of cocirculating pathogens in an area and the identification of potential ecological associations between different viruses. This array can aid in the biosurveillance of mosquito-borne viruses circulating in specific geographical areas.IMPORTANCE Viruses associated with mosquitoes have made a large impact on public and veterinary health. In the United States, several viruses, including WNV, DENV, and chikungunya virus (CHIKV), are responsible for human disease. From 2015 to 2018, imported Zika cases were reported in the United States, and in 2016 to 2017, local Zika transmission occurred in the states of Texas and Florida. With globalization and a changing climate, the frequency of outbreaks linked to arboviruses will increase, revealing a need to better detect viruses in vector populations. With the capacity of the LLMDA to detect viruses, bacteria, and fungi, this study highlights its ability to broadly screen field-collected mosquitoes and contribute to the surveillance and management of arboviral diseases.

Aedes aegypti control requires dedicated resources that are usually scarce, limiting the reach and sustainability of vector control programs. This generates a need to focus on areas at risk of disease transmission and also understand the factors that might modulate local mosquito abundance. We evaluated the eco-bio-social factors that modulate indoor and outdoor relative abundance of female Ae. aegypti in communities of South Texas. We conducted housing quality and Knowledge Attitudes and Practices surveys in households that were part of a weekly mosquito surveillance program in November of 2017 and 2018. Our results showed widespread knowledge of mosquitoes and Zika virus by our participants. However, less than 35% considered them as serious problems in this region. The presence of window-mounted air conditioning units increased the risk of female mosquito relative abundance indoors. An increase in outdoor relative abundance was associated with larger properties and a higher number of children between 6 to 17 years of age. Interestingly, we observed that an increasing number of children <5 years of age modulated both indoor and outdoor relative abundance, with a 52% increase indoors and 30% decrease outdoors. The low perception of mosquito and disease risk highlights engagement needs for vector-borne disease prevention in this region. The identified risk factors can help guide public health officials in their efforts to reduce human and vector contact.

The question of how Zika virus (ZIKV) changed from a seemingly mild virus to a human pathogen capable of microcephaly and sexual transmission remains unanswered. The unexpected emergence of ZIKV's pathogenicity and capacity for sexual transmission may be due to genetic changes, and future changes in phenotype may continue to occur as the virus expands its geographic range. Alternatively, the sheer size of the 2015-16 epidemic may have brought attention to a pre-existing virulent ZIKV phenotype in a highly susceptible population. Thus, it is important to identify patterns of genetic change that may yield a better understanding of ZIKV emergence and evolution. However, because ZIKV has an RNA genome and a polymerase incapable of proofreading, it undergoes rapid mutation which makes it difficult to identify combinations of mutations associated with viral emergence. As next generation sequencing technology has allowed whole genome consensus and variant sequence data to be generated for numerous virus samples, the task of analyzing these genomes for patterns of mutation has become more complex. However, understanding which combinations of mutations spread widely and become established in new geographic regions versus those that disappear relatively quickly is essential for defining the trajectory of an ongoing epidemic. In this study, multiscale analysis of the wealth of genomic data generated over the course of the epidemic combined with in vivo laboratory data allowed trends in mutations and outbreak trajectory to be assessed. Mutations were detected throughout the genome via deep sequencing, and many variants appeared in multiple samples and in some cases become consensus. Similarly, amino acids that were previously consensus in pre-outbreak samples were detected as low frequency variants in epidemic strains. Protein structural models indicate that most of the mutations associated with the epidemic transmission occur on the exposed surface of viral proteins. At the macroscale level, consensus data was organized into large and interactive databases to allow the spread of individual mutations and combinations of mutations to be visualized and assessed for temporal and geographical patterns. Thus, the use of multiscale modeling for identifying mutations or combinations of mutations that impact epidemic transmission and phenotypic impact can aid the formation of hypotheses which can then be tested using reverse genetics.

South Texas has experienced local transmission of Zika virus and of other mosquito-borne viruses such as chikungunya virus and dengue virus in the last decades. Using a mosquito surveillance program in the Lower Rio Grande Valley (LRGV) and San Antonio, TX, from 2016 to 2018, we detected the presence of an insect-specific virus, cell fusing agent virus (CFAV), in the Aedes aegypti mosquito population. We tested 6,326 females and 1,249 males from the LRGV and 659 females from San Antonio for CFAV by RT-PCR using specific primers. Infection rates varied from 0 to 261 per 1,000 mosquitoes in the LRGV and 115 to 208 per 1,000 in San Antonio depending on the month of collection. Infection rates per 1,000 individuals appeared higher in females collected from BG Sentinel 2 traps compared to Autocidal Gravid Ovitraps, but the ratio of the percentage of infected pools did not differ by trap type. The natural viral load in individual males ranged from 1.25 x 10² to 5.50 x 10⁶ RNA copies and in unfed females from 5.42 x 10³ to 8.70 x 10⁶ RNA copies. Gravid females were found to harbor fewer viral particles than males and unfed females.

Genetic analysis of intra-host viral populations provides unique insight into pre-emergent mutations that may contribute to the genotype of future variants. Clinical samples positive for SARS-CoV-2 collected in California during the first months of the pandemic were sequenced to define the dynamics of mutation emergence as the virus became established in the state. Deep sequencing of 90 nasopharyngeal samples showed that many mutations associated with the establishment of SARS-CoV-2 globally were present at varying frequencies in a majority of the samples, even those collected as the virus was first detected in the US. A subset of mutations that emerged months later in consensus sequences were detected as subconsensus members of intra-host populations. Spike mutations P681H, H655Y, and V1104L were detected prior to emergence in variant genotypes, mutations were detected at multiple positions within the furin cleavage site, and pre-emergent mutations were identified in the nucleocapsid and the envelope genes. Because many of the samples had a very high depth of coverage, a bioinformatics pipeline, "Mappgene", was established that uses both iVar and LoFreq variant calling to enable identification of very low-frequency variants. This enabled detection of a spike protein deletion present in many samples at low frequency and associated with a variant of concern.

Mosquito-borne viruses are emerging or re-emerging globally, afflicting millions of people around the world. Aedes aegypti, the yellow fever mosquito, is the principal vector of dengue, Zika, and chikungunya viruses, and has well-established populations across tropical and subtropical urban areas of the Americas, including the southern United States. While intense arboviral epidemics have occurred in Mexico and further south in the Americas, local transmission in the United States has been minimal. Here, we study Ae. aegypti and Culex quinquefasciatus host feeding patterns and vertebrate host communities in residential environments of South Texas to identify host-utilization relative to availability. Only 31% of Ae. aegypti blood meals were derived from humans, while 50% were from dogs and 19% from other wild and domestic animals. In Cx. quinquefasciatus, 67% of blood meals were derived from chicken, 22% came from dogs, 9% from various wild avian species, and 2% from other mammals including one human, one cat, and one pig. We developed a model for the reproductive number, R₀, for Zika virus (ZIKV) in South Texas relative to northern Mexico using human disease data from Tamaulipas, Mexico. We show that ZIKV R₀ in South Texas communities could be greater than one if the risk of human exposure to Ae. aegypti bites in these communities is at least 60% that of Northern Mexico communities. The high utilization of non-human vertebrates and low risk of human exposure in South Texas diminishes the outbreak potential for human-amplified urban arboviruses transmitted by Ae. aegypti.