Transcriptomic response of female adult moths to host and non-host plants in two closely related species

BACKGROUND

Divergent selection has been shown to promote speciation in many taxa and especially in phytophagous insects. In the Ostrinia species complex, the European corn borer (ECB) and adzuki bean borer (ABB) are two sibling species specialized to different host plants. The first is a well-known maize pest, whereas the second is a polyphagous species associated with various dicotyledons. Their specialization to host plants is driven by morphological, behavioral and physiological adaptations. In particular, previous studies have shown that ECB and ABB display marked behavior with regard to plant choice during oviposition, involving specific preference and avoidance mechanisms. In this study, our goal was to identify the mechanisms underlying this host-plant specialization in adult females through an analysis of their gene expression. We assembled and annotated a de novo reference transcriptome and measured differences in gene expression between ECB and ABB females, and between environments. We related differentially expressed genes to host preference behavior, and highlighted the functional categories involved. We also conducted a specific analysis of chemosensory genes, which are considered to be good candidates for host recognition before oviposition.

RESULTS

We recorded more differentially expressed genes in ECB than in ABB samples, and noticed that the majority of genes potentially involved in the host preference were different between the two species. At the functional level, the response to plant environment in adult females involved many processes, including the metabolism of carbohydrates, lipids, proteins, and amino acids; detoxification mechanisms and immunity; and the chemosensory repertoire (as expected). Until now, most of the olfactory receptors described in Ostrinia spp. had been tested for their putative role in pheromone recognition by males. Here we observed that one specific olfactory receptor was clearly associated with ECB's discrimination between maize and mugwort conditions, highlighting a potential new candidate involved in plant odor discrimination in adult females.

CONCLUSIONS

Our results are a first step toward the identification of candidate genes and functions involved in chemosensory processes, carbohydrate metabolism, and virus and retrovirus dynamics. These candidates provide new avenues for research into understanding the role of divergent selection between different environments in species diversification.

Larval transcriptomic response to host plants in two related phytophagous lepidopteran species: implications for host specialization and species divergence

Background

Most phytophagous insects have morphological, behavioral and physiological adaptations allowing them to specialize on one or a few plant species. Identifying the mechanisms involved in host plant specialization is crucial to understand the role of divergent selection between different environments in species diversification, and to identify sustainable targets for the management of insect pest species. In the present study, we measured larval phenotypic and transcriptomic responses to host plants in two related phytophagous lepidopteran species: the European corn borer (ECB), a worldwide pest of maize, and the adzuki bean borer (ABB), which feeds of various dicotyledons. Our aim was to identify the genes and functions underlying host specialization and/or divergence between ECB and ABB.

Results

At the phenotypic level, we observed contrasted patterns of survival, weight gain and developmental time between ECB and ABB, and within ECB and ABB reared on two different host plants. At the transcriptomic level, around 8% of the genes were differentially expressed (DE) between species and/or host plant. 70% of these DE genes displayed a divergent pattern of expression between ECB and ABB, regardless of the host, while the remaining 30% were involved in the plastic response between hosts. We further categorized plastic DE genes according to their parallel or opposite pattern between ECB and ABB to specifically identify candidate genes involved in the species divergence by host specialization. These candidates highlighted a comprehensive response, involving functions related to plant recognition, digestion, detoxification, immunity and development. Last, we detected viral, bacterial, and yeast genes whose incidence contrasted ECB and ABB samples, and maize and mugwort conditions. We suggest that these microorganism communities might influence the survival, metabolism and defense patterns observed in ECB and ABB larvae.

Conclusions

The comprehensive approach developed in the present study allowed to identify phenotypic specialization patterns and underlying candidate molecular mechanisms, and highlighted the putative role of microorganisms in the insect-host plant interaction. These findings offer the opportunity to pinpoint specific and sustainable molecular or physiological targets for the regulation of ECB pest populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4589-x) contains supplementary material, which is available to authorized users.

Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

De novo genome and transcriptome resources of the Adzuki bean borer Ostrinia scapulalis (Lepidoptera: Crambidae)

We present a draft genome assembly with a de novo prediction and automated functional annotation of coding genes, and a reference transcriptome of the Adzuki bean borer, Ostrinia scapulalis, based on RNA sequencing of various tissues and developmental stages. The genome assembly spans 419 Mb, has a GC content of 37.4% and includes 26,120 predicted coding genes. The reference transcriptome holds 33,080 unigenes and contains a high proportion of a set of genes conserved in eukaryotes and arthropods, used as quality assessment of the reconstructed transcripts. The new genomic and transcriptomic data presented here significantly enrich the public sequence databases for the Crambidae and Lepidoptera, and represent useful resources for future researches related to the evolution and the adaptation of phytophagous moths. The genome and transcriptome assemblies have been deposited and made accessible via a NCBI BioProject (id PRJNA390510) and the LepidoDB database (http://bipaa.genouest.org/sp/ostrinia_scapulalis/).

Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations

Unravelling variation among taxonomic orders regarding the rate of evolution in microsatellites is crucial for evolutionary biology and population genetics research. The mean mutation rate of microsatellites tends to be lower in arthropods than in vertebrates, but data are scarce and mostly concern accumulation of mutations in model species. Based on parent-offspring segregations and a hierarchical Bayesian model, the mean rate of mutation in the orthopteran insect Schistocerca gregaria was estimated at 2.1e(-4) per generation per untranscribed dinucleotide locus. This is close to vertebrate estimates and one order of magnitude higher than estimates from species of other arthropod orders, such as Drosophila melanogaster and Daphnia pulex. We also found evidence of a directional bias towards expansions even for long alleles and exceptionally large ranges of allele sizes. Finally, at transcribed microsatellites, the mean rate of mutation was half the rate found at untranscribed loci and the mutational model deviated from that usually considered, with most mutations involving multistep changes that avoid disrupting the reading frame. Our direct estimates of mutation rate were discussed in the light of peculiar biological and genomic features of S. gregaria, including specificities in mismatch repair and the dependence of its activity to allele length. Shedding new light on the mutational dynamics of grasshopper microsatellites is of critical importance for a number of research fields. As an illustration, we showed how our findings improve microsatellite application in population genetics, by obtaining a more precise estimation of S. gregaria effective population size from a published data set based on the same microsatellites.

Host specialization involving attraction, avoidance and performance, in two phytophagous moth species

Host specialization plays a key role in the extreme diversification of phytophagous insects. Whereas proximate mechanisms of specialization have been studied extensively, their consequences for species divergence remain unclear. Preference for, and performance on hosts are thought to be a major source of divergence in phytophagous insects. We assessed these major components of specialization in two moth species, the European corn borer (ECB) and the Adzuki bean borer (ABB), by testing their oviposition behaviour in different conditions (choice or no-choice set-ups) and their performances, by reciprocal transplant at the larval stage on the usual host and an alternative host plant. We demonstrated that both ABB and ECB have a strong preference for their host plants for oviposition, but that relative larval performances on the usual host and an alternative host differed according to the experiment and the trait considered (weight or survival). Finally, we show for the first time that the preference for maize in ECB conceals a strong avoidance of mugwort. The differences in performance, attraction and avoidance between ECB and ABB are discussed in the light of the underlying mechanisms and divergence process.

Eleven polymorphic microsatellite markers for Oedaleus decorus (Orthoptera, Acrididae), an endangered grasshopper in Central Europe

We isolated and characterized 11 microsatellite loci in the grasshopper Oedaleus decorus (Orthoptera: Acrididae), an endangered species in Central Europe. Polymorphism was studied from two populations, one out of two populations known from Switzerland (n = 20 individuals) and one site from south of France (n = 20). The number of alleles and the expected heterozygosity ranged from five to 12 and from 0.559 to 0.898, respectively, in the Swiss population, and from 14 to 23 and from 0.895 to 0.974, respectively, in the French population. These microsatellite markers are suitable for further conservation genetic studies of O. decorus.

Long microsatellites and unusually high levels of genetic diversity in the Orthoptera

Much remains to be learned about the mutational processes governing the evolution of microsatellite repeat regions and the associated levels of genetic diversity observed at microsatellite markers across populations or species. An extensive survey of microsatellite variation in 210 insect species from six major orders revealed that within Orthopterans, which are characterized by giant genomes, levels of genetic diversity were ~20% higher and microsatellite repeat arrays were longer than in any other group. Because of the mutation dependence on repeat length, this result suggests a higher microsatellite loci mutation rate in the Orthoptera. We deem it plausible that differences among insect orders, either in mismatch repair systems or in abundance of transposable element-derived microsatellites, can shape the size distribution of both genomes and microsatellite repeat regions. Our findings emphasise that observed levels of genetic diversity can greatly vary across species (orders at least) because of molecular differences in the mechanisms that determine microsatellite size, and are therefore critical to conservation and population genetics studies, where microsatellite repeat variability is primarily interpreted in terms of population demography and history.

From Russia with lobe: genetic differentiation in trilobed uncus Ostrinia spp. follows food plant, not hairy legs

Trilobed uncus taxa of the genus Ostrinia (Lepidoptera, Crambidae) illustrate the complex relationship, at early stages of speciation, between reproductive isolation and differentiation in morphology, resource use and genetic variation. On the basis of behaviour and ecology, we recently hypothesized that individuals with small mid-tibiae belong to two distinct species depending on host plant--O. nubilalis and O. scapulalis sensu Frolov et al. (2007) feeding on maize and on a number of dicotyledons, respectively. Individuals with small, medium or massive mid-tibiae would all belong to O. scapulalis as long as they feed on these dicotyledons. This contrasts with previous taxonomy, which distinguished three species by male mid-tibia morphology, regardless of host plant. Here, we test our hypothesis by examining the genetic structure of Ostrinia populations from regions with mid-tibia polymorphism--Western Russia and Kazakhstan--and comparing it with that of French populations where only small mid-tibiae occur. Results support two predictions: (1) maize- and dicotyledon-collected populations are genetically differentiated from each other like in France, and (2) dicotyledon-collected populations show no genetic evidence of consisting of more than one species. Between-species differentiation was unrelated to geographic distance, despite significant isolation by distance within species. The distinction between two and only two species differing by host plant thus holds at continental scale. Interestingly, one microsatellite locus contributed ∼10 times more than the others to differentiation between both taxa. This deserves further investigation, as it might reveal a linkage between this outlier and loci involved in host-plant adaptation and/or reproductive isolation.

Genetic structure and gene flow in French populations of two Ostrinia taxa: host races or sibling species?

Most models of ecological speciation concern phytophagous insects in which speciation is thought to be driven by host shifts and subsequent adaptations of populations. Despite the ever-increasing number of studies, the current evolutionary status of most models remains incompletely resolved, as estimates of gene flow between taxa remain extremely rare. We studied the population genetics of two taxa of the Ostrinia genus--one feeding mainly on maize and the other on mugwort and hop--occurring in sympatry throughout France. The actual level of divergence of these taxa was unknown because the genetic structure of populations had been investigated over a limited geographical area and the magnitude of gene flow between populations had not been estimated. We used 11 microsatellite markers to investigate the genetic structure of populations throughout France and the extent of gene flow between the two Ostrinia taxa at several sites at which they are sympatric. We observed clear genetic differentiation between most populations collected on the typical respective hosts of each taxon. However, populations displaying intermediate allelic frequencies were found on hop plants in southern France. Individual assignments revealed that this result could be accounted for by the presence of both taxa on the same host. Gene flow, estimated by determining the proportion of hybrids detected, was low: probably<1% per generation, regardless of site. This indicates that the two Ostrinia taxa have reached a high level of genetic divergence and should be considered sibling species rather than host races.

Embolie graisseuse pulmonaire après liposuccion

A 24-year-old woman undergoes buttock's liposuction as an outpatient procedure. As she went back home, progressive dyspnea, respiratory distress and collapse developed. At hospital admission, she was dyspneic with thoracic oppression, tachycardia and anguish. Chest X-ray and thoracic CT scan suggested a pulmonary localisation of fat emboli. Symptomatic treatment allowed complete recovery. This report discusses diagnosis of fat emboli after liposuction as well as epidemiology and physiopathology.

Introgression in natural populations of bioindicators: a case study of Carabus splendens and Carabus punctatoauratus

The evolutionary importance of hybridization in wild plants and animals has become increasingly widely recognized in the last decade. In practical terms, hybridization provides an exceptionally tough set of problems for conservation biologists. We illustrate this in a case study of two Carabidae species widely used to evaluate the impact of human activities on biodiversity. These two species live in a complex mosaic of sympatry/allopatry and are known to hybridize in controlled conditions. Hybridization has not been quantified in natural populations to date due to the lack of a simple set of phenotypic traits for identifying hybrids. We thus screened for hybrids in natural populations, by multilocus genotyping at nine microsatellite loci. A high level of genetic differentiation between these two taxa was observed, as shown by allelic frequency distributions. Two Bayesian assignment procedures without obligatory pure taxon references were used to infer different classes of hybrids (F(1), F(2) and backcrosses) and mixture proportions between the two species. A low level of hybridization (F(1) genotypes) was observed in natural populations, contrasting with results obtained in controlled conditions. A high level of introgression was, however, detected at three of 12 sites, as revealed by the detection of backcrossed genotypes. This interspecific gene flow was detected in a limited zone of the common geographical range of the two species and was not related to the pattern of sympatry/allopatry. We then considered the origin and repercussions of this introgression, based on intraspecific genetic diversity and geographical structure.

Multiple paternity in Norway lobster (Nephrops norvegicus L.) assessed with microsatellite markers

We investigated genetic diversity and the mating system of the Norway lobster (Nephrops norvegicus) in a wild population off the Portuguese coast. Approximately 100 individuals were screened for 2 microsatellite loci. For 11 ovigerous lobsters both the female and a sample of her offspring (24 eggs) were genotyped. High genetic diversity was observed for the 2 markers in the population. Paternity within broods was analyzed by comparing multilocus genotypes of each egg with the corresponding mother, and the male parent contribution was then deduced. Multiple paternity was observed in 6 of the 11 broods studied. In those cases, 2 to 3 male parents were likely to have contributed to the fertilization of the eggs. When multiple paternity was involved, the comparative reproductive success of the male parents was quite even. This is the first report of multiple paternity in the Norway lobster. Comparisons with other taxa are presented, and consequences of multiple paternity are discussed.

Cardiac transplantation in a patient with hereditary hemochromatosis: role of adjunctive phlebotomy and erythropoietin

We describe the case of a 36-year-old woman with hereditary hemochromatosis (HH) resulting in end-stage cardiomyopathy and treated successfully with orthotopic cardiac transplantation. Before and after transplantation, the patient underwent aggressive treatment with frequent phlebotomy. We used erythropoietin concomitantly to maintain adequate hematocrit to support continued phlebotomy. We believe that aggressive use of phlebotomy provided the patient hemodynamic benefit and hastened the return of endocrine function post-transplantation. We also believe that the patient's history of high-dose vitamin C usage may have accelerated iron deposition in the heart and other vital organs.

Comparison of microsatellites and amplified fragment length polymorphism markers for parentage analysis

This study compares the properties of dominant markers, such as amplified fragment length polymorphisms (AFLPs), with those of codominant multiallelic markers, such as microsatellites, in reconstructing parentage. These two types of markers were used to search for both parents of an individual without prior knowledge of their relationships, by calculating likelihood ratios based on genotypic data, including mistyping. Experimental data on 89 oak trees genotyped for six microsatellite markers and 159 polymorphic AFLP loci were used as a starting point for simulations and tests. Both sets of markers produced high exclusion probabilities, and among dominant markers those with dominant allele frequencies in the range 0.1-0.4 were more informative. Such codominant and dominant markers can be used to construct powerful statistical tests to decide whether a genotyped individual (or two individuals) can be considered as the true parent (or parent pair). Gene flow from outside the study stand (GFO), inferred from parentage analysis with microsatellites, overestimated the true GFO, whereas with AFLPs it was underestimated. As expected, dominant markers are less efficient than codominant markers for achieving this, but can still be used with good confidence, especially when loci are deliberately selected according to their allele frequencies.

Oncostatin M suppresses EGF-mediated protein tyrosine phosphorylation in breast cancer cells

The effect of oncostatin M (OM) on epidermal growth factor (EGF)-mediated protein tyrosine phosphorylation in an infiltrating ductal breast carcinoma cell line, H3922, was investigated by Western blot analysis. Pretreatment of H3922 cells with OM for 72 h suppressed EGF-stimulated protein tyrosine phosphorylation signals by 77%. Interestingly, pretreatment with OM for 6 or 48 h had little effect on these signals. EGF-mediated tyrosine phosphorylation of EGF receptor (EGFR) was suppressed by 55% in 72-h OM pretreated H3922 cells. No reduction in EGFR protein expression was detected in these cells. Flow cytometric analysis verified that OM does not suppress EGFR expression. The effect of OM could not be attributed to induction of protein tyrosine phosphatases. An H3922 subclone cell line, designated H3922-8, was found to exhibit no proliferative response to treatment with EGF. However, EGF-mediated protein tyrosine phosphorylation was detected in these cells. Radioligand EGF binding studies comparing H3922 to H3922-8 cells indicated that the clonal cells apparently lack high affinity EGF receptors. The mechanism by which OM suppresses EGF-mediated tyrosine phosphorylation has not been completely characterized. However, the suppressive effect occurs regardless of whether the cells are acutely responsive (H3922) or virtually unresponsive (H3922-8) to EGF stimulation of cell growth.

A 3-basepair deletion in repeat 1 of the LDL receptor promoter reduces transcriptional activity in a South African Pedi

We have examined a naturally occurring mutation in the promoter region of the low density lipoprotein receptor (LDLR) gene of a South African Black patient with a clinical diagnosis of familial hypercholesterolemia (FH). The mutation constitutes a 3-bp deletion at nucleotide position -92 (FH Pedi-2) in the distal Sp1 binding site in repeat 1 of the LDLR promoter. The patient carries a second mutant LDLR allele containing a 1-bp deletion in exon 2 (FH Pedi-1) that gives rise to a frameshift mutation. Consistent with low receptor activity previously observed in cultured fibroblasts from the patient (5-15%), the rate of LDL receptor synthesis was markedly reduced to less than 20% of normal. DNase I footprint analysis indicated that the -92 mutation abolished binding of Sp1 to repeat 1 in the LDLR promoter. Transcription studies in transfected cells using normal and mutant promoter fragments linked to a luciferase reporter gene demonstrated that the promoter fragment containing the -92 mutation had approximately 10% of normal promoter activity. These findings indicate that the distal Sp1 binding site is essential for maximal activity of the normal intact LDLR promoter.

Novel mechanism of transcriptional activation of hepatic LDL receptor by oncostatin M

In this paper we describe a sterol-independent regulation of low density, lipoprotein receptor (LDLR) transcription by the cytokine oncostatin M (OM) in HepG2 cells. We show that OM-induced expression is independent of cholesterol regulation and occurs at the transcriptional level. To elucidate regulatory mechanism(s), we constructed a luciferase reporter system comprising either the native LDLR promoter including repeats 1, 2, and 3, or a synthetic promoter vector containing repeats 2+3 only, allowing us to directly examine OM effects on individual elements. Specific mutants in repeats 1, 2, and 3 were made to facilitate the mapping of the OM effect on the promoter. Wildtype and mutant constructs were assayed for cholesterol and OM regulation. The results show that mutation within the core SRE-1 element of repeat 2 totally abolished cholesterol regulation but had no effect on OM inducibility. More interesting, a mutation within repeat 1 reduced basal transcription activity to 10% of the native promoter, but OM induction was unaltered. However, the identical mutation engineered in repeat 3 significantly decreased OM induction of LDLR promoter activity. These results suggest a novel regulatory role for the repeat 3 element in LDLR transcription.

Identification and characterization of (GA/CT)n-microsatellite loci from Quercus petraea

In this study a size selected genomic library from Quercus petraea was screened for (GA/CT)n-microsatellite sequences. The resulting loci were analysed by PCR for their usefulness as molecular markers in Q. petraea and Q. robur. 17 out of 52 tested primer pairs resulted in the amplification of a polymorphic single-locus pattern. The number of alleles found per locus varied from 6 to 16. Combining the genetic variation observed for the characterized loci provides a unique genotype for all the individuals tested. Using intraspecific controlled crosses of Q. robur trees Mendelian inheritance could be shown for five loci.