A Chromosome-Level Genome Assembly and Annotation for the Clouded Apollo Butterfly (Parnassius mnemosyne): A Species of Global Conservation Concern

The clouded apollo (Parnassius mnemosyne) is a palearctic butterfly distributed over a large part of western Eurasia, but population declines and fragmentation have been observed in many parts of the range. The development of genomic tools can help to shed light on the genetic consequences of the decline and to make informed decisions about direct conservation actions. Here, we present a high-contiguity, chromosome-level genome assembly of a female clouded apollo butterfly and provide detailed annotations of genes and transposable elements. We find that the large genome (1.5 Gb) of the clouded apollo is extraordinarily repeat rich (73%). Despite that, the combination of sequencing techniques allowed us to assemble all chromosomes (nc = 29) to a high degree of completeness. The annotation resulted in a relatively high number of protein-coding genes (22,854) compared with other Lepidoptera, of which a large proportion (21,635) could be assigned functions based on homology with other species. A comparative analysis indicates that overall genome structure has been largely conserved, both within the genus and compared with the ancestral lepidopteran karyotype. The high-quality genome assembly and detailed annotation presented here will constitute an important tool for forthcoming efforts aimed at understanding the genetic consequences of fragmentation and decline, as well as for assessments of genetic diversity, population structure, inbreeding, and genetic load in the clouded apollo butterfly.

Development and performance of a next generation sequencing (NGS) assay for monitoring of dd-cfDNA post solid organ transplantation

The aim of this study was to evaluate the analytical performance of a novel NGS assay, intended for monitoring of donor-derived cell-free DNA (dd-cfDNA), and describe its validity in clinical plasma samples from kidney transplanted patients. Artificial and clinical samples with increasing amounts of patient DNA were evaluated using NGS analysis of indel markers. Monitoring of dd-cfDNA with the NGS assay presented herein demonstrated a sensitivity of ≥0.1% dd-cfDNA and excellent accuracy (R2 0.99) throughout an extensive range of dd-cfDNA (0.1-30%). The precision of the test was determined for two levels (0.1% (LoD) and 1%) of dd-cfDNA. The between run precision (CV%) for the respective level was 16% and 9% and the corresponding result for the within run precision was 19% and 7%. To evaluate performance of the assay in clinical samples, 507 retrospective monitoring samples from 21 patients transplanted either with kidneys from living or deceased donors were analyzed. Monitoring samples were sampled at multiple time points from 24 h up to 90 days post-transplantation. We show that in one patient, increase of dd-cfDNA preceded increase of creatinine caused by acute cellular rejection by several days. In conclusion, the NGS assay displayed a combination of high sensitivity with good accuracy and precision in both artificial and clinical dd-cfDNA samples.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.

Migratory behaviour is positively associated with genetic diversity in butterflies

Migration is typically associated with risk and uncertainty at the population level, but little is known about its cost-benefit trade-offs at the species level. Migratory insects in particular often exhibit strong demographic fluctuations due to local bottlenecks and outbreaks. Here, we use genomic data to investigate levels of heterozygosity and long-term population size dynamics in migratory insects, as an alternative to classical local and short-term approaches such as regional field monitoring. We analyse whole-genome sequences from 97 Lepidoptera species and show that individuals of migratory species have significantly higher levels of genome-wide heterozygosity, a proxy for effective population size, than do nonmigratory species. Also, we contribute whole-genome data for one of the most emblematic insect migratory species, the painted lady butterfly (Vanessa cardui), sampled across its worldwide distributional range. This species exhibits one of the highest levels of genomic heterozygosity described in Lepidoptera (2.95 ± 0.15%). Coalescent modelling (PSMC) shows historical demographic stability in V. cardui, and high effective population size estimates of 2-20 million individuals 10,000 years ago. The study reveals that the high risks associated with migration and local environmental fluctuations do not seem to decrease overall genetic diversity and demographic stability in migratory Lepidoptera. We propose a "compensatory" demographic model for migratory r-strategist organisms in which local bottlenecks are counterbalanced by reproductive success elsewhere within their typically large distributional ranges. Our findings highlight that the boundaries of populations are substantially different for sedentary and migratory insects, and that, in the latter, local and even regional field monitoring results may not reflect whole population dynamics. Genomic diversity patterns may elucidate key aspects of an insect's migratory nature and population dynamics at large spatiotemporal scales.

Evaluating coevolution in a horizontally transmitted mutualism

Many interspecific interactions are shaped by coevolution. Transmission mode is thought to influence opportunities for coevolution within symbiotic interactions. Vertical transmission maintains partner fidelity, increasing opportunities for coevolution, but horizontal transmission may disrupt partner fidelity, potentially reducing opportunities for coevolution. Despite these predictions, the role of coevolution in the maintenance of horizontally transmitted symbioses is unclear. Leveraging a tractable insect-bacteria symbiosis, we tested for signatures of pairwise coevolution by assessing patterns of host-symbiont specialization. If pairwise coevolution defines the interaction, we expected to observe evidence of reciprocal specialization between hosts and their local symbionts. We found no evidence for local adaptation between sympatric lineages of Anasa tristis squash bugs and Caballeronia spp. symbionts across their native geographic range. We also found no evidence for specialization between three co-localized Anasa host species and their native Caballeronia symbionts. Our results demonstrate generalist dynamics underlie the interaction between Anasa insect hosts and their Caballeronia symbionts. We predict that selection from multiple host species may favor generalist symbiont traits through diffuse coevolution. Alternatively, selection for generalist traits may be a consequence of selection by hosts for fixed cooperative symbiont traits without coevolution.

Population genomics reveals variable patterns of immune gene evolution in monarch butterflies (Danaus plexippus)

Humoral and cellular immune responses provide animals with major defences against harmful pathogens. While it is often assumed that immune genes undergo rapid diversifying selection, this assumption has not been tested in many species. Moreover, it is likely that different classes of immune genes experience different levels of evolutionary constraint, resulting in varying selection patterns. We examined the evolutionary patterns for a set of 91 canonical immune genes of North American monarch butterflies (Danaus plexippus), using as an outgroup the closely related soldier butterfly (Danaus eresimus). As a comparison to these immune genes, we selected a set of control genes that were paired with each immune for approximate size and genomic location. As a whole, these immune genes had a significant but modest reduction in Tajima's D relative to paired-control genes, but otherwise did not show distinct patterns of population genetic variation or evolutionary rates. When further partitioning these immune genes into four functional classes (recognition, signalling, modulation, and effector), we found distinct differences among these groups. Relative to control genes, recognition genes exhibit increased nonsynonymous diversity and divergence, suggesting reduced constraints on evolution, and supporting the notion that coevolution with pathogens results in diversifying selection. In contrast, signalling genes showed an opposite pattern of reduced diversity and divergence, suggesting evolutionary constraints and conservation. Modulator and effector genes showed no statistical differences from controls. These results are consistent with patterns found in immune genes in fruit flies and Pieris butterflies, suggesting that consistent selective pressures on different classes of immune genes broadly govern the evolution of innate immunity among insects.

Host plant diet affects growth and induces altered gene expression and microbiome composition in the wood white (Leptidea sinapis) butterfly

In a time with decreasing biodiversity, especially among insects, a detailed understanding about specific resource utilization strategies is crucial. The physiological and behavioural responses to host switches in phytophagous insects are poorly understood. Earlier studies indicate that a host plant switch might be associated with distinctive molecular and physiological responses in different lineages. Expanding the assessment of such associations across Lepidoptera will reveal if there are general patterns in adaptive responses, or if each switch event is more of a unique character. We investigated host plant preference, fitness consequences, effects on expression profiles and gut microbiome composition in two common wood white (Leptidea sinapis) populations with different host plant preferences from the extremes of the species distribution area (Sweden and Catalonia). Our results show that female Catalonian wood whites lack preference for either host plant (Lotus corniculatus or L. dorycnium), while Swedish females laid significantly more eggs on L. corniculatus. Individuals from both populations reared on L. dorycnium had longer developmental times and smaller body size as adults. This indicates that both environmental and genetic factors determine the choice to use a specific host plant. Gene expression analysis revealed a more pronounced response to host plant in the Catalonian compared to the Swedish population. In addition, host plant treatment resulted in a significant shift in microbiome community structure in the Catalonian population. Together, this suggests that population specific plasticity associated with local conditions underlies host plant utilisation in wood whites.

Genomic evidence for gene flow between monarchs with divergent migratory phenotypes and flight performance

Monarch butterflies are known for their spectacular annual migration in eastern North America, with millions of monarchs flying up to 4,500 km to overwintering sites in central Mexico. Monarchs also live west of the Rocky Mountains, where they travel shorter distances to overwinter along the Pacific Coast. It is often assumed that eastern and western monarchs form distinct evolutionary units, but genomic studies to support this notion are lacking. We used a tethered flight mill to show that migratory eastern monarchs have greater flight performance than western monarchs, consistent with their greater migratory distances. However, analysing more than 20 million SNPs in 43 monarch genomes, we found no evidence for genomic differentiation between eastern and western monarchs. Genomic analysis also showed identical and low levels of genetic diversity, and demographic analyses indicated similar effective population sizes and ongoing gene flow between eastern and western monarchs. Gene expression analysis of a subset of candidate genes during active flight revealed differential gene expression related to nonmuscular motor activity. Our results demonstrate that eastern and western monarchs maintain migratory differences despite ongoing gene flow, and suggest that migratory differences between eastern and western monarchs are not driven by select major-effects alleles. Instead, variation in migratory distance and destination may be driven by environmentally induced differential gene expression or by many alleles of small effect.

Multilayered Tuning of Dosage Compensation and Z-Chromosome Masculinization in the Wood White (Leptidea sinapis) Butterfly

In species with genetic sex determination, dosage compensation can evolve to equal expression levels of sex-linked and autosomal genes. Current knowledge about dosage compensation has mainly been derived from male-heterogametic (XX/XY) model organisms, whereas less is understood about the process in female-heterogametic systems (ZZ/ZW). In moths and butterflies, downregulation of Z-linked expression in males (ZZ) to match the expression level in females (ZW) is often observed. However, little is known about the underlying regulatory mechanisms, or if dosage compensation patterns vary across ontogenetic stages. In this study, we assessed dynamics of Z-linked and autosomal expression levels across developmental stages in the wood white (Leptidea sinapis). We found that although expression of Z-linked genes in general was reduced compared with autosomal genes, dosage compensation was actually complete for some categories of genes, in particular sex-biased genes, but equalization in females was constrained to a narrower gene set. We also observed a noticeable convergence in Z-linked expression between males and females after correcting for sex-biased genes. Sex-biased expression increased successively across developmental stages, and male-biased genes were enriched on the Z-chromosome. Finally, all five core genes associated with the ribonucleoprotein dosage compensation complex male-specific lethal were detected in adult females, in correspondence with a reduction in the expression difference between autosomes and the single Z-chromosome. We show that tuning of gene dosage is multilayered in Lepidoptera and argue that expression balance across chromosomal classes may predominantly be driven by enrichment of male-biased genes on the Z-chromosome and cooption of available dosage regulators.

Dissecting the Effects of Selection and Mutation on Genetic Diversity in Three Wood White (Leptidea) Butterfly Species

The relative role of natural selection and genetic drift in evolution is a major topic of debate in evolutionary biology. Most knowledge spring from a small group of organisms and originate from before it was possible to generate genome-wide data on genetic variation. Hence, it is necessary to extend to a larger number of taxonomic groups, descriptive and hypothesis-based research aiming at understanding the proximate and ultimate mechanisms underlying both levels of genetic polymorphism and the efficiency of natural selection. In this study, we used data from 60 whole-genome resequenced individuals of three cryptic butterfly species (Leptidea sp.), together with novel gene annotation information and population recombination data. We characterized the overall prevalence of natural selection and investigated the effects of mutation and linked selection on regional variation in nucleotide diversity. Our analyses showed that genome-wide diversity and rate of adaptive substitutions were comparatively low, whereas nonsynonymous to synonymous polymorphism and substitution levels were comparatively high in Leptidea, suggesting small long-term effective population sizes. Still, negative selection on linked sites (background selection) has resulted in reduced nucleotide diversity in regions with relatively high gene density and low recombination rate. We also found a significant effect of mutation rate variation on levels of polymorphism. Finally, there were considerable population differences in levels of genetic diversity and pervasiveness of selection against slightly deleterious alleles, in line with expectations from differences in estimated effective population sizes.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Oxidative stress induced by fluoroquinolones on treatment for complicated urinary tract infections in Indian patients

The aim of the study is to examine the oxidative stress in patients on fluoroquinolones (ciprofloxacin, levofloxacin, gatifloxacin) therapy for complicated urinary tract infections and to correlate with plasma concentrations at different time intervals. Superoxide dismutase (SOD), glutathione, plasma antioxidant status and lipid peroxides were evaluated in 52 patients on different dosage regimens up to 5 days. There is significant and gradual elevation of lipid peroxide levels in patients on ciprofloxacin (3.6 ± 0.34 nmol/ml to 6.2 ± 0.94 nmol/ml) and levofloxacin (3.5 ± 0.84 nmol/ml to 5.1 ± 0.28 nmol/ml) dosage regimen but not with gatifloxacin (3.5 ± 0.84 nmol/ml to 3.74 ± 0.17 nmol.ml). There was substantial depletion in both SOD and glutathione levels particularly with ciprofloxacin. On the 5(th) day of treatment, plasma antioxidant status decreased by 77.6% %, 50.5%, 7.56% for ciprofloxacin, levofloxacin and gatifloxacin respectively. In conclusion ciprofloxacin and levofloxacin induce more reactive oxygen species that lead to cell damage than gatifloxacin irrespective of their concentrations in patient population.

The prognostic significance of the S-phase fraction and DNCA ploidy in node-negative estrogen receptor-positive breast cancer

The aim of this study was to determine the significance of the phase fraction S-phase fraction (SPF) and DNA ploidy as determined by DNA flow cytometry on the prognosis of patients with Stage 1 (TINOMO) and stage II A (T2NOMO) breast cancer who had estrogen receptor-positive tumors. A total of 45 patients with invasive ductal carcinoma were studied. The median follow-up was 3.8 years. The prognostic cut-off for the SPF of the 45 patients was 34% as determined by log rank test and this was used to classify patients into low (34%) and high (34%) SPF groups. 30 patients (67%) had low SPF. 36 (80%) of the tumors examined were diploid. Using the cut-off, the SPF had no prognostic value in Stage 1 breast cancer patients. There was significantly longer disease-free survival and overall survival in Stage II A patients with low SPF compared to Stage II A patients with high SPF (96% vs. 50% at 4 years, p = 0.02, and 100% vs. 62% at 4 years, p = 0.03, respectively). No significant prognostic value was noted with the DNA ploidy status in both groups of patients.