Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect.

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Close Kin Dyads Indicate Intergenerational Dispersal and Barriers.

AbstractThe movement of individuals through continuous space is typically constrained by dispersal ability and dispersal barriers. A range of approaches have been developed to investigate these. Kindisperse is a new approach that infers recent intergenerational dispersal (σ) from close kin dyads and appears particularly useful for investigating taxa that are difficult to observe individually. This study, focusing on the mosquito Aedes aegypti, shows how the same close kin data can also be used for barrier detection. We empirically demonstrate this new extension of the method using genome-wide sequence data from 266 Ae. aegypti. First, we use the spatial distribution of full-sib dyads collected within one generation to infer past movements of ovipositing female mosquitoes. These dyads indicated the relative barrier strengths of two roads and performed favorably against alternative genetic methods for detecting barriers. We then use Kindisperse to quantify recent intergenerational dispersal (σ=81.5-197.1 m generation-1/2) from the difference in variance between the sib and the first cousin spatial distributions and, from this, estimate effective population density (ρ=833-4,864 km-2). Dispersal estimates showed general agreement with those from mark-release-recapture studies. Barriers, σ, ρ, and neighborhood size (331-526) can inform forthcoming releases of dengue-suppressing Wolbachia bacteria into this mosquito population.

Testing for adaptive changes linked to range expansion following a single introduction of the fall webworm.

Adaptive evolution following colonization can affect the impact of invasive species. The fall webworm (FWW) invaded China 40 years ago through a single introduction event involving a severe bottleneck and subsequently diverged into two genetic groups. The well-recorded invasion history of FWW, coupled with a clear pattern of genetic divergence, provides an opportunity to investigate whether there is any sign of adaptive evolution following the invasion. Based on genome-wide SNPs, we identified genetically separated western and eastern groups of FWW and correlated spatial variation in SNPs with geographical and climatic factors. Geographical factors explained a similar proportion of the genetic variation across all populations compared with climatic factors. However, when the two population groups were analysed separately, environmental factors explained more variation than geographical factors. SNP outliers in populations of the western group had relatively stronger response to precipitation than temperature-related variables. Functional annotation of SNP outliers identified genes associated with insect cuticle protein potentially related to desiccation adaptation in the western group and genes associated with lipase biosynthesis potentially related to temperature adaptation in the eastern group. Our study suggests that invasive species may maintain the evolutionary potential to adapt to heterogeneous environments despite a single invasion event. The molecular data suggest that quantitative trait comparisons across environments would be worthwhile.

Local climate adaptation and gene flow in the native range of two co-occurring fruit moths with contrasting invasiveness.

Invasive species pose increasing threats to global biodiversity and ecosystems. While previous studies have characterized successful invaders based on ecological traits, characteristics related to evolutionary processes have rarely been investigated. Here we compared gene flow and local adaptation using demographic analyses and outlier tests in two co-occurring moth pests across their common native range of China, one of which (the peach fruit moth, Carposina sasakii) has maintained its native distribution, while the other (the oriental fruit moth, Grapholita molesta) has expanded its range globally during the past century. We found that both species showed a pattern of genetic differentiation and an evolutionary history consistent with a common southwestern origin and northward expansion in their native range. However, for the noninvasive species, genetic differentiation was closely aligned with the environment, and there was a relatively low level of gene flow, whereas in the invasive species, genetic differentiation was associated with geography. Genome scans indicated stronger patterns of climate-associated loci in the noninvasive species. While strong local adaptation and reduced gene flow across its native range may have decreased the invasiveness of C. sasakii, this requires further validation with additional comparisons of invasive and noninvasive species across their native range.

Comparative mitogenomics and phylogenetics of the stinging wasps (Hymenoptera: Aculeata).

The stinging wasps (Hymenoptera: Aculeata) include diverse groups such as vespid wasps, ants and bees. Phylogenetic relationships among major lineages of stinging wasps have been inferred from molecular and morphological data. However, the genomic features of the mitochondrial genomes and their phylogenetic utility remain to be explored. In this study, we determined 23 mitochondrial genomes from the Aculeata. Four Mutillidae species showed relatively low A + T content compared to other species of the Aculeata (69.7%-77.4%). Eleven out of 44 species, mainly from the Chrysididae and the Pompilidae, showed reversals of GC skews. Gene rearrangements occurred across the species. Patterns of tRNA rearrangement were conserved in some groups, including the Chrysididae, Bethylidae, Pompilidae, Scolioidea and Vespoidea. Rearrangement of protein-coding genes were found in 12 out of 44 species of the Aculeata, including all four species from the Chrysididae, both species from the Bethylidae, one species from the Dryinidae, all three Scolioidea species and two Apoidea species. Phylogenetic inference showed a long branch in species with unusual genomic features, such as in the Mutillidae and Bethylidae. By excluding these species, we found paraphyly of the Chrysidoidea and a sister group relationship between the Formicoidea and Vespoidea. These results improve our understanding of the evolution of mitochondrial genomes in the Aculeata and, in general, the evolution across this subclade.

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China.

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise FST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.

Field-evolved resistance and cross-resistance of the two-spotted spider mite, Tetranychus urticae, to bifenazate, cyenopyrafen and SYP-9625.

The acaricide bifenazate acts as complex III inhibitor whereas cyenopyrafen and SYP-9625 act as complex II inhibitors. All these acaricides are commonly used to control two-spotted spider mite (TSSM), Tetranychus urticae Koch. We examined field-evolved and laboratory-selected resistance of TSSM to these three acaricides and determined cross-resistance among them. Six field populations of TSSM showed low levels of resistance to bifenazate with resistance ratios ranging from 2.20 to 10.65 compared to a susceptible strain. SYP-9625, structurally similar to cyenopyrafen, showed slightly higher activity to TSSMs but significant cross-resistance in both field populations and a laboratory-selected strain by SYP-9625. However, low levels of resistance to these two chemicals were found in field populations even when used for short time periods. Cross-resistance was not found between bifenazate and Complex II inhibitors, cyenopyrafen and SYP-9625, in both field populations and the laboratory-selected strain. Field-evolved resistance of TSSM to the tested acaricides is still low and should be delayed by the implementation of resistance management practices. Cross-resistance between cyenopyrafen and SYP-9625 is obvious, so they should not be used together in resistance management strategies based on mode of action rotation.

Multiple refugia from penultimate glaciations in East Asia demonstrated by phylogeography and ecological modelling of an insect pest.

BACKGROUND: Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China. RESULTS: Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage. CONCLUSIONS: Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

Gene arrangement and sequence of mitochondrial genomes yield insights into the phylogeny and evolution of bees and sphecid wasps (Hymenoptera: Apoidea).

The Apoidea represent a large and common superfamily of the Hymenoptera including the bees and sphecid wasps. A robust phylogenetic tree is essential to understanding the diversity, taxonomy and evolution of the Apoidea. In this study, features of apoid mitochondrial genomes were used to reconstruct phylogenetic relationships. Twelve apoid mitochondrial genomes were newly sequenced, representing six families and nine subfamilies. Gene rearrangement events have occurred in all apoid mitochondrial genomes sequenced to date. Sphecid wasps have both tRNA and protein-coding gene rearrangements in 5 of 8 species. In bees, the only rearranged genes are tRNAs; long-tongued bees (Apidae + Megachilidae) are characterized by movement of trnA to the trnI-trnQ-trnM tRNA cluster. Phylogenetic analyses of mitochondrial gene sequences support the known paraphyly of sphecid wasps, with bees nested within this clade. The Ampulicidae is sister to the remaining Apoidea. Crabronidae is paraphyletic, split into Crabronidae s.s. and Philanthidae, with the latter group a sister clade to bees. The monophyletic bees are either classified into two clades, long-tongued bees (Apidae + Megachilidae) and short-tongued bees (Andrenidae + Halictidae + Colletidae + Melitidae), or three groups with the Melitidae sister to the other bees. Our study showed that both gene sequences and arrangements provide information on the phylogeny of apoid families.

Low vitamin D levels are associated with cognitive impairment in patients with Hashimoto thyroiditis.

BACKGROUND: Cognitive impairment is commonly observed in patients with Hashimoto thyroiditis (HT). Low levels of vitamin D have been correlated with cognitive impairment in non-HT population. We examined the association of vitamin D levels with cognitive impairment in patients with HT. METHODS: We recruited 194 patients with HT and 200 healthy volunteers. Levels of serum 25-hydroxyvitamin D (25(OH)D) were measured using a competitive protein-binding assay. Cognitive funtion was assessed using Montreal Cognitive Assessment score (MoCA). Subjects with a MoCA scores < 26 are considered as having mild cognitive impairment (MCI). Multivariate analysis was performed using logistic regression models. RESULTS: Fifty-five HT patients (28.4%) were diagnosed as having MCI. Patients with MCI had significantly lower 25(OH)D levels when compared with patients without MCI (33.9 ± 6.2 vs. 44.3 ± 9.6 nmol/L, P < 0.001). Significant differences in 25(OH)D quartiles of HT patients were observed between the patients with MCI and the patients without MCI (P < 0.001). In multivariate analyses, serum 25(OH)D levels (≤ 34.0 and ≥ 47.1 nmol/L) were significantly associated with cognitive impairment in patients with HT (OR 6.279, 95% CI 2.673-14.834, P < 0.001; OR 0.061, 95% CI 0.008-0.491, P = 0.009, respectively). CONCLUSION: Our results demonstrate an important association between serum vitamin D levels and cognitive impairment in patients with HT.

Toxicity of seven insecticides to different developmental stages of the whitefly Bemisia tabaci MED (Hemiptera: Aleyrodidae) in multiple field populations of China.

Chemical control is important in the management of the tobacco whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae). Susceptibility of B. tabaci to insecticides may vary among different developmental stages and geographical populations. In this study, we examined toxicity of seven commonly-used insecticides to B. tabaci MED in four field populations from China. Avermectin has high level of toxicity to all stages of B. tabaci MED in all four populations. Cyantraniliprole and sulfoxaflor have high toxicity to adults. Spirotetramat, cyantraniliprole and flonicamid have high toxicity to nymphs but not adults. Acetamiprid, cyantraniliprole and sulfoxaflor have high toxicity to eggs. However, the relative toxicity of B. tabaci MED to these chemicals varied across different populations, with little consistency in population differences across developmental stages. Our findings together with some instances where LC95 values were higher than field recommended dosages indicate field-evolved resistance to insecticides (such as thiamethoxam and sulfoxaflor) and stage-specific mechanisms that will influence effective control of B. tabaci MED by insecticides.

Preference and performance of the two-spotted spider mite Tetranychus urticae (Acari: Tetranychidae) on strawberry cultivars.

The two-spotted spider mite (TSSM), Tetranychus urticae Koch (Acari: Tetranychidae), is one of the most serious pests of strawberry worldwide. Understanding the preference of TSSM for particular cultivars of strawberry and performance on them helps identify host-plant resistance to this pest mite. In this study, we tested preference, developmental duration, fecundity and population levels of TSSM on 14 strawberry cultivars. TSSM showed strong preference for the Chinese cultivars of Yanxiang, Baixuegongzhu, and Jingtaoxiang. Development of TSSM on the cultivars varied from 32.32 to 36.82 days; it was longest on the cultivars Hongxiutianxiang and Baixuegongzhu, and shortest on Yanxiang, Jingzangxiang, and Darselect as well as on a wild variety (Wuye). TSSM had high fecundity on the cultivars Yanxiang, Taoxun, Hongxiutianxiang, Jingzangxiang, Albion and Baixuegongzhu as well as on Wuye, whereas egg production was lowest on Sweet Charlie, Portola, Akihime, and Benihoppe. After 28 days of plant infestation with 10 pairs of adults, the cultivars Yanxiang, Taoxun, Jingzangxiang, Jingtaoxiang, and Baixuegongzhu had the highest number of mites (> 1000 per plant), whereas mite numbers on Albion and Camarosa were low. The population size of TSSM was correlated with fecundity, but no correlation was found between other preference/performance measures. Our study suggests that a rapid increase of population size of TSSM on cultivars of strawberry is related to high fecundity, and also that there are substantial differences in preference and performance across cultivars.

Efficacy of carbon dioxide treatments for the control of the two-spotted spider mite, Tetranychus urticae, and treatment impact on plant seedlings.

To develop a new control method for the two-spotted spider mite (TSSM), Tetranychus urticae, we investigated the effect of controlled atmospheres of carbon dioxide (CO2) on TSSM mortality under different concentrations and treatment periods, and evaluated the impact of treatments on seedlings of five host plants of TSSM. Egg hatching rate of TSSM was reduced to 37.7, 5.4 or 0% after 24 h treatment involving concentrations of 16.7, 33.3 or 50%, respectively. Mobile stages (nymphs and adult) of TSSM were completely controlled after 24 h treatment at concentrations higher than 33.3%. After 4 h at concentrations of 33.3 or 50%, 1st-day survival rate for all mobile stages was 45.3 or 36.0%, respectively, whereas after 8 or 16 h treatments, all values were decreased to zero. Seedlings of four major host plants of TSSM (cucumber, eggplant, rape, green peppers) were damaged to varying degrees after 24 h at the three concentrations, but strawberry, another host plant, was not damaged. Cucumber suffered the most serious damage, resulting in wilting and death. In conclusion, controlled atmospheres of CO2 can kill TSSM, particularly at high concentrations and with long treatment times. It can be used to control TSSM on strawberry, but should be used cautiously on other host plants.

Patterns of genetic variation among geographic and host-plant associated populations of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae).

BACKGROUND: Populations of herbivorous insects may become genetically differentiated because of local adaptation to different hosts and climates as well as historical processes, and further genetic divergence may occur following the development of reproductive isolation among populations. Here we investigate the population genetic structure of the orchard pest peach fruit moth (PFM) Carposina sasakii (Lepidoptera: Carposinidae) in China, which shows distinct biological differences when characterized from different host plants. Genetic diversity and genetic structure were assessed among populations from seven plant hosts and nine regions using 19 microsatellite loci and a mitochondrial sequence. RESULTS: Strong genetic differentiation was found among geographical populations representing distinct geographical regions, but not in host-associated populations collected from the same area. Mantel tests based on microsatellite loci indicated an association between genetic differentiation and geographical distance, and to a lesser extent environmental differentiation. Approximate Bayesian Computation analyses supported the scenario that PFM likely originated from a southern area and dispersed northwards before the last glacial maximum during the Quaternary. CONCLUSIONS: Our analyses suggested a strong impact of geographical barriers and historical events rather than host plants on the genetic structure of the PFM; however, uncharacterized environmental factors and host plants may also play a role. Studies on adaptive shifts in this moth should take into account geographical and historical factors.

Development and Characterization of Novel Microsatellite Markers for the Peach Fruit Moth Carposina sasakii (Lepidoptera: Carposinidae) Using Next-Generation Sequencing.

The peach fruit moth Carposina sasakii is an economically important pest on dozens of fruits from Rosaceae and Rhamnaceae in Northeast Asia. We developed novel microsatellite markers for C. sasakii from randomly sequenced regions of the genome using next-generation sequencing. In total, 95,153 microsatellite markers were isolated from 4.70 GB genomic sequences. Thirty-five polymorphic markers were developed by assessing in 63 individuals from two geographical populations. The allele numbers ranged from 2 to 9 with an average value of 4.60 per locus, while the polymorphism information content ranged from 0.075 to 0.696 with an average value of 0.407. Furthermore, the observed and expected heterozygosity varied from 0.000 to 0.677 and 0.062 to 0.771, respectively. The microsatellites developed provide abundant molecular markers for investigating genetic structure, genetic diversity, and existence of host-plant associated biotypes of C. sasakii.

Population genetic structure and approximate Bayesian computation analyses reveal the southern origin and northward dispersal of the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae) in its native range.

The oriental fruit moth (OFM) Grapholita molesta is one of the most destructive orchard pests. Assumed to be native to China, the moth is now distributed throughout the world. However, the evolutionary history of this moth in its native range remains unknown. In this study, we explored the population genetic structure, dispersal routes and demographic history of the OFM in China and South Korea based on mitochondrial genes and microsatellite loci. The Mantel test indicated a significant correlation between genetic distance and geographical distance in the populations. Bayesian analysis of population genetic structure (baps) identified four nested clusters, while the geneland analysis inferred five genetic groups with spatial discontinuities. Based on the approximate Bayesian computation approach, we found that the OFM was originated from southern China near the Shilin area of Yunnan Province. The early divergence and dispersal of this moth was dated to the Penultimate glaciation of Pleistocene. Further dispersal from southern to northern region of China occurred before the last glacial maximum, while the expansion of population size in the derived populations in northern region of China occurred after the last glacial maximum. Our results indicated that the current distribution and structure of the OFM were complicatedly influenced by climatic and geological events and human activities of cultivation and wide dissemination of peach in ancient China. We provide an example on revealing the origin and dispersal history of an agricultural pest insect in its native range as well as the underlying factors.

The effects of fluid resuscitation according to PiCCO on the early stage of severe acute pancreatitis.

OBJECTIVES: To evaluate the therapeutic effect of early fluid resuscitation under the guidance of Pulse indicator Continuous Cardiac Output (PiCCO) on patients with severe acute pancreatitis (SAP). METHODS: Clinical data of 18 SAP patients (the study group), who had undergone fluid resuscitation under the guidance of PiCCO from October 2011 to October 2013, were analyzed prospectively. Clinical data of 25 cases (control group) who had undergone fluid resuscitation without the guidance of PiCCO from January 2009 to September 2011 were collected. Then, retrospective and prospective case-control study was carried out. RESULTS: During the first 6 h, 0-24 h, 24-48 h, and 0-72 h of admission, the study group received more volume of fluid than the control group. There were significantly faster decline of APACHE II score and the value of blood lactate in study group, as well as the length of ICU stay and the proportion of renal failure at 72 h of admission. According to the 2012 Atlanta classification, six cases in study group turned into moderate SAP (33.30%), significantly higher than the control group (8.00%) (p = 0.0049). The volume of fluid infusion and clinical parameters were linearly relative. CONCLUSIONS: The PiCCO device may be a useful adjunct for fluid resuscitation monitoring in patients with SAP in the early stage. Early fluid resuscitation under the guidance of PiCCO can improve tissue perfusion, reduce the SIRS persistence time and the length of ICU stay. This program did not increase the risk of respiratory failure and influence the mortality.

Nuclear and mitochondrial genomes of the plum fruit moth Grapholita funebrana.

The plum fruit moth Grapholita funebrana (Tortricidae, Lepidoptera) is an important pest of many wild and cultivated stone fruits and other plants in the family Rosaceae. Here, we assembled its nuclear and mitochondrial genomes using Illumina, Nanopore, and Hi-C sequencing technologies. The nuclear genome size is 570.9 Mb, with a repeat rate of 51.28%, and a BUCSO completeness of 97.7%. The karyotype for males is 2n = 56. We identified 17,979 protein-coding genes, 5,643 tRNAs, and 94 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and 2 rRNA. These genomes provide resources to understand the genetics, ecology, and genome evolution of the tortricid moths.

Nuclear and mitochondrial genomes of Polylopha cassiicola: the first assembly in Chlidanotinae (Tortricidae).

Tortricidae is one of the largest families in Lepidoptera, including subfamilies of Tortricinae, Olethreutinae, and Chlidanotinae. Here, we assembled the gap-free genome for the subfamily Chlidanotinae using Illumina, Nanopore, and Hi-C sequencing from Polylopha cassiicola, a pest of camphor trees in southern China. The nuclear genome is 302.03 Mb in size, with 36.82% of repeats and 98.4% of BUCSO completeness. The karyotype is 2n = 44 for males. We identified 15412 protein-coding genes, 1052 tRNAs, and 67 rRNAs. We also determined the mitochondrial genome of this species and annotated 13 protein-coding genes, 22 tRNAs, and one rRNA. These high-quality genomes provide valuable information for studying phylogeny, karyotypic evolution, and adaptive evolution of tortricid moths.

Chromosome-level genome assembly of the western flower thrips Frankliniella occidentalis.

The western flower thrips Frankliniella occidentalis (Thysanoptera: Thripidae) is a global invasive species that causes increasing damage by direct feeding on crops and transmission of plant viruses. Here, we assemble a previously published scaffold-level genome into a chromosomal level using Hi-C sequencing technology. The assembled genome has a size of 302.58 Mb, with a contig N50 of 1533 bp, scaffold N50 of 19.071 Mb, and BUSCO completeness of 97.8%. All contigs are anchored on 15 chromosomes. A total of 16,312 protein-coding genes are annotated in the genome with a BUSCO completeness of 95.2%. The genome contains 492 non-coding RNA, and 0.41% of interspersed repeats. In conclusion, this high-quality genome provides a convenient and high-quality resource for understanding the ecology, genetics, and evolution of thrips.