Landscape genetics of a pollinator longhorn beetle [Typocerus v. velutinus (Olivier)] on a continuous habitat surface.

Landscape connectivity, the degree to which the landscape structure facilitates or impedes organismal movement and gene flow, is increasingly important to conservationists and land managers. Metrics for describing the undulating shape of continuous habitat surfaces can expand the usefulness of continuous gradient surfaces that describe habitat and predict the flow of organisms and genes. We adopted a landscape gradient model of habitat and used surface metrics of connectivity to model the genetic continuity between populations of the banded longhorn beetle [Typocerus v. velutinus (Olivier)] collected at 17 sites across a fragmentation gradient in Indiana, USA. We tested the hypothesis that greater habitat connectivity facilitates gene flow between beetle populations against a null model of isolation by distance (IBD). We used next-generation sequencing to develop 10 polymorphic microsatellite loci and genotype the individual beetles to assess the population genetic structure. Isolation by distance did not explain the population genetic structure. The surface metrics model of habitat connectivity explained the variance in genetic dissimilarities 30 times better than the IBD model. We conclude that surface metrology of habitat maps is a powerful extension of landscape genetics in heterogeneous landscapes.

Differential expression of candidate salivary effector proteins in field collections of Hessian fly, Mayetiola destructor.

Evidence is emerging that some proteins secreted by gall-forming parasites of plants act as effectors responsible for systemic changes in the host plant, such as galling and nutrient tissue formation. A large number of secreted salivary gland proteins (SSGPs) that are the putative effectors responsible for the physiological changes elicited in susceptible seedling wheat by Hessian fly, Mayetiola destructor (Say), larvae have been documented. However, how the genes encoding these candidate effectors might respond under field conditions is unknown. The goal of this study was to use microarray analysis to investigate variation in SSGP transcript abundance amongst field collections from different geographical regions (southeastern USA, central USA, and the Middle East). Results revealed significant variation in SSGP transcript abundance amongst the field collections studied. The field collections separated into three distinct groups that corresponded to the wheat classes grown in the different geographical regions as well as to recently described Hessian fly populations. These data support previous reports correlating Hessian fly population structure with micropopulation differences owing to agro-ecosystem parameters such as cultivation of regionally adapted wheat varieties, deployment of resistance genes and variation in climatic conditions.

Biological and molecular characterization of Hessian fly (Diptera: Cecidomyiidae) from Israel.

Samples of a dipteran pest of wheat were tested to confirm identity, describe local populations and suggest the use of deploying resistance (R) genes in wheat cultivars for control of Mayetiola destructor, Hessian fly (HF). Morphological evaluation of adults and a free-choice oviposition preference test documenting that females overwhelmingly preferred to oviposit on wheat instead of barley supported they were HF. Using the cytochrome c oxidase subunit I (coxI), the Barcoding Region, nine haplotypes were revealed. Two were found only in the Israeli collections and averaged 3% sequence divergence compared to the other seven haplotypes found in the United States, Israel and Syria. In evaluations of virulence, the Israeli HF in culture was virulent to 11 of the 19 (R) genes tested, and complementation analysis documented that, for four of the R genes tested, the Israeli HF shared loci for virulence with HF from the United States. Levels of HF infestation at seven Israeli fields were at least at the 5-8% level, which historically has indicated a significant yield loss. Microsatellite genotyping of the five HF collections from Israel revealed mixed populations in Israel that are distinctly separate from the single population in Syria.

Bowman-Birk inhibitor affects pathways associated with energy metabolism in Drosophila melanogaster.

Bowman-Birk inhibitor (BBI) is toxic when fed to certain insects, including the fruit fly, Drosophila melanogaster. Dietary BBI has been demonstrated to slow growth and increase insect mortality by inhibiting the digestive enzymes trypsin and chymotrypsin, resulting in a reduced supply of amino acids. In mammals, BBI influences cellular energy metabolism. Therefore, we tested the hypothesis that dietary BBI affects energy-associated pathways in the D. melanogaster midgut. Through microarray and metabolomic analyses, we show that dietary BBI affects energy utilization pathways in the midgut cells of D. melanogaster. In addition, ultrastructure studies indicate that microvilli are significantly shortened in BBI-fed larvae. These data provide further insights into the complex cellular response of insects to dietary protease inhibitors.

Transcriptional signatures in response to wheat germ agglutinin and starvation in Drosophila melanogaster larval midgut.

One function of plant lectins such as wheat germ agglutinin is to serve as defences against herbivorous insects. The midgut is one critical site affected by dietary lectins. We observed marked cellular, structural and gene expression changes in the midguts of Drosophila melanogaster third instar larvae that were fed wheat germ agglutinin. Some of these changes were similar to those observed in the midguts of starved D. melanogaster. Dietary wheat germ agglutinin caused shortening, branching, swelling, distortion and in some cases disintegration of the midgut microvilli. Starvation was accompanied primarily by shortening of the microvilli. Microarray analyses revealed that dietary wheat germ agglutinin evoked differential expression of 61 transcripts; seven of these were also differentially expressed in starved D. melanogaster. The differentially transcribed gene clusters in wheat germ agglutinin-fed larvae were associated with (1) cytoskeleton organization; (2) digestive enzymes; (3) detoxification reactions; and (4) energy metabolism. Four possible transcription factor binding motifs were associated with the differentially expressed genes. One of these exhibited substantial similarity to MyoD, a transcription factor binding motif associated with cellular structures in mammals. These results are consistent with the hypothesis that wheat germ agglutinin caused a starvation-like effect and structural changes of midgut cells of D. melanogaster third-instar larvae.

Transcriptomic profiles of Drosophila melanogaster third instar larval midgut and responses to oxidative stress.

Oligoarray analysis was used to determine the number and nature of genes expressed in third instar Drosophila melanogaster larval midguts. The majority of transcripts were associated with protein synthesis and metabolism. Serine proteases were the main proteolytic enzymes detected. Some 40% of the cytochrome P450 genes and 74% of the glutathione S transferases (GSTs) in the genome of D. melanogaster were observed to be expressed in the midgut by oligoarray analysis. We also identified potential transcription factor binding motifs (TFBMs) of P450s, GSTs and carboxylesterases. Many of the midgut-expressed GST genes contained candidate TFBMs homologous to TFBMs in mammals that have been associated with responses to oxidative stress. We also investigated the response of GSTs in the midgut to dietary H2O2, which showed a dosage-based differential response.

Genome-wide analysis of phenobarbital-inducible genes in Drosophila melanogaster.

An oligoarray analysis was conducted to determine the differential expression of genes due to phenobarbital exposure in Drosophila melanogaster (w(1118) strain) third instar larvae. Seventeen genes were observed to be induced with increased expression by a statistical analysis of microarrays approach with a q < or = 0.05. At q < or = 0.12, four more genes (Cyp12d1, DmGstd4, and two genes with unknown function) were found to be up-regulated, and 11 genes with unknown function were found to be down-regulated. Fifteen of these genes, Cyp4d14, Cyp6a2, Cyp6a8, Cyp12d1, Cyp6d5, Cyp6w1, CG2065, DmGstd6, DmGstd7, Amy-p/Amy-d, Ugt86Dd, GC5724, Jheh1, Jheh2 and CG11893, were verified using quantitative real time polymerase chain reaction. Some of these genes have been shown to be over-transcribed in metabolically DDT-resistant Drosophila strains.

Rangewide population genetic structure of the African malaria vector Anopheles funestus.

Anopheles funestus is a primary vector of malaria in Africa south of the Sahara. We assessed its rangewide population genetic structure based on samples from 11 countries, using 10 physically mapped microsatellite loci, two per autosome arm and the X (N = 548), and 834 bp of the mitochondrial ND5 gene (N = 470). On the basis of microsatellite allele frequencies, we found three subdivisions: eastern (coastal Tanzania, Malawi, Mozambique and Madagascar), western (Burkina Faso, Mali, Nigeria and western Kenya), and central (Gabon, coastal Angola). A. funestus from the southwest of Uganda had affinities to all three subdivisions. Mitochondrial DNA (mtDNA) corroborated this structure, although mtDNA gene trees showed less resolution. The eastern subdivision had significantly lower diversity, similar to the pattern found in the codistributed malaria vector Anopheles gambiae. This suggests that both species have responded to common geographic and/or climatic constraints. The western division showed signatures of population expansion encompassing Kenya west of the Rift Valley through Burkina Faso and Mali. This pattern also bears similarity to A. gambiae, and may reflect a common response to expanding human populations following the development of agriculture. Due to the presumed recent population expansion, the correlation between genetic and geographic distance was weak. Mitochondrial DNA revealed further cryptic subdivision in A. funestus, not detected in the nuclear genome. Mozambique and Madagascar samples contained two mtDNA lineages, designated clade I and clade II, that were separated by two fixed differences and an average of 2% divergence, which implies that they have evolved independently for approximately 1 million years. Clade I was found in all 11 locations, whereas clade II was sampled only on Madagascar and Mozambique. We suggest that the latter clade may represent mtDNA capture by A. funestus, resulting from historical gene flow either among previously isolated and divergent populations or with a related species.

Molecular differentiation between chromosomally defined incipient species of Anopheles funestus.

Anopheles funestus Giles is one of the most important vectors of malaria in sub-Saharan Africa. The population structure of this mosquito in Burkina Faso, West Africa based on chromosomal inversion data led to the description of two chromosomal forms, Kiribina and Folonzo. Because both forms co-occur in the same locales yet differ significantly, both in the frequency of inverted arrangements on chromosome arms 3R and 2R and in vectorial capacity, they were hypothesized to be emerging species with at least partial barriers to gene flow. This hypothesis would be strengthened by molecular evidence of differentiation between Kiribina and Folonzo at loci outside chromosomal inversions. We surveyed molecular variation in sympatric populations of the two forms using sequences from the mitochondrial ND5 gene and genotypes at sixteen microsatellite loci distributed across the genome. Both classes of marker revealed slight but significant differentiation between the two forms (mtDNA F(ST) = 0.023, P < 0.001; microsatellite F(ST) = 0.004, P < 0.001; R(st) = 0.009, P = 0.002). Locus-by-locus analysis of the microsatellite data showed that significant differentiation was not genome-wide, but could be attributed to five loci on chromosome 3R (F(ST) = 0.010, P < 0.001; R(st) = 0.016, P = 0.002). Importantly, three of these loci are outside of, and in linkage equilibrium with, chromosomal inversions, suggesting that differentiation between chromosomal forms extends beyond the inversions themselves. The slight overall degree of differentiation indicated by both marker classes is likely an underestimate because of recent population expansion inferred for both Folonzo and Kiribina. The molecular evidence from this study is consistent with the hypothesis of incipient speciation between Kiribina and Folonzo.