First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies.

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Biallelic variants in TAMM41 are associated with low muscle cardiolipin levels, leading to neonatal mitochondrial disease.

Mitochondrial disorders are clinically and genetically heterogeneous, with variants in mitochondrial or nuclear genes leading to varied clinical phenotypes. TAMM41 encodes a mitochondrial protein with cytidine diphosphate-diacylglycerol synthase activity: an essential early step in the biosynthesis of phosphatidylglycerol and cardiolipin. Cardiolipin is a mitochondria-specific phospholipid that is important for many mitochondrial processes. We report three unrelated individuals with mitochondrial disease that share clinical features, including lethargy at birth, hypotonia, developmental delay, myopathy, and ptosis. Whole exome and genome sequencing identified compound heterozygous variants in TAMM41 in each proband. Western blot analysis in fibroblasts showed a mild oxidative phosphorylation (OXPHOS) defect in only one of the three affected individuals. In skeletal muscle samples, however, there was severe loss of subunits of complexes I-IV and a decrease in fully assembled OXPHOS complexes I-V in two subjects as well as decreased TAMM41 protein levels. Similar to the tissue-specific observations on OXPHOS, cardiolipin levels were unchanged in subject fibroblasts but significantly decreased in the skeletal muscle of affected individuals. To assess the functional impact of the TAMM41 missense variants, the equivalent mutations were modeled in yeast. All three mutants failed to rescue the growth defect of the Δtam41 strains on non-fermentable (respiratory) medium compared with wild-type TAM41, confirming the pathogenicity of the variants. We establish that TAMM41 is an additional gene involved in mitochondrial phospholipid biosynthesis and modification and that its deficiency results in a mitochondrial disorder, though unlike families with pathogenic AGK (Sengers syndrome) and TAFAZZIN (Barth syndrome) variants, there was no evidence of cardiomyopathy.

Mimicry can drive convergence in structural and light transmission features of transparent wings in Lepidoptera.

Müllerian mimicry is a positive interspecific interaction, whereby co-occurring defended prey species share a common aposematic signal. In Lepidoptera, aposematic species typically harbour conspicuous opaque wing colour patterns with convergent optical properties among co-mimetic species. Surprisingly, some aposematic mimetic species have partially transparent wings, raising the questions of whether optical properties of transparent patches are also convergent, and of how transparency is achieved. Here, we conducted a comparative study of wing optics, micro and nanostructures in neotropical mimetic clearwing Lepidoptera, using spectrophotometry and microscopy imaging. We show that transparency, as perceived by predators, is convergent among co-mimics in some mimicry rings. Underlying micro- and nanostructures are also sometimes convergent despite a large structural diversity. We reveal that while transparency is primarily produced by microstructure modifications, nanostructures largely influence light transmission, potentially enabling additional fine-tuning in transmission properties. This study shows that transparency might not only enable camouflage but can also be part of aposematic signals.

Comparative transcriptome analysis at the onset of speciation in a mimetic butterfly-The Ithomiini Melinaea marsaeus.

Ecological speciation entails divergent selection on specific traits and ultimately on the developmental pathways responsible for these traits. Selection can act on gene sequences but also on regulatory regions responsible for gene expression. Mimetic butterflies are a relevant system for speciation studies because wing colour pattern (WCP) often diverges between closely related taxa and is thought to drive speciation through assortative mating and increased predation on hybrids. Here, we generate the first transcriptomic resources for a mimetic butterfly of the tribe Ithomiini, Melinaea marsaeus, to examine patterns of differential expression between two subspecies and between tissues that express traits that likely drive reproductive isolation; WCP and chemosensory genes. We sequenced whole transcriptomes of three life stages to cover a large catalogue of transcripts, and we investigated differential expression between subspecies in pupal wing discs and antennae. Eighteen known WCP genes were expressed in wing discs and 115 chemosensory genes were expressed in antennae, with a remarkable diversity of chemosensory protein genes. Many transcripts were differentially expressed between subspecies, including two WCP genes and one odorant receptor. Our results suggest that in M. marsaeus the same genes as in other mimetic butterflies are involved in traits causing reproductive isolation, and point at possible candidates for the differences in those traits between subspecies. Differential expression analyses of other developmental stages and body organs and functional studies are needed to confirm and expand these results. Our work provides key resources for comparative genomics in mimetic butterflies, and more generally in Lepidoptera.

Developmental Consequences of Defective ATG7-Mediated Autophagy in Humans.

BACKGROUND: Autophagy is the major intracellular degradation route in mammalian cells. Systemic ablation of core autophagy-related (ATG) genes in mice leads to embryonic or perinatal lethality, and conditional models show neurodegeneration. Impaired autophagy has been associated with a range of complex human diseases, yet congenital autophagy disorders are rare. METHODS: We performed a genetic, clinical, and neuroimaging analysis involving five families. Mechanistic investigations were conducted with the use of patient-derived fibroblasts, skeletal muscle-biopsy specimens, mouse embryonic fibroblasts, and yeast. RESULTS: We found deleterious, recessive variants in human ATG7, a core autophagy-related gene encoding a protein that is indispensable to classical degradative autophagy. Twelve patients from five families with distinct ATG7 variants had complex neurodevelopmental disorders with brain, muscle, and endocrine involvement. Patients had abnormalities of the cerebellum and corpus callosum and various degrees of facial dysmorphism. These patients have survived with impaired autophagic flux arising from a diminishment or absence of ATG7 protein. Although autophagic sequestration was markedly reduced, evidence of basal autophagy was readily identified in fibroblasts and skeletal muscle with loss of ATG7. Complementation of different model systems by deleterious ATG7 variants resulted in poor or absent autophagic function as compared with the reintroduction of wild-type ATG7. CONCLUSIONS: We identified several patients with a neurodevelopmental disorder who have survived with a severe loss or complete absence of ATG7, an essential effector enzyme for autophagy without a known functional paralogue. (Funded by the Wellcome Centre for Mitochondrial Research and others.).

Involvement of Arabidopsis BIG protein in cell death mediated by Myo-inositol homeostasis.

Programmed cell death (PCD) is essential for several aspects of plant life. We previously identified the mips1 mutant of Arabidopsis thaliana, which is deficient for the enzyme catalysing myo-inositol synthesis, and that displays light-dependent formation of lesions on leaves due to Salicylic Acid (SA) over-accumulation. Rationale of this work was to identify novel regulators of plant PCD using a genetic approach. A screen for secondary mutations that abolish the mips1 PCD phenotype identified a mutation in the BIG gene, encoding a factor of unknown molecular function that was previously shown to play pleiotropic roles in plant development and defence. Physiological analyses showed that BIG is required for lesion formation in mips1 via SA-dependant signalling. big mutations partly rescued transcriptomic and metabolomics perturbations as stress-related phytohormones homeostasis. In addition, since loss of function of the ceramide synthase LOH2 was not able to abolish cell death induction in mips1, we show that PCD induction is not fully dependent of sphingolipid accumulation as previously suggested. Our results provide further insights into the role of the BIG protein in the control of MIPS1-dependent cell death and also into the impact of sphingolipid homeostasis in this pathway.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata.

BACKGROUND: Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata. RESULTS: We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns. CONCLUSIONS: This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures.

Explaining the maintenance of adaptive diversity within populations is a long-standing goal in evolutionary biology, with important implications for conservation, medicine, and agriculture. Adaptation often leads to the fixation of beneficial alleles, and therefore it erodes local diversity so that understanding the coexistence of multiple adaptive phenotypes requires deciphering the ecological mechanisms that determine their respective benefits. Here, we show how antagonistic frequency-dependent selection (FDS), generated by natural and sexual selection acting on the same trait, maintains mimicry polymorphism in the toxic butterfly Heliconius numata Positive FDS imposed by predators on mimetic signals favors the fixation of the most abundant and best-protected wing-pattern morph, thereby limiting polymorphism. However, by using mate-choice experiments, we reveal disassortative mate preferences of the different wing-pattern morphs. The resulting negative FDS on wing-pattern alleles is consistent with the excess of heterozygote genotypes at the supergene locus controlling wing-pattern variation in natural populations of H. numata The combined effect of positive and negative FDS on visual signals is sufficient to maintain a diversity of morphs displaying accurate mimicry with other local prey, although some of the forms only provide moderate protection against predators. Our findings help understand how alternative adaptive phenotypes can be maintained within populations and emphasize the need to investigate interactions between selective pressures in other cases of puzzling adaptive polymorphism.

Selection of Valid Reference Genes for Reverse Transcription Quantitative PCR Analysis in Heliconius numata (Lepidoptera: Nymphalidae).

Identifying the genetic basis of adaptive variation is challenging in non-model organisms and quantitative real time PCR. is a useful tool for validating predictions regarding the expression of candidate genes. However, comparing expression levels in different conditions requires rigorous experimental design and statistical analyses. Here, we focused on the neotropical passion-vine butterflies Heliconius, non-model species studied in evolutionary biology for their adaptive variation in wing color patterns involved in mimicry and in the signaling of their toxicity to predators. We aimed at selecting stable reference genes to be used for normalization of gene expression data in RT-qPCR analyses from developing wing discs according to the minimal guidelines described in Minimum Information for publication of Quantitative Real-Time PCR Experiments (MIQE). To design internal RT-qPCR controls, we studied the stability of expression of nine candidate reference genes (actin, annexin, eF1α, FK506BP, PolyABP, PolyUBQ, RpL3, RPS3A, and tubulin) at two developmental stages (prepupal and pupal) using three widely used programs (GeNorm, NormFinder and BestKeeper). Results showed that, despite differences in statistical methods, genes RpL3, eF1α, polyABP, and annexin were stably expressed in wing discs in late larval and pupal stages of Heliconius numata This combination of genes may be used as a reference for a reliable study of differential expression in wings for instance for genes involved in important phenotypic variation, such as wing color pattern variation. Through this example, we provide general useful technical recommendations as well as relevant statistical strategies for evolutionary biologists aiming to identify candidate-genes involved adaptive variation in non-model organisms.

Characterisation of alleles of tomato light signalling genes generated by TILLING.

Targeting Induced Local Lesions IN Genomes (TILLING) combines chemical mutagenesis with high throughput screening to allow the generation of alleles of selected genes. In this study, TILLING has been applied to produce a series of mutations in genes encoding essential components of the tomato light signal transduction pathway in an attempt to enhance fruit nutritional quality. Point mutations to DEETIOLATED1 (DET1), which is responsible for the high pigment2 (hp2) tomato mutant, resulted in elevated levels of both carotenoid and phenylpropanoid phytonutrients in ripe fruit, whilst immature fruit showed increased chlorophyll content, photosynthetic capacity and altered fruit morphology. Furthermore, genotypes with mutations to the UV-DAMAGED DNA BINDING PROTEIN 1 (DDB1), COP1 and COP1like were also characterised. These genotypes largely did not display phenotypes characteristic of mutation to light signalling components but their characterisation has enabled interrogation of structure function relationships of the mutated genes.