Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex.

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

What is "many-body" dispersion and should I worry about it?

Inclusion of dispersion effects in density-functional calculations is now standard practice in computational chemistry. In many dispersion models, the dispersion energy is written as a sum of pairwise atomic interactions consisting of a damped asymptotic expansion from perturbation theory. There has been much recent attention drawn to the importance of "many-body" dispersion effects, which by their name imply limitations with a pairwise atomic expansion. In this perspective, we clarify what is meant by many-body dispersion, as this term has previously referred to two very different physical phenomena, here classified as electronic and atomic many-body effects. Atomic many-body effects refer to the terms in the perturbation-theory expansion of the dispersion energy involving more than two atoms, the leading contribution being the Axilrod-Teller-Muto three-body term. Conversely, electronic many-body effects refer to changes in the dispersion coefficients of the pairwise terms induced by the atomic environment. Regardless of their nature, many-body effects cause pairwise non-additivity in the dispersion energy, such that the dispersion energy of a system does not equal the sum of the dispersion energies of its atomic pairs taken in isolation. A series of examples using the exchange-hole dipole moment (XDM) method are presented to assess the relative importance of electronic and atomic many-body effects on the dispersion energy. Electronic many-body effects can result in variation in the leading-order C6 dispersion coefficients by as much as 50%; hence, their inclusion is critical for good performance of a pairwise asymptotic dispersion correction. Conversely, atomic many-body effects represent less than 1% of the total dispersion energy and are much less significant than higher-order (C8 and C10) pairwise terms. Their importance has been previously overestimated through empirical fitting, where they can offset underlying errors stemming either from neglect of higher-order pairwise terms or from the base density functional.

Incongruence between molecules and morphology: A seven-gene phylogeny of Dacini fruit flies paves the way for reclassification (Diptera: Tephritidae).

Molecular and morphological research often suggest conflicting results. Selective pressure on certain morphologies can confound understanding of evolutionary relationships. Dacini is one of the most diverse tribes of tephritid flies and contains many economically important pest species. Their black and yellow patterned body markings are presumed to act as wasp mimicry, and the characters separating species and groups are limited and in some cases phenotypically plastic. The traditional taxonomy of the tribe is controversial because groupings are based on unique combinations of morphological characters without the use of cladistic methods, though recent phylogenetic and taxonomic analyses have resulted in significant changes to their taxonomy. The monophyly of the three largest genera in the tribe has been tested with only small numbers of representatives per genus and a limited number of genes. To further understand the taxonomy and evolution of Dacini we sequenced seven genes from 167 Dacini species and five dipteran outgroups to construct a robust phylogeny and test phylogenetic relationships between genera, subgenera, and species complexes. Our phylogeny confirms the monophyly of Dacus, Bactrocera, and Zeugodacus. However, most groups below the genus level are not monophyletic, and only through further revision will we be able to understand their evolution and clarify the taxonomy within this tribe.

Pervasive Delocalisation Error Causes Spurious Proton Transfer in Organic Acid-Base Co-Crystals.

Dispersion-corrected density-functional theory (DFT-D) methods have become the workhorse of many computational protocols for molecular crystal structure prediction due to their efficiency and convenience. However, certain limitations of DFT, such as delocalisation error, are often overlooked or are too expensive to remedy in solid-state applications. This error can lead to artificial stabilisation of charge transfer and, in this work, it is found to affect the correct identification of the protonation site in multicomponent acid-base crystals. As such, commonly used DFT-D methods cannot be applied with any reliability to the study of acid-base co-crystals or salts, while hybrid functionals remain too restrictive for routine use. This presents an impetus for the development of new functionals with reduced delocalisation error for solid-state applications; the structures studied herein constitute an excellent benchmark for this purpose.

Competing fragmentation processes of ß-substituted propanoate ions upon collision-induced dissociation.

RATIONALE: When subjected to collisional activation, gas-phase carboxylate ions typically undergo decarboxylation. However, alternative fragmentation processes dominate when the carboxylate group is located within certain structural motifs. In this work, the fragmentation processes of ß-substituted carboxylate ions are characterized to improve correlations between reactivity and structure. METHODS: Mass spectra were collected using both ion trap and triple quadrupole mass spectrometers operating in the negative ion mode; collision-induced dissociation (CID) of ions was used to study the relationship between product ions and the structures of their precursor ions. Quantum mechanical computations were performed on a full range of reaction geometries at the MP2/6-311++G(2d,p)//B3LYP/6-31++G(2d,p) level of theory. RESULTS: For a series of ß-substituted carboxylate ions, a product ion corresponding to the anion of the ß-substituent was obtained upon CID. Detailed computations indicated that decarboxylative elimination and at least one other fragmentation mechanism had feasible energetics for the formation of substituent anions differing in their gas-phase basicities. Predicted energetics for anti- and synperiplanar alignments in the transition structures for decarboxylative elimination correlated with the positions of crossover points in breakdown curves acquired for conformationally constrained ions. CONCLUSIONS: The feasibility of more than one mechanism was established for the fragmentation of ß-substituted propanoates. The contribution of each mechanistic pathway to the formation of the substituent anion was influenced by structural variations and conformational constraints, but mostly depended on the nature of the substituent. Copyright © 2016 John Wiley & Sons, Ltd.

Torquoselectivity in the Nazarov reactions of allenyl vinyl ketones.

Nazarov reactions mediated by BF3-etherate of a series of carbon-substituted allenyl vinyl ketones provided intermediates in which substituents on the termini of the allenes had rotated away from the vinyl moieties, and these intermediates were trapped by (4 + 3)-cyclizations. A computational examination of the torquoselectivity of these Nazarov reactions confirmed a kinetic preference for the observed isomers and pointed to steric interactions and the degree of allene deformation as significant factors in determining the torquoselectivity. The study also suggested that the high proportion of one geometrical isomer in the Nazarov products might also be due to some preferential trapping of the major Nazarov intermediate.

Phenyl group participation in rearrangements during collision-induced dissociation of deprotonated phenoxyacetic acid.

RATIONALE: The identification of trace constituents in biological and environmental samples is frequently based on the fragmentation patterns resulting from the collision-induced dissociation (CID) of gas-phase ions. Credible mechanistic characterization of fragmentation processes, including rearrangements, is required to make reliable assignments for structures of precursor and product ions. METHODS: Mass spectra were collected using both ion trap and triple quadrupole mass spectrometers operating in the negative ion mode. Precursor ion scans and CID of ions generated in-source were used to establish precursor-product ion relationships. Density functional theory (DFT) computations were performed at the MP2/6-311++G(2d,p)//B3LYP/6-31++G(2d,p) level of theory. RESULTS: Product ions at m/z 93 and 107 obtained upon CID of phenoxyacetate were attributed to phenoxide and o-methylphenoxide, respectively. An isotopic labeling experiment and computations showed that the phenoxide ion was formed by intramolecular displacement with formation of an α-lactone and also by a Smiles rearrangement. Rearrangement of phenoxyacetate via the ion-neutral complex formed in the α-lactone displacement pathway gave the isomeric o-hydroxyphenylacetate ion which yielded o-methylphenoxide upon decarboxylation. Computations provided feasible energetics for these pathways. CONCLUSIONS: Previously unrecognized and energetically favorable rearrangements during the collision-induced fragmentation of phenoxyacetate have been characterized using isotopic labeling and DFT computations. Notably, the phenyl substituent plays an indispensable role in each rearrangement process resulting in multiple pathways for the fragmentation of phenoxyacetate.

Density Functional Theory Study of BF3-Mediated Additions of Enols and [(Trimethylsilyl)oxy]alkenes to an Oxyallyl Cation: Homologous Mukaiyama Reactions.

The addition of enols and [(trimethylsilyl)oxy]alkenes, bearing methyl substituents at various positions, to a cyclic, BF3-complexed oxyallyl cation has been studied at the M06/6-311G(d)//B3LYP/6-31G(d) level of theory. The reactions with the [(trimethylsilyl)oxy]alkenes are homologous Mukaiyama reactions, which have not been examined computationally previously. In most instances a number of transition states were located, and the difference in energy between these transition states was not large, which pointed to low levels of diastereoselectivity in the reactions of the oxyallyl cation model compound. The lowest energy transition states were those with a synclinal geometry in which the alkene was positioned over the cyclic oxyallyl cation, and the relative orientation of the alkene and the oxyallyl cation was rationalized in terms of stabilizing intermolecular interactions, revealed by NBO analysis, between one or more fluorines of the complexed BF3 and hydrogens on the alkene moiety, and between the oxygen on the alkene and the π-system of the oxyallyl cation. Because, in most instances with these simple models, two or more transition states of relatively low energy were located, predictions of diastereoselectivity in more complex examples that are based on simple models cannot be recommended.

Genetic Diversity of Bactrocera dorsalis (Diptera: Tephritidae) on the Hawaiian Islands: Implications for an Introduction Pathway Into California.

Population genetic diversity of the oriental fruit fly, Bactrocera dorsalis (Hendel), on the Hawaiian islands of Oahu, Maui, Kauai, and Hawaii (the Big Island) was estimated using DNA sequences of the mitochondrial cytochrome c oxidase subunit I gene. In total, 932 flies representing 36 sampled sites across the four islands were sequenced for a 1,500-bp fragment of the gene named the C1500 marker. Genetic variation was low on the Hawaiian Islands with >96% of flies having just two haplotypes: C1500-Haplotype 1 (63.2%) or C1500-Haplotype 2 (33.3%). The other 33 flies (3.5%) had haplotypes similar to the two dominant haplotypes. No population structure was detected among the islands or within islands. The two haplotypes were present at similar frequencies at each sample site, suggesting that flies on the various islands can be considered one population. Comparison of the Hawaiian data set to DNA sequences of 165 flies from outbreaks in California between 2006 and 2012 indicates that a single-source introduction pathway of Hawaiian origin cannot explain many of the flies in California. Hawaii, however, could not be excluded as a maternal source for 69 flies. There was no clear geographic association for Hawaiian or non-Hawaiian haplotypes in the Bay Area or Los Angeles Basin over time. This suggests that California experienced multiple, independent introductions from different sources.

Towards a better future for DNA barcoding: Evaluating monophyly- and distance-based species identification using COI gene fragments of Dacini fruit flies.

The utility of a universal DNA 'barcode' fragment (658 base pairs of the Cytochrome C Oxidase I [COI] gene) has been established as a useful tool for species identification, and widely criticized as one for understanding the evolutionary history of a group. Large amounts of COI sequence data have been produced that hold promise for rapid species identification, for example, for biosecurity. The fruit fly tribe Dacini holds about a thousand species, of which 80 are pests of economic concern. We generated a COI reference library for 265 species of Dacini containing 5601 sequences that span most of the COI gene using circular consensus sequencing. We compared distance metrics versus monophyly assessments for species identification and although we found a 'soft' barcode gap around 2% pairwise distance, the exceptions to this rule dictate that a monophyly assessment is the only reliable method for species identification. We found that all fragments regularly used for Dacini fruit fly identification >450 base pairs long provide similar resolution. 11.3% of the species in our dataset were non-monophyletic in a COI tree, which is mostly due to species complexes. We conclude with recommendations for the future generation and use of COI libraries. We revise the generic assignment of Dacus transversus stat. rev. Hardy 1982, and Dacus perpusillus stat. rev. Drew 1971 and we establish Dacus maculipterus White 1998 syn. nov. as a junior synonym of Dacus satanas Liang et al. 1993.

NMR metabolomics analysis of the effects of 5-lipoxygenase inhibitors on metabolism in glioblastomas.

Changes across metabolic networks are emerging as an integral part of cancer development and progression. Increasing comprehension of the importance of metabolic processes as well as metabolites in cancer is stimulating exploration of novel, targeted treatment options. Arachidonic acid (AA) is a major component of phospholipids. Through the cascade catalyzed by cyclooxygenases and lipoxygenases, AA is also a precursor to cellular signaling molecules as well as molecules associated with a variety of diseases including cancer. 5-Lipoxygenase catalyzes the transformation of AA into leukotrienes (LT), important mediators of inflammation. High-throughput analysis of metabolic profiles was used to investigate the response of glioblastoma cell lines to treatment with 5-lipoxygenase inhibitors. Metabolic profiling of cells following drug treatment provides valuable information about the response and metabolic alterations induced by the drug action and give an indication of both on-target and off-target effects of drugs. Four different 5-lipoxygenase inhibitors and antioxidants were tested including zileuton, caffeic acid, and its analogues caffeic acid phenethyl ester and caffeic acid cyclohexethyl ester. A NMR approach identified metabolic signatures resulting from application of these compounds to glioblastoma cell lines, and metabolic data were used to develop a better understanding of the mode of action of these inhibitors.

Synthesis and antiradical/antioxidant activities of caffeic acid phenethyl ester and its related propionic, acetic, and benzoic acid analogues.

Caffeic acid phenethyl ester (CAPE) is a bioactive component isolated from propolis. A series of CAPE analogues was synthesized and their antiradical/antioxidant effects analyzed. The effect of the presence of the double bond and of the conjugated system on the antioxidant effect is evaluated with the analogues obtained from 3-(3,4-dihydroxyphenyl) propanoic acid. Those obtained from 2-(3,4-dihydroxyphenyl) acetic acid and 3,4-dihydroxybenzoic acid allow the evaluation of the effect of the presence of two carbons between the carbonyl and aromatic system.

Controlling anisotropic properties by manipulating the orientation of chiral small molecules.

Chiral π-conjugated molecules bring new functionality to technological applications and represent an exciting, rapidly expanding area of research. Their functional properties, such as the absorption and emission of circularly polarized light or the transport of spin-polarized electrons, are highly anisotropic. As a result, the orientation of chiral molecules critically determines the functionality and efficiency of chiral devices. Here we present a strategy to control the orientation of a small chiral molecule (2,2'-dicyano[6]helicene) by the use of organic and inorganic templating layers. Such templating layers can either force 2,2'-dicyano[6]helicene to adopt a face-on orientation and self-assemble into upright supramolecular columns oriented with their helical axis perpendicular to the substrate, or an edge-on orientation with parallel-lying supramolecular columns. Through such control, we show that low- and high-energy chiroptical responses can be independently 'turned on' or 'turned off'. The templating methodologies described here provide a simple way to engineer orientational control and, by association, anisotropic functional properties of chiral molecular systems for a range of emerging technologies.

Novel lures and COI sequences reveal cryptic new species of Bactrocera fruit flies in the Solomon Islands (Diptera, Tephritidae, Dacini).

Results from a snap-shot survey of Dacine fruit flies carried out on three of the Solomon Islands in April 2018 are reported. Using traps baited with the male lures cue-lure, methyl eugenol, and zingerone, 30 of the 48 species previously known to occur in the Solomon Islands were collected. Six species are newly described here: Bactroceraallodistincta sp. nov., B.geminosimulata sp. nov., B.kolombangarae sp. nov., B.quasienochra sp. nov., B.tsatsiai sp. nov., and B.vargasi sp. nov., all authored by Leblanc & Doorenweerd. An illustrated key to the 54 species now known to be present in the country is provided.

A Universal Force Field for Materials, Periodic GFN-FF: Implementation and Examination.

In this study, the adaption of the recently published molecular GFN-FF for periodic boundary conditions (pGFN-FF) is described through the use of neighbor lists combined with appropriate charge sums to handle any dimensionality from 1D polymers to 2D surfaces and 3D solids. Numerical integration over the Brillouin zone for the calculation of π bond orders of periodic fragments is also included. Aside from adapting the GFN-FF method to handle periodicity, improvements to the method are proposed in regard to the calculation of topological charges through the inclusion of a screened Coulomb term that leads to more physical charges and avoids a number of pathological cases. Short-range damping of three-body dispersion is also included to avoid collapse of some structures. Analytic second derivatives are also formulated with respect to both Cartesian and strain variables, including prescreening of terms to accelerate the dispersion/coordination number contribution to the Hessian. The modified pGFN-FF scheme is then applied to a wide range of different materials in order to examine how well this universal model performs.

Evaluation of methyl eugenol and cue-lure traps with solid lure and insecticide dispensers for fruit fly monitoring and male annihilation in the Hawaii Areawide Pest Management Program.

Methyl eugenol (ME) and cue-lure (C-L) traps with solid lure dispensers were deployed in areas with low and high populations of oriental fruit fly, Bactrocera dorsalis (Hendel), and melon fly, Bactrocera cucurbitae (Coquillett), respectively. In low-density areas, standard Jackson traps or Hawaii Fruit Fly Areawide Pest Management (AWPM) traps with FT Mallet ME wafers impregnated with dimethyl dichloro-vinyl phosphate (DDVP) or AWPM traps with Scentry ME cones and vapor tape performed equally as well as standard Jackson traps with liquid ME/C-L and naled. Standard Jackson traps or AWPM traps with FT Mallet C-L wafers impregnated with DDVP or AWPM traps with Scentry C-L plugs with vapor tape performed equally as well as standard Jackson traps with a lure-naled solution. In high density areas, captures with traps containing FT Mallet wafers (ME and C-L) outperformed AWPM traps with Scentry cones and plugs (ME and C-L) with DDVP insecticidal strips over a 6-mo period. Captures of B. dorsalis and B. cucurbitae with wafers containing both ME and raspberry ketone (FT Mallet MC) were equivalent to those containing separate lures. From a worker safety and convenience standpoint, FT Mallet ME and C-L wafers with DDVP or Scentry plugs, with or without DDVP vapor tape, are more convenient and safer to handle than standard liquid insecticide formulations used for monitoring and male annihilation programs in Hawaii, and for detections traps used on the U.S. mainland. Furthermore, the FT Mallet MC wafer might be used in a single trap in place of two separate traps for detection of both ME and C-L responding fruit flies.

Highly variable COI haplotype diversity between three species of invasive pest fruit fly reflects remarkably incongruent demographic histories.

Distance decay principles predict that species with larger geographic ranges would have greater intraspecific genetic diversity than more restricted species. However, invasive pest species may not follow this prediction, with confounding implications for tracking phenomena including original ranges, invasion pathways and source populations. We sequenced an 815 base-pair section of the COI gene for 441 specimens of Bactrocera correcta, 214 B. zonata and 372 Zeugodacus cucurbitae; three invasive pest fruit fly species with overlapping hostplants. For each species, we explored how many individuals would need to be included in a study to sample the majority of their haplotype diversity. We also tested for phylogeographic signal and used demographic estimators as a proxy for invasion potency. We find contrasting patterns of haplotype diversity amongst the species, where B. zonata has the highest diversity but most haplotypes were represented by singletons; B. correcta has ~7 dominant haplotypes more evenly distributed; Z. cucurbitae has a single dominant haplotype with closely related singletons in a 'star-shape' surrounding it. We discuss how these differing patterns relate to their invasion histories. None of the species showed meaningful phylogeographic patterns, possibly due to gene-flow between areas across their distributions, obscuring or eliminating substructure.

Asymptotic Pairwise Dispersion Corrections Can Describe Layered Materials Accurately.

A recent study by Tawfik et al. [ Phys. Rev. Mater. 2018, 2, 034005] found that few density functionals, none of which are asymptotic pairwise dispersion methods, describe the geometry and binding of layered materials accurately. Here, we show that the exchange-hole dipole moment (XDM) dispersion model attains excellent results for graphite, hexagonal BN, and transition-metal dichalcogenides. Contrary to what has been argued, successful modeling of layered materials does not necessitate meta-GGA exchange, nonlocal correlation functionals, or the inclusion of three-body dispersion terms. Rather, a GGA functional, combined with a simple asymptotic pairwise dispersion correction, can be reliably used, provided that it properly accounts for the geometric dependence of the dispersion coefficients. The overwhelming contribution to the variation of the pairwise dispersion coefficients comes from the immediate vicinity of an atom and is already present for single layers. Longer-range and interlayer effects are examined in detail for graphite.

Dispersion XDM with Hybrid Functionals: Delocalization Error and Halogen Bonding in Molecular Crystals.

The accurate calculation of relative lattice energies of molecular crystals is important in polymorph ranking and crystal structure prediction. Delocalization error has been shown to affect calculated intermolecular binding energies in DFT and is similarly expected to affect the lattice energies of some classes of molecular crystals. In this work, we explore the use of dispersion-corrected hybrid functionals in the planewave-pseudopotentials approach to reduce delocalization error. We combine several hybrid functionals with the exchange-hole dipole moment (XDM) model for dispersion and show that they generally outperform GGA functionals in the calculation of both gas-phase binding energies and molecular crystal lattice energies. We apply the resulting XDM-corrected functionals to four halogen-bonded crystals: Cl2, Br2, I2, and ICl. GGA functionals severely overestimate their lattice energies, while hybrid functionals give accurate values. The preference of GGA functionals for monatomic structures in the Br2 and Cl2 crystals is also explained. Finally, we apply a recently developed method to calculate Bader's delocalization indices to examine the extent of intermolecular delocalization in the halogen molecular crystals. It is shown that intermolecular delocalization indices can be used to measure the strength of halogen bonds within the crystal, as well as detect the presence of delocalization error.

Six years of fruit fly surveys in Bangladesh: a new species, 33 new country records and discovery of the highly invasive Bactrocera carambolae (Diptera, Tephritidae).

We engaged in six years of snap-shot surveys for fruit flies in rural environments and ten protected forest areas of Bangladesh, using traps baited with male lures (cue-lure, methyl eugenol, zingerone). Our work has increased the recorded number of species of Tephritidae in the country from seven to 37. We summarize these surveys and report eight new country occurrence records, and a new species (Zeugodacus madhupuri Leblanc & Doorenweerd, sp. nov.) is described. The highlight among the new records is the discovery, and significant westward range extension, of Bactrocera carambolae Drew & Hancock, a major fruit pest detected in the Chattogram and Sylhet Divisions. We rectify the previously published erroneous record of Bactrocera bogorensis (Hardy), which was based on a misidentification of Zeugodacus diaphorus (Hendel). We also report the occurrence in Bangladesh of nine other Tephritidae, the rearing of three primary fruit fly parasitoids from Zeugodacus, and records of non-target attraction to fruit fly lures.