Understanding disease burden and transmission dynamics in resource-limited, developing countries like Nepal is often challenging due to a lack of adequate surveillance systems. These issues are exacerbated by limited access to diagnostic and research facilities throughout the country. Nepal has one of the highest COVID-19 case rates (915 cases per 100,000 people) in South Asia, with densely-populated Kathmandu experiencing the highest number of cases. Swiftly identifying case clusters and introducing effective intervention programs is crucial to mounting an effective containment strategy. The rapid identification of circulating SARS-CoV-2 variants can also provide important information on viral evolution and epidemiology. Genomic-based environmental surveillance can help in the early detection of outbreaks before clinical cases are recognized, and identify viral micro-diversity that can be used for designing real-time risk-based interventions. This research aimed to develop a genomic-based environmental surveillance system by detecting and characterizing SARS-CoV-2 in sewage samples of Kathmandu using portable next-generation DNA sequencing devices. Out of 20 selected sites in the Kathmandu Valley, sewage samples from 16 (80%) sites had detectable SARS-CoV-2. A heat-map was created to visualize transmission activity in the community based on viral load intensity and corresponding geospatial data. Further, 41 mutations were observed in the SARS-CoV-2 genome. Some detected mutations (n=9, 2%) were novel and yet to be reported in the global database, with one indicating a frameshift deletion in the spike gene. We also observed more transition than transversion on detected mutations, indicating rapid viral evolution in the host. Our study has demonstrated the feasibility of rapidly obtaining vital information on community transmission and disease dynamics of SARS-CoV-2 using genomic-based environmental surveillance.

Tiger (Panthera tigris) populations are in danger across their entire range due to habitat loss, poaching and the demand for tiger parts. The Bengal tiger (Panthera tigris tigris) is an endangered apex predator with a population size estimated to be less than 200 in Nepal. In spite of strict wildlife protection laws, illegal trade of tiger parts is increasing; and Nepal has become one of the major sources and transit routes for poached wildlife parts. Identification of wildlife parts is often challenging for law enforcement officials due to inadequate training and lack of available tools. Here, we describe a molecular forensic approach to gain insight into illegally trafficked tiger parts seized across Nepal. We created Nepal’s first comprehensive reference genetic database of wild tigers through the Nepal Tiger Genome Project (2011–2013). This database has nuclear DNA microsatellite genotype and sex profiles, including geo-spatial information, of over 60% (n = 120) of the wild tigers of Nepal. We analyzed 15 putative cases of confiscated poached tiger parts and all were confirmed to be of tiger. Ten samples were identified as male and five were female. We determined probable geo-source location for 9 of the 14 samples with 6–8 nuclear DNA microsatellite loci using inferences from four different statistical assignment methods. Six samples were assigned to Bardia National Park and one of these was an exact match to a female tiger previously profiled in our fecal DNA reference database. Two tiger samples were assigned to Shuklaphanta Wildlife Reserve and one to Chitwan National Park. We are unable to definitively assign five tiger samples which could be offspring dispersers or might have come from tiger population outside of Nepal. Our study revealed that the western region, particularly Bardia National Park, is a poaching hotspot for illegal tiger trade in Nepal. We present feasibility of using molecular forensic based evidence to incriminate criminals in a court of law in the fight against wildlife crime.

BackgroundCarbonic anhydrases (CAs) are ubiquitous, essential enzymes which catalyze the conversion of carbon dioxide and water to bicarbonate and H+ ions. Vertebrate genomes generally contain gene loci for 15–21 different CA isoforms, three of which are enzymatically inactive. CA VI is the only secretory protein of the enzymatically active isoforms. We discovered that non-mammalian CA VI contains a C-terminal pentraxin (PTX) domain, a novel combination for both CAs and PTXs.MethodsWe isolated and sequenced zebrafish (Danio rerio) CA VI cDNA, complete with the sequence coding for the PTX domain, and produced the recombinant CA VI–PTX protein. Enzymatic activity and kinetic parameters were measured with a stopped-flow instrument. Mass spectrometry, analytical gel filtration and dynamic light scattering were used for biophysical characterization. Sequence analyses and Bayesian phylogenetics were used in generating hypotheses of protein structure and CA VI gene evolution. A CA VI–PTX antiserum was produced, and the expression of CA VI protein was studied by immunohistochemistry. A knock-down zebrafish model was constructed, and larvae were observed up to five days post-fertilization (dpf). The expression of ca6 mRNA was quantitated by qRT-PCR in different developmental times in morphant and wild-type larvae and in different adult fish tissues. Finally, the swimming behavior of the morphant fish was compared to that of wild-type fish.ResultsThe recombinant enzyme has a very high carbonate dehydratase activity. Sequencing confirms a 530-residue protein identical to one of the predicted proteins in the Ensembl database (ensembl.org). The protein is pentameric in solution, as studied by gel filtration and light scattering, presumably joined by the PTX domains. Mass spectrometry confirms the predicted signal peptide cleavage and disulfides, and N-glycosylation in two of the four observed glycosylation motifs. Molecular modeling of the pentamer is consistent with the modifications observed in mass spectrometry. Phylogenetics and sequence analyses provide a consistent hypothesis of the evolutionary history of domains associated with CA VI in mammals and non-mammals. Briefly, the evidence suggests that ancestral CA VI was a transmembrane protein, the exon coding for the cytoplasmic domain was replaced by one coding for PTX domain, and finally, in the therian lineage, the PTX-coding exon was lost. We knocked down CA VI expression in zebrafish embryos with antisense morpholino oligonucleotides, resulting in phenotype features of decreased buoyancy and swim bladder deflation in 4 dpf larvae.DiscussionThese findings provide novel insights into the evolution, structure, and function of this unique CA form.

Bengal tigers ( Panthera tigris tigris ) serve a pivotal role as an apex predator in forest ecosystems. To increase our knowledge on factors impacting the viability and health of this endangered species, we studied the gut microbiota in 32 individual Bengal tigers from three geographically separated areas (Chitwan National Park (CNP), Bardia National Park (BNP) and Suklaphanta Wildlife Reserve (SWR)) in Nepal, using noninvasive genetic sampling methods. Gut microbiota influence the immune system, impact various physiological functions, and modulates metabolic reactions, that ultimately impact the host health, behavior and development. Across the tiger populations in Nepal, we found significant differences in the composition of microbial communities based on their geographic locations. Specifically, we detected significant differences between CNP and the other two protected areas (CNP vs BNP: pseudo t = 1.944, P = 0.006; CNP vs SWR: pseudo t = 1.9942, P = 0.0071), but no differences between BNP and SWR. This mirrors what has been found for tiger gene flow in the same populations, suggesting gut microbiota composition and host gene flow may be linked. Furthermore, predictive metagenome functional content analysis (PICRUSt) revealed a higher functional enrichment and diversity for significant gut microbiota in the Chitwan tiger population and the lowest enrichment and diversity in Suklaphanta. The CNP tiger population contained higher proportions of microbiota that are associated with predicted functions relevant for metabolism of amino acid, lipid, xenobiotics biodegradation, terpenoides and polyketides than the SWR population. We conclude the tiger population structure, gut microbiota profile and associated functional metabolic categories are correlated, with geographically most separated CNP and SWR tiger population having the most distinct and different host genotype and microbiota profiles. Our work dramatically expands the understanding of tiger microbiota in wild populations and provides a valuable case study on how to investigate genetic diversity at different hierarchical levels, including hosts as well as their microbial communities.