With the advent of high-throughput technologies for measuring genome-wide expression profiles, a large number of methods have been proposed for discovering diagnostic markers that can accurately discriminate between different classes of a disease. However, factors such as the small sample size of typical clinical data, the inherent noise in high-throughput measurements, and the heterogeneity across different samples, often make it difficult to find reliable gene markers. To overcome this problem, several studies have proposed the use of pathway-based markers, instead of individual gene markers, for building the classifier. Given a set of known pathways, these methods estimate the activity level of each pathway by summarizing the expression values of its member genes, and use the pathway activities for classification. It has been shown that pathway-based classifiers typically yield more reliable results compared to traditional gene-based classifiers. In this paper, we propose a new classification method based on probabilistic inference of pathway activities. For a given sample, we compute the log-likelihood ratio between different disease phenotypes based on the expression level of each gene. The activity of a given pathway is then inferred by combining the log-likelihood ratios of the constituent genes. We apply the proposed method to the classification of breast cancer metastasis, and show that it achieves higher accuracy and identifies more reproducible pathway markers compared to several existing pathway activity inference methods.
The versatility of DNA copy number amplifications for profiling and categorization of various tissue samples has been widely acknowledged in the biomedical literature. For instance, this type of measurement techniques provides possibilities for exploring sets of cancerous tissues to identify novel subtypes. The previously utilized statistical approaches to various kinds of analyses include traditional algorithmic techniques for clustering and dimension reduction, such as independent and principal component analyses, hierarchical clustering, as well as model-based clustering using maximum likelihood estimation for latent class models.
While purely algorithmic methods are usually easily applicable, their suboptimal performance and limitations in making formal inference have been thoroughly discussed in the statistical literature. Here we introduce a Bayesian model-based approach to simultaneous identification of underlying tissue groups and the informative amplifications. The model-based approach provides the possibility of using formal inference to determine the number of groups from the data, in contrast to the ad hoc methods often exploited for similar purposes. The model also automatically recognizes the chromosomal areas that are relevant for the clustering.
Validatory analyses of simulated data and a large database of DNA copy number amplifications in human neoplasms are used to illustrate the potential of our approach. Our software implementation BASTA for performing Bayesian statistical tissue profiling is freely available for academic purposes at
Margin-based classifiers have been popular in both machine learning and statistics for classification problems. Among numerous classifiers, some are hard classifiers while some are soft ones. Soft classifiers explicitly estimate the class conditional probabilities and then perform classification based on estimated probabilities. In contrast, hard classifiers directly target on the classification decision boundary without producing the probability estimation. These two types of classifiers are based on different philosophies and each has its own merits. In this paper, we propose a novel family of large-margin classifiers, namely large-margin unified machines (LUMs), which covers a broad range of margin-based classifiers including both hard and soft ones. By offering a natural bridge from soft to hard classification, the LUM provides a unified algorithm to fit various classifiers and hence a convenient platform to compare hard and soft classification. Both theoretical consistency and numerical performance of LUMs are explored. Our numerical study sheds some light on the choice between hard and soft classifiers in various classification problems.
Class Probability Estimation; DWD; Fisher Consistency; Regularization; SVM
Under-age drinking is a long-standing public health problem in the USA and the identification of underage drinkers suffering alcohol-related problems has been difficult by using diagnostic criteria that were developed in adult populations. For this reason, it is important to characterize patterns of drinking in adolescents that are associated with alcohol-related problems. Latent class analysis is a statistical technique for explaining heterogeneity in individual response patterns in terms of a smaller number of classes. However, the latent class analysis assumption of local independence may not be appropriate when examining behavioural profiles and could have implications for statistical inference. In addition, if covariates are included in the model, non-differential measurement is also assumed. We propose a flexible set of models for local dependence and differential measurement that use easily interpretable odds ratio parameterizations while simultaneously fitting a marginal regression model for the latent class prevalences. Estimation is based on solving a set of second-order estimating equations. This approach requires only specification of the first two moments and allows for the choice of simple ‘working’ covariance structures. The method is illustrated by using data from a large-scale survey of under-age drinking. This new approach indicates the effectiveness of introducing local dependence and differential measurement into latent class models for selecting substantively interpretable models over more complex models that are deemed empirically superior.
Differential measurement; Latent class; Local dependence; Marginal regression; Odds ratio; Second-order estimating equations
Latent class models (LCMs) are used increasingly for addressing a broad variety of problems, including sparse modeling of multivariate and longitudinal data, model-based clustering, and flexible inferences on predictor effects. Typical frequentist LCMs require estimation of a single finite number of classes, which does not increase with the sample size, and have a well-known sensitivity to parametric assumptions on the distributions within a class. Bayesian nonparametric methods have been developed to allow an infinite number of classes in the general population, with the number represented in a sample increasing with sample size. In this article, we propose a new nonparametric Bayes model that allows predictors to flexibly impact the allocation to latent classes, while limiting sensitivity to parametric assumptions by allowing class-specific distributions to be unknown subject to a stochastic ordering constraint. An efficient MCMC algorithm is developed for posterior computation. The methods are validated using simulation studies and applied to the problem of ranking medical procedures in terms of the distribution of patient morbidity.
Factor analysis; Latent variables; Mixture model; Model-based clustering; Nested Dirichlet process; Order restriction; Random probability measure; Stick breaking
Obsessive-compulsive disorder (OCD) is phenomenologically heterogeneous, and findings of underlying structure classification based on symptom grouping have been ambiguous to date. Variable-centered approaches, primarily factor analysis, have been used to identify homogeneous groups of symptoms, but person-centered latent methods have seen little use. This study was designed to uncover sets of homogeneous groupings within 1611 individuals with OCD, based on symptoms.
Latent class analysis (LCA) models using 61 obsessive-compulsive symptoms (OCS) collected from the Yale-Brown Obsessive-Compulsive Scale were fit. Relationships between latent class membership and treatment response, gender, symptom severity and comorbid tic disorders were tested for relationship to class membership.
LCA models of best fit yielded three classes. Classes differed only in frequency of symptom endorsement. Classes with higher symptom endorsement were associated with earlier age of onset, being male, higher YBOCS symptom severity scores, and comorbid tic disorders. There were no differences in treatment response between classes.
These results provide support for the validity of a single underlying latent OCD construct, in addition to the distinct symptom factors identified previously via factor analyses.
obsessions; compulsions; latent class
Classification is a very useful statistical tool for information extraction. In particular, multicategory classification is commonly seen in various applications. Although binary classification problems are heavily studied, extensions to the multicategory case are much less so. In view of the increased complexity and volume of modern statistical problems, it is desirable to have multicategory classifiers that are able to handle problems with high dimensions and with a large number of classes. Moreover, it is necessary to have sound theoretical properties for the multicategory classifiers. In the literature, there exist several different versions of simultaneous multicategory Support Vector Machines (SVMs). However, the computation of the SVM can be difficult for large scale problems, especially for problems with large number of classes. Furthermore, the SVM cannot produce class probability estimation directly. In this article, we propose a novel efficient multicategory composite least squares classifier (CLS classifier), which utilizes a new composite squared loss function. The proposed CLS classifier has several important merits: efficient computation for problems with large number of classes, asymptotic consistency, ability to handle high dimensional data, and simple conditional class probability estimation. Our simulated and real examples demonstrate competitive performance of the proposed approach.
Fisher Consistency; Multicategory Classification; Regression; Squared loss; Support Vector Machine
We propose a novel learning method that combines multiple experimental modalities to improve the MHC Class-I binding prediction. Multiple experimental modalities are often accessible in the context of a binding problem. Such modalities can provide different labels of data, such as binary classifications, affinity measurements, or direct estimations of the binding profile. Current machine learning algorithms usually focus on a given label type. We here present a novel Multi-Label Vector Optimization (MLVO) formalism to produce classifiers based on the simultaneous optimization of multiple labels. Within this methodology, all label types are combined into a single constrained quadratic dual optimization problem.
We apply the MLVO to MHC class-I epitope prediction. We combine affinity measurements (IC50/EC50), binary classifications of epitopes as T cell activators and existing algorithms. The multi-label vector optimization algorithms produce classifiers significantly better than the ones resulting from any of its components. These matrix based classifier are better or equivalent to the existing state of the art MHC-I epitope prediction tools in the studied alleles.
Multiclass classification of microarray data samples with a reduced number of genes is a rich and challenging problem in Bioinformatics research. The problem gets harder as the number of classes is increased. In addition, the performance of most classifiers is tightly linked to the effectiveness of mandatory gene selection methods. Critical to gene selection is the availability of estimates about the maximum number of genes that can be handled by any classification algorithm. Lack of such estimates may lead to either computationally demanding explorations of a search space with thousands of dimensions or classification models based on gene sets of unrestricted size. In the former case, unbiased but possibly overfitted classification models may arise. In the latter case, biased classification models unable to support statistically significant findings may be obtained.
A novel bound on the maximum number of genes that can be handled by binary classifiers in binary mediated multiclass classification algorithms of microarray data samples is presented. The bound suggests that high-dimensional binary output domains might favor the existence of accurate and sparse binary mediated multiclass classifiers for microarray data samples.
A comprehensive experimental work shows that the bound is indeed useful to induce accurate and sparse multiclass classifiers for microarray data samples.
A range of psychiatric symptoms and cognitive deficits occur in Parkinson’s disease (PD), and symptom overlap and co-morbidity complicate the classification of non-motor symptoms. The objective of this study was to use analytic-based approaches to classify psychiatric and cognitive symptoms in PD.
Cross-sectional evaluation of a convenience sample of patients in specialty care.
Two outpatient movement disorders centers at the University of Pennsylvania and Philadelphia Veterans Affairs Medical Center.
177 patients with mild-moderate idiopathic PD and without significant global cognitive impairment.
Subjects were assessed with an extensive psychiatric, neuropsychological, and neurological battery. Latent class analysis (LCA) was used to statistically delineate group-level symptom profiles across measures of psychiatric and cognitive functioning. Predictors of class membership were also examined.
Results from the LCA indicated that a four-class solution best fit the data. 32.3% of the sample had good psychiatric and normal cognitive functioning, 17.5% had significant psychiatric co-morbidity but normal cognition, 26.0% had few psychiatric symptoms but had poorer cognitive functioning across a range of cognitive domains, and 24.3% had both significant psychiatric co-morbidity and poorer cognitive functioning. Age, disease severity, and medication use predicted class membership.
LCA delineates four classes of patients in mild-moderate PD, three of which experience significant non-motor impairments and comprise over two-thirds of patients. Neuropsychiatric symptoms and cognitive deficits follow distinct patterns in PD, and further study is needed to determine if these classes are generalizable, stable, predict function, quality of life and long-term outcomes, and are amenable to treatment at a class level.
Parkinson’s disease; neuropsychology; psychiatry; cognition; depression
Bioinformatics data analysis often deals with additive mixtures of signals for which only class labels are known. Then, the overall goal is to estimate class related signals for data mining purposes. A convenient application is metabolic monitoring of patients using infrared spectroscopy. Within an infrared spectrum each single compound contributes quantitatively to the measurement.
In this work, we propose a novel factorization technique for additive signal factorization that allows learning from classified samples. We define a composed loss function for this task and analytically derive a closed form equation such that training a model reduces to searching for an optimal threshold vector. Our experiments, carried out on synthetic and clinical data, show a sensitivity of up to 0.958 and specificity of up to 0.841 for a 15-class problem of disease classification. Using class and regression information in parallel, our algorithm outperforms linear SVM for training cases having many classes and few data.
The presented factorization method provides a simple and generative model and, therefore, represents a first step towards predictive factorization methods.
Uncertainty often affects molecular biology experiments and data for different reasons. Heterogeneity of gene or protein expression within the same tumor tissue is an example of biological uncertainty which should be taken into account when molecular markers are used in decision making. Tissue Microarray (TMA) experiments allow for large scale profiling of tissue biopsies, investigating protein patterns characterizing specific disease states. TMA studies deal with multiple sampling of the same patient, and therefore with multiple measurements of same protein target, to account for possible biological heterogeneity. The aim of this paper is to provide and validate a classification model taking into consideration the uncertainty associated with measuring replicate samples.
We propose an extension of the well-known Naïve Bayes classifier, which accounts for biological heterogeneity in a probabilistic framework, relying on Bayesian hierarchical models. The model, which can be efficiently learned from the training dataset, exploits a closed-form of classification equation, thus providing no additional computational cost with respect to the standard Naïve Bayes classifier. We validated the approach on several simulated datasets comparing its performances with the Naïve Bayes classifier. Moreover, we demonstrated that explicitly dealing with heterogeneity can improve classification accuracy on a TMA prostate cancer dataset.
The proposed Hierarchical Naïve Bayes classifier can be conveniently applied in problems where within sample heterogeneity must be taken into account, such as TMA experiments and biological contexts where several measurements (replicates) are available for the same biological sample. The performance of the new approach is better than the standard Naïve Bayes model, in particular when the within sample heterogeneity is different in the different classes.
Study objective: To examine social inequalities in minor psychiatric morbidity as measured by the GHQ-12 using lagged models of psychiatric morbidity and changing job status.
Design: GHQ scores were modelled using two level hierarchical regression models with measurement occasions nested within individuals. The paper compares and contrasts three different ways of describing social position: income, social advantage and lifestyle (the Cambridge scale), and social class (the new National Statistics Socio-Economic Classification), and adjusts for attrition.
Setting: Survey interviews for a nationally representative sample of adults of working age living in Britain.
Participants: 8091 original adult respondents in 1991 who remain of working age during 1991–1998 from the British Household Panel Survey (BHPS).
Main results: There was a relation of GHQ-12 to social position when social position was combined with employment status. This relation itself varied according to a person's psychological health in the previous year.
Conclusions: The relation between social position and minor psychiatric morbidity depended on whether or not a person was employed, unemployed, or economically inactive. It was stronger in those with previously less good psychological health. Among employed men and women in good health, GHQ-12 varied little according to social class, status, or income. There was a "classic" social gradient in psychiatric morbidity, with worse health in less advantaged groups, among the economically inactive. Among the unemployed, a "reverse" gradient was found: the impact of unemployment on minor psychiatric morbidity was higher for those who were previously in a more advantaged social class position.
The advent of microarray technology has made it possible to classify disease states based on gene expression profiles of patients. Typically, marker genes are selected by measuring the power of their expression profiles to discriminate among patients of different disease states. However, expression-based classification can be challenging in complex diseases due to factors such as cellular heterogeneity within a tissue sample and genetic heterogeneity across patients. A promising technique for coping with these challenges is to incorporate pathway information into the disease classification procedure in order to classify disease based on the activity of entire signaling pathways or protein complexes rather than on the expression levels of individual genes or proteins. We propose a new classification method based on pathway activities inferred for each patient. For each pathway, an activity level is summarized from the gene expression levels of its condition-responsive genes (CORGs), defined as the subset of genes in the pathway whose combined expression delivers optimal discriminative power for the disease phenotype. We show that classifiers using pathway activity achieve better performance than classifiers based on individual gene expression, for both simple and complex case-control studies including differentiation of perturbed from non-perturbed cells and subtyping of several different kinds of cancer. Moreover, the new method outperforms several previous approaches that use a static (i.e., non-conditional) definition of pathways. Within a pathway, the identified CORGs may facilitate the development of better diagnostic markers and the discovery of core alterations in human disease.
The advent of microarray technology has drawn immense interest to identify gene expression levels that can serve as biomarkers for disease. Marker genes are selected by examining each individual gene to see how well its expression level discriminates different disease types. In complex diseases such as cancer, good marker genes can be hard to find due to cellular heterogeneity within the tissue and genetic heterogeneity across patients. A promising technique for addressing these challenges is to incorporate biological pathway information into the marker identification procedure, permitting disease classification based on the activity of entire pathways rather than simply on the expression levels of individual genes. However, previous pathway-based methods have not significantly outperformed gene-based methods. Here, we propose a new pathway-based classification procedure in which markers are encoded not as individual genes, nor as the set of genes making up a known pathway, but as subsets of “condition-responsive genes (CORGs)” within those pathways. Using expression profiles from seven different microarray studies, we show that the accuracy of this method is significantly better than both the conventional gene- and pathway- based diagnostics. Furthermore, the identified CORGs may facilitate the development of effective diagnostic markers and the discovery of molecular mechanisms underlying disease.
Mitochondria exist as a network of interconnected organelles undergoing constant fission and fusion. Current approaches to study mitochondrial morphology are limited by low data sampling coupled with manual identification and classification of complex morphological phenotypes. Here we propose an integrated mechanistic and data-driven modeling approach to analyze heterogeneous, quantified datasets and infer relations between mitochondrial morphology and apoptotic events. We initially performed high-content, multi-parametric measurements of mitochondrial morphological, apoptotic, and energetic states by high-resolution imaging of human breast carcinoma MCF-7 cells. Subsequently, decision tree-based analysis was used to automatically classify networked, fragmented, and swollen mitochondrial subpopulations, at the single-cell level and within cell populations. Our results revealed subtle but significant differences in morphology class distributions in response to various apoptotic stimuli. Furthermore, key mitochondrial functional parameters including mitochondrial membrane potential and Bax activation, were measured under matched conditions. Data-driven fuzzy logic modeling was used to explore the non-linear relationships between mitochondrial morphology and apoptotic signaling, combining morphological and functional data as a single model. Modeling results are in accordance with previous studies, where Bax regulates mitochondrial fragmentation, and mitochondrial morphology influences mitochondrial membrane potential. In summary, we established and validated a platform for mitochondrial morphological and functional analysis that can be readily extended with additional datasets. We further discuss the benefits of a flexible systematic approach for elucidating specific and general relationships between mitochondrial morphology and apoptosis.
The purpose of the study was to explore heterogeneity and differential treatment outcome among a sample of patients with binge eating disorder (BED).
A latent class analysis was conducted with 205 treatment-seeking, overweight or obese individuals with BED randomized to Interpersonal Psychotherapy (IPT), Behavioral Weight Loss (BWL), or guided self-help based on Cognitive Behavioral Therapy (CBTgsh). A latent transition analysis tested the predictive validity of the latent class analysis model.
A 4-class model yielded the best overall fit to the data. Class 1 was characterized by a lower mean body mass index (BMI) and increased physical activity. Individuals in class 2 reported the most binge eating, shape and weight concerns, compensatory behaviors, and negative affect. Class 3 patients reported similar binge eating frequencies to class 2 with lower levels of exercise or compensation. Class 4 was characterized by the highest average BMI, the most overeating episodes, fewer binge episodes, and an absence of compensatory behaviors. Classes 1 and 3 had the highest and lowest percentage of individuals with a past eating disorder diagnosis, respectively. The latent transition analysis found a higher probability of remission from binge eating among those receiving IPT in Class 2 and CBTgsh in Class 3.
The latent class analysis identified four distinct classes using baseline measures of eating disorder and depressive symptoms, body weight, and physical activity. Implications of the observed differential treatment response are discussed.
binge eating disorder; latent class analysis; latent transition analysis; treatment specificity
Non ignorable missing data is a common problem in longitudinal studies. Latent class models are attractive for simplifying the modeling of missing data when the data are subject to either a monotone or intermittent missing data pattern. In our study, we propose a new two-latent-class model for categorical data with informative dropouts, dividing the observed data into two latent classes; one class in which the outcomes are deterministic and a second one in which the outcomes can be modeled using logistic regression. In the model, the latent classes connect the longitudinal responses and the missingness process under the assumption of conditional independence. Parameters are estimated by the method of maximum likelihood estimation based on the above assumptions and the tetrachoric correlation between responses within the same subject. We compare the proposed method with the shared parameter model and the weighted GEE model using the areas under the ROC curves in the simulations and the application to the smoking cessation data set. The simulation results indicate that the proposed two-latent-class model performs well under different missing procedures. The application results show that our proposed method is better than the shared parameter model and the weighted GEE model.
Area under ROC curve; Informative dropout; Latent class; Tetrachoric correlation
Very few analytical approaches have been reported to resolve the variability in microarray measurements stemming from sample heterogeneity. For example, tissue samples used in cancer studies are usually contaminated with the surrounding or infiltrating cell types. This heterogeneity in the sample preparation hinders further statistical analysis, significantly so if different samples contain different proportions of these cell types. Thus, sample heterogeneity can result in the identification of differentially expressed genes that may be unrelated to the biological question being studied. Similarly, irrelevant gene combinations can be discovered in the case of gene expression based classification.
We propose a computational framework for removing the effects of sample heterogeneity by "microdissecting" microarray data in silico. The computational method provides estimates of the expression values of the pure (non-heterogeneous) cell samples. The inversion of the sample heterogeneity can be facilitated by providing accurate estimates of the mixing percentages of different cell types in each measurement. For those cases where no such information is available, we develop an optimization-based method for joint estimation of the mixing percentages and the expression values of the pure cell samples. We also consider the problem of selecting the correct number of cell types.
The efficiency of the proposed methods is illustrated by applying them to a carefully controlled cDNA microarray data obtained from heterogeneous samples. The results demonstrate that the methods are capable of reconstructing both the sample and cell type specific expression values from heterogeneous mixtures and that the mixing percentages of different cell types can also be estimated. Furthermore, a general purpose model selection method can be used to select the correct number of cell types.
The cumulative incidence function is commonly reported in studies with competing risks. The aim of this paper is to compute the treatment-specific cumulative incidence functions, adjusting for potentially imbalanced prognostic factors among treatment groups. The underlying regression model considered in this study is the proportional hazards model for a subdistribution function . We propose estimating the direct adjusted cumulative incidences for each treatment using the pooled samples as the reference population. We develop two SAS macros for estimating the direct adjusted cumulative incidence function for each treatment based on two regression models. One model assumes the constant subdistribution hazard ratios between the treatments and the alternative model allows each treatment to have its own baseline subdistribution hazard function. The macros compute the standard errors for the direct adjusted cumulative incidence estimates, as well as the standard errors for the differences of adjusted cumulative incidence functions between any two treatments. Based on the macros’ output, one can assess treatment effects at predetermined time points. A real bone marrow transplant data example illustrates the practical utility of the SAS macros.
Competing risks; Cumulative incidence function; Cox model; Subdistribution
This study aimed to investigate empirically how and in what way individuals with symptoms of functional somatic syndromes should be classified. We also aimed to look into genetic and environmental influences on the classification.
A total of 28531 twins aged 41–64 underwent screening interviews via a computer-assisted data collection system from 1998 to 2002. Nine functional somatic symptoms (abnormal tiredness, general muscular pain, recurrent abdominal discomfort, back pain, gastroesophageal reflux, recurrent headache, recurrent urinary problem, dizziness, breathlessness at rest) were assessed using structured questions in a blinded manner. Latent class analysis was applied to the data. Structural equation modeling was further performed in order to estimate the relative importance of genetic and environmental influences on class probability.
Latent class analysis resulted in a 5-class solution. Individuals in the first class did not show any health problems. Those assigned to the second, third, and fourth classes tended to have abnormal tiredness, gastrointestinal problems, and pain-related symptoms, respectively. Individuals in the fifth class had multiple symptoms to a greater extent than the other classes. All the five classes showed modest genetic influences (7 – 29% of the total variation) with gender differences except Class 3; however, the majority of influences on the class membership derived from unique environmental effects.
The findings suggested the necessity of re-defining the existing classification criteria for functional somatic syndromes in terms of single (uncomplicated) or multiple (complicated) syndromes. Environmental influences are important for the aetiology of functional somatic syndromes.
functional somatic syndromes; chronic fatigue syndrome; chronic widespread pain; irritable bowel syndrome; comorbidity; latent class analysis
In attention/hyperactivity and aggressive behavior problems were measured in 335 children from school entry throughout adolescence, at 3-year intervals. Children were participants in a high-risk prospective study of substance use disorders and comorbid problems. A parallel process latent growth model found aggressive behavior decreasing throughout childhood and adolescence, whereas inattentive/hyperactive behavior levels were constant. Growth mixture modeling, in which developmental trajectories are statistically classified, found two classes for inattention/hyperactivity and two for aggressive behavior, resulting in a total of four trajectory classes. Different influences of the family environment predicted development of the two types of behavior problems when the other behavior problem was held constant. Lower emotional support and lower intellectual stimulation by the parents in early childhood predicted membership in the high problem class of inattention/hyperactivity when the trajectory of aggression was held constant. Conversely, conflict and lack of cohesiveness in the family environment predicted membership in a worse developmental trajectory of aggressive behavior when the inattention/hyperactivity trajectories were held constant. The implications of these findings for the development of inattention/hyperactivity and for the development of risk for the emergence of substance use disorders are discussed.
The braconid subfamily Rogadinae is a large, cosmopolitan group of endoparasitoid wasps characterised by 'mummifying' their lepidopteran host larvae, from which the adult subsequently emerges. Rogadines attack a variety of both macro- and microlepidopteran taxa, although the speciose genus Aleiodes almost exclusively attacks macrolepidopterans. Here, we investigate the phylogenetic history of the Rogadinae, revise their higher-level classification and assess the evolution of their host ranges and mummy types. We also assess the divergence times within the subfamily and discuss the reasons for the extraordinary evolutionary diversification of Aleiodes.
Our Bayesian analyses weakly support the monophyly of the subfamily. A clade comprising all Aleiodes species and some other taxa is not nested within the tribe Rogadini as previously supposed, but instead is recovered as sister to the Yeliconini, with the remaining Rogadini genera being recovered as sister to the Stiropiini. The Rogadinae is estimated to have originated during the mid to late Eocene, 36.1–51.62 MYA. Molecular dating gives a more recent origin for the Aleiodes clade (17.98–41.76 MYA) compared to the origins proposed for two of its principal lepidopteran host groups (Noctuidae: 60.7–113.4 MYA; Geometridae 48–62 MYA). The Bayesian ancestral reconstruction of the emergence habits from the mummified hosts weakly recovered an anterior emergence as the ancestral condition for the subfamily. Producing a hard mummy has evolved at various times independently, though most of the species with this biology belong to the Aleiodes clade.
Based on our results, we erect the tribe Aleiodini nov. to include Aleiodes and Heterogamus stat. rev. Cordylorhogas, Pholichora and Hemigyroneuron are synonymised with Aleiodes. The molecular dating of clades and the ancestral reconstruction of host ranges support the hypothesis that radiation within Aleiodes s. s. was due to host recruitment leading to host range expansion followed by speciation, and not to parasitoid-host coevolution. Within the Rogadinae, variation in the site of emergence from the mummified host probably evolved as a consequence of the mummy's site and mode of formation, and the extent of mummy tanning/hardness to the degree of protection needed in relation to the cost of providing it.
An important statistical problem associated with tumor classification is the identification of new tumor classes using gene-expression profiles. One aspect of this clustering problem is estimating the number of clusters in a dataset. A new prediction-based resampling method, Clest, was developed to address this problem.
Microarray technology is increasingly being applied in biological and medical research to address a wide range of problems, such as the classification of tumors. An important statistical problem associated with tumor classification is the identification of new tumor classes using gene-expression profiles. Two essential aspects of this clustering problem are: to estimate the number of clusters, if any, in a dataset; and to allocate tumor samples to these clusters, and assess the confidence of cluster assignments for individual samples. Here we address the first of these problems.
We have developed a new prediction-based resampling method, Clest, to estimate the number of clusters in a dataset. The performance of the new and existing methods were compared using simulated data and gene-expression data from four recently published cancer microarray studies. Clest was generally found to be more accurate and robust than the six existing methods considered in the study.
Focusing on prediction accuracy in conjunction with resampling produces accurate and robust estimates of the number of clusters.
Macro-hormones and macro-enzymes are high molecular weight conjugates of hormones or enzymes, respectively, often with immunoglobulins. These are referred to as macromolecular complexes, and may cause artefactually elevated biochemical tests results. Macro enzymes of the most commonly measured serum enzymes have been identified and are recognised as a source of elevated measurements that may cause diagnostic confusion; macro-creatine kinase and macro-amylase are the two better known macro-enzymes in clinical practice. Literature on macro-hormones is largely restricted to macro-prolactin. We present a case of a clinically euthyroid patient, who had persistently elevated thyroid stimulating hormone (TSH) but free thyroxine within the reference limits. She underwent repeated thyroid investigations and thyroid hormone interference studies, until macro-TSH was identified as the most likely cause of unexplained elevated TSH. Following the identification and characterisation of this biochemical abnormality, she is no longer subject to repeated blood tests for assessment of thyroid function; the patient currently remains clinically euthyroid.
An essential step of a metagenomic study is the taxonomic classification, that is, the identification of the taxonomic lineage of the organisms in a given sample. The taxonomic classification process involves a series of decisions. Currently, in the context of metagenomics, such decisions are usually based on empirical studies that consider one specific type of classifier. In this study we propose a general framework for analyzing the impact that several decisions can have on the classification problem. Instead of focusing on any specific classifier, we define a generic score function that provides a measure of the difficulty of the classification task. Using this framework, we analyze the impact of the following parameters on the taxonomic classification problem: (i) the length of n-mers used to encode the metagenomic sequences, (ii) the similarity measure used to compare sequences, and (iii) the type of taxonomic classification, which can be conventional or hierarchical, depending on whether the classification process occurs in a single shot or in several steps according to the taxonomic tree.
We defined a score function that measures the degree of separability of the taxonomic classes under a given configuration induced by the parameters above. We conducted an extensive computational experiment and found out that reasonable values for the parameters of interest could be (i) intermediate values of n, the length of the n-mers; (ii) any similarity measure, because all of them resulted in similar scores; and (iii) the hierarchical strategy, which performed better in all of the cases.
As expected, short n-mers generate lower configuration scores because they give rise to frequency vectors that represent distinct sequences in a similar way. On the other hand, large values for n result in sparse frequency vectors that represent differently metagenomic fragments that are in fact similar, also leading to low configuration scores. Regarding the similarity measure, in contrast to our expectations, the variation of the measures did not change the configuration scores significantly. Finally, the hierarchical strategy was more effective than the conventional strategy, which suggests that, instead of using a single classifier, one should adopt multiple classifiers organized as a hierarchy.
Metagenomics; classification problem; taxonomic classification