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.

Plasticity and genetic effects contribute to different axes of neural divergence in a community of mimetic Heliconius butterflies.

Changes in ecological preference, often driven by spatial and temporal variation in resource distribution, can expose populations to environments with divergent information content. This can lead to adaptive changes in the degree to which individuals invest in sensory systems and downstream processes, to optimize behavioural performance in different contexts. At the same time, environmental conditions can produce plastic responses in nervous system development and maturation, providing an alternative route to integrating neural and ecological variation. Here, we explore how these two processes play out across a community of Heliconius butterflies. Heliconius communities exhibit multiple Mullerian mimicry rings, associated with habitat partitioning across environmental gradients. These environmental differences have previously been linked to heritable divergence in brain morphology in parapatric species pairs. They also exhibit a unique dietary adaptation, known as pollen feeding, that relies heavily on learning foraging routes, or trap-lines, between resources, which implies an important environmental influence on behavioural development. By comparing brain morphology across 133 wild-caught and insectary-reared individuals from seven Heliconius species, we find strong evidence for interspecific variation in patterns of neural investment. These largely fall into two distinct patterns of variation; first, we find consistent patterns of divergence in the size of visual brain components across both wild and insectary-reared individuals, suggesting genetically encoded divergence in the visual pathway. Second, we find interspecific differences in mushroom body size, a central component of learning and memory systems, but only among wild caught individuals. The lack of this effect in common-garden individuals suggests an extensive role for developmental plasticity in interspecific variation in the wild. Finally, we illustrate the impact of relatively small-scale spatial effects on mushroom body plasticity by performing experiments altering the cage size and structure experienced by individual H. hecale. Our data provide a comprehensive survey of community level variation in brain structure, and demonstrate that genetic effects and developmental plasticity contribute to different axes of interspecific neural variation.

Decision tree for accurate infection timing in individuals newly diagnosed with HIV-1 infection.

BACKGROUND: There is today no gold standard method to accurately define the time passed since infection at HIV diagnosis. Infection timing and incidence measurement is however essential to better monitor the dynamics of local epidemics and the effect of prevention initiatives. METHODS: Three methods for infection timing were evaluated using 237 serial samples from documented seroconversions and 566 cross sectional samples from newly diagnosed patients: identification of antibodies against the HIV p31 protein in INNO-LIA, SediaTM BED CEIA and SediaTM LAg-Avidity EIA. A multi-assay decision tree for infection timing was developed. RESULTS: Clear differences in recency window between BED CEIA, LAg-Avidity EIA and p31 antibody presence were observed with a switch from recent to long term infection a median of 169.5, 108.0 and 64.5 days after collection of the pre-seroconversion sample respectively. BED showed high reliability for identification of long term infections while LAg-Avidity is highly accurate for identification of recent infections. Using BED as initial assay to identify the long term infections and LAg-Avidity as a confirmatory assay for those classified as recent infection by BED, explores the strengths of both while reduces the workload. The short recency window of p31 antibodies allows to discriminate very early from early infections based on this marker. BED recent infection results not confirmed by LAg-Avidity are considered to reflect a period more distant from the infection time. False recency predictions in this group can be minimized by elimination of patients with a CD4 count of less than 100 cells/mm3 or without no p31 antibodies. For 566 cross sectional sample the outcome of the decision tree confirmed the infection timing based on the results of all 3 markers but reduced the overall cost from 13.2 USD to 5.2 USD per sample. CONCLUSIONS: A step-wise multi assay decision tree allows accurate timing of the HIV infection at diagnosis at affordable effort and cost and can be an important new tool in studies analyzing the dynamics of local epidemics or the effects of prevention strategies.

A modified method to analyse cell proliferation using EdU labelling in large insect brains.

The study of neurogenesis is critical to understanding of the evolution of nervous systems. Within invertebrates, this process has been extensively studied in Drosophila melanogaster, which is the predominant model thanks to the availability of advanced genetic tools. However, insect nervous systems are extremely diverse, and by studying a range of taxa we can gain additional information about how nervous systems and their development evolve. One example of the high diversity of insect nervous system diversity is provided by the mushroom bodies. Mushroom bodies have critical roles in learning and memory and vary dramatically across species in relative size and the type(s) of sensory information they process. Heliconiini butterflies provide a useful snapshot of this diversity within a closely related clade. Within Heliconiini, the genus Heliconius contains species where mushroom bodies are 3-4 times larger than other closely related genera, relative to the rest of the brain. This variation in size is largely explained by increases in the number of Kenyon cells, the intrinsic neurons which form the mushroom body. Hence, variation in mushroom body size is the product of changes in cell proliferation during Kenyon cell neurogenesis. Studying this variation requires adapting labelling techniques for use in less commonly studied organisms, as methods developed for common laboratory insects often do not work. Here, we present a modified protocol for EdU staining to examine neurogenesis in large-brained insects, using Heliconiini butterflies as our primary case, but also demonstrating applicability to cockroaches, another large-brained insect.

Rapid expansion and visual specialisation of learning and memory centres in the brains of Heliconiini butterflies.

Changes in the abundance and diversity of neural cell types, and their connectivity, shape brain composition and provide the substrate for behavioral evolution. Although investment in sensory brain regions is understood to be largely driven by the relative ecological importance of particular sensory modalities, how selective pressures impact the elaboration of integrative brain centers has been more difficult to pinpoint. Here, we provide evidence of extensive, mosaic expansion of an integration brain center among closely related species, which is not explained by changes in sites of primary sensory input. By building new datasets of neural traits among a tribe of diverse Neotropical butterflies, the Heliconiini, we detected several major evolutionary expansions of the mushroom bodies, central brain structures pivotal for insect learning and memory. The genus Heliconius, which exhibits a unique dietary innovation, pollen-feeding, and derived foraging behaviors reliant on spatial memory, shows the most extreme enlargement. This expansion is primarily associated with increased visual processing areas and coincides with increased precision of visual processing, and enhanced long term memory. These results demonstrate that selection for behavioral innovation and enhanced cognitive ability occurred through expansion and localized specialization in integrative brain centers.

High level of novelty under the hood of convergent evolution.

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

The evolution of adult pollen feeding did not alter postembryonic growth in Heliconius butterflies.

For many animals, the availability and provision of dietary resources can vary markedly between juvenile and adult stages, often leading to a temporal separation of nutrient acquisition and use. Juvenile developmental programs are likely limited by the energetic demands of many adult tissues and processes with early developmental origins. Enhanced dietary quality in the adult stage may, therefore, alter selection on life history and growth patterns in juvenile stages. Heliconius are unique among butterflies in actively collecting and digesting pollen grains, which provide an adult source of essential amino acids. The origin of pollen feeding has therefore previously been hypothesized to lift constraints on larval growth rates, allowing Heliconius to spend less time as larvae when they are most vulnerable to predation. By measuring larval and pupal life-history traits across three pollen-feeding and three nonpollen-feeding Heliconiini, we provide the first test of this hypothesis. Although we detect significant interspecific variation in larval and pupal development, we do not find any consistent shift associated with pollen feeding. We discuss how this result may fit with patterns of nitrogen allocation, the benefits of nitrogenous stores, and developmental limitations on growth. Our results provide a framework for studies aiming to link innovations in adult Heliconius to altered selection regimes and developmental programs in early life stages.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius.

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

Heliconius butterflies have bright patterns on their wings that tell potential predators that they are toxic. As a result, predators learn to avoid eating them. Over time, unrelated species of butterflies have evolved similar patterns to avoid predation through a process known as Müllerian mimicry. Worldwide, there are over 180,000 species of butterflies and moths, most of which have different wing patterns. How do genes create this pattern diversity? And do butterflies use similar genes to create similar wing patterns? One of the genes involved in creating wing patterns is called cortex. This gene has a large region of DNA around it that does not code for proteins, but instead, controls whether cortex is on or off in different parts of the wing. Changes in this non-coding region can act like switches, turning regions of the wing into different colours and creating complex patterns, but it is unclear how these switches have evolved. Butterfly wings get their colour from tiny structures called scales, which each have their own unique set of pigments. In Heliconius butterflies, there are three types of scales: yellow/white scales, black scales, and red/orange/brown scales. Livraghi et al. used a DNA editing technique called CRISPR to find out whether the cortex gene affects scale type. First, Livraghi et al. confirmed that deleting cortex turned black and red scales yellow. Next, they used the same technique to manipulate the non-coding DNA around the cortex gene to see the effect on the wing pattern. This manipulation turned a black-winged butterfly into a butterfly with a yellow wing band, a pattern that occurs naturally in Heliconius butterflies. The next step was to find the mutation responsible for the appearance of yellow wing bands in nature. It turns out that a bit of extra genetic code, derived from so-called 'jumping genes', had inserted itself into the non-coding DNA around the cortex gene, 'flipping' the switch and leading to the appearance of the yellow scales. Genetic information contains the instructions to generate shape and form in most organisms. These instructions evolve over millions of years, creating everything from bacteria to blue whales. Butterfly wings are visual evidence of evolution, but the way their genes create new patterns isn't specific to butterflies. Understanding wing patterns can help researchers to learn how genetic switches control diversity across other species too.

High frequency of new recombinant forms in HIV-1 transmission networks demonstrated by full genome sequencing.

The HIV-1 epidemic in Belgium is primarily driven by MSM. In this patient population subtype B predominates but an increasing presence of non-B subtypes has been reported. We aimed to define to what extent the increasing subtype heterogeneity in a high at risk population induces the formation and spread of new recombinant forms. The study focused on transmission networks that reflect the local transmission to an important extent. One hundred and five HIV-1 transmission clusters were identified after phylogenetic analysis of 2849 HIV-1 pol sequences generated for the purpose of baseline drug resistance testing between 2013 and 2017. Of these 105 clusters, 62 extended in size during the last two years and were therefore considered as representing ongoing transmission. These 62 clusters included 774 patients in total. From each cluster between 1 and 3 representative patients were selected for near full-length viral genome sequencing. In total, the full genome sequence of 101 patients was generated. Indications for the presence of a new recombinant form were found for 10 clusters. These 10 clusters represented 105 patients or 13.6% of the patients covered by the study. The findings clearly show that new recombinant strains highly contribute to local transmission, even in an epidemic that is largely MSM and subtype B driven. This is an evolution that needs to be monitored as reshuffling of genome fragments through recombination may influence the transmissibility of the virus and the pathology of the infection. In addition, important changes in the sequence of the viral genome may challenge the performance of tests used for diagnosis, patient monitoring and drug resistance analysis.

Single genome sequencing of near full-length HIV-1 RNA using a limiting dilution approach.

Sequencing very long stretches of the HIV-1 genome can advance studies on virus evolution and in vivo recombination but remains technically challenging. We developed an efficient procedure to sequence near full-length HIV-1 RNA using a two-amplicon approach. The whole genome was successfully amplified for 107 (88%) of 121 plasma samples including samples from patients infected with HIV-1 subtype A1, B, C, D, F1, G, H, CRF01_AE and CRF02_AG. For the 17 samples with a viral load below 1000 c/ml and the 104 samples with a viral load above 1000 c/ml, the amplification efficiency was respectively 53% and 94%. The sensitivity of the method was further evaluated using limiting dilution of RNA extracted from a plasma pool containing an equimolar mixture of three HIV-1 subtypes (B, C and CRF02_AG) and diluted before and after cDNA generation. Both RNA and cDNA dilution showed comparable sensitivity and equal accuracy in reflecting the subtype distribution of the plasma pool. One single event of in vitro recombination was detected amongst the 41 sequences obtained after cDNA dilution but no indications for in vitro recombination were found after RNA dilution. In conclusion, a two-amplicon strategy and limiting dilution of viral RNA followed by reverse transcription, nested PCR and Sanger sequencing, allows near full genome sequencing of individual HIV-1 RNA molecules. This method will be a valuable tool in the study of virus evolution and recombination.

Quantification of total HIV-1 DNA in buffy coat cells, feasibility and potential added value for clinical follow-up of HIV-1 infected patients on ART.

BACKGROUND: Successfully treated HIV-1 infected patients have a sustained undetectable viral RNA load. In these cases the total HIV-1 DNA load may constitute a valuable tool to further follow the overall viral burden. The value of this marker outside of cure research has been rarely studied. OBJECTIVES: To develop a quantitative (q)PCR for total HIV-1 DNA quantification in buffy coat cells and to evaluate the value of this parameter in clinical follow-up. STUDY DESIGN: A qPCR using primers and a probe in the conserved HIV-1 LTR region was adapted for use on DNA extracted from buffy coat cells. Sensitivity, accuracy and reproducibility were evaluated using 8E5 cells and samples from naive and treatment experienced patients. The clinical value of DNA load analysis was assessed by testing 119 longitudinal samples from 9 patients before and after ART initiation and 249 cross sectional samples from therapy-experienced patients. RESULTS: Inter- and intra-assay coefficients of variability were 5.56 and 5.94 (%CV). HIV-1 DNA was detected in 249 of the 263 (94.7%) patients on ART for at least 5 months (median: 53 months; IQR: 28-84 months). The HIV-1 DNA load varied between 0.60 and 3.37 copies/106 blood cells and showed significant correlation with the pre-ART CD4+ T-cell count nadir and peak viral RNA load. ART initiation resulted in a slow and limited decline of the total HIV-1 DNA concentration. CONCLUSIONS: Quantification of total HIV-1 DNA from buffy coat cells is feasible, sensitive and reliable. Although determination of the on-therapy HIV-1 DNA load may be informative, regular testing has limited clinical value because of the very slow evolution.

Frequency of occurrence of HIV-1 dual infection in a Belgian MSM population.

INTRODUCTION: HIV-1 dual infection is a condition that results from infection with at least two HIV-1 variants from different sources. The scarceness of information on this condition is partly due to the fact that its detection is technically challenging. Using next-generation sequencing we defined the extent of HIV-1 dual infection in a cohort of men who have sex with men (MSM). MATERIAL & METHODS: Eighty-six MSM, diagnosed with HIV-1 subtype B infection between 2008 and 2013 were selected for next-generation sequencing of the HIV-1 envelope V3. Sequencing was performed on 2 plasma samples collected with an interval of > 6 months before the initiation of antiretroviral therapy. Maximum likelihood phylogenetic trees were inspected for dual infection, defined as the presence of two or more monophyletic clusters with ≥ 90% bootstrap support and a mean between-cluster genetic distance of ≥ 10%. To confirm dual infection, deep V3 sequencing of intermediate samples was performed as well as clonal sequencing of the HIV-1 protease-reverse transcriptase gene. RESULTS: Five of the 74 patients (6.8%) for whom deep sequencing was successful, showed clear evidence of dual infection. In 4 of them, the second strain was absent in the first sample but occurred in subsequent samples. This was highly suggestive for superinfection. In 3 patients both virus variants were of subtype B, in 2 patients at least one of the variants was a subtype B/non-B recombinant virus. CONCLUSIONS: Dual infection was confirmed in 6.8% of MSM diagnosed with HIV-1 in Belgium. This prevalence is probably an underestimation, because stringent criteria were used to classify viral variants as originating from a new infection event.

Experimental heatwaves compromise sperm function and cause transgenerational damage in a model insect.

Climate change is affecting biodiversity, but proximate drivers remain poorly understood. Here, we examine how experimental heatwaves impact on reproduction in an insect system. Male sensitivity to heat is recognised in endotherms, but ectotherms have received limited attention, despite comprising most of biodiversity and being more influenced by temperature variation. Using a flour beetle model system, we find that heatwave conditions (5 to 7 °C above optimum for 5 days) damaged male, but not female, reproduction. Heatwaves reduce male fertility and sperm competitiveness, and successive heatwaves almost sterilise males. Heatwaves reduce sperm production, viability, and migration through the female. Inseminated sperm in female storage are also damaged by heatwaves. Finally, we discover transgenerational impacts, with reduced reproductive potential and lifespan of offspring when fathered by males, or sperm, that had experienced heatwaves. This male reproductive damage under heatwave conditions provides one potential driver behind biodiversity declines and contractions through global warming.

Longitudinal sequencing of HIV-1 infected patients with low-level viremia for years while on ART shows no indications for genetic evolution of the virus.

HIV-infected patients on antiretroviral therapy (ART) may present low-level viremia (LLV) above the detection level of current viral load assays. In many cases LLV is persistent but does not result in overt treatment failure or selection of drug resistant viral variants. To elucidate whether LLV reflects active virus replication, we extensively sequenced pol and env genes of the viral populations present before and during LLV in 18 patients and searched for indications of genetic evolution. Maximum likelihood phylogenetic trees were inspected for temporal structure both visually and by linear regression analysis of root-to-tip and pairwise distances. Viral coreceptor tropism was assessed at different time points before and during LLV. In none of the patients consistent indications for genetic evolution were found over a median period of 4.8 years of LLV. As such these findings could not provide evidence that active virus replication is the main driver of LLV.