Brad S. Coates

Two differentially expressed ommochrome-binding protein-like genes (obp1 and obp2) in larval fat body of the European corn borer, Ostrinia nubilalis

Abstract Ommochrome-binding proteins function in coloration and detoxification pathways by transporting tryptophan metabolites, and increase in hemolymph concentration prior to diapause. Two ommochrome-binding protein genes from the European corn borer Ostrinia nubilalis (Hübner) (Onobp1 and Onobp2; GenBank accession nos. AY819651 to AY819655 and AY862870) were isolated. Relatedness to OBP-encoding genes was suggested by peptide similarity, phylogenetic reconstruction, and expression data. 21 single nucleotide polymorphisms between obp1 and 23 polymorphisms between obp2 alleles were identified, and resultant genomic markers were inherited in a Mendelian fashion. RT-PCR showed fat body specific Onobp1 and Onobp2 transcription. The Onobp1 transcript was RT-PCR amplified from fat body of 5th instars, whereas Onobp2 was expressed in fat body of 4th and 5th instars, and peaked in 5th instar wandering and 1 week old diapausing larvae. Expression suggests gene duplicates are maintained by change in temporal expression. The significance of Onobp1 and 2 gene products to O. nubilalis diapause physiology requires additional investigation. DAP diapause-associated polypeptide OBP ommochrome binding protein Onobp1 Ostrinia nubilalis ommochrome binding protein 1 Onobp2 Ostrinia nubilalis ommochrome binding protein 2

Distribution of Genes and Repetitive Elements in the Diabrotica virgifera virgifera Genome Estimated Using BAC Sequencing

Feeding damage caused by the western corn rootworm, Diabrotica virgifera virgifera, is destructive to corn plants in North America and Europe where control remains challenging due to evolution of resistance to chemical and transgenic toxins. A BAC library, DvvBAC1, containing 109,486 clones with 104 ± 34.5 kb inserts was created, which has an ~4.56X genome coverage based upon a 2.58 Gb (2.80 pg) flow cytometry-estimated haploid genome size. Paired end sequencing of 1037 BAC inserts produced 1.17 Mb of data (~0.05% genome coverage) and indicated ~9.4 and 16.0% of reads encode, respectively, endogenous genes and transposable elements (TEs). Sequencing genes within BAC full inserts demonstrated that TE densities are high within intergenic and intron regions and contribute to the increased gene size. Comparison of homologous genome regions cloned within different BAC clones indicated that TE movement may cause haplotype variation within the inbred strain. The data presented here indicate that the D. virgifera virgifera genome is large in size and contains a high proportion of repetitive sequence. These BAC sequencing methods that are applicable for characterization of genomes prior to sequencing may likely be valuable resources for genome annotation as well as scaffolding.

Two sex-chromosome-linked microsatellite loci show geographic variance among North American Ostrinia nubilalis

Abstract PCR-based O. nubilalis population and pedigree analysis indicated female specificity of a (GAAAAT)n microsatellite, and male specificity of a CAYCARCGTCACTAA repeat unit marker. These loci were respectively named Ostrinia nubilalis W-chromosome 1 (ONW1) and O. nubilalis Z-chromosome 1 (ONZ1). Intact repeats of three, four, or five GAAAAT units are present among ONW1 alleles, and biallelic variation exists at the ONZ1 locus. Screening of 493 male at ONZ1 and 448 heterogametic females at ONZ1 and ONW1 loci from eleven North American sample sites was used to construct genotypic data. Analysis of molecular variance (AMOVA) and F-statistics indicated no female haplotype or male ONZ1 allele frequency differentiation between voltinism ecotypes. Four subpopulations from northern latitudes, Minnesota and South Dakota, showed the absence of a single female haplotype, a significant deviation of ONZ1 data from Hardy-Weinberg expectation, and low-level geographic divergence from other subpopulations. Low ONZ1 and ONW1 allele diversity could be attributed either to large repeat unit sizes, low repeat number, reduced effective population (Ne) size of sex chromosomes, or the result of recent O. nubilalis introduction and population expansion, but likely could not be due to inbreeding. These sequences have been deposited in GenBank AF442958, and AY102618 to AY102620. ONW1   Ostrinia nubilalis W-chromosome marker number 1 ONZ1   Ostrinia nubilalis Z-chromosome marker number 1

Comparative Performance of Single Nucleotide Polymorphism and Microsatellite Markers for Population Genetic Analysis

Microsatellite loci are standard genetic markers for population genetic analysis, whereas single nucleotide polymorphisms (SNPs) are more recent tools that require assessment of neutrality and appropriate use in population genetics. Twelve SNP markers were used to describe the genetic structure of Diabrotica virgifera virgifera (LeConte; Coleoptera: Chrysomelidae) in the United States of America and revealed a high mean observed heterozygosity (0.40 +/- 0.059) and low global F(ST) (0.029). Pairwise F(ST) estimates ranged from 0.007 to 0.045, and all but 2 populations showed significant levels of genetic differentiation (P < or = 0.008). Population parameters and conclusions based on SNP markers were analogous to that obtained by use of microsatellite markers from the identical population samples. SNP-based F(ST) estimates were 3-fold higher than corresponding estimates from microsatellites, wherein lower microsatellite F(ST) estimates likely resulted from an overestimate of migration rates between subpopulations due to convergence of allele size (homoplasy). No significant difference was observed in the proportion of SNP or microsatellite markers loci that were nonneutral within populations. SNP markers provided estimates of population genetic parameters consistent with those from microsatellite data, and their low back mutation rates may result in reduced propensity for error in estimation of population parameters.

Genome sequencing of the sweetpotato whitefly Bemisia tabaci MED/Q

Abstract The sweetpotato whitefly Bemisia tabaci is a highly destructive agricultural and ornamental crop pest. It damages host plants through both phloem feeding and vectoring plant pathogens. Introductions of B. tabaci are difficult to quarantine and eradicate because of its high reproductive rates, broad host plant range, and insecticide resistance. A total of 791 Gb of raw DNA sequence from whole genome shotgun sequencing, and 13 BAC pooling libraries were generated by Illumina sequencing using different combinations of mate-pair and pair-end libraries. Assembly gave a final genome with a scaffold N50 of 437 kb, and a total length of 658 Mb. Annotation of repetitive elements and coding regions resulted in 265.0 Mb TEs (40.3%) and 20 786 protein-coding genes with putative gene family expansions, respectively. Phylogenetic analysis based on orthologs across 14 arthropod taxa suggested that MED/Q is clustered into a hemipteran clade containing A. pisum and is a sister lineage to a clade containing both R. prolixus and N. lugens. Genome completeness, as estimated using the CEGMA and Benchmarking Universal Single-Copy Orthologs pipelines, reached 96% and 79%. These MED/Q genomic resources lay a foundation for future ‘pan-genomic’ comparisons of invasive vs. noninvasive, invasive vs. invasive, and native vs. exotic Bemisia, which, in return, will open up new avenues of investigation into whitefly biology, evolution, and management.

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.

Allelic variation of a Beauveria bassiana (Ascomycota: Hypocreales) minisatellite is independent of host range and geographic origin

The minisatellite locus, BbMin1, was isolated from a partial Beauveria bassiana genomic library that consisted of poly(GA) flanked inserts. Polymerase chain reaction (PCR) of the BbMin1 repeat demonstrated allele size variation among 95 B. bassiana isolates. Amplification was also observed from single isolates of Beauveria amorpha, Beauveria brongniartii, and Beauveria caledonica. Eight alleles were identified at the haploid locus, where repeat number fluctuated between one and fourteen. AMOVA and theta (Fst) indicated that fixation of repeat number has not occurred within pathogenic ecotypes or geographically isolated samples of B. bassiana. Selective neutrality of allele size, the rate of BbMin1 mutation, and the age of the species may contribute to host and geographic independence of the marker. Presence of alleles with a large number of repeat units may be attributed to the rare occurrence of somatic recombination or DNA replication error. The molecular genetic marker was useful for the identification of genetic types of B. bassiana and related species.

Mining an Ostrinia nubilalis midgut expressed sequence tag (EST) library for candidate genes and single nucleotide polymorphisms (SNPs)

Genes expressed in lepidopteran midgut tissues are involved in digestion and Bacillus thuringiensis (Bt) toxin resistance traits. Five hundred and thirty five unique transcripts were annotated from 1745 high quality O. nubilalis larval midgut expressed sequence tags (ESTs). Full‐length cDNA sequence of 12 putative serine proteinase genes and 3 partial O. nubilalis aminopeptidase N protein genes, apn1, apn3, and apn4, were obtained, and genes may have roles in plant feeding and Bt toxin resistance traits of Ostrinia larvae. The EST library was not normalized and insert frequencies reflect transcript levels under the initial treatment conditions and redundancy of inserts from highly expressed transcripts allowed prediction of putative single nucleotide polymorphisms (SNPs). Ten di‐, tri‐ or tetranucleotide repeat unit microsatellite loci were identified, and minisatellite repeats were observed within the C‐termini of two encoded serine proteinases. Molecular markers showed polymorphism at 28 SNP loci and one microsatellite locus, and Mendelian inheritance indicated that markers were applicable to genome mapping applications. This O. nubilalis larval midgut EST collection is a resource for gene discovery, expression information, and allelic variation for use in genetic marker development.

Beauveria bassiana haplotype determination based on nuclear rDNA internal transcribed spacer PCR±RFLP

DNA sequence alignment of the nuclear 5.8S ribosomal RNA (rRNA) gene and internal transcribed spacers (ITS) from Beauveria bassiana demonstrated that 6.62% sequence variation existed between nine isolates. A higher level of mutation was observed within the ITS regions, where 8.39% divergence occurred. Polymerase chain reaction restriction fragment length polymorphism, PCR-RFLP, and DNA sequence alignment of the ITS1 and ITS2 regions identified seven polymorphic restriction endonuclease sites, Alu I, Hha I, Hinf I, Sin I, Tru 9AI and two Tha I restriction sites. The allelic frequency of each genetic marker was determined from 96 isolates. PCR-RFLP defined 24 B. bassiana genotypes within the sample set, from which eight phylogenetic clusters were predicted to exist. AMOVA and Fst (θ) indicated that no significant correlation existed between B. bassiana haplotype and insect host range as defined by insect order from which each isolate was derived.

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.