HOW SIMILAR ARE GENETIC CORRELATION STRUCTURES? DATA FROM MICE AND RATS

An integral assumption of many models of morphometric evolution is the equality of the genetic variance-covariance structure across evolutionary time. To examine this assumption, the quantitative-genetic aspects of morphometric form are examined for eight pelvic traits in laboratory rats (Rattus norvegicus) and random-bred ICR mice (Mus musculus). In both species, all traits are significantly heritable, and there are significant phenotypic and genetic correlations among traits, although environmental correlations among the eight traits are low. The size relations among the pelvic variables are isometric. Three matrix-permutation tests are used to examine similarity of phenotypic, genetic, and environmental covariance and correlation matrices within and between species. Independent patterns of morphometric covariation and correlation arise from genetic and environmental effects within each species and from environmental effects between species. The patterns of phenotypic and genetic covariation and correlation are similar within each species, and the phenotypic and genetic correlations are also similar between these species. However, genetic covariance matrices show no significant statistical association between species. It is suggested that the assumption of equality of genetic variance-covariance structures across divergent taxa should be approached with caution.     

The stability of P in coral reef fishes

The constancy of phenotypic variation and covariation is an assumption that underlies most recent investigations of past selective regimes and attempts to predict future responses to selection. Few studies have tested this assumption of constancy despite good reasons to expect that the pattern of phenotypic variation and covariation may vary in space and time. We compared phenotypic variance-covariance matrices (P) estimated for populations of six species of distantly related coral reef fishes sampled at two locations on Australia's Great Barrier Reef separated by more than 1000 km. The intraspecific similarity between these matrices was estimated using two methods: matrix correlation and common principal component analysis. Although there was no evidence of equality between pairs of P, both statistical approaches indicated a high degree of similarity in morphology between the two populations for each species. In general, the hierarchical decomposition of the variance-covariance structure of these populations indicated that all principal components of phenotypic variance-covariance were shared but that they differed in the degree of variation associated with each of these components. The consistency of this pattern is remarkable given the diversity of morphologies and life histories encompassed by these species. Although some phenotypic instability was indicated, these results were consistent with a generally conserved pattern of multivariate selection between populations.     

The presence–absence matrix reloaded: the use and interpretation of range–diversity plots

Aim A great deal of information on distribution and diversity can be extracted from presence–absence matrices (PAMs), the basic analytical tool of many biogeographic studies. This paper presents numerical procedures that allow the analysis of such information by taking advantage of mathematical relationships within PAMs. In particular, we show how range–diversity (RD) plots summarize much of the information contained in the matrices by the simultaneous depiction of data on distribution and diversity. Innovation We use matrix algebra to extract and process data from PAMs. Information on the distribution of species and on species richness of sites is computed using the traditional R (by rows) and Q (by columns) procedures, as well as the new Rq (by rows, considering the structure of columns) and Qr (by columns, considering the structure by rows) methods. Matrix notation is particularly suitable for summarizing complex calculations using PAMs, and the associated algebra allows the implementation of efficient computational programs. We show how information on distribution and species richness can be depicted simultaneously in RD plots, allowing a direct examination of the relationship between those two aspects of diversity. We explore the properties of RD plots with a simple example, and use null models to show that while parameters of central tendency are not affected by randomization, the dispersion of points in RD plots does change, showing the significance of patterns of co‐occurrence of species and of similarity among sites. Main conclusion Species richness and range size are both valid measures of diversity that can be analysed simultaneously with RD plots. A full analysis of a system requires measures of central tendency and dispersion for both distribution and species richness.     

Goodness-of-fit in generalized nonlinear mixed-effects models

In recent years, generalized linear and nonlinear mixed-effects models have proved to be powerful tools for the analysis of unbalanced longitudinal data. To date, much of the work has focused on various methods for estimating and comparing the parameters of mixed-effects models. Very little work has been done in the area of model selection and goodness-of-fit, particularly with respect to the assumed variance-covariance structure. In this paper, we present a goodness-of-fit statistic which can be used in a manner similar to the R2 criterion in linear regression for assessing the adequacy of an assumed mean and variance-covariance structure. In addition, we introduce an approximate pseudo-likelihood ratio test for testing the adequacy of the hypothesized convariance structure. These methods are illustrated and compared to the usual normal theory likelihood methods (Akaike's information criterion and the likelihood ratio test) using three examples. Simulation results indicate the pseudo-likelihood ratio test compares favorably with the standard normal theory likelihood ratio test, but both procedures are sensitive to departures from normality.     

HOMOGENEITY OF THE GENETIC VARIANCE-COVARIANCE MATRIX FOR ANTIPREDATOR TRAITS IN TWO NATURAL POPULATIONS OF THE GARTER SNAKE THAMNOPHIS ORDINOIDES

Quantitative genetic models of evolution rely on the genetic variance-covariance matrix to predict the phenotypic response to selection. Both prospective and retrospective studies of phenotypic evolution across generations rely on assumptions about the constancy of patterns of genetic covariance through time. In the absence of robust theoretical predictions about the stability of genetic covariances, this assumption must be tested with empirical comparisons of genetic parameters among populations and species. Genetic variance-covariance matrices were estimated for a suite of antipredator traits in two populations of the northwestern garter snake, Thamnophis ordinoides. The characters studied include color pattern and antipredator behaviors that interact to facilitate escape from predators. Significant heritabilities for all traits were detected in both populations. Genetic correlations and covariances were found among behaviors in both populations and between color pattern and behavior in one of the populations. Phenotypic means differed among populations, but pairwise comparisons revealed no heterogeneity of genetic parameters between the populations. The structure of the genetic variance-covariance matrix has apparently not changed significantly during the divergence of these two populations.     

VARIATION IN GENETIC ARCHITECTURE OF CALLING SONG AMONG POPULATIONS OF ALLONEMOBIUS SOCIUS,A. FASCIATUS,AND A HYBRID POPULATION: DRIFT OR SELECTION?

Predictions using quantitative genetic models generally assume that the variance-covariance matrices remain constant over time. This assumption is based on the supposition that selection is generally weak and hence variation lost through selection can be replaced by new mutations. Whether this is generally true can only be ascertained from empirical studies. Ideally for such a study we should be able to make a prediction concerning the relative strength of selection versus genetic drift. If the latter force is prevalent then the variance-covariances matrices should be proportional to each other. Previous studies have indicated that females in the two sibling cricket species Allonemobius socius and A. fasciatus do not discriminate between males of the two species by their calling song. Therefore, differences between the calling song of the two males most likely result from drift rather than sexual selection. We test this hypothesis by comparing the genetic architecture of calling song of three populations of A. fasciatus with two populations of A. socius. We found no differences among populations within species, but significant differences in the G (genetic) and P (phenotypic) matrices between species, with the matrices being proportional as predicted under the hypothesis of genetic drift. Because of the proportional change in the (co)variances no differences between species are evident in the heritabilities or genetic correlations. Comparison of the two species with a hybrid population from a zone of overlap showed highly significant nonproportional variation in genetic architecture. This variation is consistent with a general mixture of two separate genomes or selection. Qualitative conclusions reached using the phenotypic matrices are the same as those reached using the genetic matrices supporting the hypothesis that the former may be used as surrogate measures of the latter.     

Purification of cadmium-binding proteins from related species of terrestrial helicidae (gastropoda,mollusca): A comparative study

Summary Three species of terrestrial Helicidae (Helix pomatia, Cepaea hortensis andArianta arbustorum) were fed cadmium-rich diet in the laboratory. The snails accumulated high amounts of the metal in their hepatopancreas. Most cadmium and some zinc were found, after centrifugation, in the soluble fractions from which a cadmium-binding protein was isolated for each species by ion exchange and gel chromatography. The proteins contained different amounts of cadmium, but little or no zinc, and showed high absorption at 254 nm indicating the presence of cadmium-mercaptide bonds. After gel filtration, a molecular weight of 12000 was found for cadmium-binding proteins fromHelix pomatia andArianta arbustorum, whereas a molecular weight of 10 000 was found for a cadmium-binding protein fromCepaea hortensis. SDS-polyacrylamide gel electrophoresis showed one single band for each protein fromHelix pomatia andArianta arbustorum and suggested a molecular weight of 11000 for both species. Amino acid analysis revealed, for each protein, high amounts of cysteine (12–20%), glycine (15–19%), and serine (12–14%), and moderately elevated contents of lysine (9–13%) and alanine (4–8%), but no methionine and only traces, if any, of aromatic amino acids. The ratios of cadmium to cysteine were 1:5, 1:10 and 1:3 in the proteins fromHelix pomatia,Cepaea hortensis andArianta arbustorum, respectively. Some features of the isolated proteins resembled mammalian metallothioneins. Most characteristics, however, differed from true metallothioneins and were similar to cadmium-binding proteins found in some marine molluscs.     

Morphometrics of the callitrichid forelimb: A case study in size and shape

Mosimann and colleagues formulated a technique that distinguishes between size and shape, based on the concept of geometric similarity and the distinction between log size-and-shape and log shape variables. We extend these formulations in an examination of the forelimb of three callitrichid species (adultSaguinus oedipus, Saguinus fuscicollis, andCallithrix jacchus). We employ principal components analysis to explore the relationship between variance explained by size-and-shape versus shape alone. Independence of shape vectors is examined via correlation analysis. Then we use log shape data to construct intersample (species means) and total sample (between all paris of individuals) matrices of average taxonomic distances. These distance matrices are subjected to cluster analysis and principal coordinate ordinations. Results of principal components analysis suggest that after isometric size is removed, there remains sufficient shape information to discriminate among the three taxa. Careful examination and quantification of the relationships between shape and size suggest that size information (e.g., geometric mean) is fundamental for understanding shape differences within and among callitrichid species; in other words, most aspects of forelimb shape are significantly correlated with size. Contrary to conventional wisdom, we also demonstrate that such correlations are not spurious. Ordinations and clustering of log shape distance matrices (based on means and individuals) support the notion that, despite differences in size, the two tamarins are more similar in shape than either is toC. jacchus (despite size similarity betweenS. fuscicollis andC. jacchus). Although shape variation in the forelimb of calliirichids may have a functional component, the phylogenetic signal remains strong and serves to group individuals accordingly.     

Identification of bacterial strains by laboratories participating in the Deutsches Krebsforschungszentrum quality assurance programme

The quality assurance programme (QAP) of the Deutsches Krebsforschungszentrum (DKFZ) is a proficiency testing system developed to service the laboratory animal discipline. QAP comprises the quarterly distribution of two bacterial strains originating from various species of animals for identification to the species level and antibiotic susceptibility testing. We compared identification results reported by QAP participants over the years 1996-2004 with those obtained by the Dutch Bacterial Diagnostics reference laboratory on 68 samples comprising 71 bacterial strains and a fungus. Significant differences were found in the frequency of reported and correct identifications when bacteria were assigned to different groups based on morphology by Gram stain and on origin (animal versus environmental, rodent and rabbit versus other animal species, pathogen versus non-pathogens). Rodent and rabbit pathogens yielded 73% correct identifications, and with all bacterial strains only 60% of the identifications were correct. We assume that most QAP participants were from laboratory animal diagnostic laboratories. If this is true, the capabilities of laboratories in the laboratory animal discipline to correctly identify bacterial species are well below what are considered acceptable limits for human diagnostic laboratories. The distribution of cultured bacteria circumvents the most difficult step in the microbiological monitoring of animals, namely primary culture from clinical samples. We propose to set up a QAP that comprises the distribution of specimens mimicking clinical samples normally submitted to laboratory animal diagnostic laboratories.     

Embryology of the Turbellaria and its phylogenetic significance

Developmental characters — including oocyte and yolk cell structure, patterns of cleavage, and modes of gastrulation — are presented and examined in relation to the phylogeny of the Turbellaria. Eggshell granules, which have been demonstrated to occur in the oocytes of entolecithal eggs and the yolk cells of ectolecithal eggs, are compared among species, and their potential value as a taxonomic character is discussed. The quartet 4d spiral cleavage of the entolecithal egg of polyclads is described as reminiscent of the primitive pattern of early development for the Turbellaria. This is compared to duet spiral cleavage of acoels, and possible phylogenetic schemes involving the two types of spiral cleavage are reviewed. The link between the precise spiral cleavage, which characterizes development of most archoophorans, and blastomere separation (Blastomeren-Anarchie), which occurs in several neoophoran orders, is established by the occurrence of quartet 4d spiral cleavage in one neoophoran order, and of both quartet spiral cleavage and Blastomeren-Anarchie in different species of a second neoophoran order. The epibolic gastrulation of polyclads is described as primitive for the Turbellaria because of its similarity to that of other members of the Spiralia. Although no identical process occurs in neoophoran development, the earlier event of formation of the hull membrane in some neoophorans, and the later event of formation of the definitive epidermis in all neoophorans studied are presented as processes of possible homology to the epibolic gastrulation of polyclads. The lack of correspondence between polyclads and neoophorans in the relationship of the definitive body axes to the egg axis is discussed, and an hypothesis is advanced to account for the differences. The phylogenetic relationships indicated by known developmental phenomena differ only slightly from the scheme presented by Karling in 1974.     

QUANTITATIVE GENETICS OF MORPHOLOGICAL DIFFERENTIATION IN PEROMYSCUS. I. TESTS OF THE HOMOGENEITY OF GENETIC COVARIANCE STRUCTURE AMONG SPECIES AND SUBSPECIES

Phenotypic and additive genetic covariance matrices were estimated for 15 morphometric characters in three species and subspecies of Peromyscus. Univariate and multivariate ANOVAs indicate these groups are highly diverged in all characters, P. leucopus having the largest body size, P. maniculatus bairdii the smallest, and P. maniculatus nebrascensis being intermediate. Comparing the structure of P and G within each taxon revealed significant similarities in all three cases. This proportionality was strong enough to justify using P in the place of G to analyze evolutionary processes using quantitative genetic models when G can not be estimated, as in fossil material. However, the similarity between genetic and phenotypic covariance structures is sufficiently low that estimates of the genetic parameters should be used when possible. The additive genetic covariance matrices were compared to examine the assumption that they remain constant during evolution, an assumption which underlies many applications of quantitative-genetic models. While matrix permutation tests indicated statistically significant proportionality between the genetic covariance structures of the two P. maniculatus subspecies, there is no evidence of significant genetic structural similarity between species. This result suggests that the assumption of constant genetic covariance structure may be valid only within species. (It does not, however, necessarily imply a causal relationship between speciation and heterogeneity of genetic covariance structures.) The low matrix correlation for the two P. maniculatus subspecies' genetic covariance matrices indicates G may not be functionally constant, even within species. The lack of similarity observed here may be due partly to sampling variation.     

Numerical taxonomy of some taxa of the genus Fucus in the Iberian Peninsula

A study of morphological variation patterns in three Iberian Fucus species has been carried out. Data on 36 features from 113 specimens have been studied by numerical methods (PCA, discriminant analysis and ANOVA).The PCA ordination displayed a continuous trend of variation in the following Fucus species: F. ceranoides, F. vesiculosus and F. spiralis. However, no obvious morphological discontinuities have been observed. In addition, three varieties of F. spiralis and five of F. vesiculosus have been identified.At the varietal level discriminant functions could be used for identifying 93–100% of the individuals within F. vesiculosus and 72–91% of F. spiralis.The most discriminating diagnostic character appeared to be the midrib size, as detected by ANOVA procedures.Some qualitative characters, such as presence of vesicles, dioecious-monoecious character and receptacle edge shape, frond and receptacle shape, have been also important for discriminating specific and infraspecific taxa.A nomenclatural discussion of varietal taxa is also included.     

Ototyphlonemertes correae sp. nov. and a redescription of O. duplex (Nemertea: Monostilifera: Ototyphlonemertidae), with a phylogenetic consideration of the genus

The species Ototyphlonemertes correae sp. nov. and O. duplex Bürger, 1895 (Ototyphlonemertidae) are described. The new species possesses a smooth stylet, oligogranular statolith and nephridia. Phylogenetic relationships are reconstructed for 21 Ototyphlonemertes species, based on an analysis of 18 morphological characters, totalling 43 character states. The most informative characters are stylet structure, basis shape, statolith shape and nephridia. Based on the phylogeny, the subgenus Otohelicophora (Greek, helicophora meaning the spiral carriers) is separated, which is defined with reference to the synapomorphy 'spiral stylet'.     

Comparative sperm ultrastructure in Nemertea

Although the monophyly of Nemertea is strongly supported by unique morphological characters and results of molecular phylogenetic studies, their ingroup relationships are largely unresolved. To contribute solving this problem we studied sperm ultrastructure of 12 nemertean species that belong to different subtaxa representing the commonly recognized major monophyletic groups. The study yielded a set of 26 characters with an unexpected variation among species of the same genus (Tubulanus and Procephalothrix species), whereas other species varied in metric values or only one character state (Ramphogordius). In some species, the sperm nucleus has grooves (Zygonemertes virescens, Amphiporus imparispinosus) that may be twisted and give a spiral shape to the sperm head (Paranemertes peregrina, Emplectonema gracile). To make the characters from sperm ultrastructure accessible for further phylogenetic analyses, they were coded in a character matrix. Published data for eight species turned out to be sufficiently detailed to be included. Comparative evaluation of available information on the sperm ultrastructure suggests that subtaxa of Heteronemertea and Hoplonemertea are supported as monophyletic by sperm morphology. However, the data do not provide information on the existing contradictions regarding the internal relationships of “Palaeonemertea.” Nevertheless, our study provides evidence that sperm ultrastructure yields numerous potentially informative characters that will be included in upcoming phylogenetic analyses. J. Morphol. 2010. © 2010 Wiley‐Liss, Inc.     

Multiple regression on distance matrices: a multivariate spatial analysis tool

I explore the use of multiple regression on distance matrices (MRM), an extension of partial Mantel analysis, in spatial analysis of ecological data. MRM involves a multiple regression of a response matrix on any number of explanatory matrices, where each matrix contains distances or similarities (in terms of ecological, spatial, or other attributes) between all pair-wise combinations of n objects (sample units); tests of statistical significance are performed by permutation. The method is flexible in terms of the types of data that may be analyzed (counts, presence–absence, continuous, categorical) and the shapes of response curves. MRM offers several advantages over traditional partial Mantel analysis: (1) separating environmental distances into distinct distance matrices allows inferences to be made at the level of individual variables; (2) nonparametric or nonlinear multiple regression methods may be employed; and (3) spatial autocorrelation may be quantified and tested at different spatial scales using a series of lag matrices, each representing a geographic distance class. The MRM lag matrices model may be parameterized to yield very similar inferences regarding spatial autocorrelation as the Mantel correlogram. Unlike the correlogram, however, the lag matrices model may also include environmental distance matrices, so that spatial patterns in species abundance distances (community similarity) may be quantified while controlling for the environmental similarity between sites. Examples of spatial analyses with MRM are presented.     

Can tooth differentiation help to understand species coexistence? The case of wood mice in China

Five wood mice Apodemus species occur across China, in allopatry but also in sympatry up to cases of syntopy. They all share a similar external appearance, similar habitats of grasslands and forests and a generalist feeding behaviour. This overall similarity raises questions about the mechanisms insuring competition avoidance and allowing the coexistence of the species. In this context, a morphometric analysis of two characters related to feeding (mandible and molar) addressed the following issues: (1) Were the species actually different in size and/or shape of these characters, supporting their role in resource partitioning? (2) Did this pattern of phenotypic divergence match the neutral genetic differentiation, suggesting that differentiation might have occurred in a former phase of allopatry as a result of stochastic processes? (3) Did the species provide evidence of character displacement when occurring in sympatry, supporting an ongoing role of competition in the interspecific divergence? Results evidenced first that different traits, here mandibles and molars, provided discrepant pictures of the evolution of the Apodemus group in China. Mandible shape appeared as prone to vary in response to local conditions, blurring any phylogenetic or ecological pattern, whereas molar shape evolution appeared to be primarily driven by the degree of genetic differentiation. Molar size and shape segregated the different species in the morphospace, suggesting that these features may be involved in a resource partitioning between Apodemus species. The morphological segregation of the species, likely achieved by processes of differentiation in isolation promoted by the complex landscape of China, could contribute to competition avoidance and hence explain why no evidence was found of character displacement.     

HIERARCHICAL COMPARISON OF GENETIC VARIANCE-COVARIANCE MATRICES. I. USING THE FLURY HIERARCHY

The comparison of additive genetic variance-covariance matrices (G-matrices) is an increasingly popular exercise in evolutionary biology because the evolution of the G-matrix is central to the issue of persistence of genetic constraints and to the use of dynamic models in an evolutionary time frame. The comparison of G-matrices is a nontrivial statistical problem because family structure induces nonindependence among the elements in each matrix. Past solutions to the problem of G-matrix comparison have dealt with this problem, with varying success, but have tested a single null hypothesis (matrix equality or matrix dissimilarity). Because matrices can differ in many ways, several hypotheses are of interest in matrix comparisons. Flury (1988) has provided an approach to matrix comparison in which a variety of hypotheses are tested, including the two extreme hypotheses prevalent in the evolutionary literature. The hypotheses are arranged in a hierarchy and involve comparisons of both the principal components (eigenvectors) and eigenvalues of the matrix. We adapt Flury's hierarchy of tests to the problem of comparing G-matrices by using randomization testing to account for nonindependence induced by family structure. Software has been developed for carrying out this analysis for both genetic and phenotypic data. The method is illustrated with a garter snake test case.     

Geometric and traditional morphometrics for the assessment of character state identity: multivariate statistical analyses of character variation in the genus Arrenurus (Acari,Hydrachnidia, Arrenuridae)

Geometric and traditional morphometric approaches are tested to describe and reveal taxonomic characters and character states in variation of shape and size of idiosoma morphology in the Arrenuridae. Patterns of variation of idiosoma and glandularia features of males from 11 Mexican species of Arrenurus (Megaluracarus) and two species of subgenus Dadayella were explored with five landmark configurations and three sets of interlandmark distances. Separate principal component analyses (PCA) and canonical variate analyses (CVA) were performed for each data set. The eight multivariate analyses of variance among 13 a priori groups (species) detected significantly different morphometric variants, which were interpreted as different taxonomic character states. Patterns of character state similarity among species were examined with unweighted‐pair grouping method using averages (UPGMA) cluster analyses on Mahalanobis distances. Analyses of five landmark configurations revealed important taxonomical variation in the anterior idiosoma outline (10 character states), the outline of the posterior region or cauda (13 states), the distribution of postocularia, and the second and third pairs of dorsoglandularia (nine states), the fourth pair of dorsoglandularia (three states), and ventroglandularia on the posterior side of idiosoma (nine character states). Multivariate analyses of three sets of distance measurements also resulted in the detection of potential taxonomic characters related to idiosoma size (12 character states), postocularia and dorsoglandularia (13 states), and ventroglandularia (nine character states). Morphometric analyses of distances and shapes provide a formal basis for the interpretation of taxonomic characters, and for the discovery of character states. These characters should be investigated further in a wider sample of species for the phylogenetic systematics of these water mites. In the meantime, idiosoma regions and structures were tested for congruence in a phylogenetic analysis, and were proposed as homologous among the species sampled.     

The comparative anatomy of the abdominal gastrointestinal tract of six species of African mole‐rats (Rodentia,Bathyergidae)

The gastrointestinal tracts (GITs) of six species of African mole‐rats (Bathyergidae) were compared. The aim was to provide a comprehensive anatomical comparison between the different species. The relative shape, length, and surface areas were taken into account to determine whether the GITs are phylogenetically constrained or exhibit anatomical adaptations in response to diets. In all six species the stomach was simple and glandular. With the exception of Heterocephalus glaber, the caecum was coiled in a flat spiral, the ascending colon was arranged in a loop of varying lengths, and a mucosal colonic papillary‐lined groove was present in the ascending colon in all species. By contrast, the caecum in H. glaber was uncoiled, the ascending colon was not looped, and the groove was not papillated. A caeco‐appendix was observed only in Bathyergus suillus and Georychus capensis. Hierarchical multivariate cluster analysis on the presence/absence of nine anatomical structures associated with the GIT of mole‐rats revealed that H. glaber was anatomically the least similar of the six species (77.6% similarity) with respect to the nine GIT variables included. All Cryptomys species were the same (100% similarity), and two species B. suillus and G. capensis grouped together and were more similar to the Cryptomys genus (95% similarity) than they were to H. glaber. These findings support previous phylogenetic classifications. The voluminous caeco‐colon in B. suillus may be explained by its ingestion of grasses in addition to below‐ground storage organs of plants. We conclude that phylogeny and diet affect the GIT anatomy of the African mole rats studied here. J. Morphol., 2010. © 2009 Wiley‐Liss, Inc.