l e t t e r sHow an insect evolves to become a successful herbivore is of profound biological and practical importance. Herbivores are often adapted to feed on a specific group of evolutionarily and biochemically related host plants 1 , but the genetic and molecular bases for adaptation to plant defense compounds remain poorly understood 2 . We report the first whole-genome sequence of a basal lepidopteran species, Plutella xylostella, which contains 18,071 protein-coding and 1,412 unique genes with an expansion of gene families associated with perception and the detoxification of plant defense compounds. A recent expansion of retrotransposons near detoxification-related genes and a wider system used in the metabolism of plant defense compounds are shown to also be involved in the development of insecticide resistance. This work shows the genetic and molecular bases for the evolutionary success of this worldwide herbivore and offers wider insights into insect adaptation to plant feeding, as well as opening avenues for more sustainable pest management.The global pest P. xylostella (Lepidoptera: Yponomeutidae) is thought to have coevolved with the crucifer plant family 3 ( Supplementary Fig. 1) and has become the most destructive pest of economically important food crops, including rapeseed, cauliflower and cabbage 4 . Recently, the total cost of damage and management worldwide was estimated at $4-5 billion per annum 5,6 . This insect is the first species to have evolved resistance to dichlorodiphenyltrichloroethane (DDT) in the 1950s 7 and to Bacillus thuringiensis (Bt) toxins in the 1990s 8 and has developed resistance to all classes of insecticide, making it increasingly difficult to control 9,10 . P. xylostella provides an exceptional system for understanding the genetic and molecular bases of how insect herbivores cope with the broad range of plant defenses and chemicals encountered in the environment (Supplementary Fig. 2).We used a P. xylostella strain (Fuzhou-S) collected from a field in Fuzhou in southeastern China (26.08 °N, 119.28 °E) for sequencing ( Supplementary Fig. 1). Whole-genome shotgun-based Illumina sequencing of single individuals (Supplementary Table 1), even after ten generations of laboratory inbreeding, resulted in a poor initial assembly (N50 = 2.4 kb, representing the size above which 50% of the total length of the sequences is included), owing to high levels of heterozygosity ( Supplementary Figs. 3 and 4 and Supplementary Table 2). Subsequently, we sequenced 100,800 fosmid clones (comprising ~10× the genome length) to a depth of 200× ( Supplementary Fig. 5 and Supplementary Tables 3-5), assembling the resulting sequence data into 1,819 scaffolds, with an N50 of 737 kb, spanning ~394 Mb of the genome sequence (version 1; Supplementary Fig. 6 and Supplementary Table 6). The assembly covered 85.5% of a set of protein-coding ESTs (Supplementary Tables 7 and 8) generated by transcriptome sequencing 11 . Alignment of a subject scaffold against a 126-kb BAC (GenBank GU058050) from an altern...

Herbivore specialists adapt to feed on a specific group of host plants by evolving various mechanisms to respond to plant defenses. Insects also possess complex gut microbiotas but their potential role in adaptation is poorly understood. Our previous study of the genome of diamondback moth, Plutella xylostella, revealed an intrinsic capacity to detoxify plant defense compounds, which is an important factor in its success as a pest. Here we expand on that work with a complete taxonomic and functional profile of the P. xylostella gut microbiota obtained by metagenomic sequencing. Gene enrichment in the metagenome, accompanied by functional identification, revealed an important role of specific gut bacteria in the breakdown of plant cell walls, detoxification of plant phenolics, and synthesis of amino acids. Microbes participating in these pathways mainly belonged to three highly abundant bacteria: Enterobacter cloacae, Enterobacter asburiae, and Carnobacterium maltaromaticum. Results show that while the gut microbial community may be complex, a small number of functionally active species can be disproportionally important. The presence of specific enzymes in the microbiota community, such as supporting amino acid synthesis, digestion and detoxification functions, demonstrates the beneficial interactions between P. xylostella and its gut microbiota. These interactions can be potential targets for manipulation to provide novel pest management approaches.

Cytochrome P450 monooxygenases are present in almost all organisms and can play vital roles in hormone regulation, metabolism of xenobiotics and in biosynthesis or inactivation of endogenous compounds. In the present study, a genome-wide approach was used to identify and analyze the P450 gene family of diamondback moth, Plutella xylostella, a destructive worldwide pest of cruciferous crops. We identified 85 putative cytochrome P450 genes from the P. xylostella genome, including 84 functional genes and 1 pseudogene. These genes were classified into 26 families and 52 subfamilies. A phylogenetic tree constructed with three additional insect species shows extensive gene expansions of P. xylostella P450 genes from clans 3 and 4. Gene expression of cytochrome P450s was quantified across multiple developmental stages (egg, larva, pupa and adult) and tissues (head and midgut) using P. xylostella strains susceptible or resistant to insecticides chlorpyrifos and fiprinol. Expression of the lepidopteran specific CYP367s predominantly occurred in head tissue suggesting a role in either olfaction or detoxification. CYP340s with abundant transposable elements and relatively high expression in the midgut probably contribute to the detoxification of insecticides or plant toxins in P. xylostella. This study will facilitate future functional studies of the P. xylostella P450s in detoxification.

Aflatoxins (AFs) are a group of highly oxygenated polyketidese-derived toxins mainly produced by Aspergillus flavus and A. parasiticus, whose biosynthesis mechanisms are extremely sophisticated. Methylation is known as the major form of epigenetic regulation, which is correlated with gene expression. As the DNA methylation inhibitor 5-azacytidine (5-AC) blocks AF production, we studied AFB1 metabolism and morphological changes of A. flavus by treatment with 5-AC in liquid culture. The results show that 5-AC caused a decrease in AF production and concurrent changes in morphology. In addition, we isolated a non-aflatoxigenic mutant of A. flavus, showing a significant reduction in pigment production, after 5-AC treatment. This mutant showed significant reduction in the expression of genes in the AF biosynthesis pathway, and conidia formation. Furthermore, as AF biosynthesis and oxidative stress are intimately related events, we assessed the viability of A. flavus to oxidative stress after treatment with 5-AC, which showed that the mutant was more sensitive to the strong oxidant hydrogen peroxide. We found that the non-aflatoxigenic mutant showed a decrease in reactive oxygen species (ROS) and metabolites indicative of oxidative stress, which may be caused by the disruption of the defence system against excessive ROS formation after 5-AC treatment. These data indicate that 5-AC, as an inactivator of DNA methyltransferase, plays a very important role in AFB1 metabolism and the development of A. flavus, which might provide an effective strategy to pre- or post-harvest control of AFs.

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D-genome) of the allopolyploid (AD-genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin-immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon-related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A-genome and related diploid species (B-, F- and G-genomes), indicating that they colonized the centromeres of D-genome lineage after the divergence of the A- and D- ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A- and D-subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D- and AD-genome species, yet localized to just the NORs in A-, B-, F-, and G-genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.

Aflatoxins (AFs) have always been regarded as the most effective carcinogens, posing a great threat to agriculture, food safety, and human health. Aspergillus flavus is the major producer of aflatoxin contamination in crops. The prevention and control of A. flavus and aflatoxin continues to be a global problem. In this study, we demonstrated that the cell-free culture filtrate of Aspergillus oryzae and a non-aflatoxigenic A. flavus can effectively inhibit the production of AFB1 and the growth and reproduction of A. flavus, indicating that both of the non-aflatoxigenic Aspergillus strains secrete inhibitory compounds. Further transcriptome sequencing was performed to analyze the inhibitory mechanism of A. flavus treated with fermenting cultures, and the results revealed that genes involved in the AF biosynthesis pathway and other biosynthetic gene clusters were significantly downregulated, which might be caused by the reduced expression of specific regulators, such as AflS, FarB, and MtfA. The WGCNA results further revealed that genes involved in the TCA cycle and glycolysis were potentially involved in aflatoxin biosynthesis. Our comparative transcriptomics also revealed that two conidia transcriptional factors, brlA and abaA, were found to be significantly downregulated, which might lead to the downregulation of conidiation-specific genes, such as the conidial hydrophobins genes rodA and rodB. In summary, our research provides new insights for the molecular mechanism of controlling AF synthesis to control the proliferation of A. flavus and AF pollution.

Aspergillus flavus is one of the most opportunistic pathogens invading many important oilseed crops and foodstuffs with such toxic secondary metabolites as aflatoxin (AF) and Cyclopiazonic acid. We previously used the DNA methylation inhibitor 5-azacytidine to treat with an AF-producing A. flavus A133 strain, and isolated a mutant (NT) of A. flavus, which displayed impaired abilities of AF biosynthesis and fungal development. In this study, gas chromatography–mass spectrometry (GC-MS) analysis was used to reveal the metabolic changes between these two strains. A total of 1181 volatiles were identified in these two strains, among which 490 volatiles were found in these two strains in vitro and 332 volatiles were found in vivo. The NT mutant was found to produce decreasing volatile compounds, among which most of the fatty acid-derived volatiles were significantly downregulated in the NT mutant compared to the A133 strain, which are important precursors for AF biosynthesis. Two antioxidants and most of the amino acids derived volatiles were found significantly upregulated in the NT mutant. Overall, our results reveal the difference of metabolic profiles in two different A. flavus isolates, which may provide valuable information for controlling infections of this fungal pathogen.