Prediction of species composition ratios in pooled specimens of the Anopheles Hyrcanus group using quantitative sequencing.

BACKGROUND: Plasmodium vivax is transmitted by members of the Anopheles Hyrcanus Group that includes six species in the Republic of Korea: Anopheles sinensis sensu stricto (s.s.), Anopheles pullus, Anopheles kleini, Anopheles belenrae, Anopheles lesteri, and Anopheles sineroides. Individual Anopheles species within the Hyrcanus Group demonstrate differences in their geographical distributions, vector competence and insecticide resistance, making it crucial for accurate species identification. Conventional species identification conducted using individual genotyping (or barcoding) based on species-specific molecular markers requires extensive time commitment and financial resources. RESULTS: A population-based quantitative sequencing (QS) protocol developed in this study provided a rapid estimate of species composition ratios among pooled mosquitoes as a cost-effective alternative to individual genotyping. This can be accomplished by using species- or group-specific nucleotide sequences of the mitochondrial cytochrome C oxidase subunit I (COI) and the ribosomal RNA internal transcribed spacer 2 (ITS2) region as species identification alleles in a two-step prediction protocol. Standard genomic DNA fragments of COI and ITS2 genes were amplified from each Anopheles species using group-specific universal primer sets. Following sequencing of the COI or ITS2 amplicons generated from sets of standard DNA mixtures, equations were generated via linear regression to predict species-specific nucleotide sequence frequencies at different positions. Species composition ratios between An. sineroides, An. pullus and An. lesteri were estimated from QS of the COI amplicons based on the mC.260A, mC.122C and mC.525C alleles at the first step, followed by the prediction of species composition ratios between An. sinensis, An. kleini and An. belenrae based on QS of the ITS2 amplicons using the rI.370G and rI.389T alleles. The COI copy number was not significantly different between species, suggesting the reliability of COI-based prediction. In contrast, ITS2 showed a slightly but significantly higher copy number in An. belenrae, requiring an adjustment of its predicted composition ratio. A blind test proved that predicted species composition ratios either from pooled DNA specimens or pooled mosquito specimens were not statistically different from the actual values, demonstrating that the QS-based prediction is accurate and reliable. CONCLUSIONS: This two-step prediction protocol will facilitate rapid estimation of the species composition ratios in field-collected Anopheles Hyrcanus Group populations and is particularly useful for studying the vector ecology of Anopheles population and epidemiology of malaria.

Changes in intestinal flora in patients with extrahepatic cholangiocarcinoma / 临床肝胆病杂志

Objective To investigate the changes in intestinal flora in patients with extrahepatic cholangiocarcinoma (ECC) and related influencing factors. Methods Fecal samples were collected from 16 patients with ECC who were hospitalized and treated in Department of Hepatobiliary and Pancreatic Surgery, The First Affiliated Hospital of Zhengzhou University, from January to December 2020, and absolute quantitative bacterial 16S rRNA was used for sequencing. A comparison was made with 20 patients with common bile duct stones (CBDS group) and 10 healthy controls (normal group), and the three groups were compared in terms of the differences in intestinal flora and the association with clinical indices. A one-way analysis of variance was used for comparison of normally distributed data with homogeneity of variance between the three groups, the t -test was used for comparison between two groups; the Kruskal-Wallis H test was used for comparison of non-normally distributed data between the three groups, and the Nemenyi test was used for comparison between groups. The chi-square test or the Fisher's exact test were used for comparison of categorical data between the three groups, and a Spearman correlation analysis was used to investigate correlation. Results The ECC group had significantly higher levels of total bilirubin (TBil) and direct bilirubin (DBil) than the CBDS group and the normal group. The α diversity analysis showed that there were no significant differences in observed species, Chao1 Index, and Shannon Index between the three groups (all P > 0.05), while there were significant differences in Shannon Index and Simpson Index between the three groups. The ECC group had a similar species diversity to the normal group and a significantly greater species diversity than the CBDS group ( P < 0.05), and the CBDS group had a significantly greater species diversity than the normal group ( P < 0.05). The β diversity analysis showed that the structure of intestinal flora in the ECC group was significantly different from that in the normal group and the CBDS group ( P < 0.05). The analysis of the difference in bacterial composition showed that Prevotella , Lactobacillus , Megasphaera , and Sutterella were significantly enriched in the ECC group. The correlation analysis showed that Prevotella was negatively correlated with the use of antibiotics, acid inhibitors, and liver-protecting drugs, and Lactobacillus , Megasphaera , and Sutterella were positively correlated with TBil and DBil. Conclusion There is a significant change in intestinal flora in patients with ECC, which is closely associated with liver function and the use of drugs.

Regional and seasonal detection of resistance mutation frequencies in field populations of Anopheles Hyrcanus Group and Culex pipiens complex in Korea.

Pyrethroid (PYR) and organophosphate (OP) insecticides have been extensively used for mosquito control for several decades in South Korea, and has resulted in the rapid development of resistance in the field. In this study, quantitative sequencing (QS) protocols were developed for the frequency prediction of insecticide resistance alleles [e.g., the L1014F/C mutation on the voltage sensitive sodium channel as a PYR resistance allele and the G119S mutation on the acetylcholinesterase 1 as OP resistance alleles] in four regional populations of Anopheles Hyrcanus Group and Culex pipiens complex. Both of the L1014F/C and G119S mutations were observed in all examined regional populations of An. Hyrcanus Group, suggesting a wide distribution of both PYR and OP resistance. In contrast, populations of the Cx. pipiens complex were determined to possess almost no G119S mutation, but relatively higher frequencies of the L1014F mutation, showing differential resistance patterns between different mosquito groups. The mutation frequencies were also monitored throughout a mosquito season (May-October) at one collection site to determine the seasonal changes of resistance mutation frequency in mosquito populations. Dramatic decreases of both L1014F/C and G119S mutation frequencies were observed in the An. Hyrcanus Group toward the fall, with no mutations observed in the early spring, suggesting a connection between the fitness costs of overwintering and insecticide resistance. However, no apparent trends were detectable in the Cx. pipiens complex populations due to low or zero mutation frequencies.

Monitoring of Pyrethroid Resistance Allele Frequency in the Common Bed Bug (Cimex lectularius) in the Republic of Korea.

Two-point mutations (V419L and L925I) on the voltage-sensitive sodium channel of bed bugs (Cimex lectularius) are known to confer pyrethroid resistance. To determine the status of pyrethroid resistance in bed bugs in Korea, resistance allele frequencies of bed bug strains collected from several US military installations in Korea and Mokpo, Jeollanamdo, from 2009-2019 were monitored using a quantitative sequencing. Most bed bugs were determined to have both of the point mutations except a few specimens, collected in 2009, 2012 and 2014, having only a single point mutation (L925I). No susceptible allele was observed in any of the bed bugs examined, suggesting that pyrethroid resistance in bed bug populations in Korea has reached a serious level. Large scale monitoring is required to increase our knowledge on the distribution and prevalence of pyrethroid resistance in bed bug populations in Korea. Based on present study, it is urgent to restrict the use of pyrethroids and to introduce effective alternative insecticides. A nation-wide monitoring program to determine the pyrethroid resistance level in bed bugs and to select alternative insecticides should be implemented.

Pollination Drop Proteome and Reproductive Organ Transcriptome Comparison in Gnetum Reveals Entomophilous Adaptation.

Gnetum possesses morphologically bisexual but functionally unisexual reproductive structures that exude sugary pollination drops to attract insects. Previous studies have revealed that the arborescent species (G. gnemon L.) and the lianoid species (G. luofuense C.Y.Cheng) possess different pollination syndromes. This study compared the proteome in the pollination drops of these two species using label-free quantitative techniques. The transcriptomes of fertile reproductive units (FRUs) and sterile reproductive units (SRUs) for each species were furthermore compared using Illumina Hiseq sequencing, and integrated proteomic and transcriptomic analyses were subsequently performed. Our results show that the differentially expressed proteins between FRUs and SRUs were involved in carbohydrate metabolism, the biosynthesis of amino acids and ovule defense. In addition, the differentially expressed genes between the FRUs and SRUs (e.g., MADS-box genes) were engaged in reproductive development and the formation of pollination drops. The integrated protein-transcript analyses revealed that FRUs and their exudates were relatively conservative while the SRUs and their exudates were more diverse, probably functioning as pollinator attractants. The evolution of reproductive organs appears to be synchronized with changes in the pollination drop proteome of Gnetum, suggesting that insect-pollinated adaptations are not restricted to angiosperms but also occur in gymnosperms.

Estimation of the influence of sequencing errors and distribution of random-sequence tags on quantitative sequencing.

To simultaneously sequence and quantify target DNA, quantitative sequencing (qSeq) employs stochastic labeling of target DNA molecules with random-sequence tags (RSTs). This recently developed approach allows parallel quantification of hundreds of microorganisms in natural habitats in a single sequencing run. Yet, no study has addressed to what extent sequencing errors affect quantification and how many sequence reads are needed for quantification. Here, we addressed those issues by using numerical simulations and experimental data from second-generation sequencing of various RSTs. We found that heterogeneous distribution of observed RSTs affected the number of sequence reads required to quantitate target genes, whereas the effect of sequencing errors is smaller than of the RSTs distribution. Because of the heterogeneous RSTs distribution, 15-fold more sequence reads than the number of observed RSTs should be obtained to retrieve almost all RSTs needed for quantification; in that case, quantification error is estimated to be within 5%.

Expansion of the Knockdown Resistance Frequency Map for Human Head Lice (Phthiraptera: Pediculidae) in the United States Using Quantitative Sequencing.

Pediculosis is a prevalent parasitic infestation of humans, which is increasing due, in part, to the selection of lice resistant to either the pyrethrins or pyrethroid insecticides by the knockdown resistance (kdr) mechanism. To determine the extent and magnitude of the kdr-type mutations responsible for this resistance, lice were collected from 138 collection sites in 48 U.S. states from 22 July 2013 to 11 May 2015 and analyzed by quantitative sequencing. Previously published data were used for comparisons of the changes in the frequency of the kdr-type mutations over time. Mean percent resistance allele frequency (mean % RAF) values across the three mutation loci were determined from each collection site. The overall mean % RAF (±SD) for all analyzed lice was 98.3 ± 10%. 132/138 sites (95.6%) had a mean % RAF of 100%, five sites (3.7%) had intermediate values, and only a single site had no mutations (0.0%). Forty-two states (88%) had a mean % RAF of 100%. The frequencies of kdr-type mutations did not differ regardless of the human population size that the lice were collected from, indicating a uniformly high level of resistant alleles. The loss of efficacy of the Nix formulation (Prestige Brand, Tarrytown, NY) from 1998 to 2013 was correlated to the increase in kdr-type mutations. These data provide a plausible reason for the decrease in the effectiveness of permethrin in the Nix formulation, which is the parallel increase of kdr-type mutations in lice over time.

Toxicodynamic mechanisms and monitoring of acaricide resistance in the two-spotted spider mite.

The two-spotted spider (Tetranychus urticae) is one of the most serious pests world-wide and has developed resistance to many types of acaricides. Various mutations on acaricide target site genes have been determined to be responsible for toxicodynamic resistance, and the genotyping and frequency prediction of these mutations can be employed as an alternative resistance monitoring strategy. A quantitative sequencing (QS) protocol was reported as a population-based genotyping technique, and applied for the determination of resistance allele frequencies in T. urticae field populations. In addition, a modified glass vial bioassay method (residual contact vial bioassay, RCV) was implemented as a rapid on-site resistance monitoring tool. The QS protocol, together with the RCV, would greatly facilitate monitoring of T. urticae resistance. Recent completion of T. urticae genome analysis should facilitate the identification of additional resistance genetic markers that can be employed for molecular resistance monitoring.

Allele-specific qRT-PCR demonstrates superior detection of single nucleotide polymorphisms as genetic markers for West Nile virus compared to Luminex® and quantitative sequencing.

To enable in vivo and in vitro competitive fitness comparisons among West Nile viruses (WNV), three reference viruses were marked genetically by site-directed mutagenesis with five synonymous nucleotide substitutions in the envelope gene region of the genome. Phenotypic neutrality of the mutants was assessed experimentally by competitive replication in cell culture and genetic stability of the substituted nucleotides was confirmed by direct sequencing. Luminex(®) technology, quantitative sequencing and quantitative RT-PCR (qRT-PCR) were compared in regard to specificity, sensitivity and accuracy for quantitation of wildtype and genetically marked viruses in mixed samples based on RNA obtained from samples of known viral titers. Although Luminex(®) technology and quantitative sequencing provided semi-quantitative or qualitative measurements, a sequence-specific primer extension approach using a specific reverse primer set in singleplex qRT-PCR demonstrated the best quantitation and specificity in the detection of RNA from wildtype and mutant viruses.