Phylogenomic resolution of scorpions reveals multilevel discordance with morphological phylogenetic signal

Scorpions represent an iconic lineage of arthropods, historically renowned for their unique bauplan, ancient fossil record and venom potency. Yet, higher level relationships of scorpions, based exclusively on morphology, remain virtually untested, and no multilocus molecular phylogeny has been deployed heretofore towards assessing the basal tree topology. We applied a phylogenomic assessment to resolve scorpion phylogeny, for the first time, to our knowledge, sampling extensive molecular sequence data from all superfamilies and examining basal relationships with up to 5025 genes. Analyses of supermatrices as well as species tree approaches converged upon a robust basal topology of scorpions that is entirely at odds with traditional systematics and controverts previous understanding of scorpion evolutionary history. All analyses unanimously support a single origin of katoikogenic development, a form of parental investment wherein embryos are nurtured by direct connections to the parent''s digestive system. Based on the phylogeny obtained herein, we propose the following systematic emendations: Caraboctonidae is transferred to Chactoidea new superfamilial assignment; superfamily Bothriuroidea revalidated is resurrected and Bothriuridae transferred therein; and Chaerilida and Pseudochactida are synonymized with Buthida new parvordinal synonymies.     

On the performance of adaptive preprocessing technique in analyzing high‐dimensional censored data

Preprocessing for high‐dimensional censored datasets, such as the microarray data, is generally considered as an important technique to gain further stability by reducing potential noise from the data. When variable selection including inference is carried out with high‐dimensional censored data the objective is to obtain a smaller subset of variables and then perform the inferential analysis using model estimates based on the selected subset of variables. This two stage inferential analysis is prone to circularity bias because of the noise that might still remain in the dataset. In this work, I propose an adaptive preprocessing technique that uses sure independence screening (SIS) idea to accomplish variable selection and reduces the circularity bias by some popularly known refined high‐dimensional methods such as the elastic net, adaptive elastic net, weighted elastic net, elastic net‐AFT, and two greedy variable selection methods known as TCS, PC‐simple all implemented with the accelerated lifetime models. The proposed technique addresses several features including the issue of collinearity between important and some unimportant covariates, which is often the case in high‐dimensional setting under variable selection framework, and different level of censoring. Simulation studies along with an empirical analysis with a real microarray data, mantle cell lymphoma, is carried out to demonstrate the performance of the adaptive pre‐processing technique.     

How does dune morphology shape coastal EC habitats occurrence? A remote sensing approach using airborne LiDAR on the Mediterranean coast

We examined the relationship between coastal habitats (sensu European Union Habitats Directive) and local dune morphology along a Mediterranean coastal dune system by integrating field collected vegetation data and remotely sensed imagery. Specifically, we described the morphological profile of each EC habitat based on the morphological variables that are most likely to affect their occurrence, including elevation, slope, curvature, northness, eastness and sea distance. In addition, we assessed the role and strength of each morphological variable in determining the occurrence of EC habitats.We used 394 random vegetation plots representative of six EC habitats (Habitat 1210: “Annual vegetation of drift lines”; Habitat 2110: “Embryonic shifting dunes”; Habitat 2120: “Shifting dunes along the shoreline with Ammophila arenaria”; Habitat 2210 and 2230: “Crucianellion maritimae fixed beach dunes” and “Malcolmietalia dune grasslands”; Habitat 2250: “Coastal dunes with Juniperus spp.”; Habitat 2260: “Cisto-Lavanduletalia dune sclerophyllous scrubs”) found along the Tyrrhenian coast of central Italy. We derived each morphological variable from a DTM (Digital Terrain Model) obtained from 2-m resolution LiDAR (Light Detection And Range) images. The mean value of each variable was calculated at different spatial scales using buffer areas of increasing radius (2 m, 4 m, 8 m) around each vegetation plot. Mean morphological values for each EC habitat were compared using Kruskal-Wallis rank test. The role and strength of the relationship between habitat type and the morphological variables were assessed using Generalized Linear Models.EC habitats occur differentially across dune morphology, and the role and strength of each morphological variable define habitat specificity. Dune elevation and sea distance were determined to be the key factors in shaping EC habitat occurrence along this section of the Mediterranean coast. Identification of the close relationship between habitat type and morphological variables deriving from airborne LiDAR imagery points to the high potential of such remote sensing tool for analyzing and monitoring the integrity of coastal dune ecosystems. As airborne LiDAR enables the rapid collection of extremely accurate topographic data over large areas, it also offers useful information for the management of these threatened and fragile ecosystems.     

Market basket analysis of grasshopper (Orthoptera: Acrididae) assemblages in eastern Wyoming: a 17‐year case study using associative analysis for ecological insights into grasshopper outbreaks

1. This study utilised an associative analysis (AA) technique named market basket analysis (MBA) to investigate whether particular grasshopper (Orthoptera: Acrididae) species associations are common during outbreaks (>9.6 grasshoppers m^?2) that last >3 years. This study is the first of its kind to use MBA on animal communities. 2. A subset of the 17 years of grasshopper density data from the Wyoming Grasshopper Survey Dataset was used to explore associations among grasshopper species. 3. Associations of certain species were found with over 80% confidence. Life‐history traits of those species commonly found together were examined and compared (a posteriori), creating opportunities to hypothesise certain ecological relationships (e.g. interspecific competition, indirect mutualism) for future studies. 4. This case study shows that further MBA analysis of grasshopper assemblages should prove useful in discovering ecological relationships of grasshopper species during outbreaks. Preliminary examples are demonstrated.     

Integrating Linguistics,Social Structure,and Geography to Model Genetic Diversity within India

India represents an intricate tapestry of population substructure shaped by geography, language, culture, and social stratification. Although geography closely correlates with genetic structure in other parts of the world, the strict endogamy imposed by the Indian caste system and the large number of spoken languages add further levels of complexity to understand Indian population structure. To date, no study has attempted to model and evaluate how these factors have interacted to shape the patterns of genetic diversity within India. We merged all publicly available data from the Indian subcontinent into a data set of 891 individuals from 90 well-defined groups. Bringing together geography, genetics, and demographic factors, we developed Correlation Optimization of Genetics and Geodemographics to build a model that explains the observed population genetic substructure. We show that shared language along with social structure have been the most powerful forces in creating paths of gene flow in the subcontinent. Furthermore, we discover the ethnic groups that best capture the diverse genetic substructure using a ridge leverage score statistic. Integrating data from India with a data set of additional 1,323 individuals from 50 Eurasian populations, we find that Indo-European and Dravidian speakers of India show shared genetic drift with Europeans, whereas the Tibeto-Burman speaking tribal groups have maximum shared genetic drift with East Asians.     

Mechanical biomarkers of medial compartment knee osteoarthritis diagnosis and severity grading: Discovery phase

ObjectiveTo investigate, as a discovery phase, if 3D knee kinematics assessment parameters can serve as mechanical biomarkers, more specifically as diagnostic biomarker and burden of disease biomarkers, as defined in the Burden of Disease, Investigative, Prognostic, Efficacy of Intervention and Diagnostic classification scheme for osteoarthritis (OA) (Altman et al., 1986). These biomarkers consist of a set of biomechanical parameters discerned from 3D knee kinematic patterns, namely, flexion/extension, abduction/adduction, and tibial internal/external rotation measurements, during gait recording.Methods100 medial compartment knee OA patients and 40 asymptomatic control subjects participated in this study. OA patients were categorized according to disease severity, by the Kellgren and Lawrence grading system. The proposed biomarkers were identified by incremental parameter selection in a regression tree of cross-sectional data. Biomarker effectiveness was evaluated by receiver operating characteristic curve analysis, namely, the area under the curve (AUC), sensitivity and specificity.ResultsDiagnostic biomarkers were defined by a set of 3 abduction/adduction kinematics parameters. The performance of these biomarkers reached 85% for the AUC, 80% for sensitivity and 90% for specificity; the likelihood ratio was 8%. Burden of disease biomarkers were defined by a 3-decision tree, with sets of kinematics parameters selected from all 3 movement planes.ConclusionThe results demonstrate, as part of a discovery phase, that sets of 3D knee kinematic parameters have the potential to serve as diagnostic and burden of disease biomarkers of medial compartment knee OA.     

PFA: Program for the Quantitative Assessment of Cell Metabolism by Spectral Data Analysis

Assessment of mitochondrial oxidative metabolism has wide-ranging importance, from pharmacokinetic analysis to studies in cell viability and apoptosis. Here we present the Perfusion File Analyzer (PFA) application for the real-time analysis of spectral data to measure cytochrome c reduction, cytochrome a3 reduction, and other parameters important to cellular metabolism, which are collected during tissue perfusion experiments. Our current efforts are focused on quantitating changes in mitochondrial function by normalizing baseline drift of spectral data while addressing two major challenges: (1) a lack of real-time feedback from the system when aiming is compromised, and (2) an inability to adjust calculated data in the event of spectral shift. PFA has been developed to address these issues, and is currently used for quality assessment of human islets prior to clinical transplantation.     

Simulating tick distributions over sub-Saharan Africa: the use of observed and simulated climate surfaces

Aim A broad suit of climate data sets is becoming available for use in predictive species modelling. We compare the efficacy of using interpolated climate surfaces [Center for Resource and Environmental Studies (CRES) and Climate Research Unit (CRU)] or high‐resolution model‐derived climate data [Division of Atmospheric Research limited‐area model (DARLAM)] for predictive species modelling, using tick distributions from sub‐Saharan Africa. Location The analysis is restricted to sub‐Saharan Africa. The study area was subdivided into 3000 grids cells with a resolution of 60 × 60 km. Methods Species distributions were predicted using an established multivariate climate envelope modelling approach and three very different climate data sets. The recorded variance in the climate data sets was quantified by employing omnidirectional variograms. To further compare the interpolated tick distributions that flowed from using three climate data sets, we calculated true positive (TP) predictions, false negative (FN) predictions as well as the proportional overlaps between observed and modelled tick distributions. In addition, the effect of tick data set size on the performance of the climate data sets was evaluated by performing random draws of known tick distribution records without replacement. Results The predicted distributions were consistently wider ranging than the known records when using any of the three climate data sets. However, the proportional overlap between predicted and known distributions varied as follows: for Rhipicephalus appendiculatus Neumann (Acari: Ixodidae), these were 60%, 60% and 70%; for Rhipicephalus longus Neumann (Acari: Ixodidae) 60%, 57% and 75%; for Rhipicephalus zambeziensis Walker, Norval & Corwin (Acari: Ixodidae) 57%, 51% and 62%, and for Rhipicephalus capensis Koch (Acari: Ixodidae) 70%, 60% and 60% using the CRES, CRU and DARLAM climate data sets, respectively. All data sets were sensitive to data size but DARLAM performed better when using smaller species data sets. At a 20% data subsample level, DARLAM was able to capture more than 50% of the known records and captured more than 60% of known records at higher subsample levels. Main conclusions The use of data derived from high‐resolution nested climate models (e.g. DARLAM) provided equal or even better species distribution modelling performance. As the model is dynamic and process based, the output data are available at the modelled resolution, and are not hamstrung by the sampling intensity of observed climate data sets (c. one sample per 30,000 km² for Africa). In addition, when exploring the biodiversity consequences of climate change, these modelled outputs form a more useful basis for comparison with modelled future climate scenarios.     

Predictability of the terrestrial carbon cycle

Terrestrial ecosystems sequester roughly 30% of anthropogenic carbon emission. However this estimate has not been directly deduced from studies of terrestrial ecosystems themselves, but inferred from atmospheric and oceanic data. This raises a question: to what extent is the terrestrial carbon cycle intrinsically predictable? In this paper, we investigated fundamental properties of the terrestrial carbon cycle, examined its intrinsic predictability, and proposed a suite of future research directions to improve empirical understanding and model predictive ability. Specifically, we isolated endogenous internal processes of the terrestrial carbon cycle from exogenous forcing variables. The internal processes share five fundamental properties (i.e., compartmentalization, carbon input through photosynthesis, partitioning among pools, donor pool‐dominant transfers, and the first‐order decay) among all types of ecosystems on the Earth. The five properties together result in an emergent constraint on predictability of various carbon cycle components in response to five classes of exogenous forcing. Future observational and experimental research should be focused on those less predictive components while modeling research needs to improve model predictive ability for those highly predictive components. We argue that an understanding of predictability should provide guidance on future observational, experimental and modeling research.