Conflicting phylogenetic signals in genomic data of the coffee family (Rubiaceae)

Reconstructions of phylogenetic relationships in the flowering plant family Rubiaceae have up until now relied heavily on single‐ or multi‐gene data, primarily from the plastid compartment. With the availability of cost‐ and time‐efficient techniques for generating complete genome sequences, the opportunity arises to resolve some of the relationships that, up until now, have proven problematic. Here, we contribute new data from complete 58 plastid genome sequences, representing 55 of the currently 65 recognized tribes of the Rubiaceae. Also contributed are new data from the nuclear rDNA cistrons for corresponding taxa. Phylogenetic analyses are conducted on two plastid data sets, one including data from the protein coding genes only, and a second where protein coding data are combined with non‐coding regions, and on a nuclear rDNA data set. Our results clearly show that simply adopting a “more characters” approach does not resolve the relationships in the Rubiaceae. More importantly, we identify conflicting phylogenetic signals in the data. Analyses of the same plastid data, treated as nucleotides or as codon‐degenerated data, resolve and support conflicting topologies in the subfamily Cinchonoideae. As these analyses use the same data, we interpret the conflict to result from erroneous assumptions in the models used to reconstruct our phylogenies. Conflicting signals are also identified in the analyses of the plastid versus the nuclear rDNA data sets. These analyses use data from different genomic compartments, with different inheritance patterns, and we interpret the conflicts as representing “real” conflicts, reflecting biological processes of the past.     

Expert Knowledge Influences Decision-Making for Couples Receiving Positive Prenatal Chromosomal Microarray Testing Results

To assess how participants receiving abnormal prenatal genetic testing results seek information and understand the implications of results, 27 US female patients and 12 of their male partners receiving positive prenatal microarray testing results completed semi-structured phone interviews. These interviews documented participant experiences with chromosomal microarray testing, understanding of and emotional response to receiving results, factors affecting decision-making about testing and pregnancy termination, and psychosocial needs throughout the testing process. Interview data were analyzed using a modified grounded theory approach. In the absence of certainty about the implications of results, understanding of results is shaped by biomedical expert knowledge (BEK) and cultural expert knowledge (CEK). When there is a dearth of BEK, as in the case of receiving results of uncertain significance, participants rely on CEK, including religious/spiritual beliefs, “gut instinct,” embodied knowledge, and social network informants. CEK is a powerful platform to guide understanding of prenatal genetic testing results. The utility of culturally situated expert knowledge during testing uncertainty emphasizes that decision-making occurs within discourses beyond the biomedical domain. These forms of “knowing” may be integrated into clinical consideration of efficacious patient assessment and counseling.     

Comprehensive Realism's Weight-of-Evidence Based Distributional Dose-Response Characterization

Challenges to low-dose linearity and other default assumptions in cancer risk assessment and the limitations associated with NOAELs, LOAELs, and constant uncertainty factor values in the evaluation of noncancer health effects have stimulated the continued evolution of risk assessment methodologies. The increasing need for more realistic estimates of the dose-response relationship, better uncertainty characterization, and greater utilization of cost-benefit analyses have also contributed to this evolution. “Comprehensive Realism” is an emerging quantitative weight-of-evidence based risk assessment methodology for both cancer and noncancer health effects which utilizes probability distributions and decision analysis techniques to reflect more of the relevant human exposure data, more of the available and pertinent human and animal dose-response data, and the current state of knowledge about the relative plausibility of alternative dose-response analyses. A tree (like a decision tree and a probability tree) is used to decompose the dose-response assessment into component factors, to provide a structure for explicitly considering multiple alternatives for each factor, and to explicitly incorporate the current state of knowledge about the relative plausibility of these alternatives. Groundbreaking work has demonstrated the feasibility of weight-of-evidence based distributional characterizations, and provided initial examples. Computer software implementations are also available.     

Evolution of testing strategies for genetic toxicity

Shortly following the inception of genetic toxicology as a distinct discipline within toxicology, questions arose regarding the type and number of tests needed to classify a chemical as a mutagenic hazard or as a potential carcinogen. To some degree the discipline separated into two sub-specialties, (1) genetic risk assessment and (2) cancer prediction since data from experimental oncology also supports the existence of a genotoxic step in tumor initiation. The issue of which and how many tests continued to be debated, but is now focused more tightly around two independent phenomena. Tier or sequential testing was initially proposed as a logical and cost-effective method, but was discarded on the basis that the lower tier tests appeared to have too many false responses to force or exclude further testing of the test agent. Matrix (battery) testing was proposed for screening on the hypothesis that combinations of endpoints and multiple phylogenetic target organisms were needed to achieve satisfactory predictability. As the results from short-term test 'validation' studies for carcinogen prediction and evaluations of EPA's Gene-Tox data accumulated, it became obvious that qualitative differences remained between predictive and definitive tests and by assembling different combinations of short-term assays investigators did not appear to resolve the lack of concordance. Recent trends in genetic toxicology testing have focused on mathematical models for test selection, and standardized systems for multi-test data assessment.     

A Combined Approach to Cartographic Displacement for Buildings Based on Skeleton and Improved Elastic Beam Algorithm

Scale reduction from source to target maps inevitably leads to conflicts of map symbols in cartography and geographic information systems (GIS). Displacement is one of the most important map generalization operators and it can be used to resolve the problems that arise from conflict among two or more map objects. In this paper, we propose a combined approach based on constraint Delaunay triangulation (CDT) skeleton and improved elastic beam algorithm for automated building displacement. In this approach, map data sets are first partitioned. Then the displacement operation is conducted in each partition as a cyclic and iterative process of conflict detection and resolution. In the iteration, the skeleton of the gap spaces is extracted using CDT. It then serves as an enhanced data model to detect conflicts and construct the proximity graph. Then, the proximity graph is adjusted using local grouping information. Under the action of forces derived from the detected conflicts, the proximity graph is deformed using the improved elastic beam algorithm. In this way, buildings are displaced to find an optimal compromise between related cartographic constraints. To validate this approach, two topographic map data sets (i.e., urban and suburban areas) were tested. The results were reasonable with respect to each constraint when the density of the map was not extremely high. In summary, the improvements include (1) an automated parameter-setting method for elastic beams, (2) explicit enforcement regarding the positional accuracy constraint, added by introducing drag forces, (3) preservation of local building groups through displacement over an adjusted proximity graph, and (4) an iterative strategy that is more likely to resolve the proximity conflicts than the one used in the existing elastic beam algorithm.     

COMBINED AND SEPARATE ANALYSES OF MORPHOLOGICAL AND MOLECULAR DATA IN THE PLANT FAMILY RUBIACEAE

Abstract— The Rubiaceae are one of the largest of the families of angiosperms, with over 10000 species. The tribal and subfamilial classification is provisional due to the lack of phylogenetic hypotheses. The present study of the Rubiaceae is based on 33 genera and three data sets, one morphological and two molecular from chloroplast DNA, restriction sites andrbcL sequences. There is much congruence between the morphological and the molecular data sets, but also conflict. For parsimony reasons, the best phylogenetic hypothesis is a tree based on an analysis of the combined data sets. The so-called “total evidence” criterion for the combined analysis is simply a reiteration of the principle of parsimony. In this particular study, the classification would be almost the same even if based on the separate analyses instead of the combined. Despite the inapplicability of consensus trees or trees from separate analyses for phylogenetic hypotheses and classification, separate analyses may provide important information. It is the best way to reveal conflicts between different data sets. Knowledge of the conflicts can promote further detailed investigation in order to improve understanding of characters and phylogenetic hypotheses. In this study, the tribe Vanguerieae provides such an example; morphological data support a position in the subfamily Cinchonoideae, but DNA and a tree based on the combined data support a position in subfamily Ixoroideae. The tribe's position in the morphological tree is probably due to missing information concerning the correct pollen presentation system. Bootstrap fractions and K. Bremer's branch support values are used to evaluate the stability of particular nodes in the trees. Interestingly these values are not always correlated, e.g. in the morphological tree, the node with the highest branch support value has very low bootstrap fraction. The reasons for these differences are unclear, but large differences are presumably more likely to occur on short branches.     

Why 'suitable' in vitro methods, as defined in the final EU REACH legislation, are an inappropriate basis for risk assessment

The final text of the REACH legislation, that has recently been published by the European Parliament, and which will come into force in the middle of 2007, implies that results obtained from so-called 'suitable' in vitro methods could be used for testing. It is also stated that a 'suitable' test is one that has been deemed to be ready to enter a validation study - for example, according to criteria specified by ECVAM. It is argued that, because the wording of the legal text is too vague and imprecise, it is totally unsatisfactory, mainly because it seems to endorse the application of non-validated methods for regulatory risk assessment. While their use might be suitable for providing information for classification in a weight-of-evidence approach, setting safe dose levels should only ever be based on the use of data from a properly validated test method. This concern, and other problems raised by the final legislation, is discussed, and recommendations are made which could avoid the premature use of potentially unsuitable methods for risk assessment, whilst ensuring that potentially suitable tests can be used for this purpose, once they have been formally validated.     

Exploring Component Stability Using Life-Stage Concordance In Sabethine Mosquitoes (Diptera: Culicidae)

Abstract — Morphological characters from sabethine mosquitoes were coded from larvae, pupae and adults, and life-stage partitions were evaluated to determine the contribution of each to the topology of a combined cladogram. Initial tests failed to find congruence between characters partitioned by life stage. However, when components from the combined analysis were tested using reduced taxon sets, a high degree of concordance between partitions was observed. A procedure for assessing individual life-stage contribution is employed, in which exhaustive searches are used to explore all possible arrangements for each of the selected components. Seven of the 10 components examined were able to recover the combined topology with a reduced taxon set. Congruent arrangements of taxa were typically observed for two or more life stages, although partitioned data were less resolved and frequently included aberrant topologies (those not supported by other partitioned or combined reduced taxon tree sets). In addition, none of the partitioned data sets gave robust results for all tests, suggesting that studies which emphasize character data from single life stages may support misleading arrangements of taxa. One component on the combined cladogram was not supported by any of the life-stage partitions when analysed separately. These results are complementary to “total evidence” approach, and demonstrate that partitions of data are useful for examining suits of characters which may cause some components of the “total     

A Pooled Response Strategy for Combining Multiple Lines of Evidence to Quantitatively Estimate Impact

The impacts of sediment contaminants can be evaluated by different lines of evidence, including toxicity tests and ecological community studies. Responses from 10 different toxicity assays/tests were combined to arrive at a “site score.” We employed a relatively simple summary measure, pooled P-values where we quantify a potential decrement in response in a contaminated site relative to nominally clean reference sites. The response-specific P-values were defined relative to a “null” distribution of responses in reference sites, and were then pooled using standard meta-analytic methods. Ecological community data were also evaluated using an analogous strategy. A distribution of distances of the reference sites from thecentroid of the reference sites was obtained. The distance from each of the test sites from the centroid of the reference sites was then calculated, and the proportion of reference distances that exceed the test site difference was used to define an empirical P-value for that test site. A plot of the toxicity P-value versus the community P-value was used to identify sites based on both alteration in community structure and toxicity, that is, by weight-of-evidence. This approach provides a useful strategy for examining multiple lines of evidence that should be accessible to the broader scientific community. The use of a large collection of reference sites to empirically define P-values is appealing in that parametric distribution assumptions are avoided, although this does come at the cost of assuming the reference sites provide an appropriate comparison group for test sites.     

Progress in achieving quantitative classification of psychopathology

Shortcomings of approaches to classifying psychopathology based on expert consensus have given rise to contemporary efforts to classify psychopathology quantitatively. In this paper, we review progress in achieving a quantitative and empirical classification of psychopathology. A substantial empirical literature indicates that psychopathology is generally more dimensional than categorical. When the discreteness versus continuity of psychopathology is treated as a research question, as opposed to being decided as a matter of tradition, the evidence clearly supports the hypothesis of continuity. In addition, a related body of literature shows how psychopathology dimensions can be arranged in a hierarchy, ranging from very broad “spectrum level” dimensions, to specific and narrow clusters of symptoms. In this way, a quantitative approach solves the “problem of comorbidity” by explicitly modeling patterns of co‐occurrence among signs and symptoms within a detailed and variegated hierarchy of dimensional concepts with direct clinical utility. Indeed, extensive evidence pertaining to the dimensional and hierarchical structure of psychopathology has led to the formation of the Hierarchical Taxonomy of Psychopathology (HiTOP) Consortium. This is a group of 70 investigators working together to study empirical classification of psychopathology. In this paper, we describe the aims and current foci of the HiTOP Consortium. These aims pertain to continued research on the empirical organization of psychopathology; the connection between personality and psychopathology; the utility of empirically based psychopathology constructs in both research and the clinic; and the development of novel and comprehensive models and corresponding assessment instruments for psychopathology constructs derived from an empirical approach.     

Is the “Ecological Fallacy” a Fallacy?

Ecological studies of health effects due to agent exposure are generally considered to be a blunt instrument of scientific investigation, unfit to determine the “true” exposure-effect relationship for an agent. Based on this widely accepted tenet, ecological studies of the correlation between the local air concentration of radon and the local lung cancer mortality as measured by Cohen have been criticized as being subject to the “Ecological Fallacy” and thus producing invalid risk data. Here we discuss the data that a risk assessment needs as a minimum requirement for making a valid risk estimate. The examination of these data and a “thought experiment” show that it is Cohen's raw ecological data, uncorrected for population characteristic factors, which are the proper data for a risk assessment. Consequently, the “true” exposure-effect relationship is less and less important the more population characteristic factors are identified and the larger they are. This reversal of the usual argument is due to our approach: Here, the prediction of the health effects in an exposed population is of primary importance, not the shape of the “true” exposure-effect relationship. The results derived in this paper hold for ecological studies of any agent causing any health or other effect.     

TOPIARY PRUNING AND WEIGHTING REINFORCE AN AFRICAN ORIGIN FOR THE HUMAN MITOCHONDRIAL DNA TREE

A method is presented for removing recent homoplastic events from a phylogenetic tree. This “topiary pruning” method produces a series of progressively modified duplicates of the original set of data, from which more and more of the most recent substitutions have been removed. The edited sets of data have increased amounts of information per remaining taxon, while similar but randomized data sets subjected to topiary pruning do not. The ability of topiary pruning to “unscramble” artificial data sets that have high levels of homoplasy is demonstrated, and is shown to be similar in its effects to the weighting method of Kluge and Farris (1969), although with the additional advantage of reducing the number of taxa to the point where bootstrapping is feasible. Pruning and weighting used together produce closer approximations to the “true” tree than either method used separately. It is further shown that in these artificial data sets midpoint rooting is more likely to be accurate than outgroup rooting. When pruning and weighting are applied to the extensive sets of mitochondrial DNA data of Cann et al. (1987) and Vigilant et al. (1991), trees result that have deep branch points, some of which lead to entirely African branches. In the case of the Vigilant et al. data, the three African branches have bootstrap values between 0.94 and 1.0, and the consensus and bootstrap midpoint roots also have high bootstrap values and occur on these African branches near their junction. An African origin of the human mitochondrial tree is not proved by this approach, particularly since sequences from non-African groups are underrepresented in current data sets, but it is rendered more likely.     

The Apache Longbow–Hellfire Missile Test at Yuma Proving Ground: Ecological Risk Assessment for Missile Firing

A multiple stressor risk assessment was conducted at Yuma Proving Ground, Arizona, as a demonstration of the Military Ecological Risk Assessment Framework. The focus was a testing program at Cibola Range, which involved an Apache Longbow helicopter firing Hellfire missiles at moving targets, that is, M60-A1 tanks. This article describes the ecological risk assessment for the missile launch and detonation. The primary stressor associated with this activity was sound. Other minor stressors included the detonation impact, shrapnel, and fire. Exposure to desert mule deer (Odocoileus hemionus crooki) was quantified using the Army sound contour program BNOISE2, as well as distances from the explosion to deer. Few effects data were available from related studies. Exposure-response models for the characterization of effects consisted of human “disturbance” and hearing damage thresholds in units of C-weighted decibels (sound exposure level) and a distance-based No-Observed-Adverse-Effects Level for moose and cannonfire. The risk characterization used a weight-of-evidence approach and concluded that risk to mule deer behavior from the missile firing was likely for a negligible number of deer, but that no risk to mule deer abundance and reproduction is expected.     

Integrated testing strategies for toxicity employing new and existing technologies

We have developed individual, integrated testing strategies (ITS) for predicting the toxicity of general chemicals, cosmetics, pharmaceuticals, inhaled chemicals, and nanoparticles. These ITS are based on published schemes developed previously for the risk assessment of chemicals to fulfil the requirements of REACH, which have been updated to take account of the latest developments in advanced in chemico modelling and in vitro technologies. In addition, we propose an ITS for neurotoxicity, based on the same principles, for incorporation in the other ITS. The technologies are deployed in a step-wise manner, as a basis for decision-tree approaches, incorporating weight-of-evidence stages. This means that testing can be stopped at the point where a risk assessment and/or classification can be performed, with labelling in accordance with the requirements of the regulatory authority concerned, rather than following a checklist approach to hazard identification. In addition, the strategies are intelligent, in that they are based on the fundamental premise that there is no hazard in the absence of exposure - which is why pharmacokinetic modelling plays a key role in each ITS. The new technologies include the use of complex, three-dimensional human cell tissue culture systems with in vivo-like structural, physiological and biochemical features, as well as dosing conditions. In this way, problems of inter-species extrapolation and in vitro/in vivo extrapolation are minimised. This is reflected in the ITS placing more emphasis on the use of volunteers at the whole organism testing stage, rather than on existing animal testing, which is the current situation.     

Assessment of a large number of empirical plant species niche models by elicitation of knowledge from two national experts

Quantitative models play an increasing role in exploring the impact of global change on biodiversity. To win credibility and trust, they need validating. We show how expert knowledge can be used to assess a large number of empirical species niche models constructed for the British vascular plant and bryophyte flora. Key outcomes were (a) scored assessments of each modeled species and niche axis combination, (b) guidance on models needing further development, (c) exploration of the trade‐off between presenting more complex model summaries, which could lead to more thorough validation, versus the longer time these take to evaluate, (d) quantification of the internal consistency of expert opinion based on comparison of assessment scores made on a random subset of models evaluated by both experts. Overall, the experts assessed 39% of species and niche axis combinations to be “poor” and 61% to show a degree of reliability split between “moderate” (30%), “good” (25%), and “excellent” (6%). The two experts agreed in only 43% of cases, reaching greater consensus about poorer models and disagreeing most about models rated as better by either expert. This low agreement rate suggests that a greater number of experts is required to produce reliable assessments and to more fully understand the reasons underlying lack of consensus. While area under curve (AUC) statistics showed generally very good ability of the models to predict random hold‐out samples of the data, there was no correspondence between these and the scores given by the experts and no apparent correlation between AUC and species prevalence. Crowd‐sourcing further assessments by allowing web‐based access to model fits is an obvious next step. To this end, we developed an online application for inspecting and evaluating the fit of each niche surface to its training data.     

An integrated environment and cost assessment method based on LCA and LCC for mechanical product manufacturing

Purpose

The purpose of this study is to provide an integrated method to identify the resource consumption, environmental emission, and economic cost for mechanical product manufacturing from economic and ecological dimensions and ultimately to provide theoretical and data support of energy conservation and emission reduction for mechanical product manufacturing.

Methods

The applied research methods include environmental life cycle assessment (LCA) and life cycle cost (LCC). In life cycle environmental assessment, the inventory data are referred from Chinese Life Cycle Database and midpoint approach and EDIP2003 and CML2001 models of life cycle impact assessment (LCIA) are selected. In life cycle cost assessment, three cost categories are considered. The proposed environment and cost assessment method is based on the theory of social willingness to pay for potential environmental impacts. With the WD615 Steyr engine as a case, life cycle environment and cost are analyzed and evaluated.

Results and discussion

The case study indicates that, in different life cycle phases, the trend of cost result is generally similar to the environmental impacts; the largest proportion of cost and environmental impact happened in the two phases of “material production” and “component manufacturing” and the smallest proportion in “material transport” and “product assembly.” The environmental impact category of Chinese resource depletion potential (CRDP) accounted for the largest proportion, followed by global warming potential (GWP) and photochemical ozone creation potential (POCP), whereas the impacts of eutrophication potential (EP) and acidification potential (AP) are the smallest. The life cycle “conventional cost” accounted for almost all the highest percentage in each phase (except “material transport” phase), which is more than 80% of the total cost. The “environmental cost” and “possible cost” in each phase are relatively close, and the proportion of which is far below the “conventional cost.”

Conclusions

The proposed method enhanced the conventional LCA. The case results indicate that, in a life cycle framework, the environment and cost analysis results could support each other, and focusing on the environment and cost analysis for mechanical product manufacturing will contribute to a more comprehensive eco-efficiency assessment. Further research on the life cycle can be extended to phases of “early design,” “product use,” and “final disposal.” Other LCIA models and endpoint indicators are advocated for this environmental assessment. Environmental cost can also be further investigated, and the relevant social willingness to pay for more environmental emissions is advocated to be increased.

     

Morphometric, biochemical and molecular traits in Caucasian wood mice (podemus/Sylvaemus), with remarks on species divergence

We analysed Caucasian wood mice from Georgia (n = 60) and supplementary reference material of theApodemus/Sylvaemus species group to evaluate the reliability of taxon identification. Traditional “expert knowledge” plus three different methodological approaches were employed and combined to perceive their discriminatory power for a reliable taxon assignment. Graphs of principal component scores derived from the analysis of 14 skull metrics displayed taxon membership of individuals. Individual multi--locus (L = 18) electrophoretic profiles were used to re-assess specimens to a specific genepool by an assignment test based on allele frequencies indicative of populational taxon samples of the respective sampling locations. Genotyped individuals were re-allocated to those taxa, for which they yielded the highest probability score. Genetic distances among the taxa were computed and clustered in a neighbour-joining tree. PCR-fragments of 1074bp amplified from the mitochondrial cytochromeb gene were cut with 2 six- and 4 four-cutter restriction enzymes, and resulting RFLP patterns were analysed phenetically to classify the specimens according to their molecular similarity. Partial cytochromeb sequences were used to construct a phylogenetic tree by computing neighbour-joining clusters from a matrix of percent nucleotide differences. The power of the combined classification approaches and their congruence is discussed. It is concluded that the joint application of traditional, morphometric and biochemical or genetic techniques for taxon allocation of specimens of wood mice encountered problems in species delimitation. The mtDNA topology obtained was not congruent with protein polymorphism that indicated differential historical and/or recent introgression and incomplete lineage sorting in substructured populations. Cytochromeb sequence DNA data analysed were not as adequate as expected to resolve phylogenetic relationships among Caucasian and European members of theApodemus-Sylvaemus complex. Altogether, morphometric, biochemical and sequence data sets did not support the hypothesis of the evolutionary independence of European and Caucasian lineages of wood mice. Nonetheless, extended combined morphological and genetic analyses are considered necessary prerequisites to an in-depth study of the evolutionary lineages of theApodemus/Sylvaemus group. More sequence data of a variety of genes (and plenty of nuclear markers) are needed to resolve the various levels of differentiation of the extant lineages.     

Relevance and follow-up of positive results in in vitro genetic toxicity assays: an ILSI-HESI initiative

In vitro genotoxicity assays are often used to screen and predict whether chemicals might represent mutagenic and carcinogenic risks for humans. Recent discussions have focused on the high rate of positive results in in vitro tests, especially in those assays performed in mammalian cells that are not confirmed in vivo. Currently, there is no general consensus in the scientific community on the interpretation of the significance of positive results from the in vitro genotoxicity assays. To address this issue, the Health and Environmental Sciences Institute (HESI), held an international workshop in June 2006 to discuss the relevance and follow-up of positive results in in vitro genetic toxicity assays. The goals of the meeting were to examine ways to advance the scientific basis for the interpretation of positive findings in in vitro assays, to facilitate the development of follow-up testing strategies and to define criteria for determining the relevance to human health. The workshop identified specific needs in two general categories, i.e., improved testing and improved data interpretation and risk assessment. Recommendations to improve testing included: (1) re-examine the maximum level of cytotoxicity currently required for in vitro tests; (2) re-examine the upper limit concentration for in vitro mammalian studies; (3) develop improved testing strategies using current in vitro assays; (4) define criteria to guide selection of the appropriate follow-up in vivo studies; (5) develop new and more predictive in vitro and in vivo tests. Recommendations for improving interpretation and assessment included: (1) examine the suitability of applying the threshold of toxicological concern concepts to genotoxicity data; (2) develop a structured weight of evidence approach for assessing genotoxic/carcinogenic hazard; and (3) re-examine in vitro and in vivo correlations qualitatively and quantitatively. Conclusions from the workshop highlighted a willingness of scientists from various sectors to change and improve the current paradigm and move from a hazard identification approach to a "realistic" risk-based approach that incorporates information on mechanism of action, kinetics, and human exposure..     

Differences among Yanomama Indian villages: Do the patterns of allele frequencies,anthropometrics and map locations correspond?

In order to determine the degree of correspondence between sets of multivariate observations based on different kinds of traits, two new methods, derived from fundamentally different notions of “correspondence,” are adopted here and compared. Using networks or trees to represent contemporary relationships, the first method tests the similarity of the cluster or hierarchic structures implicit in two sets of data. The second approach tests the departure from perfect geometric congruence or superimposability. Computer simulation was used to generate the distributions needed for significance tests under the null hypothesis. By the first technique, we find significant correspondence among the cluster structures for geographic, allele frequency, and anthropometric data on 19 Yanomama Indian villages. The results are similar and more precise for a subset consisting of seven villages. Some of these results differ from the conclusions which would be reached with the conventional correlations based upon entries in distance tables. The direct test of congruence, used only for the data on the subset of seven villages, gives results which differ substantially from those based on cluster-structure. There are, however, similarities between the measure of congruence and the simple correlations based on entries in the distance tables. The significant correspondences observed call for some explanation. Cultural and demographic features determine the particular non-random allocation of individuals to village fragments when a village splits. These social phenomena are invoked in tentative explanation of the agreement among historical, biological, and geographic relationships of villages.