Molecular Characterization and Phylogenetic Analysis of Culex quinquefasciatus by DNA Barcoding ()
1. Introduction
Culex quinquefasciatus, a mosquito that belongs to the phylum Arthopoda, is an important factor causing filariasis [1] , St. Louis encephalitis [2] , West Nile virus [3] and Avion malaria [4] . Culex quinquefasciatus, also called the southern house mosquito, is extensively studied as it transmits crucial diseases [5] . Culex mosquitoes are found majorly in the urban cities due to swift urbanization [6] , because of uncanny city growth, which eventually results in the progress of the disease transmission [7] - [9] . Lymphatic filariasis is spread worldwide; affecting 120 million people [10] . They generally bite on the lower limbs [11] and causes skin allergies [12] . Since 1979, there are reports on Japanese encephalitis in India [13] . This is also known as cosmopolitan mosquito, because, it acts as vector for protozoan parasites [14] , filarial worms [15] , and for arboviruses [16] [17] . Culex quinquefasciatus is the prime vector in India causing filariasis [18] ; 91% of them are caused by Wuchereria bancrofti Cobbold [19] . For the human welfare, it is very important to eradicate this species. The genome of mosquito species is generally different in different geometrical areas. Although, the genome sequence of C. quinquefasciatus [20] triggered a new hope for detection, however, the researchers are still involved in identifying and characterizing of these species at molecular level [21] . First, we reported DNA barcodes of Aedes mosquito species [22] . In fact, molecular detection involves identification of species specific sequence. In prokaryotes, particularly in bacteria 16S rRNA based molecular characterization is in practice [23] . Contrastingly, in Molecular nomenclature of the mosquitoes, in mosquito species, the method of DNA barcoding is widely used [24] [25] . The DNA barcoding concept, among all mitochondrial genome genes, uses a 648-base pare region of a mitochondria gene (cytochrome c oxidase I, or COX1) in a study of DNA barcoding of mosquito species.
In this respect, the objective of the present study was focused on the molecular detection of Culex species by DNA barcoding. In addition, from different geographical regions, sequence comparisons and evolutionary relationship among Culex species were studied.
2. Materials and Methods
2.1. Larvae Collection and Characterization
The wild type Culex larva of C. quinquefasciatus was collected from different parts of the Hyderabad city, India, during the breeding season: stagnant water, coconut shells, tires etc. and reared in the Laboratory under specific conditions according (the room temperature maintained 27˚C to 29˚C and 65% room humidity) manual for Mosquito rearing and experimental techniques published by American Mosquito Control Association (AMCA bulletin no. 5, January 1979) nude mice were used for blood meal to mosquitoes to obtain the pure cultures. This pure culture was reared again to get the F1 generation and the 4th instar larvae of F1 generation were used for genomic DNA extraction [26] . The C. quinquefasciatus was identified morphologically with the help of Entomologist deposited at The National History of Museum, Department of Zoology, Osmania University with Accession No (OUNHM.ART.2013.3). Collected larvae were reared and identified morphologically for Culex quinquefasciatus under the microscope. Microscopic analysis revealed the following characteristics, short and stout head, yellow long mouth brushes. It was also observed that the abdomen has eight segments, the saddle and the siphon, which is four times longer that its width, securing the larvae to be Culex quinquefasciatus. The larvae were allowed to grow into an adult. After morphological characterization, further study was performed from the F1 generation of the pure culture.
2.2. DNA Extraction
Genomic DNA was extracted following the procedure of Co.en, and Glover, G. (1982) [27] . From the 4th instar larvae, the total DNA was extracted and was then finely grounded in 50 µl of homogenate buffer. The homogenate was left in the thermo cycler for 30 minutes at 60˚C with a quick addition of 7 µl of 8M Potassium acetate (CH3COOK). Then, incubated the tubes in ice for 30 minutes and centrifuged at maximum speed for 15 minutes; transferred the supernatant into fresh tubes. To precipitate the DNA, incubated the tubes for 15 minutes and after adding 100 µl of 95% ethanol. DNA pellet was collected after rapid centrifugation for 15 minutes at maximum speed. Pellet was air dried and solubilized in tris buffer. DNA was stored at −20˚C till further experimental procedure had been carried.
2.3. COI Partial Sequence Amplification and Sequencing
The isolated DNA was further amplified by Polymerase Chain Reaction (PCR) by mitochondrial Cytochrome c Oxidase I (COI) gene, which can differentiate between the diversity of taxa. The mitochondrial COI gene of 501 bp was amplified by primers, the forward primer (5'- ATT GGA TTA TTA GGA TTT ATT G -3') which was designed to start approximately 76 bp into the CO1 5' region and the reverse primer (5'- GCA GGA GGA AGA GTA TGA TAT C -3'), that are available in literature [22] . The reaction mixture of 25 µl consisting of 10 - 50 ng of DNA template, 200 µM of dNTPs, 1 µ of Taq DNA polymerase, 1x assay buffer and 5 pmol concentration of each primer. PCR conditions: initial denaturation at 94˚C for 2 min, cycle denaturation at 95˚C for 30 sec, primer annealing at 55˚C for 30 sec, extension at 70˚C for 30 sec, and final extension at 70˚C for 10 min. The amplicon was confirmed in Agarose gel electrophoresis. The product was given for commercial sequencing (Bio serve Biotechnologies Pvt. Ltd.).
2.4. Multiple Sequence Alignment
Multiple sequence alignment (MSA) was performed in Clustal W [28] [29] for partial COI gene sequence of 501 bp with the default parameters. Sequence alignment studies elucidated the similarity and differences among different species in India along with other parts of the world.
2.5. Phylogenetic Analysis
The pattern of the evolution of Culex species was studied by phylogenetic analysis. The results of the DNA barcoding were made increasingly vivid with the phylogenetic analysis by the construction of a phylogenetic tree. For phylogeny, the available nucleotide sequences from National Center for Biotechnology Information (NCBI) were collected: Maharashtra, Pondicherry, Japan and UK. In MEGA tool [30] [31] , the analysis of the phylogeny was attained by maximum likelihood method with the deletion of gaps and missing data. Furthermore, Bootstrap replication (500 numbers) was used to validate the tree.
3. Results and Discussions
The Southern House Mosquito C. quinquefasciatus is a well-known vector, which played a prime role in filariasis and the incidences have seen in different parts of the world. Hence, there is an utmost need to address this problem with priority. In order to snub the prognosis of this vector, it is very essential to take certain control measures. To attain this, the mosquito larvae were collected and grown into an adult. The DNA was isolated and amplified and run on the agarose gel electrophoresis. The agarose wells comprised of sample from different geographical locations along with 200 bp marker DNA is show in Figure 1. The 501 bp amplicons were seen between 400 bp and 600 bp of the DNA marker and were then subjected to sequencing. The sequence is shown in Figure 2; it was noted that the sequence was identified to be novel and was submitted to NCBI-gene bank nucleotide database and the accession number assigned was JX088701.
Figure 1. Gel electrophoresis of mitochondrial COI gene isolated from Hyderabad Culex quinquefasciatus mosquito. Lane 1: 200 bp DNA ladder; lane 2-7: PCR products of CO I gene (the samples isolated from different parts of the Hyderabad region). The amplified product of 501bp was corresponded between 400 bp and 600 bp.
Further, Partial mitochondrial COI genes were subjected to multiple sequence alignment on clustal W to speculate the sequence conservations among Culex. species from different geographical locations. From the MSA results, similarities, variations, insertions and deletions were observed. However, Hyderabad Culex species sequence seems more close to UK species as shown in Figure 3 and other species exhibited certain extent of sequence conservation.
Evolutionary pattern of partial mitochondrial COI gene sequence of Hyderabad was checked with Maharashtra, Pondicherry, Japan and UK that were collected from NCBI-Nucleotide database. Phylogenetic analysis reveals that Hyderabad species is more evolutionarily more close to UK species; both the sequences were branched as out group in the tree as shown in Figure 4. All Indian species clustered except one Maharashtra species; this is because of high sequence diversity. In addition to that all Japan species clustered in one group; indicates that they are different from the rest of the species, however, close to Indian species.
Similarly the Evolutionary pattern of partial mitochondrial COI gene sequence of Hyderabad was checked within Indian Culex genus by MSA by bootstrap method in MEGA6 in Figure 5. The critical observation was
Figure 2. Partial mitochondrial sequence of COI gene in NCBI-Nucleotide database with accession number JX088701 shown in FASTA format.
Figure 3. Multiple sequence alignment of partial mitochondrial COI sequences of Culex quinquefasciatus from Hyderabad, Maharashtra, Japan and UK. Sequence alignment was carried out in Clustal W. The sequence of Hyderabad species is more similar to UK species.
Figure 4. Phylogenetic analysis of partial mitochondrial COI sequences of Culex quinquefasciatus from Hyderabad, Maharashtra, Japan and UK. Noticeably, Hyderabad Culex species is closely related to UK specie. And both emerged as out group in the phylogenetic tree. For phylogenetic analysis Maximum likelihood method with Tamura 3 Parameter model and uniform rates among sites were used. Boot strap method with 500 boot strap replications was used to validate the tree.
Figure 5. Phylogenetic analysis of partial mitochondrial COI sequences of Culex quinquefasciatus, from Hyderabad within Indian Culex genera. For phylogenetic analysis Maximum likelihood method with Tamura 3 Parameter model and uniform rates among sites were used. Boot strap method with 500 boot strap replications was used to validate the tree.
found that the Hyderabad species have uniqueness in CO I gene which is evolutionary close to Culexvishnui.
4. Conclusion
The medium sized brown mosquito, Culex quinquefasciatus, predominantly exists all through the tropics acts as a vector causing several parasitic diseases. Usually active at night, it is an opportunistic blood feeder which allows the parasites to use the humans as hosts. It is hence important to address their control with priority. In the present study, partial mitochondrial COI DNA sequence was successfully isolated, sequenced and was submitted to NCBI-Nucleotide database. To check the sequence similarity with other species from different geographical regions, multiple sequence alignment of Culex species was accomplished. Sequence alignment suggesting that, Hyderabad species exhibited decent alignment with UK species. To confirm that further phylogeny was constructed. It was noticed that the DNA sequence of Hyderabad Culex species was different from other species and can be used as DNA barcode to identify the organism.
Acknowledgements
The authors thank the DST-PURSE-OU, DBT-ISLARE and Department Of Zoology for providing infrastructure, RGNF-UGC, New Delhi for financial assistance.
NOTES
*Corresponding author.