Population genomics of two closely related anhydrobiotic midges reveals differences in adaptation to extreme desiccation

The sleeping chironomid Polypedilum vanderplanki is capable of anhydrobiosis, a striking example of adaptation to extreme desiccation. Tolerance to complete desiccation in this species is associated with the emergence of multiple paralogs of protective genes. One of the gene families highly expressed under anhydrobiosis and involved in this process are protein-L-isoaspartate (D-aspartate) O-methyltransferases (PIMTs). Recently, a closely related anhydrobiotic midge from Malawi, P. pembai, showing the ability to tolerate complete desiccation similar to that of P. vanderplanki, but experiences more frequent desiccation-rehydration cycles due to differences in ecology, was discovered. Here, we sequenced and assembled the genome of P. pembai and performed a population genomics analysis of several populations of P. vanderplanki and a population of P. pembai. We observe positive selection and radical changes in the genetic architecture of the PIMT locus between the two species, including multiple duplication events in the P. pembai lineage. In particular, PIMT-4, the most highly expressed of these PIMTs, is present in six copies in the P. pembai; these copies differ in expression profiles, suggesting possible sub- or neofunctionalization. The nucleotide diversity (π) of the genomic region carrying these new genes is decreased in P. pembai, but not in the orthologous region carrying the ancestral gene in P. vanderplanki, providing evidence for a selective sweep associated with post-duplication adaptation in the former. Overall, our results suggest an extensive recent and likely ongoing, adaptation of the mechanisms of anhydrobiosis.

Intracellular localization and gene expression analysis provide new insights on LEA proteins’ diversity in anhydrobiotic midge

Anhydrobiosis, an adaptive ability to withstand complete desiccation, is in insects limited to a single species: the nonbiting midge Polypedilum vanderplanki (the sleeping chironomid). Evolution of anhydrobiosis in a single representative of a large genus is associated with drastic changes in genome structure, including the emergence of new multimember gene families directly involved in desiccation tolerance. Among them, Late Embryogenesis Abundant (LEA) proteins, which protect other proteins from aggregation caused by desiccation, are believed to originate via horizontal gene transfer from a bacterial donor. To obtain new insights on the biological background of the expanded 27-member LEA protein group in P. vanderplanki, we investigated the expression of corresponding genes in a P. vanderplanki-derived cell line, capable of anhydrobiosis, in a normal state and during induction of desiccation tolerance. We found that all LEA proteins genes identified in P. vanderplanki’s genome except PvLea16 and PvLea17 are also expressed in Pv11 cells. Their expression was elevated in response to anhydrobiosis-inducing trehalose treatment. Expression patterns of PvLea genes were well preserved in Pv11 cells in comparison to P. vanderplanki’s larvae both in the control group and during the anhydrobiosis cycle. We also investigated localization of LEA proteins in Pv11 cells and Sf9 cells and found a different level of conservation in intracellular localization of the protein expressed in mammalian and insect cells.

Cooption of heat shock regulatory system for anhydrobiosis in the sleeping chironomid Polypedilum vanderplanki

Polypedilum vanderplanki is a striking and unique example of an insect that can survive almost complete desiccation. Its genome and a set of dehydration–rehydration transcriptomes, together with the genome of Polypedilum nubifer (a congeneric desiccation-sensitive midge), were recently released. Here, using published and newly generated datasets reflecting detailed transcriptome changes during anhydrobiosis, as well as a developmental series, we show that the TCTAGAA DNA motif, which closely resembles the binding motif of the Drosophila melanogaster heat shock transcription activator (Hsf), is significantly enriched in the promoter regions of desiccation-induced genes in P. vanderplanki, such as genes encoding late embryogenesis abundant (LEA) proteins, thioredoxins, or trehalose metabolism-related genes, but not in P. nubifer. Unlike P. nubifer, P. vanderplanki has double TCTAGAA sites upstream of the Hsf gene itself, which is probably responsible for the stronger activation of Hsf in P. vanderplanki during desiccation compared with P. nubifer. To confirm the role of Hsf in desiccation-induced gene activation, we used the Pv11 cell line, derived from P. vanderplanki embryo. After preincubation with trehalose, Pv11 cells can enter anhydrobiosis and survive desiccation. We showed that Hsf knockdown suppresses trehalose-induced activation of multiple predicted Hsf targets (including P. vanderplanki-specific LEA protein genes) and reduces the desiccation survival rate of Pv11 cells fivefold. Thus, cooption of the heat shock regulatory system has been an important evolutionary mechanism for adaptation to desiccation in P. vanderplanki.