BackgroundAlthough they are important disease vectors mosquito biodiversity in Pakistan is poorly known. Recent epidemics of dengue fever have revealed the need for more detailed understanding of the diversity and distributions of mosquito species in this region. DNA barcoding improves the accuracy of mosquito inventories because morphological differences between many species are subtle, leading to misidentifications.Methodology/Principal FindingsSequence variation in the barcode region of the mitochondrial COI gene was used to identify mosquito species, reveal genetic diversity, and map the distribution of the dengue-vector species in Pakistan. Analysis of 1684 mosquitoes from 491 sites in Punjab and Khyber Pakhtunkhwa during 2010–2013 revealed 32 species with the assemblage dominated by Culex quinquefasciatus (61% of the collection). The genus Aedes (Stegomyia) comprised 15% of the specimens, and was represented by six taxa with the two dengue vector species, Ae. albopictus and Ae. aegypti, dominant and broadly distributed. Anopheles made up another 6% of the catch with An. subpictus dominating. Barcode sequence divergence in conspecific specimens ranged from 0–2.4%, while congeneric species showed from 2.3–17.8% divergence. A global haplotype analysis of disease-vectors showed the presence of multiple haplotypes, although a single haplotype of each dengue-vector species was dominant in most countries. Geographic distribution of Ae. aegypti and Ae. albopictus showed the later species was dominant and found in both rural and urban environments.ConclusionsAs the first DNA-based analysis of mosquitoes in Pakistan, this study has begun the construction of a barcode reference library for the mosquitoes of this region. Levels of genetic diversity varied among species. Because of its capacity to differentiate species, even those with subtle morphological differences, DNA barcoding aids accurate tracking of vector populations.

BackgroundAlthough whiteflies (Bemisia tabaci complex) are an important pest of cotton in Pakistan, its taxonomic diversity is poorly understood. As DNA barcoding is an effective tool for resolving species complexes and analyzing species distributions, we used this approach to analyze genetic diversity in the B. tabaci complex and map the distribution of B. tabaci lineages in cotton growing areas of Pakistan.Methods/Principal FindingsSequence diversity in the DNA barcode region (mtCOI-5′) was examined in 593 whiteflies from Pakistan to determine the number of whitefly species and their distributions in the cotton-growing areas of Punjab and Sindh provinces. These new records were integrated with another 173 barcode sequences for B. tabaci, most from India, to better understand regional whitefly diversity. The Barcode Index Number (BIN) System assigned the 766 sequences to 15 BINs, including nine from Pakistan. Representative specimens of each Pakistan BIN were analyzed for mtCOI-3′ to allow their assignment to one of the putative species in the B. tabaci complex recognized on the basis of sequence variation in this gene region. This analysis revealed the presence of Asia II 1, Middle East-Asia Minor 1, Asia 1, Asia II 5, Asia II 7, and a new lineage “Pakistan”. The first two taxa were found in both Punjab and Sindh, but Asia 1 was only detected in Sindh, while Asia II 5, Asia II 7 and “Pakistan” were only present in Punjab. The haplotype networks showed that most haplotypes of Asia II 1, a species implicated in transmission of the cotton leaf curl virus, occurred in both India and Pakistan.ConclusionsDNA barcodes successfully discriminated cryptic species in B. tabaci complex. The dominant haplotypes in the B. tabaci complex were shared by India and Pakistan. Asia II 1 was previously restricted to Punjab, but is now the dominant lineage in southern Sindh; its southward spread may have serious implications for cotton plantations in this region.

Commiphora gileadensis and C. foliacea (family Burseraceae) are pantropical in nature and known for producing fragrant resin (myrrh). Both the tree species are economically and medicinally important however, least genomic understanding is available for this genus. Herein, we report the complete chloroplast genome sequences of C. gileadensis and C. foliacea and comparative analysis with related species (C. wightii and Boswellia sacra). A modified chloroplast DNA extraction method was adopted, followed with next generation sequencing, detailed bioinformatics and PCR analyses. The results revealed that the cp genome sizes of C. gileadensis and C. foliacea, are 160,268 and 160,249 bp, respectively, with classic quadripartite structures that comprises of inverted repeat’s pair. Overall, the organization of these cp genomes, GC contents, gene order, and codon usage were comparable to other cp genomes in angiosperm. Approximately, 198 and 175 perfect simple sequence repeats were detected in C. gileadensis and C. foliacea genomes, respectively. Similarly, 30 and 25 palindromic, 15 and 25 forward, and 20 and 25 tandem repeats were determined in both the cp genomes, respectively. Comparison of these complete cp genomes with C. wightii and B. sacra revealed significant sequence resemblance and comparatively highest deviation in intergenic spacers. The phylo-genomic comparison showed that C. gileadensis and C. foliacea form a single clade with previously reported C. wightii and B. sacra from family Burseraceae. Current study reports for the first time the cp genomics of species from Commiphora, which could be helpful in understanding genetic diversity and phylogeny of this myrrh producing species.

Plantago ovata (plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from plantaginaceae. the results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. the divergence time estimation as employed in BeASt revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other plantaginaceae species.

Tuberculosis (TB) is an important public health issue around the globe which is a chronic infectious disease and is still one of the major challenges for developing countries. The emergence of drug-resistant TB makes the condition worse and there is an urgent need of fast, highly sensitive diagnostic methods. This study was undertaken to evaluate the performance of GeneXpert® MTB/RIF assay and MTB culture for the detection of Mycobacterium tuberculosis (MTB) in sputum smear-negative pulmonary TB/drug-resistant tuberculosis (DR-TB) suspects. A total of 168 sputum smear-negative TB suspects were recruited for the study. Among the suspected TB cases, 52.98% were male and 47.02% were females with the mean age of 42 ± 17.6 years. All the sputum specimens collected from the study population were subjected to Ziehl–Neelsen (ZN) smear microscopy, GeneXpert MTB/RIF assay, and MTB culture. The results revealed that, out of 168 acid-fast bacilli (AFB)/ZN smear microscopy–negative sputum specimens, 48 (28.57%) and 58 (34.52%) were detected MTB positive by GeneXpert MTB/RIF assay and MTB culture, respectively, while 120 (71.43%) and 110 (65.48%) suspected TB cases were confirmed negative by GeneXpert MTB/RIF assay and MTB culture, respectively. The study concluded that GeneXpert assay was found to be a rapid and accurate tool for MTB detection in smear-negative sputum specimens. GeneXpert has advantage over ZN smear microscopy and MTB culture as it detects MTB and rifampicin resistance simultaneously within 2 h with minimal biohazards.

The citrus mealybug, Planococcus citri, is an important plant pest with a very broad plant host range. P. citri is a phloem feeder and loss of plant vigor and stunting are characteristic symptoms induced on a range of host plants, but P. citri also reduces fruit quality and causes fruit drop leading to significant yield reductions. Better strategies for managing this pest are greatly needed. RNA interference (RNAi) is an emerging tool for functional genomics studies and is being investigated as a practical tool for highly targeted insect control. Here we investigated whether RNAi effects can be induced in P. citri and whether candidate mRNAs could be identified as possible targets for RNAi-based P. citri control. RNAi effects were induced in P. citri, as demonstrated by specific target reductions of P. citri actin, chitin synthase 1 and V-ATPase mRNAs after injection of the corresponding specific double-stranded RNA inducers. We also used recombinant Tobacco mosaic virus (TMV) to express these RNAi effectors in Nicotiana benthamiana plants. We found that P. citri showed lower fecundity and pronounced death of crawlers after feeding on recombinant TMV-infected plants. Taken together, our data show that actin, chitin synthase 1 and V-ATPase mRNAs are potential targets for RNAi against P. citri, and that recombinant TMV is an effective tool for evaluating candidate RNAi effectors in plants.

DNA barcodes were obtained for 81 butterfly species belonging to 52 genera from sites in north-central Pakistan to test the utility of barcoding for their identification and to gain a better understanding of regional barcode variation. These species represent 25% of the butterfly fauna of Pakistan and belong to five families, although the Nymphalidae were dominant, comprising 38% of the total specimens. Barcode analysis showed that maximum conspecific divergence was 1.6%, while there was 1.7–14.3% divergence from the nearest neighbour species. Barcode records for 55 species showed <2% sequence divergence to records in the Barcode of Life Data Systems (BOLD), but only 26 of these cases involved specimens from neighbouring India and Central Asia. Analysis revealed that most species showed little incremental sequence variation when specimens from other regions were considered, but a threefold increase was noted in a few cases. There was a clear gap between maximum intraspecific and minimum nearest neighbour distance for all 81 species. Neighbour-joining cluster analysis showed that members of each species formed a monophyletic cluster with strong bootstrap support. The barcode results revealed two provisional species that could not be clearly linked to known taxa, while 24 other species gained their first coverage. Future work should extend the barcode reference library to include all butterfly species from Pakistan as well as neighbouring countries to gain a better understanding of regional variation in barcode sequences in this topographically and climatically complex region.

Although insects dominate the terrestrial fauna, sampling constraints and the poor taxonomic knowledge of many groups have limited assessments of their diversity. Passive sampling techniques and DNA-based species assignments now make it possible to overcome these barriers. For example, Malaise traps collect specimens with minimal intervention while the Barcode Index Number (BIN) system automates taxonomic assignments. The present study employs Malaise traps and DNA barcoding to extend understanding of insect diversity in one of the least known zoogeographic regions, the Saharo-Arabian. Insects were collected at four sites in three countries (Egypt, Pakistan, Saudi Arabia) by deploying Malaise traps. The collected specimens were analyzed by sequencing 658 bp of cytochrome oxidase I (DNA barcode) and assigning BINs on the Barcode of Life Data Systems. The year-long deployment of a Malaise trap in Pakistan and briefer placements at two Egyptian sites and at one in Saudi Arabia collected 53,092 specimens. They belonged to 17 insect orders with Diptera and Hymenoptera dominating the catch. Barcode sequences were recovered from 44,432 (84%) of the specimens, revealing the occurrence of 3,682 BINs belonging to 254 families. Many of these taxa were uncommon as 25% of the families and 50% of the BINs from Pakistan were only present in one sample. Family and BIN counts varied significantly through the year, but diversity indices did not. Although more than 10,000 specimens were analyzed from each nation, just 2% of BINs were shared by Pakistan and Saudi Arabia, 4% by Egypt and Pakistan, and 7% by Egypt and Saudi Arabia. The present study demonstrates how the BIN system can circumvent the barriers imposed by limited access to taxonomic specialists and by the fact that many insect species in the Saharo-Arabian region are undescribed.