Tolulope A. Agunbiade

Mitochondrial genome sequence and expression profiling for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae).

We report the assembly of the 14,054 bp near complete sequencing of the mitochondrial genome of the legume pod borer (LPB), Maruca vitrata (Lepidoptera: Crambidae), which we subsequently used to estimate divergence and relationships within the lepidopteran lineage. The arrangement and orientation of the 13 protein-coding, 2 rRNA, and 19 tRNA genes sequenced was typical of insect mitochondrial DNA sequences described to date. The sequence contained a high A+T content of 80.1% and a bias for the use of codons with A or T nucleotides in the 3rd position. Transcript mapping with midgut and salivary gland ESTs for mitochondrial genome annotation showed that translation from protein-coding genes initiates and terminates at standard mitochondrial codons, except for the coxI gene, which may start from an arginine CGA codon. The genomic copy of coxII terminates at a T nucleotide, and a proposed polyadenylation mechanism for completion of the TAA stop codon was confirmed by comparisons to EST data. EST contig data further showed that mature M. vitrata mitochondrial transcripts are monocistronic, except for bicistronic transcripts for overlapping genes nd4/nd4L and nd6/cytb, and a tricistronic transcript for atp8/atp6/coxIII. This processing of polycistronic mitochondrial transcripts adheres to the tRNA punctuated cleavage mechanism, whereby mature transcripts are cleaved only at intervening tRNA gene sequences. In contrast, the tricistronic atp8/atp6/coxIII in Drosophila is present as separate atp8/atp6 and coxIII transcripts despite the lack of an intervening tRNA. Our results indicate that mitochondrial processing mechanisms vary between arthropod species, and that it is crucial to use transcriptional information to obtain full annotation of mitochondrial genomes.

Comparison of the transcriptional profiles of head and body lice

Head and body lice are both blood‐feeding parasites of humans although only the body louse is a potent disease vector. In spite of numerous morphological and life history differences, head and body lice have recently been hypothesized to be ecotypes of the same species. We took a comparative genomics approach to measure nucleotide diversity by comparing expressed sequence tag data sets from head and body lice. A total of 10 771 body louse and 10 770 head louse transcripts were predicted from a combined assembly of Roche 454 and Illumina sequenced cDNAs from whole body tissues collected at all life stages and during pesticide exposure and bacterial infection treatments. Illumina reads mapped to the 10 775 draft body louse gene models from the whole genome assembly predicted nine presence/absence differences, but PCR confirmation resulted in a single gene difference. Read per million base pair estimates indicated that 14 genes showed significant differential expression between head and body lice under our treatment conditions. One novel microRNA was predicted in both lice species and 99% of the 544 transcripts from Candidatus riesia indicate that they share the same endosymbiont. Overall, few differences exist, which supports the hypothesis that these two organisms are ecotypes of the same species.

Transcriptome Sequencing, and Rapid Development and Application of SNP Markers for the Legume Pod Borer Maruca vitrata (Lepidoptera: Crambidae)

The legume pod borer, Maruca vitrata (Lepidoptera: Crambidae), is an insect pest species of crops grown by subsistence farmers in tropical regions of Africa. We present the de novo assembly of 3729 contigs from 454- and Sanger-derived sequencing reads for midgut, salivary, and whole adult tissues of this non-model species. Functional annotation predicted that 1320 M. vitrata protein coding genes are present, of which 631 have orthologs within the Bombyx mori gene model. A homology-based analysis assigned M. vitrata genes into a group of paralogs, but these were subsequently partitioned into putative orthologs following phylogenetic analyses. Following sequence quality filtering, a total of 1542 putative single nucleotide polymorphisms (SNPs) were predicted within M. vitrata contig assemblies. Seventy one of 1078 designed molecular genetic markers were used to screen M. vitrata samples from five collection sites in West Africa. Population substructure may be present with significant implications in the insect resistance management recommendations pertaining to the release of biological control agents or transgenic cowpea that express Bacillus thuringiensis crystal toxins. Mutation data derived from transcriptome sequencing is an expeditious and economical source for genetic markers that allow evaluation of ecological differentiation.

Development of Reference Transcriptomes for the Major Field Insect Pests of Cowpea: A Toolbox for Insect Pest Management Approaches in West Africa

Cowpea is a widely cultivated and major nutritional source of protein for many people that live in West Africa. Annual yields and longevity of grain storage is greatly reduced by feeding damage caused by a complex of insect pests that include the pod sucking bugs, Anoplocnemis curvipes Fabricius (Hemiptera: Coreidae) and Clavigralla tomentosicollis Stål (Hemiptera: Coreidae); as well as phloem-feeding cowpea aphids, Aphis craccivora Koch (Hemiptera: Aphididae) and flower thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae). Efforts to control these pests remain a challenge and there is a need to understand the structure and movement of these pest populations in order to facilitate the development of integrated pest management strategies (IPM). Molecular tools have the potential to help facilitate a better understanding of pest populations. Towards this goal, we used 454 pyrosequencing technology to generate 319,126, 176,262, 320,722 and 227,882 raw reads from A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. The reads were de novo assembled into 11,687, 7,647, 10,652 and 7,348 transcripts for A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. Functional annotation of the resulting transcripts identified genes putatively involved in insecticide resistance, pathogen defense and immunity. Additionally, sequences that matched the primary aphid endosymbiont, Buchnera aphidicola, were identified among A. craccivora transcripts. Furthermore, 742, 97, 607 and 180 single nucleotide polymorphisms (SNPs) were respectively predicted among A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti transcripts, and will likely be valuable tools for future molecular genetic marker development. These results demonstrate that Roche 454-based transcriptome sequencing could be useful for the development of genomic resources for cowpea pest insects in West Africa.

Genetic differentiation among Maruca vitrata F. (Lepidoptera: Crambidae) populations on cultivated cowpea and wild host plants: implications for insect resistance management and biological control strategies.

Maruca vitrata Fabricius (Lepidoptera: Crambidae) is a polyphagous insect pest that feeds on a variety of leguminous plants in the tropics and subtropics. The contribution of host-associated genetic variation on population structure was investigated using analysis of mitochondrial cytochrome oxidase 1 (cox1) sequence and microsatellite marker data from M. vitrata collected from cultivated cowpea (Vigna unguiculata L. Walp.), and alternative host plants Pueraria phaseoloides (Roxb.) Benth. var. javanica (Benth.) Baker, Loncocarpus sericeus (Poir), and Tephrosia candida (Roxb.). Analyses of microsatellite data revealed a significant global FST estimate of 0.05 (P≤0.001). The program STRUCTURE estimated 2 genotypic clusters (co-ancestries) on the four host plants across 3 geographic locations, but little geographic variation was predicted among genotypes from different geographic locations using analysis of molecular variance (AMOVA; among group variation −0.68%) or F-statistics (F ST Loc = −0.01; P = 0.62). These results were corroborated by mitochondrial haplotype data (φSTLoc = 0.05; P = 0.92). In contrast, genotypes obtained from different host plants showed low but significant levels of genetic variation (F ST Host = 0.04; P = 0.01), which accounted for 4.08% of the total genetic variation, but was not congruent with mitochondrial haplotype analyses (φSTHost = 0.06; P = 0.27). Variation among host plants at a location and host plants among locations showed no consistent evidence for M. vitrata population subdivision. These results suggest that host plants do not significantly influence the genetic structure of M. vitrata, and this has implications for biocontrol agent releases as well as insecticide resistance management (IRM) for M. vitrata in West Africa.

The spatial genetic differentiation of the legume pod borer, Maruca vitrata F. (Lepidoptera: Crambidae) populations in West Africa

The legume pod borer, Maruca vitrata, is an endemic insect pest that causes significant yield loss to the cowpea crop in West Africa. The application of population genetic tools is important in the management of insect pests but such data on M. vitrata is lacking. We applied a set of six microsatellite markers to assess the population structure of M. vitrata collected at five sites from Burkina Faso, Niger and Nigeria. Observed polymorphisms ranged from one (marker 3393) to eight (marker 32008) alleles per locus. Observed and expected heterozygosities ranged from 0.0 to 0.8 and 0.0 to 0.6, respectively. Three of the loci in samples from Nigeria and Burkina Faso deviated significantly from Hardy-Weinberg Equilibrium (HWE), whereas no loci deviated significantly in samples from Niger. Analysis of molecular variance (AMOVA) indicated that 67.3% level of the genetic variation was within individuals compared to 17.3% among populations. A global estimate of F ST=0.1 (ENA corrected F ST=0.1) was significant (P⩽0.05) and corroborated by pairwise F ST values that were significant among all possible comparisons. A significant correlation was predicted between genetic divergence and geographic distance between subpopulations (R2=0.6, P=0.04), and cluster analysis by the program STRUCTURE predicted that co-ancestry of genotypes were indicative of three distinct populations. The spatial genetic variance among M. vitrata in West Africa may be due to limited gene flow, south-north seasonal movement pattern or other reproductive barriers. This information is important for the cultural, chemical and biological control strategies for managing M. vitrata.

Potential use of Sesbania pachycarpa (Fabace ae: Papilionoideae) as a refugia for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

Deployment of cowpea with Bacillus thuringiensis (Bt-cowpea), to control Maruca vitrata Fab., must be preceded by the development of an insect resistance management (IRM) plan to ensure a sustainable use of the in-plant protection offered by the transgenic variety. One of the components of a resistance management plan involves the use of wild or cultivated host plants as refugia. In West Africa, wild refugia have the potential to be a major component of such an IRM strategy. In the current study, we examined the occurrence of M. vitrata on three cultivated cowpea varieties and one wild alternative host, Sesbania pachycarpa D.C. Our results indicate that M. vitrata population overlapped on the wild host and cowpea. The overall insect population on S. pachycarpa represents 9–13% of the population on cultivated cowpea. Based on these findings, we suggest that S. pachycarpa may contribute as a refuge for M. vitrata population in case of deployment of Bt-cowpea.

Effect of fungi–termite interaction on the chlorophyll content of three rice varieties grown on ultisol soil of Ikenne, south-west Nigeria

The effect of fungi–termite interaction on three rice varieties was conducted in a screen house at the Africa Rice Center (AfricaRice) Ibadan, Nigeria. Of the 10 fungi species (Fusarium verticilloides, Trichoderma sp., Aspergillus niger, Aspergillus flavus, Macrophoma sp., Neurospora sp., Botryodiplodia theobromae, Penicillum sp., Rhizopus sp. and Sclerotium rolfsii) isolated from termites, soil and rice plants, F. verticilloides, Trichoderma sp. and B. theobromae were used for the interaction study. Each fungus was inoculated singly and in combinations with termite into the root of each rice variety in potted soil. Leaf samples were taken to measure the chlorophyll content which is a major parameter to estimate effect of termite–fungi interaction. The chlorophyll content of the inoculated rice plants was significantly reduced when compared with the control. Fungi interaction with termite had significant reduction on the chlorophyll content. The synergistic relationship between the fungi and the termite was discussed.

Comparison of the mitochondrial genomes of the Old and New World strains of the legume pod borer, Maruca vitrata (Lepidoptera: Crambidae)

Maruca vitrata (Fabricius, 1787) is a cryptic pantropical species of Lepidoptera that are comprised of two unique strains that inhabit the American continents (New World strain) and regions spanning from Africa through to Southeast Asia and Northern Australia (Old World strain). In this study, we de novo assembled the complete mitochondrial genome sequence of the New World legume pod borer, M. vitrata, from shotgun sequence data generated on an Illumina HiSeq 2000. Phylogenomic comparisons were made with other previously published mitochondrial genome sequences from crambid moths, including the Old World strain of M. vitrata. The 15,385 bp M. vitrata (New World) sequence has an 80.7% A+T content and encodes the 13 protein-coding, 2 ribosomal RNA and 22 transfer RNA genes in the typical orientation and arrangement of lepidopteran mitochondrial DNAs. Mitochondrial genome-wide comparison between the New and Old World strains of M. vitrata detected 476 polymorphic sites (4.23% nucleotide divergence) with an excess of synonymous substitution as a result of purifying selection. Furthermore, this level of sequence variation suggests that these strains diverged from ~1.83 to 2.12 million years ago, assuming a linear rate of short-term substitution. The de novo assemblies of mitochondrial genomes from next-generation sequencing (NGS) reads provide readily available data for similar comparative studies.