Frizzled2 receives WntA signaling during butterfly wing pattern formation.

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. Although the secreted ligand WntA has been shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologs of the Frizzled-family of Wnt receptors. Here, we show that CRISPR mosaic knockouts of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss of function in multiple nymphalids. Whereas WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning and shed light on the functional diversity of insect Frizzled receptors.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility.

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Pyrin Inflammasome Activation Defines Colchicine-Responsive SURF Patients from FMF and Other Recurrent Fevers.

Syndrome of undifferentiated recurrent fever (SURF) is characterized by recurrent fevers, a lack of confirmed molecular diagnosis, and a complete or partial response to colchicine. Despite the clinical similarities to familial Mediterranean fever (FMF), the underlying inflammatory mechanisms of SURF are not yet understood. We here analyzed the in vitro activation of the pyrin inflammasome in a cohort of SURF patients compared to FMF and PFAPA patients. Peripheral blood mononuclear cells (PBMC) were collected from SURF (both colchicine-treated and untreated), FMF, PFAPA patients, and healthy donors. PBMC were stimulated ex vivo with Clostridium difficile toxin A (TcdA) and a PKC inhibitor (UCN-01), in the presence or absence of colchicine. The assembly of the pyrin inflammasome was evaluated by measuring the presence of apoptosis-associated Speck-like protein containing caspase recruitment domain (ASC) specks in monocytes using flow cytometry. IL-1ß secretion was quantified using an ELISA assay. No differences in TcdA-induced activation of pyrin inflammasome were observed among FMF, PFAPA, and healthy donors. Untreated SURF patients showed a reduced response to TcdA, which was normalized after colchicine treatment. In contrast to FMF, SURF patients, similar to PFAPA patients and healthy donors, did not exhibit pyrin inflammasome activation in response to UCN-01-mediated pyrin dephosphorylation. These data demonstrate that in vitro functional analysis of pyrin inflammasome activation can differentiate SURF from FMF and PFAPA patients, suggesting the involvement of the pyrin inflammasome in the pathophysiology of SURF.

Recolonization of secondary forests by a locally extinct Caribbean anole through the lens of range expansion theory.

Disturbance and recovery dynamics are characteristic features of many ecosystems. Disturbance dynamics are widely studied in ecology and conservation biology. Still, we know less about the ecological processes that drive ecosystem recovery. The ecological processes that mediate ecosystem recovery stand at the intersection of many theoretical frameworks. Range expansion theory is one of these complementary frameworks that can provide unique insights into the population-level processes that mediate ecosystem recovery, particularly fauna recolonization. Although the biodiversity patterns that follow the fauna recolonization of recovering forests have been well described in the literature, the ecological processes at the population level that drive these patterns remain conspicuously unknown. In this study, we tested three fundamental predictions of range expansion theory during the recolonization of recovering forests in Puerto Rico by a shade specialist anole, Anolis gundlachi. Range expansion theory predicts that individuals at the early stages of recolonization (i.e., younger forests) would have a high prevalence of dispersive traits, experience less density dependence, and suffer less parasitism. To test these predictions, we conducted a chronosequence study applying space-for-time substitution where we compared phenotypic traits (i.e., body size, body condition, and relative limb size), population density, population growth rates, and Plasmodium parasitism rates among lizard populations living in young (<30 years), mid (~40-70 years), and old-growth forests (>75 years). Lizard populations in younger forests had lower densities, higher population growth rates, and lower rates of Plasmodium parasitism compared with old-growth forests. Still, while we found that individuals had larger body sizes, and longer forelimbs in young forests in one site, this result was not consistent among sites. This suggests a potential trade-off between the traits that provide a dispersal advantage during the initial stages of recolonization and those that are advantageous to establish in novel environmental conditions. Overall, our study emphasizes the suitability of range expansion theory to describe fauna recolonization but also highlights that the ecological processes that drive recolonization are time-dependent, complex, and nuanced.

SARS-CoV-2 Omicron Replacement of Delta as Predominant Variant, Puerto Rico.

We reconstructed the SARS-CoV-2 epidemic caused by Omicron variant in Puerto Rico by sampling genomes collected during October 2021-May 2022. Our study revealed that Omicron BA.1 emerged and replaced Delta as the predominant variant in December 2021. Increased transmission rates and a dynamic landscape of Omicron sublineage infections followed.

Widespread Gene Expression Divergence in Butterfly Sensory Tissues Plays a Fundamental Role During Reproductive Isolation and Speciation.

Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation.

Divergence of chemosensing during the early stages of speciation.

Chemosensory communication is essential to insect biology, playing indispensable roles during mate-finding, foraging, and oviposition behaviors. These traits are particularly important during speciation, where chemical perception may serve to establish species barriers. However, identifying genes associated with such complex behavioral traits remains a significant challenge. Through a combination of transcriptomic and genomic approaches, we characterize the genetic architecture of chemoperception and the role of chemosensing during speciation for a young species pair of Heliconius butterflies, Heliconius melpomene and Heliconius cydno We provide a detailed description of chemosensory gene-expression profiles as they relate to sensory tissue (antennae, legs, and mouthparts), sex (male and female), and life stage (unmated and mated female butterflies). Our results untangle the potential role of chemical communication in establishing barriers during speciation and identify strong candidate genes for mate and host plant choice behaviors. Of the 252 chemosensory genes, HmOBP20 (involved in volatile detection) and HmGr56 (a putative synephrine-related receptor) emerge as strong candidates for divergence in pheromone detection and host plant discrimination, respectively. These two genes are not physically linked to wing-color pattern loci or other genomic regions associated with visual mate preference. Altogether, our results provide evidence for chemosensory divergence between H. melpomene and H. cydno, two rarely hybridizing butterflies with distinct mate and host plant preferences, a finding that supports a polygenic architecture of species boundaries.

An Update on Familial Mediterranean Fever.

(1) Background: Familial Mediterranean Fever (FMF) is the prototypal autoinflammatory disease, characterized by recurrent bursts of neutrophilic inflammation. (2) Methods: In this study we look at the most recent literature on this condition and integrate it with novel information on treatment resistance and compliance. (3) Results: The canonical clinical presentation of FMF is in children with self-limited episodes of fever and polyserositis, associated with severe long-term complications, such as renal amyloidosis. It has been described anecdotally since ancient times, however only recently it has been characterized more accurately. We propose an updated overview on the main aspects of pathophysiology, genetics, diagnosis and treatment of this intriguing disease. (4) Conclusions: Overall, this review presents the all the main aspects, including real life outcome of the latest recommendation on treatment resistance of FMF, a disease, that not only helped understanding the pathophysiology of the auto inflammatory process but also the functioning of the innate immune system itself.

Novel Doublesex Duplication Associated with Sexually Dimorphic Development of Dogface Butterfly Wings.

Sexually dimorphic development is responsible for some of the most remarkable phenotypic variation found in nature. Alternative splicing of the transcription factor gene doublesex (dsx) is a highly conserved developmental switch controlling the expression of sex-specific pathways. Here, we leverage sex-specific differences in butterfly wing color pattern to characterize the genetic basis of sexually dimorphic development. We use RNA-seq, immunolocalization, and motif binding site analysis to test specific predictions about the role of dsx in the development of structurally based ultraviolet (UV) wing patterns in Zerene cesonia (Southern Dogface). Unexpectedly, we discover a novel duplication of dsx that shows a sex-specific burst of expression associated with the sexually dimorphic UV coloration. The derived copy consists of a single exon that encodes a DNA binding but no protein-binding domain and has experienced rapid amino-acid divergence. We propose the novel dsx paralog may suppress UV scale differentiation in females, which is supported by an excess of Dsx-binding sites at cytoskeletal and chitin-related genes with sex-biased expression. These findings illustrate the molecular flexibility of the dsx gene in mediating the differentiation of secondary sexual characteristics.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry.

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Type I interferon activation in RAS-associated autoimmune leukoproliferative disease (RALD).

RAS-associated autoimmune leukoproliferative disease (RALD) is a rare immune dysregulation syndrome caused by somatic gain-of-function mutations of either NRAS or KRAS gene in hematopoietic cells. We describe a 27-year-old patient presenting at 5 months of age with recurrent infections and generalized lymphadenopathy who developed a complex multi-organ autoimmune syndrome with hypogammaglobulinemia, partially controlled with oral steroids, hydroxichloroquine, mofetil mycophenolate and IVIG prophylaxis. Activation of type I interferon pathway was observed in peripheral blood. Since 18 years of age, the patient developed regenerative nodular hyperplasia of the liver evolving into hepatopulmonary syndrome. Whole exome sequencing analysis of the peripheral blood DNA showed the NRAS p.Gly13Asp mutation validated as somatic. Our report highlights the possibility of detecting somatic NRAS gene mutations in patients with inflammatory immune dysregulation and type I interferon activation.

Perfect mimicry between Heliconius butterflies is constrained by genetics and development.

Müllerian mimicry strongly exemplifies the power of natural selection. However, the exact measure of such adaptive phenotypic convergence and the possible causes of its imperfection often remain unidentified. Here, we first quantify wing colour pattern differences in the forewing region of 14 co-mimetic colour pattern morphs of the butterfly species Heliconius erato and Heliconius melpomene and measure the extent to which mimicking colour pattern morphs are not perfectly identical. Next, using gene-editing CRISPR/Cas9 KO experiments of the gene WntA, which has been mapped to colour pattern diversity in these butterflies, we explore the exact areas of the wings in which WntA affects colour pattern formation differently in H. erato and H. melpomene. We find that, while the relative size of the forewing pattern is generally nearly identical between co-mimics, the CRISPR/Cas9 KO results highlight divergent boundaries in the wing that prevent the co-mimics from achieving perfect mimicry. We suggest that this mismatch may be explained by divergence in the gene regulatory network that defines wing colour patterning in both species, thus constraining morphological evolution even between closely related species.

Next generation sequencing panel in undifferentiated autoinflammatory diseases identifies patients with colchicine-responder recurrent fevers.

OBJECTIVES: The number of innate immune system disorders classified as systemic autoinflammatory diseases (SAID) has increased in recent years. More than 70% of patients with clinical manifestations of SAID did not receive a molecular diagnosis, thus being classed as so-called undifferentiated or undefined SAID (uSAID). The aim of the present study was to evaluate a next-generation sequencing (NGS)-based clinically oriented protocol in patients with uSAID. METHODS: We designed a NGS panel that included 41 genes clustered in seven subpanels. Patients with uSAID were classified into different groups according to their clinical features and sequenced for the coding portions of the 41 genes. RESULTS: Fifty patients were enrolled in the study. Thirty-four patients (72%) displayed recurrent fevers not consistent with a PFAPA phenotype. Sixteen patients displayed a chronic inflammatory disease course. A total of 100 gene variants were found (mean 2 per patient; range 0-6), a quarter of which affected suspected genes. Mutations with a definitive diagnostic impact were detected in two patients. Patients with genetically negative recurrent fevers displayed a prevalent gastrointestinal, skin and articular involvement. Patients responded to steroids on demands (94%) and colchicine, with a response rate of 78%. CONCLUSION: Even with a low molecular diagnostic rate, a NGS-based approach is able to provide a final diagnosis in a proportion of uSAID patients with evident cost-effectiveness. It also allows the identification of a subgroup of genetically negative patients with recurrent fever responding to steroid on demand and colchicine.

Clinical characteristics and genetic analyses of 187 patients with undefined autoinflammatory diseases.

OBJECTIVES: To describe the clinical characteristics, treatment response and genetic findings in a large cohort of patients with undefined systemic autoinflammatory diseases (SAIDs). METHODS: Clinical and genetic data from patients with undefined SAIDs were extracted from the Eurofever registry, an international web-based registry that retrospectively collects clinical information on patients with autoinflammatory diseases. RESULTS: This study included 187 patients. Seven patients had a chronic disease course, 180 patients had a recurrent disease course. The median age at disease onset was 4.3 years. Patients had a median of 12 episodes per year, with a median duration of 4 days. Most commonly reported symptoms were arthralgia (n=113), myalgia (n=86), abdominal pain (n=89), fatigue (n=111), malaise (n=104) and mucocutaneous manifestations (n=128). In 24 patients, relatives were affected as well. In 15 patients, genetic variants were found in autoinflammatory genes. Patients with genetic variants more often had affected relatives compared with patients without genetic variants (p=0.005). Most patients responded well to non-steroidal anti-inflammatory drugs (NSAIDs), corticosteroids, colchicine and anakinra. Complete remission was rarely achieved with NSAIDs alone. Notable patterns were found in patients with distinctive symptoms. Patients with pericarditis (n=11) were older at disease onset (33.8 years) and had fewer episodes per year (3.0/year) compared with other patients. Patients with an intellectual impairment (n=8) were younger at disease onset (2.2 years) and often had relatives affected (28.6%). CONCLUSION: This study describes the clinical characteristics of a large cohort of patients with undefined SAIDs. Among these, patients with pericarditis and intellectual impairment appear to comprise distinct subsets.

Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography.

Sex chromosomes are disproportionately involved in reproductive isolation and adaptation. In support of such a "large-X" effect, genome scans between recently diverged populations and species pairs often identify distinct patterns of divergence on the sex chromosome compared to autosomes. When measures of divergence between populations are higher on the sex chromosome compared to autosomes, such patterns could be interpreted as evidence for faster divergence on the sex chromosome, that is "faster-X", barriers to gene flow on the sex chromosome. However, demographic changes can strongly skew divergence estimates and are not always taken into consideration. We used 224 whole-genome sequences representing 36 populations from two Heliconius butterfly clades (H. erato and H. melpomene) to explore patterns of Z chromosome divergence. We show that increased divergence compared to equilibrium expectations can in many cases be explained by demographic change. Among Heliconius erato populations, for instance, population size increase in the ancestral population can explain increased absolute divergence measures on the Z chromosome compared to the autosomes, as a result of increased ancestral Z chromosome genetic diversity. Nonetheless, we do identify increased divergence on the Z chromosome relative to the autosomes in parapatric or sympatric species comparisons that imply postzygotic reproductive barriers. Using simulations, we show that this is consistent with reduced gene flow on the Z chromosome, perhaps due to greater accumulation of incompatibilities. Our work demonstrates the importance of taking demography into account to interpret patterns of divergence on the Z chromosome, but nonetheless provides evidence to support the Z chromosome as a strong barrier to gene flow in incipient Heliconius butterfly species.

Population genomics of parallel hybrid zones in the mimetic butterflies, H. melpomene and H. erato.

Hybrid zones can be valuable tools for studying evolution and identifying genomic regions responsible for adaptive divergence and underlying phenotypic variation. Hybrid zones between subspecies of Heliconius butterflies can be very narrow and are maintained by strong selection acting on color pattern. The comimetic species, H. erato and H. melpomene, have parallel hybrid zones in which both species undergo a change from one color pattern form to another. We use restriction-associated DNA sequencing to obtain several thousand genome-wide sequence markers and use these to analyze patterns of population divergence across two pairs of parallel hybrid zones in Peru and Ecuador. We compare two approaches for analysis of this type of data-alignment to a reference genome and de novo assembly-and find that alignment gives the best results for species both closely (H. melpomene) and distantly (H. erato, ∼15% divergent) related to the reference sequence. Our results confirm that the color pattern controlling loci account for the majority of divergent regions across the genome, but we also detect other divergent regions apparently unlinked to color pattern differences. We also use association mapping to identify previously unmapped color pattern loci, in particular the Ro locus. Finally, we identify a new cryptic population of H. timareta in Ecuador, which occurs at relatively low altitude and is mimetic with H. melpomene malleti.