The genetic structure of Asian corn borer, Ostrinia furnacalis, populations in China: haplotype variance in northern populations and potential impact on management of resistance to transgenic maize

Multilocus evidence provides insight into the demographic history and asymmetrical gene flow between Ostrinia furnacalis and Ostrinia nubilalis (Lepidoptera: Crambidae) in the Yili area, Xinjiang, China

Tianshan Mountains provide a model for studying biological evolution and speciation. Here we assess the evolutionary history of Ostrinia furnacalis (ACB) and Ostrinia nubilalis (ECB), which are sympatric in the Yili River Valley in Xinjiang, China.Our study is based on the historical gene flow analyses of two species by using three mitochondrial DNA (mtDNA, COI, COII, Cytb) and four nuclear DNA (nuDNA, EF-1α, Wingless, RPS5, CAD) markers obtained from representatives of HC (Huocheng), YN (Yining), XY (Xinyuan), and MNS (Manasi).Our results reveal that there is an asymmetrical gene flow pattern between the four populations. The population migratory pathways between these different populations show inflow into HC and YN, outflow from XY, and that MNS maintained a flow balance. Bayesian divergence time dating based on the COI gene suggests that the genetic divergence between the two species in this area may have occurred in Holocene at 0.008 Mya. Neutrality tests (Tajima's D, Fu's F s), and mismatch distribution test results suggest that population expansion events may not have occurred in the recent past. The demographic history and gene flow pattern between ACB and ECB may follow the "mountain isolation" hypothesis. The ML and BI trees of the mtDNA haplotype dataset show that ECB haplotypes are grouped together in a distinct clade and are clearly separate from ACB haplotypes. However, the geographical pattern of haplotype distribution is less clear for both ACB and ECB, supporting that there has been frequent gene flow among the geographic populations in the Tianshan Mountains.These findings indicate that the Tianshan Mountains are less likely a barrier to gene flow of the two species.

Knockout of ABC Transporter ABCG4 Gene Confers Resistance to Cry1 Proteins in Ostrinia furnacalis

Ostrinia furnacalis is an important borer on maize. Long-term and large-scale planting of transgenic corn has led O. furnacalis evolving resistance and reducing the control effect. Recently, high levels of resistance to Bt Cry1 toxins have been reported to be genetically linked to the mutation or down-regulation of ABC transporter subfamily G gene ABCG4 in O. furnacalis. In order to further determine the relationship between ABCG4 gene and the resistance to Cry1 toxins in O. furnacalis, the novel CRISPR/Cas9 genome engineering system was utilized to successfully construct ABCG4-KO knockout homozygous strain. Bioassay results indicated that an ABCG4-KO strain had a higher resistance to Cry1 proteins compared with a susceptible strain (ACB-BtS). The result indicates that the ABCG4 gene may act as a receptor of the Bt Cry1 toxin in O. furnacalis. Furthermore, the development time was significantly changed in the early stage ABCG4-KO larvae, and the population parameters were also significantly changed. In summary, our CRISPR/Cas9-mediated genome editing study presents evidence that ABCG4 gene is a functional receptor for Bt Cry1 toxins, laying the foundation for further clarification of the Bt resistance mechanism.

Landscape genetic analyses reveal host association of mitochondrial haplotypes in the Asian corn borer, Ostrinia furnacalis

Crop expansion often leads to high pest pressure. These pests may have fitness trade-offs related to host use, and some host-associated genotypes may benefit and increase in frequency. However, evidence concerning the effect of host availability on spatial distribution and frequency of mitochondrial haplotypes is scarce. We studied genetic variation of the Asian corn borer, Ostrinia furnacalis (Guenée), across a large area during 2 years (2016 and 2017). Mitochondrial sequence data were obtained from 530 individuals collected from 79 locations in Shandong Province, China. In total, 155 haplotypes were found based on the combined cytochrome oxidase subunit 1 (COI) and COII genes. Three haplotypes (H2, H12, and H23) were dominant, whereas most of the other haplotypes occurred in low frequency. A haplotype network showed that the 155 haplotypes can be grouped into three clusters. Haplotype clusters seemed to be randomly distributed. The frequency of H12 (in Cluster 1) was positively correlated with maize crop proportion, but negatively correlated with other crops (primarily vegetables, oilseed crops, and cotton) at all spatial scales (1-, 3-, and 5-km radius). Cluster 2 had haplotype H23, and this cluster was negatively correlated with semi-natural habitats. Cluster 3 had no dominant haplotype and was not affected by landscape factors. We conclude that H12 may be a maize-associated haplotype. Further study is needed to verify the possibility that the carriers of this haplotype may possess some fitness trade-offs. Our study highlights the importance of host availability in O. furnacalis haplotype distribution and frequency.

Effects of NPF on larval taste responses and feeding behaviors in Ostrinia furnacalis

The insect taste system regulates insect feeding behavior and patterns of food consumption. In this study, we showed that the medial and lateral sensilla styloconica in the mouthparts of 5th-instar Asian corn borer larvae are sensitive to fructose and sucrose in a concentration-dependent way. The two sensilla produced significant electrophysiological responses (greater than100 spikes/s) by exposure to 10 mM fructose or sucrose. However, electrophysiological responses and feeding preferences to fructose or sucrose were inhibited by neuropeptide F double-stranded RNA (dsNPF). Additionally, the medial sensilla styloconica are sensitive to low concentrations of the deterrents caffeine and nicotine. However, starvation, followed by increases in larval npf expression plus feeding, led to increases in spike frequencies of related sensilla to fructose, sucrose, and deterrents. In contrast, these responses were reduced on the dsNPF treatment. Our results suggest that NPF plays an important role influencing caterpillar feeding behavior through regulating the taste neurons of the sensilla styloconica.

Influence of voltine ecotype and geographic distance on genetic and haplotype variation in the Asian corn borer

Diapause is an adaptive dormancy strategy by which arthropods endure extended periods of adverse climatic conditions. Seasonal variation in larval diapause initiation and duration in Ostrinia furnacalis may influence adult mating generation number (voltinism) across different local environments. The degree to which voltine ecotype, geographic distance, or other ecological factors influence O. furnacalis population genetic structure remains uncertain. Genetic differentiation was estimated between voltine ecotypes collected from 8 locations. Mitochondrial haplotypes were significantly different between historically allopatric univoltine and bivoltine locations, but confounded by a strong correlation with geographic distance. In contrast, single nucleotide polymorphism (SNP) genotypes show low but significant levels of variation and a lack of influence of geographic distance between allopatric voltine locations. Regardless, 11 of 257 SNP loci were predicted to be under selection, suggesting population genetic homogenization except at loci proximal to factors putatively under selection. These findings provide evidence of haplotype divergent voltine ecotypes that may be maintained in allopatric and sympatric areas despite relatively high rates of nuclear gene flow, yet influence of voltinism on maintenance of observed haplotype divergence remains unresolved.

Baseline Susceptibility and Laboratory Selection of Resistance to Bt Cry1Ab Protein of Chinese Populations of Yellow Peach Moth, Conogethes punctiferalis (Guenée)

Yellow Peach Moth (YPM), Conogethes punctiferalis (Guenée), is one of the most destructive maize pests in the Huang-Huai-Hai summer maize region of China. Transgenic Bacillus thuringiensis (Bt) maize provides an effective means to control this insect pest in field trials. However, the establishment of Bt resistance to target pests is endangering the continued success of Bt crops. To use Bt maize against YPM, the baseline susceptibility of the local populations in the targeted areas needs to be verified. Diet-overlay bioassay results showed that all the fourteen YPM populations in China are highly susceptible to Cry1Ab. The LC50 values ranged from 0.35 to 2.38 ng/cm2 over the two years of the collection, and the difference between the most susceptible and most tolerant populations was sevenfold. The upper limit of the LC99 estimates of six pooled populations produced >99% larval mortality for representative eight populations collected in 2020 and was designated as diagnostic concentrations for monitoring susceptibility in YPM populations in China. Hence, we evaluated the laboratory selection of resistance in YPM to Cry1Ab using the diet-overlay bioassay method. Although the resistant ratio was generally low, YPM potentially could evolve resistance to Cry1Ab. The potential developmentof resistance by target pests points out the necessity to implement resistance management strategies for delaying the establishment of pest resistance to Bt crops.

Prediction of mitochondrial genome-wide variation through sequencing of mitochondrion-enriched extracts

Although mitochondrial DNA (mtDNA) haplotype variation is often applied for estimating population dynamics and phylogenetic relationships, economical and generalized methods for entire mtDNA genome enrichment prior to high-throughput sequencing are not readily available. This study demonstrates the utility of differential centrifugation to enrich for mitochondrion within cell extracts prior to DNA extraction, short-read sequencing, and assembly using exemplars from eight maternal lineages of the insect species, Ostrinia nubilalis. Compared to controls, enriched extracts showed a significant mean increase of 48.2- and 86.1-fold in mtDNA based on quantitative PCR, and proportion of subsequent short sequence reads that aligned to the O. nubilalis reference mitochondrial genome, respectively. Compared to the reference genome, our de novo assembled O. nubilalis mitochondrial genomes contained 82 intraspecific substitution and insertion/deletion mutations, and provided evidence for correction of mis-annotated 28 C-terminal residues within the NADH dehydrogenase subunit 4. Comparison to a more recent O. nubilalis mtDNA assembly from unenriched short-read data analogously showed 77 variant sites. Twenty-eight variant positions, and a triplet ATT codon (Ile) insertion within ATP synthase subunit 8, were unique within our assemblies. This study provides a generalizable pipeline for whole mitochondrial genome sequence acquisition adaptable to applications across a range of taxa.

Computational Screening of Potential Inhibitors of β-N-Acetyl-D-Hesosaminidases Using Combined Core-Fragment Growth and Pharmacophore Restraints

As a type of β-N-acetyl-D-hexosaminidase enzyme purified from the Ostriniafurnacalis (Asian corn borer), OfHex1 has been previously reported to participate in chitin degradation, indicating that it may be an ideal target for designing new environmentally friendly pesticides. Besides, a natural product, TMG-chitotriomycin, has been found to be an effective inhibitor of OfHex1, and some studies have shown that the interactions between TMG unit and residues in - 1 subsite of OfHex1 are very conservative and important, inspiring us to design new inhibitors of β-N-acetyl-D-hexosaminidase with a new strategy. In the present study, the virtual screening of TMG unit as the core fragment was conducted using the combined computational methods, such as docking, molecular dynamics, pharmacophore model, and pesticide-likeness rule. Nine compounds with the binding free energy lower than TMG-β-(GlcNAc)₂ were obtained. According to the decomposition energy and the interactions analysis, compounds 2, 3, 6 and 8, forming the hydrogen bonds with Val327 and Trp490, were considered as the possible lead structures for the further study. Our findings indicated that fragment-based lead discovery strategy might provide valuable insights into designing novel potential OfHex1 inhibitors.

Evolution of Asian Corn Borer Resistance to Bt Toxins Used Singly or in Pairs

Transgenic crops producing insecticidal proteins from Bacillus thuringiensis (Bt) have revolutionized pest control, but the benefits of this approach have been reduced by the evolution of resistance in pests. The widely adopted 'pyramid strategy' for delaying resistance entails transgenic crops producing two or more distinct toxins that kill the same pest. The limited experimental evidence supporting this strategy comes primarily from a model system under ideal conditions. Here we tested the pyramid strategy under nearly worst-case conditions, including some cross-resistance between the toxins in the pyramid. In a laboratory selection experiment with an artificial diet, we used Bt toxins Cry1Ab, Cry1F, and Cry1Ie singly or in pairs against Ostrinia furnacalis, one of the most destructive pests of corn in Asia. Under the conditions evaluated, pairs of toxins did not consistently delay the evolution of resistance relative to single toxins.

Introgression between divergent corn borer species in a region of sympatry: Implications on the evolution and adaptation of pest arthropods

The Asian corn borer, Ostrinia furnacalis, and European corn borer, O. nubilalis (Lepidoptera: Crambidae), cause damage to cultivated maize in spatially distinct geographies and have evolved divergent hydrocarbons as the basis of sexual communication. The Yili area of Xinjiang Uyghur Autonomous Region in China represents the only known region where O. furnacalis has invaded a native O. nubilalis range, and these two corn borer species have made secondary contact. Genetic differentiation was estimated between Ostrinia larvae collected from maize plants at 11 locations in Xinjiang and genotyped using high-throughput SNP and microsatellite markers. Maternal lineages were assessed by direct sequencing of mitochondrial cytochrome c oxidase subunit I and II haplotypes, and a high degree of genotypic diversity was demonstrated between lineages based on SNP genotypes. Furthermore, historical introgression was predicted among SNP genotypes only at sympatric locations in the Yili area, whereas in Xinjiang populations only O. furnacalis haplotypes were detected and no analogous introgressed genotypes were predicted. Our detection of putative hybrids and historical evidence of introgression defines Yili area as a hybrid zone between the species in normal ecological interactions and furthermore, might indicate that adaptive traits could spread even between seemingly divergent species through horizontal transmission. Results of this study indicate there may be a continuum in the degree of reproductive isolation between Ostrinia species and that the elegance of distinct and complete speciation based on modifications to the pheromone communication might need to be reconsidered.

Feeding regulation by neuropeptide Y on Asian corn borer Ostrinia furnacalis

The Asian Corn Borer Ostrinia furnacalis is a major agricultural pest. In this study, a full-length neuropeptide Y (npy) gene in O. furnacalis was sequenced and cloned from cDNA library, which contains an ORF of 273 bp by encoding 90 amino acid residues. The mature OfurNPY is composed of 29 amino acids with amidation in C-terminal. The spatiotemporal expression analysis showed that npy highest expression level was in the midgut of the fifth instar larvae (the gluttony period). When the expression of npy was knocked down by feeding or injecting dsNPY, larval food consumption, body size, and body weight were significantly inhibited compared to controls. These results indicate that NPY is an important regulator in the control of feeding of O. furnacalis.

A potential and novel type transgenic corn plant for control of the Corn Borer

The corn borer is a world-wide agricultural pest. In this study, a full-length neuropeptide F (npf) gene in Ostrinia furnacalis was sequenced and cloned from a cDNA library, in which the npf gene produces two splicing mRNA variants - npf1 and npf2 (with a 120 bp segment inserted into the npf1 sequence to generate npf2). A spatio-temporal expression analysis showed that the highest expression level of npf was in the midgut of 5^th instar larvae (the gluttony period), and their npf expression and food consumption were significantly promoted after food deprivation for 6 h. When npf was knocked down by double-stranded RNA for NPF, larval food intake, weight and body size were effectively inhibited through changes of a biosynthesis and metabolism pathway; i.e. gene silencing of NPF causes decreases of total lipid and glycogen and increases of trehalose production. Moreover, we produced transgenic corn plants with stably expressed dsNPF. Results showed that O. furnacalis larvae fed on these transgenic leaves had lower food consumption and smaller body size compared to controls. These results indicate that NPF is important in the feeding control of O. furnacalis and valuable for production of potential transgenic corn.

Variation of life-history traits of the Asian corn borer, Ostrinia furnacalis in relation to temperature and geographical latitude

Life-history traits from four geographical populations (tropical Ledong population [LD], subtropical Guangzhou [GZ] and Yongxiu populations, and temperate Langfang population [LF]) of the Asian corn borer, Ostrinia furnacalis were investigated at a wide range of temperatures (20-32°C). The larval and pupal times were significantly decreased with increasing rearing temperature, and growth rate was positively correlated with temperature. The relationship between body weight and rearing temperature in O. furnacalis did not follow the temperature-size rule (TSR); all populations exhibited the highest pupal and adult weights at high temperatures or intermediate temperatures. However, development time, growth rate, and body weight did not show a constant latitudinal gradient. Across all populations at each temperature, female were significantly bigger than males, showing a female-biased sexual size dimorphism (SSD). Contrary to Rensch's rule, the SSD tended to increase with rising temperature. The subtropical GZ population exhibited the largest degree of dimorphism while the temperate LF exhibited the smallest. Male pupae lose significantly more weight at metamorphosis compared to females. The proportionate weight losses of different populations were significantly different. Adult longevity was significantly decreased with increasing temperature. Between sexes, all populations exhibit a rather female-biased adult longevity. Finally, we discuss the adaptive significance of higher temperature-inducing high body weight in the moth's life history and why the moth exhibits the reverse TSR.