Interpreting colonization of the Calathus (Coleoptera: Carabidae) on the Canary Islands and Madeira through the application of the parametric bootstrap

The Canary Islands have proven to be an interesting archipelago for the phylogeographic study of colonization and diversification with a number of recent studies reporting evolutionary patterns and processes across a diversity of floral and faunal groups. The Canary Islands differ from the Hawaiian and Galapagos Islands by their close proximity to a continental land mass, being 110 km from the northwestern coast of Africa. This close proximity to a continent obviously increases the potential for colonization, and it can be expected that at the level of the genus some groups will be the result of more than one colonization. In this study we investigate the phylogeography of a group of carabid beetles from the genus Calathus on the Canary Islands and Madeira, located 450 km to the north of the Canaries and 650 km from the continent. The Calathus are well represented on these islands with a total of 29 species, and on the continent there are many more. Mitochondrial cytochrome oxidase I and II sequence data has been used to identify the phylogenetic relationships among the island species and a selection of continental species. Specific hypotheses of monophyly for the island fauna are tested with parametric bootstrap analysis. Data suggest that the Canary Islands have been colonized three times and Madeira twice. Four of these colonizations are of continental origin, but it is possible that one Madeiran clade may be monophyletic with a Canarian clade. The Calathus faunas of Tenerife and Madeira are recent in origin, similar to patterns previously reported for La Gomera, El Hierro, and Gran Canaria.     

VESSEL COLLISIONS WITH WHALES: THE PROBABILITY OF LETHAL INJURY BASED ON VESSEL SPEED

Historical records demonstrate that the most numerous, per capita, ocean-going-vessel strikes recorded among large-whale species accrue to the North Atlantic right whale ( Eubalaena glacialis ). As vessel speed restrictions are being considered to reduce the likelihood and severity of vessel collisions with right whales, we present an analysis of the published historical records of vessels striking large whales. We examine the influence of vessel speed in contributing to either a lethal injury (defined as killed or severely injured) or a nonlethal injury (defined as minor or no apparent injury) to a large whale when struck. A logistic regression model fitted to the observations, and consistent with a bootstrap model, demonstrates that the greatest rate of change in the probability of a lethal injury ( P _lethal) to a large whale occurs between vessel speeds of 8.6 and 15 knots where P _lethal increases from 0.21 to 0.79. The probability of a lethal injury drops below 0.5 at 11.8 knots. Above 15 knots, P _lethal asymptotically approaches 1. The uncertainties in the logistic regression estimates are relatively large at relatively low speeds ( e.g. , at 8 knots the probability is 0.17 with a 95% CI of 0.03–0.6). The results we provide can be used to assess the utility of vessel speed limits that are being considered to reduce the lethality of vessels striking the critically endangered North Atlantic right whale and other large whales that are frequent victims of vessel strikes.     

Ant Diversity in Dominant Vegetation Types of Southern Cameroon

Ants have been shown as particularly affected by land disturbance through deforestation and conversion of forest to agriculture. The effect of land use change on ant diversity in the Congo Basin is not well known. We conducted intensive sampling along a gradient of increasing vegetation disturbance to test the effect of habitat disturbance on ant diversity and Functional Groups composition. Sampling was conducted in 30 plots (5 study sites × 3 habitat × 2 plots/habitat), replicated six times in 1 year. In each plot, ants were monitored with pitfall traps, quadrats and baits. We recorded 237 ant morphospecies grouped in 10 subfamilies and 43 genera. Myrmicaria opaciventris was the most abundant species followed by Anoplolepis tenella. Forest had greater ant diversity compared with fallows and mixed‐crop fields. Functional groups were dominated by Opportunists, followed by Omnivorous Arboreal Dominants and Generalized Mymicinae. Their composition was not affected by the disturbance, but occurrence of Specialist Predators decreased with increasing disturbance. Occurrence of Generalized Myrmicinae, Opportunists and Subordinate Camponotini increased with disturbance. These results indicate that forest conversion into mixed‐crop fields reduce ant diversity. It can also increase abundance of species with generalized diet that predominates where stress and disturbance limits other ants.     

The comparison of the confidence regions in phylogeny

In this paper, several different procedures for constructing confidence regions for the true evolutionary tree are evaluated both in terms of coverage and size without considering model misspecification. The regions are constructed on the basis of tests of hypothesis using six existing tests: Shimodaira Hasegawa (SH), SOWH, star form of SOWH (SSOWH), approximately unbiased (AU), likelihood weight (LW), generalized least squares, plus two new tests proposed in this paper: single distribution nonparametric bootstrap (SDNB) and single distribution parametric bootstrap (SDPB). The procedures are evaluated on simulated trees both with small and large number of taxa. Overall, the SH, SSOWH, AU, and LW tests led to regions with higher coverage than the nominal level at the price of including large numbers of trees. Under the specified model, the SOWH test gives accurate coverage and relatively small regions. The SDNB and SDPB tests led to the small regions with occasional undercoverage. These two procedures have a substantial computational advantage over the SOWH test. Finally, the cutoff levels for the SDNB test are shown to be more variable than those for the SDPB test.     

Identifying the basal angiosperm node in chloroplast genome phylogenies: sampling one's way out of the Felsenstein zone

While there has been strong support for Amborella and Nymphaeales (water lilies) as branching from basal-most nodes in the angiosperm phylogeny, this hypothesis has recently been challenged by phylogenetic analyses of 61 protein-coding genes extracted from the chloroplast genome sequences of Amborella, Nymphaea, and 12 other available land plant chloroplast genomes. These character-rich analyses placed the monocots, represented by three grasses (Poaceae), as sister to all other extant angiosperm lineages. We have extracted protein-coding regions from draft sequences for six additional chloroplast genomes to test whether this surprising result could be an artifact of long-branch attraction due to limited taxon sampling. The added taxa include three monocots (Acorus, Yucca, and Typha), a water lily (Nuphar), a ranunculid (Ranunculus), and a gymnosperm (Ginkgo). Phylogenetic analyses of the expanded DNA and protein data sets together with microstructural characters (indels) provided unambiguous support for Amborella and the Nymphaeales as branching from the basal-most nodes in the angiosperm phylogeny. However, their relative positions proved to be dependent on the method of analysis, with parsimony favoring Amborella as sister to all other angiosperms and maximum likelihood (ML) and neighbor-joining methods favoring an Amborella + Nymphaeales clade as sister. The ML phylogeny supported the later hypothesis, but the likelihood for the former hypothesis was not significantly different. Parametric bootstrap analysis, single-gene phylogenies, estimated divergence dates, and conflicting indel characters all help to illuminate the nature of the conflict in resolution of the most basal nodes in the angiosperm phylogeny. Molecular dating analyses provided median age estimates of 161 MYA for the most recent common ancestor (MRCA) of all extant angiosperms and 145 MYA for the MRCA of monocots, magnoliids, and eudicots. Whereas long sequences reduce variance in branch lengths and molecular dating estimates, the impact of improved taxon sampling on the rooting of the angiosperm phylogeny together with the results of parametric bootstrap analyses demonstrate how long-branch attraction might mislead genome-scale phylogenetic analyses.     

Bias-corrected maximum likelihood estimator of the negative binomial dispersion parameter

We derive a first-order bias-corrected maximum likelihood estimator for the negative binomial dispersion parameter. This estimator is compared, in terms of bias and efficiency, with the maximum likelihood estimator investigated by Piegorsch (1990, Biometrics46, 863-867), the moment and the maximum extended quasi-likelihood estimators investigated by Clark and Perry (1989, Biometrics45, 309-316), and a double-extended quasi-likelihood estimator. The bias-corrected maximum likelihood estimator has superior bias and efficiency properties in most instances. For ease of comparison we give results for the two-parameter negative binomial model. However, an example involving negative binomial regression is given.     

DNA marker mining of ILSTS035 microsatellite locus on chromosome 6 of Hanwoo cattle

We describe tests for detecting and locating quantitative trait loci (QTL) for traits in Hanwoo cattle. From results of a permutation test to detect QTL for marbling, we selected the microsatellite locus ILSTS035 on chromosome 6 for further analysis. K-means clustering analysis applied to five traits and nine DNA markers in ILSTS035 resulted in three cluster groups. Finally we employed the bootstrap test method to calculate confidence intervals using the resampling method to find major DNA markers. We conclude that the major markers of ILSTS035 locus on chromosome 6 of Hanwoo cattle are markers 235 bp and 266 bp.     

The influence of temperature and precipitation on spring dispersal of Frankliniella fusca changes as the season progresses

Effects of temperature and precipitation on spring dispersal patterns of Frankliniella fusca (Hinds) (Thysanoptera: Thripidae) caught on yellow sticky traps were estimated in central and eastern North Carolina and eastern Virginia, USA, from 1997 to 2001, and in 2004 and 2007. Data were collected from 44 sites within 14 locations over 7 years, resulting in 30 location‐year data combinations. The following independent variables were examined to determine their relationship to the number of F. fusca caught on sticky traps during specified time intervals: cumulative degree days (base 10.5 °C) from 1 November to the start of each trapping interval (DD), number of days with temperatures favorable for flight during each trapping interval (DTFF), and an index of rainfall during specific intervals prior to and during the trapping interval (RI). Regression models that contained various combinations of these variables explained 62, 79, 74, and 68% of the variation in the number of dispersing F. fusca captured during 1–15 April, 16–30 April, 1–15 May, and 16–31 May, respectively. The results provide strong evidence that the suppressive effects of precipitation on growth of local populations developing during late winter and early spring are subsequently manifest at the landscape scale as reductions in the populations of dispersing adults that may persist for as long as 5–6 weeks after the precipitation occurs.     

A semiparametric response surface model for assessing drug interaction

Summary .   When multiple drugs are administered simultaneously, investigators are often interested in assessing whether the drug combinations are synergistic, additive, or antagonistic. Existing response surface models are not adequate to capture the complex patterns of drug interactions. We propose a two-component semiparametric response surface model with a parametric function to describe the additive effect of a combination dose and a nonparametric function to capture the departure from the additive effect. The nonparametric function is estimated using the technique developed in thin plate splines, and the pointwise bootstrap confidence interval for this function is constructed. The proposed semiparametric model offers an effective way of formulating the additive effect while allowing the flexibility of modeling a departure from additivity. Example and simulations are given to illustrate that the proposed model provides an excellent estimation for different patterns of interactions between two drugs.