Total infectome characterization of respiratory infections in pre-COVID-19 Wuhan, China.

At the end of 2019 Wuhan witnessed an outbreak of "atypical pneumonia" that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their "total infectome", including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen-Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016-2017.

SvSTL1 in the large subunit family of ribonucleotide reductases plays a major role in chloroplast development of Setaria viridis.

Ribonucleotide reductases (RNRs) are essential enzymes in DNA synthesis. However, little is known about the RNRs in plants. Here, we identified a svstl1 mutant from the self-created ethyl methanesulfonate (EMS) mutant library of Setaria viridis. The mutant leaves exhibited a bleaching phenotype at the heading stage. Paraffin section analysis showed the destruction of the C4 Kranz anatomy. Transmission electron microscopy results further demonstrated the severely disturbed development of some chloroplasts. MutMap analysis revealed that the SvSTL1 gene is the primary candidate, encoding a large subunit of RNRs. Complementation experiments confirmed that SvSTL1 is responsible for the phenotype of svstl1. There are two additional RNR large subunit homologs in S. viridis, SvSTL2 and SvSTL3. To further understand the functions of these three RNR large subunit genes, a series of mutants were generated via CRISPR/Cas9 technology. In striking contrast to the finding that all three SvSTLs interact with the RNR small subunit, the phenotype varied along with the copies of chloroplast genome among different svstl single mutants: the svstl1 mutant exhibited pronounced chloroplast development and significantly fewer copies of the chloroplast genome than the svstl2 or svstl3 single mutants. These results suggested that SvSTL1 plays a major role in the optimal function of RNRs and is essential for chloroplast development. Furthermore, through the analysis of double and triple mutants, the study provides new insights into the finely tuned coordination among SvSTLs to maintain normal chloroplast development in the emerging C4 model plant S. viridis.

Optical design mode based on fast automatic design process for freeform reflective imaging systems with modest FOV.

Traditional optical design methods require designer intervention in the system's evolution from the starting point to the final design. Trial-and-error during design optimization improves system performance step by step but requires much time and effort. A new optical design framework, end-to-end fast automatic design, is proposed and achieved for the freeform reflective optics in this paper, which promotes a new optical design mode. Compared with the traditional mode through improving performance after each trial, an optical system with good image quality can be directly obtained in the end-to-end design process with simple input and no human involvement within a short time. If there is still the possibility for performance improvement of the obtained system, the designer can vary the input parameters repeatedly to obtain multiple systems with good image quality. Finally, the desired system is selected from these systems. Compared with the step-by-step trials in traditional optimization, this new optical design mode involves high-speed trials of the end-to-end automatic design process, reducing the dependence on experience and skill. In this paper, an end-to-end fast automatic design method for freeform imaging systems is developed based on a new design route. Using an initial plane system as an input, a freeform system with excellent image quality can be designed automatically within 1-2 min. After several trials of the end-to-end fast design process, three high-performance freeform systems are designed successfully that consider volume control, beam obscuration, and mirror interference.

Freeform imaging system with resolution that varies with the field angle in two dimensions.

The human eye's resolution varies with the field angle and has high center resolution and low edge resolution characteristics. In this paper, a freeform imaging system is presented that has resolution distribution characteristics similar to those of the human eye. Field-dependent parameters are used to describe the system's optical properties and a direct design method is proposed to realize the novel functionality. An off-axis reflective freeform imaging system with high center resolution and low edge resolution within a square 30°×30° field of view (FOV) is designed using this method. The maximum instantaneous field of view (IFOV) ratio of center field resolution to edge field resolution is 0.47. Only three freeform surfaces are used to attain good image quality. Simultaneous improvements are observed in both resolution and FOV while the detector remains fixed.

Fast automatic design method for freeform imaging systems through system construction and correction.

In traditional optical design, a starting point is selected and coefficients optimization is then performed using software. The process requires considerable time and the involvement of a human with design skills and experience. In this Letter, a fast automatic method for freeform imaging systems design is proposed. Using a plane system as the input, a freeform optical system with high image quality can be designed automatically at high speed. The method consists of system construction and system correction, combining the advantages of the direct design method and the methods based on aberration analysis. After system construction generates a system with fundamental optical parameters, system correction is an iterative process that alternates between image plane correction and surfaces correction to improve the image quality to a high level. Two examples required 5 min 56 s and 6 min 10 s to design freeform systems with near-diffraction-limit image quality.

Role of systemic injection of rabies immunoglobulin in rabies vaccination.

To determine the role of systemic injection of rabies immunoglobulin (RIG) in rabies vaccination, we analyzed the level of antibody against rabies virus in the serum of mice that received various doses of RIG combined with rabies vaccine. Our results indicate that systemic injection of RIG does not contribute detectably to passive or adaptive immunization, suggesting that the main function of RIG in individuals with category III exposure is to neutralize rabies virus via immediate local infiltration of the wound.

Probing biological nanotopology via diffusion of weakly constrained plasmonic nanorods with optical coherence tomography.

Biological materials exhibit complex nanotopology, i.e., a composite liquid and solid phase structure that is heterogeneous on the nanoscale. The diffusion of nanoparticles in nanotopological environments can elucidate biophysical changes associated with pathogenesis and disease progression. However, there is a lack of methods that characterize nanoprobe diffusion and translate easily to in vivo studies. Here, we demonstrate a method based on optical coherence tomography (OCT) to depth-resolve diffusion of plasmon-resonant gold nanorods (GNRs) that are weakly constrained by the biological tissue. By using GNRs that are on the size scale of the polymeric mesh, their Brownian motion is minimally hindered by intermittent collisions with local macromolecules. OCT depth-resolves the particle-averaged translational diffusion coefficient (DT) of GNRs within each coherence volume, which is separable from the nonequilibrium motile activities of cells based on the unique polarized light-scattering properties of GNRs. We show how this enables minimally invasive imaging and monitoring of nanotopological changes in a variety of biological models, including extracellular matrix (ECM) remodeling as relevant to carcinogenesis, and dehydration of pulmonary mucus as relevant to cystic fibrosis. In 3D ECM models, DT of GNRs decreases with both increasing collagen concentration and cell density. Similarly, DT of GNRs is sensitive to human bronchial-epithelial mucus concentration over a physiologically relevant range. This novel method comprises a broad-based platform for studying heterogeneous nanotopology, as distinct from bulk viscoelasticity, in biological milieu.

Transcriptomic Investigation of the Virus Spectrum Carried by Midges in Border Areas of Yunnan Province.

Yunnan province in China shares its borders with three neighboring countries: Myanmar, Vietnam, and Laos. The region is characterized by a diverse climate and is known to be a suitable habitat for various arthropods, including midges which are notorious for transmitting diseases which pose significant health burdens affecting both human and animal health. A total of 431,100 midges were collected from 15 different locations in the border region of Yunnan province from 2015 to 2020. These midges were divided into 37 groups according to the collection year and sampling site. These 37 groups of midges were then homogenized to extract nucleic acid. Metatranscriptomics were used to analyze their viromes. Based on the obtained cytochrome C oxidase I gene (COI) sequences, three genera were identified, including one species of Forcipomyia, one species of Dasyhelea, and twenty-five species of Culicoides. We identified a total of 3199 viruses in five orders and 12 families, including 1305 single-stranded positive-stranded RNA viruses (+ssRNA) in two orders and seven families, 175 single-stranded negative-stranded RNA viruses (-ssRNA) in two orders and one family, and 1719 double-stranded RNA viruses in five families. Six arboviruses of economic importance were identified, namely Banna virus (BAV), Japanese encephalitis virus (JEV), Akabane virus (AKV), Bluetongue virus (BTV), Tibetan circovirus (TIBOV), and Epizootic hemorrhagic disease virus (EHDV), all of which are capable, to varying extents, of causing disease in humans and/or animals. The survey sites in this study basically covered the current distribution area of midges in Yunnan province, which helps to predict the geographic expansion of midge species. The complexity and diversity of the viral spectrum carried by midges identified in the study calls for more in-depth research, which can be utilized to monitor arthropod vectors and to predict the emergence and spread of zoonoses and animal epidemics, which is of great significance for the control of vector-borne diseases.

Total Infectome Characterization of Respiratory Infections during the 2022-23 COVID-19 Outbreak in China Revealed Extensive Coinfections with Links to SARS-CoV-2 Status, Age, and Disease Severity.

Between 7 December 2022 and 28 February 2023, China experienced a new wave of COVID-19 that swept across the entire country and resulted in an increasing amount of respiratory infections and hospitalizations. The purpose of this study is to reveal the intensity and composition of coinfecting microbial agents. In total, 196 inpatients were recruited from The Third People's Hospital of Shenzhen, and 169 respiratory and 73 blood samples were collected for metagenomic next-generation sequencing. The total "Infectome" was characterized and compared across different groups defined by the SARS-CoV-2 detection status, age groups, and severity of disease. Our results revealed a total of 22 species of pathogenic microbes (4 viruses, 13 bacteria, and 5 fungi), and more were discovered in the respiratory tract than in blood. The diversity of the total infectome was highly distinguished between respiratory and blood samples, and it was generally higher in patients that were SARS-CoV-2-positive, older in age, and with more severe disease. At the individual pathogen level, HSV-1 seemed to be the major contributor to these differences observed in the overall comparisons. Collectively, this study reveals the highly complex respiratory infectome and high-intensity coinfection in patients admitted to the hospital during the period of the 2023 COVID-19 pandemic in China.

Meta-transcriptomic analysis of companion animal infectomes reveals their diversity and potential roles in animal and human disease.

Companion animals such as cats and dogs harbor diverse microbial communities that can potentially impact human health due to close and frequent contact. To better characterize their total infectomes and assess zoonotic risks, we characterized the overall infectomes of companion animals (cats and dogs) and evaluated their potential zoonotic risks. Meta-transcriptomic analyses were performed on 239 samples from cats and dogs collected across China, identifying 24 viral species, 270 bacterial genera, and two fungal genera. Differences in the overall microbiome and infectome composition were compared across different animal species (cats or dogs), sampling sites (rectal or oropharyngeal), and health status (healthy or diseased). Diversity analyses revealed that viral abundance was generally higher in diseased animals compared to healthy ones, while differences in microbial composition were mainly driven by sampling site, followed by animal species and health status. Disease association analyses validated the pathogenicity of known pathogens and suggested potential pathogenic roles of previously undescribed bacteria and newly discovered viruses. Cross-species transmission analyses identified seven pathogens shared between cats and dogs, such as alphacoronavirus 1, which was detected in both oropharyngeal and rectal swabs albeit with differential pathogenicity. Further analyses showed that some viruses, like alphacoronavirus 1, harbored multiple lineages exhibiting distinct pathogenicity, tissue, or host preferences. Ultimately, a systematic evolutionary screening identified 27 potential zoonotic pathogens in this sample set, with far more bacterial than viral species, implying potential health threats to humans. Overall, our meta-transcriptomic analysis reveals a landscape of actively transcribing microorganisms in major companion animals, highlighting key pathogens, those with the potential for cross-species transmission, and possible zoonotic threats. IMPORTANCE: This study provides a comprehensive characterization of the entire community of infectious microbes (viruses, bacteria, and fungi) in companion animals like cats and dogs, termed the "infectome." By analyzing hundreds of samples from across China, the researchers identified numerous known and novel pathogens, including 27 potential zoonotic agents that could pose health risks to both animals and humans. Notably, some of these zoonotic pathogens were detected even in apparently healthy pets, highlighting the importance of surveillance. The study also revealed key microbial factors associated with respiratory and gastrointestinal diseases in pets, as well as potential cross-species transmission events between cats and dogs. Overall, this work sheds light on the complex microbial landscapes of companion animals and their potential impacts on animal and human health, underscoring the need for monitoring and management of these infectious agents.

The clinical outcome of COVID-19 is strongly associated with microbiome dynamics in the upper respiratory tract.

OBJECTIVES: The respiratory tract is the portal of entry for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although a variety of respiratory pathogens other than SARS-CoV-2 have been associated with severe cases of COVID-19 disease, the dynamics of the upper respiratory microbiota during disease the course of disease, and how they impact disease manifestation, remain uncertain. METHODS: We collected 349 longitudinal upper respiratory samples from a cohort of 65 COVID-19 patients (cohort 1), 28 samples from 28 recovered COVID-19 patients (cohort 2), and 59 samples from 59 healthy controls (cohort 3). All COVID-19 patients originated from the earliest stage of the epidemic in Wuhan. Based on a modified clinical scale, the disease course was divided into five clinical disease phases (pseudotimes): "Healthy" (pseudotime 0), "Incremental" (pseudotime 1), "Critical" (pseudotime 2), "Complicated" (pseudotime 3), "Convalescent" (pseudotime 4), and "Long-term follow-up" (pseudotime 5). Using meta-transcriptomics, we investigated the features and dynamics of transcriptionally active microbes in the upper respiratory tract (URT) over the course of COVID-19 disease, as well as its association with disease progression and clinical outcomes. RESULTS: Our results revealed that the URT microbiome exhibits substantial heterogeneity during disease course. Two clusters of microbial communities characterized by low alpha diversity and enrichment for multiple pathogens or potential pathobionts (including Acinetobacter and Candida) were associated with disease progression and a worse clinical outcome. We also identified a series of microbial indicators that classified disease progression into more severe stages. Longitudinal analysis revealed that although the microbiome exhibited complex and changing patterns during COVID-19, a restoration of URT microbiomes from early dysbiosis toward more diverse status in later disease stages was observed in most patients. In addition, a group of potential pathobionts were strongly associated with the concentration of inflammatory indicators and mortality. CONCLUSION: This study revealed strong links between URT microbiome dynamics and disease progression and clinical outcomes in COVID-19, implying that the treatment of severe disease should consider the full spectrum of microbial pathogens present.

Metagenomic analysis of individual mosquito viromes reveals the geographical patterns and drivers of viral diversity.

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Here, using a meta-transcriptomic approach, we determined the viromes of 2,438 individual mosquitoes (81 species), spanning ~4,000 km along latitudes and longitudes in China. From these data we identified 393 viral species associated with mosquitoes, including 7 (putative) species of arthropod-borne viruses (that is, arboviruses). We identified potential mosquito species and geographic hotspots of viral diversity and arbovirus occurrence, and demonstrated that the composition of individual mosquito viromes was strongly associated with host phylogeny. Our data revealed a large number of viruses shared among mosquito species or genera, enhancing our understanding of the host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, perhaps reflecting long-distance mosquito dispersal. Together, these results greatly expand the known mosquito virome, linked viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the biogeography and diversity of viruses in insect vectors.

Motility-, autocorrelation-, and polarization-sensitive optical coherence tomography discriminates cells and gold nanorods within 3D tissue cultures.

We propose a method for differentiating classes of light scatterers based upon their temporal and polarization properties computed from time series of polarization-sensitive optical coherence tomography (PS-OCT) images. The amplitude (motility) and time scale (autocorrelation decay time) of the speckle fluctuations are combined with the cross-polarization pixel-wise to render Motility-, autocorrelation-, and polarization-sensitive (MAPS) OCT contrast images. This combination of metrics provides high specificity for discriminating diffusive gold nanorods and mammary epithelial cell spheroids within 3D tissue culture, based on their unique MAPS signature. This has implications toward highly specific contrast in molecular (nanoparticle-based) and functional (cellular activity) imaging using standard PS-OCT hardware.

Genetic Diversity and Evolution of Viruses Infecting Felis catus: A Global Perspective.

Cats harbor many important viral pathogens, and the knowledge of their diversity has been greatly expanded thanks to increasingly popular molecular sequencing techniques. While the diversity is mostly described in numerous regionally defined studies, there lacks a global overview of the diversity for the majority of cat viruses, and therefore our understanding of the evolution and epidemiology of these viruses was generally inadequate. In this study, we analyzed 12,377 genetic sequences from 25 cat virus species and conducted comprehensive phylodynamic analyses. It revealed, for the first time, the global diversity for all cat viruses known to date, taking into account highly virulent strains and vaccine strains. From there, we further characterized and compared the geographic expansion patterns, temporal dynamics and recombination frequencies of these viruses. While respiratory pathogens such as feline calicivirus showed some degree of geographical panmixes, the other viral species are more geographically defined. Furthermore, recombination rates were much higher in feline parvovirus, feline coronavirus, feline calicivirus and feline foamy virus than the other feline virus species. Collectively, our findings deepen the understanding of the evolutionary and epidemiological features of cat viruses, which in turn provide important insight into the prevention and control of cat pathogens.

Nelson Bay Reovirus Isolated from Bats and Blood-Sucking Arthropods Collected in Yunnan Province, China.

Nelson Bay reovirus (NBV) is an emerging zoonotic virus that can cause acute respiratory disease in humans. These viruses are mainly discovered in Oceania, Africa, and Asia, and bats have been identified as their main animal reservoir. However, despite recent expansion of diversity for NBVs, the transmission dynamics and evolutionary history of NBVs are still unclear. This study successfully isolated two NBV strains (MLBC1302 and MLBC1313) from blood-sucking bat fly specimens (Eucampsipoda sundaica) and one (WDBP1716) from the spleen specimen of a fruit bat (Rousettus leschenaultii), which were collected at the China-Myanmar border area of Yunnan Province. Syncytia cytopathic effects (CPE) were observed in BHK-21 and Vero E6 cells infected with the three strains at 48 h postinfection. Electron micrographs of ultrathin sections showed numerous spherical virions with a diameter of approximately 70 nm in the cytoplasm of infected cells. The complete genome nucleotide sequence of the viruses was determined by metatranscriptomic sequencing of infected cells. Phylogenetic analysis demonstrated that the novel strains were closely related to Cangyuan orthoreovirus, Melaka orthoreovirus, and human-infecting Pteropine orthoreovirus HK23629/07. Simplot analysis revealed the strains originated from complex genomic reassortment among different NBVs, suggesting the viruses experienced a high reassortment rate. In addition, strains successfully isolated from bat flies also implied that blood-sucking arthropods might serve as potential transmission vectors. IMPORTANCE Bats are the reservoir of many viral pathogens with strong pathogenicity, including NBVs. Nevertheless, it is unclear whether arthropod vectors are involved in transmitting NBVs. In this study, we successfully isolated two NBV strains from bat flies collected from the body surface of bats, which implies that they may be vectors for virus transmission between bats. While the potential threat to humans remains to be determined, evolutionary analyses involving different segments revealed that the novel strains had complex reassortment histories, with S1, S2, and M1 segments highly similar to human pathogens. Further experiments are required to determine whether more NBVs are vectored by bat flies, their potential threat to humans, and transmission dynamics.

Meta-transcriptomics for the diversity of tick-borne virus in Nujiang, Yunnan Province.

Ticks, an arthropod known for transmitting various pathogens such as viruses, bacteria, and fungi, pose a perpetual public health concern. A total of 2,570 ticks collected from Nujiang Prefecture in Yunnan Province between 2017 and 2022 were included in the study. Through the meta-transcriptomic sequencing of four locally distributed tick species, we identified 13 RNA viruses belonging to eight viral families, namely, Phenuiviridae, Nairoviridae, Peribunyaviridae, Flaviviridae, Chuviridae, Rhabdoviridae, Orthomyxoviridae, and Totiviridae. The most prevalent viruses were members of the order Bunyavirales, including three of Phenuiviridae, two were classified as Peribunyaviridae, and one was associated with Nairoviridae. However, whether they pose a threat to human health still remains unclear. Indeed, this study revealed the genetic diversity of tick species and tick-borne viruses in Nujiang Prefecture based on COI gene and tick-borne virus research. These data clarified the genetic evolution of some RNA viruses and furthered our understanding of the distribution pattern of tick-borne pathogens, highlighting the importance and necessity of monitoring tick-borne pathogens.

Origin, genomic diversity and evolution of African swine fever virus in East Asia.

Since 2018, the outbreaks of genotype II African swine fever virus (ASFV) in China and several eastern Asian countries have caused a huge impact on the local swine industry, resulting in huge economic losses. However, little is known about the origin, genomic diversity, evolutionary features, and epidemiological history of the genotype II ASFV. Here, 14 high-quality complete genomes of ASFVs were generated via sequencing of samples collected from China over the course of 3 years, followed by phylogenetic and phylodynamic analyses. The strains identified were relatively homogeneous, with a total of 52 SNPs and 11 indels compared with the prototype strain HLJ/2018, among which there were four exceptionally large deletions (620-18,023 nt). Evolutionary analyses revealed that ASFV strains distributed in eastern Asia formed a monophyly and a 'star-like' structure centered around the prototype strain, suggesting a single origin. Additionally, phylogenetic network analysis and ancestral reconstruction of geographic state indicated that genotype II ASFV strains in eastern Asia likely originated from Western Europe. Overall, these results contribute to the understanding of the history and current status of genotype II ASFV strains in eastern Asian, which could be of considerable importance in disease control and prevention.

Individual bat virome analysis reveals co-infection and spillover among bats and virus zoonotic potential.

Bats are reservoir hosts for many zoonotic viruses. Despite this, relatively little is known about the diversity and abundance of viruses within individual bats, and hence the frequency of virus co-infection and spillover among them. We characterize the mammal-associated viruses in 149 individual bats sampled from Yunnan province, China, using an unbiased meta-transcriptomics approach. This reveals a high frequency of virus co-infection (simultaneous infection of bat individuals by multiple viral species) and spillover among the animals studied, which may in turn facilitate virus recombination and reassortment. Of note, we identify five viral species that are likely to be pathogenic to humans or livestock, based on phylogenetic relatedness to known pathogens or in vitro receptor binding assays. This includes a novel recombinant SARS-like coronavirus that is closely related to both SARS-CoV and SARS-CoV-2. In vitro assays indicate that this recombinant virus can utilize the human ACE2 receptor such that it is likely to be of increased emergence risk. Our study highlights the common occurrence of co-infection and spillover of bat viruses and their implications for virus emergence.

Metagenomic analysis of individual mosquitos reveals the ecology of insect viruses.

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Using a meta-transcriptomic approach, we analysed the virome of 2,438 individual mosquitos (79 species), spanning ~4000 km along latitudes and longitudes in China. From these data we identified 393 core viral species associated with mosquitos, including seven (putative) arbovirus species. We identified potential species and geographic hotspots of viral richness and arbovirus occurrence, and demonstrated that host phylogeny had a strong impact on the composition of individual mosquito viromes. Our data revealed a large number of viruses shared among mosquito species or genera, expanding our knowledge of host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, possibly facilitated by long-distance mosquito migrations. Together, our results greatly expand the known mosquito virome, linked the viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the ecology of viruses of insect vectors.

A time-series meta-transcriptomic analysis reveals the seasonal, host, and gender structure of mosquito viromes.

Although metagenomic sequencing has revealed high numbers of viruses in mosquitoes sampled globally, our understanding of how their diversity and abundance varies in time and space as well as by host species and gender remains unclear. To address this, we collected 23,109 mosquitoes over the course of 12 months from a bat-dwelling cave and a nearby village in Yunnan province, China. These samples were organized by mosquito species, mosquito gender, and sampling time for meta-transcriptomic sequencing. A total of 162 eukaryotic virus species were identified, of which 101 were novel, including representatives of seventeen RNA virus multi-family supergroups and four species of DNA virus from the families Parvoviridae, Circoviridae, and Nudiviridae. In addition, two known vector-borne viruses-Japanese encephalitis virus and Banna virus-were found. Analyses of the entire virome revealed strikingly different viral compositions and abundance levels in warmer compared to colder months, a strong host structure at the level of mosquito species, and no substantial differences between those viruses harbored by male and female mosquitoes. At the scale of individual viruses, some were found to be ubiquitous throughout the year and across four mosquito species, while most of the other viruses were season and/or host specific. Collectively, this study reveals the diversity, dynamics, and evolution of the mosquito virome at a single location and sheds new lights on the ecology of these important vector animals.