A viral histone H4 (CpBV-H4) is encoded in a polydnavirus, Cotesia plutellae bracovirus (CpBV). It plays a crucial role in parasitism of an endoparasitoid wasp, C. plutellae, against diamondback moth, Plutella xylostella, by altering host gene expression in an epigenetic mode by its N-terminal tail after joining host nucleosomes. Comparative transcriptomic analysis between parasitized and nonparasitized P. xylostella by RNA-Seq indicated that 1,858 genes were altered at more than two folds in expression levels at late parasitic stage, including 877 up-regulated genes and 981 down-regulated genes. Among parasitic factors altering host gene expression, CpBV-H4 alone explained 16.3% of these expressional changes. To characterize the joining sites of CpBV-H4 on host chromosomes, ChIP-Seq (chromatin immunoprecipitation followed by deep sequencing) was applied to chromatins extracted from parasitized larvae. It identified specific 538 ChIP targets. Joining sites were rich (60.2%) in AT sequence. Almost 40% of ChIP targets included short nucleotide repeat sequences presumably recognizable by transcriptional factors and chromatin remodeling factors. To further validate these CpBV-H4 targets, CpBV-H4 was transiently expressed in nonparasitized host at late larval stage and subjected to ChIP-Seq. Two kinds of ChIP-Seqs shared 51 core joining sites. Common targets were close (within 1 kb) to genes regulated at expression levels by CpBV-H4. However, other host genes not close to CpBV-H4 joining sites were also regulated by CpBV-H4. These results indicate that CpBV-H4 joins specific chromatin regions of P. xylostella and controls about one sixth of the total host genes that were regulated by C. plutellae parasitism in an epigenetic mode.

The extent of the conservation of synteny and gene order in aphids has been previously investigated only by comparing a small subset of linkage groups between the pea aphid Acyrthosiphon pisum and a few other aphid species. Here we compared the localization of eight A. pisum scaffolds (covering more than 5 Mb and 83 genes) in respect to the Drosophila melanogaster Muller elements identifying orthologous loci spanning all the four A. pisum chromosomes. Comparison of the genetic maps revealed a conserved synteny across different loci suggesting that the study of the fruit fly Muller elements could favour the identification of chromosomal markers useful for the study of chromosomal rearrangements in aphids. A. pisum is the first aphid species to have its genome sequenced and the finding that there are several chromosomal regions in synteny between Diptera and Hemiptera indicates that the genomic tools developed in A. pisum will be broadly useful not only for the study of other aphids but also for other insect species.

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Telomeric sequences are conserved across species. The most common sequence reported among insects is (TTAGG)n, but its universal occurrence is not a consensus because other canonical motifs have been reported. In the present study, we used fluorescence in situ hybridization (FISH) using telomeric probes with (TTAGG)6 repeats to describe the telomere composition of leafcutter ants. We performed the molecular cytogenetic characterization of six Acromyrmex Mayr, 1865 and one Atta Fabricius, 1804 species (Acromyrmex ambiguus (Emery, 1888), Ac. crassispinus (Forel, 1909), Ac. lundii (Guérin-Mèneville, 1838), Ac. nigrosetosus (Forel, 1908), Ac. rugosus (Smith, 1858), Ac. subterraneus subterraneus (Forel, 1893), and Atta sexdens (Linnaeus, 1758)) and described it using a karyomorphometric approach on their chromosomes. The diploid chromosome number 2n = 38 was found in all Acromyrmex species, and the karyotypic formulas were as follows: Ac. ambiguus 2K = 14M + 12SM + 8ST + 4A, Ac. crassispinus 2K = 12M + 20SM + 4ST + 2A, Ac. lundii 2K = 10M + 14SM + 10ST + 4A, Ac. nigrosetosus 2K = 12M + 14SM + 10ST + 2A, and Ac. subterraneus subterraneus 2K = 14M + 18SM + 4ST + 2A. The exact karyotypic formula was not established for Ac. rugosus. FISH analyses revealed the telomeric regions in all the chromosomes of the species studied in the present work were marked by the (TTAGG)6 sequence. These results reinforce the premise that Formicidae presents high homology between their genera for the presence of the canonical sequence (TTAGG)n.

Anthracnose diseases caused by Colletotrichum species are among the most common fungal diseases. These symptoms typically manifest as dark, sunken lesions on leaves, stems, and fruit. In China, mango anthracnose seriously affects fruit yield and quality. Genome sequencing of several species shows the presence of mini-chromosomes. These are thought to contribute to virulence, but their formation and activity remain to be fully elucidated. Here, we assembled 17 Colletotrichum genomes (16 isolated from mango plus one from persimmon) through PacBio long-read sequencing. Half of the assembled scaffolds had telomeric repeats at both ends indicating full-length chromosomes. Based on comparative genomics analysis at interspecies and intraspecies levels, we identified extensive chromosomal rearrangements events. We analyzed mini-chromosomes of Colletotrichum spp. and found large variation among close relatives. In C. fructicola, homology between core chromosomes and mini-chromosomes suggested that some mini-chromosomes were generated by recombination of core chromosomes. In C. musae GZ23-3, we found 26 horizontally transferred genes arranged in clusters on mini-chromosomes. In C. asianum FJ11-1, several potential pathogenesis-related genes on mini-chromosomes were upregulated, especially in strains with highly pathogenic phenotypes. Mutants of these upregulated genes showed obvious defects in virulence. Our findings provide insights into the evolution and potential relationships to virulence associated with mini-chromosomes. IMPORTANCE Colletotrichum is a cosmopolitan fungal genus that seriously affects fruit yield and quality of many plant species. Mini-chromosomes have been found to be related to virulence in Colletotrichum. Further examination of mini-chromosomes can help us elucidate some pathogenic mechanisms of Colletotrichum. In this study, we generated novel assemblies of several Colletotrichum strains. Comparative genomic analyses within and between Colletotrichum species were conducted. We then identified mini-chromosomes in our sequenced strains systematically. The characteristics and generation of mini-chromosomes were investigated. Transcriptome analysis and gene knockout revealed pathogenesis-related genes located on mini-chromosomes of C. asianum FJ11-1. This study represents the most comprehensive investigation of chromosome evolution and potential pathogenicity of mini-chromosomes in the Colletotrichum genus.

A very common protein modification in multicellular organisms is protein glycosylation or the addition of carbohydrate structures to the peptide backbone. Although the Class of the Insecta is the largest animal taxon on Earth, almost all information concerning glycosylation in insects is derived from studies with only one species, namely the fruit fly Drosophila melanogaster.

QTLs for insect resistance parameters, trichome type IV development, and more than 200 non-volatile metabolites, including 76 acyl sugars, all co-locate at the end of Chromosome 2 of Solanum galapagense. Host plant resistance is gaining importance as more and more insecticides are being banned due to environmental concerns. In tomato, resistance towards insects is found in wild relatives and has been attributed to the presence of glandular trichomes and their specific phytochemical composition. In this paper, we describe the results from a large-scale QTL mapping of data from whitefly resistance tests, trichome phenotyping and a comprehensive metabolomics analysis in a recombinant inbred line population derived from a cross between the cultivated Solanum lycopersicum and the wild relative S. galapagense, which is resistant to a range of pest insects. One major QTL (Wf-1) was found to govern the resistance against two different whitefly species. This QTL co-localizes with QTLs for the presence of trichomes type IV and V, as well as all 76 acyl sugars detected and about 150 other non-volatile phytochemicals, including methyl esters of the flavonols myricetin and quercetin. Based on these results, we hypothesize that Wf-1 is regulating the formation of glandular trichome type IV on the leaf epidermis, enabling the production and accumulation of bioactive metabolites in this type of trichomes.

Sex determination and the regulation of sexual dimorphism are among the most fascinating topics in modern biology. As the most species-rich group of sexually reproducing organisms on Earth, insects have multiple sex determination systems. Though sex chromosomes and sex-biased genes are well-studied in dozens of insects, their gene sequences are scattered in various databases. Moreover, a shortage of annotation hinders the deep mining of these data. Here, we collected the chromosome-level sex chromosome data of 49 insect species, including 34 X chromosomes, 15 Z chromosomes, 5 W chromosomes and 2 Y chromosomes. We also obtained Y-linked contigs of four insects species-Anopheles gambiae, Drosophila innubila, Drosophila yakuba and Tribolium castaneum. The unannotated chromosome-level sex chromosomes were annotated using a standard pipeline, yielding a total of 123 030 protein-coding genes, 2 159 427 repeat sequences, 894 miRNAs, 1574 rRNAs, 5105 tRNAs, 395 snoRNAs (small nucleolar RNA), 54 snRNAs (small nuclear RNA) and 5959 other ncRNAs (non-coding RNA). In addition, 36 781 sex-biased genes were identified by analyzing 62 RNA-seq (RNA sequencing) datasets. Together with 5707 sex-biased genes from the Drosophila genus collected from the Sex-Associated Gene Database, we obtained a total of 42 488 sex-biased genes from 13 insect species. All these data were deposited into InSexBase, a new user-friendly database of insect sex chromosomes and sex-biased genes. Database URL: http://www.insect-genome.com/Sexdb/.

Mechanisms and evolutionary dynamics of sex-determination systems are of particular interest in insect vectors of human pathogens like mosquitoes because novel control strategies aim to convert pathogen-transmitting females into nonbiting males, or rely on accurate sexing for the release of sterile males. In Aedes aegypti, the main vector of dengue and Zika viruses, sex determination is governed by a dominant male-determining locus, previously thought to reside within a small, nonrecombining, sex-determining region (SDR) of an otherwise homomorphic sex chromosome. Here, we provide evidence that sex chromosomes in Ae. aegypti are genetically differentiated between males and females over a region much larger than the SDR. Our linkage mapping intercrosses failed to detect recombination between X and Y chromosomes over a 123-Mbp region (40% of their physical length) containing the SDR. This region of reduced male recombination overlapped with a smaller 63-Mbp region (20% of the physical length of the sex chromosomes) displaying high male-female genetic differentiation in unrelated wild populations from Brazil and Australia and in a reference laboratory strain originating from Africa. In addition, the sex-differentiated genomic region was associated with a significant excess of male-to-female heterozygosity and contained a small cluster of loci consistent with Y-specific null alleles. We demonstrate that genetic differentiation between sex chromosomes is sufficient to assign individuals to their correct sex with high accuracy. We also show how data on allele frequency differences between sexes can be used to estimate linkage disequilibrium between loci and the sex-determining locus. Our discovery of large-scale genetic differentiation between sex chromosomes in Ae. aegypti lays a new foundation for mapping and population genomic studies, as well as for mosquito control strategies targeting the sex-determination pathway.

InSatDb presents an interactive interface to query information regarding microsatellite characteristics per se of five fully sequenced insect genomes (fruit-fly, honeybee, malarial mosquito, red-flour beetle and silkworm). InSatDb allows users to obtain microsatellites annotated with size (in base pairs and repeat units); genomic location (exon, intron, up-stream or transposon); nature (perfect or imperfect); and sequence composition (repeat motif and GC%). One can access microsatellite cluster (compound repeats) information and a list of microsatellites with conserved flanking sequences (microsatellite family or paralogs). InSatDb is complete with the insects information, web links to find details, methodology and a tutorial. A separate 'Analysis' section illustrates the comparative genomic analysis that can be carried out using the output. InSatDb is available at www.cdfd.org.in/insatdb.

Germline-restricted DNA has evolved in diverse animal taxa and is found in several vertebrate clades, nematodes, and flies. In these lineages, either portions of chromosomes or entire chromosomes are eliminated from somatic cells early in development, restricting portions of the genome to the germline. Little is known about why germline-restricted DNA has evolved, especially in flies, in which 3 diverse families, Chironomidae, Cecidomyiidae, and Sciaridae, carry germline-restricted chromosomes (GRCs). We conducted a genomic analysis of GRCs in the fungus gnat Bradysia (Sciara) coprophila (Diptera: Sciaridae), which has 2 large germline-restricted "L" chromosomes. We sequenced and assembled the genome of B. coprophila and used differences in sequence coverage and k-mer frequency between somatic and germline tissues to identify GRC sequence and compare it to the other chromosomes in the genome. We found that the GRCs in B. coprophila are large, gene rich, and have many genes with divergent homologs on other chromosomes in the genome. We also found that 2 divergent GRCs exist in the population we sequenced. GRC genes are more similar in sequence to genes from another Dipteran family (Cecidomyiidae) than to homologous genes from Sciaridae. This unexpected finding suggests that these chromosomes likely arose in Sciaridae through hybridization with a related lineage. These results provide a foundation from which to answer many questions about the evolution of GRCs in Sciaridae, such as how this hybridization event resulted in GRCs and what features on these chromosomes cause them to be restricted to the germline.

The field performance of Sterile Insect Technique (SIT) is improved by sex-sorting and releasing only sterile males. This can be accomplished by resource-intensive separation of males from females by morphology. Alternatively, sex-ratio biasing genetic constructs can be used to selectively remove one sex without the need for manual or automated sorting, but the resulting genetically engineered (GE) control agents would be subject to additional governmental regulation. Here we describe and demonstrate a genetic method for the batch production of non-GE males. This method could be applied to generate the heterogametic sex (XY, or WZ) in any organism with chromosomal sex determination. We observed up to 100% sex-selection with batch cultures of more than 103 individuals. Using a stringent transgene detection assay, we demonstrate the potential of mass production of transgene free males.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

The Hessian fly (Mayetiola destructor) is an important insect pest of wheat. It has tractable genetics, polytene chromosomes, and a small genome (158 Mb). Investigation of the Hessian fly presents excellent opportunities to study plant-insect interactions and the molecular mechanisms underlying genome imprinting and chromosome elimination. A physical map is needed to improve the ability to perform both positional cloning and comparative genomic analyses with the fully sequenced genomes of other dipteran species.

R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.

Tissue culture is performed to maintain isolated portions of multicellular organisms in an artificial milieu that is outside the individual organism and for considerable periods of time; cells derived from cultured explants are, in general, different from cells of the corresponding tissue in a living organism. The changes in cultured tissues that precede and often explain the subsequent cell proliferation of explant-derived cells have been partially studied, but little is known about the molecular and genomic basis of these changes. Comparative transcriptomics of intact and cultured (90 hours in MGM-450 insect medium) Bombyx mori tissues revealed that fewer genes represented a larger portion of the transcriptome of intact fat body tissues than of cultured fat body tissues. This analysis also indicated that expression of genes encoding sugar transporters and immune response proteins increased during culture and that expression of genes encoding lipoproteins and cuticle proteins decreased during culture. These results provide support for hypotheses that cultured tissues respond immunologically to surgery, adapt to the medium by accelerating sugar uptake, and terminate their identity as part of an intact organism by becoming independent of that organism.

The evolution of the diverse insect lineages is one of the most fascinating issues in evolutionary biology. Despite extensive research in this area, the resolution of insect phylogeny especially of interordinal relationships has turned out to be still a great challenge. One of the challenges for insect systematics is the radiation of the polyneopteran lineages with several contradictory and/or unresolved relationships. Here, we provide the first transcriptomic data for three enigmatic polyneopteran orders (Dermaptera, Plecoptera, and Zoraptera) to clarify one of the most debated issues among higher insect systematics. We applied different approaches to generate 3 data sets comprising 78 species and 1,579 clusters of orthologous genes. Using these three matrices, we explored several key mechanistic problems of phylogenetic reconstruction including missing data, matrix selection, gene and taxa number/choice, and the biological function of the genes. Based on the first phylogenomic approach including these three ambiguous polyneopteran orders, we provide here conclusive support for monophyletic Polyneoptera, contesting the hypothesis of Zoraptera + Paraneoptera and Plecoptera + remaining Neoptera. In addition, we employ various approaches to evaluate data quality and highlight problematic nodes within the Insect Tree that still exist despite our phylogenomic approach. We further show how the support for these nodes or alternative hypotheses might depend on the taxon- and/or gene-sampling.

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

This study focused on analyzing the distribution of microsatellites in holocentric chromosomes of the Triatominae subfamily, insect vectors of Chagas disease. We employed a non-denaturing FISH technique to determine the chromosomal distribution of sixteen microsatellites across twenty-five triatomine species, involving five genera from the two principal tribes: Triatomini and Rhodniini. Three main hybridization patterns were identified: strong signals in specific chromosomal regions, dispersed signals dependent on microsatellite abundance and the absence of signals in certain chromosomal regions or entire chromosomes. Significant variations in hybridization patterns were observed between Rhodniini and Triatomini species. Rhodniini species displayed weak and scattered hybridization signals, indicating a low abundance of microsatellites in their genomes. In contrast, Triatomini species exhibited diverse and abundant hybridization patterns, suggesting that microsatellites are a significant repetitive component in their genomes. One particularly interesting finding was the high abundance of GATA repeats, and to a lesser extent AG repeats, in the Y chromosome of all analyzed Triatomini species. In contrast, the Y chromosome of Rhodniini species did not show enrichment in GATA and AG repeats. This suggests that the richness of GATA repeats on the Y chromosome likely represents an ancestral trait specific to the Triatomini tribe. Furthermore, this information can be used to elucidate the evolutionary relationships between Triatomini and other groups of reduviids, contributing to the understanding of the subfamily's origin. Overall, this study provides a comprehensive understanding of the composition and distribution of microsatellites within Triatominae genomes, shedding light on their significance in the evolutionary processes of these species.

Eukaryotic genomes are packaged by Histone proteins in a structure called chromatin. There are different chromatin types. Euchromatin is typically associated with decondensed, transcriptionally active regions and heterochromatin to more condensed regions of the chromosomes. Methylation of Lysine 9 of Histone H3 (H3K9me) is a conserved biochemical marker of heterochromatin. In many organisms, heterochromatin is usually localized at telomeric as well as pericentromeric regions but can also be found at interstitial chromosomal loci. This distribution may vary in different species depending on their general chromosomal organization. Holocentric species such as Spodoptera frugiperda (Lepidoptera: Noctuidae) possess dispersed centromeres instead of a monocentric one and thus no observable pericentromeric compartment. To identify the localization of heterochromatin in such species we performed ChIP-Seq experiments and analyzed the distribution of the heterochromatin marker H3K9me2 in the Sf9 cell line and whole 4th instar larvae (L4) in relation to RNA-Seq data.