A Monte-Carlo algorithm for maximum likelihood estimation of variance components

International audienc

Primate modularity and evolution: first anatomical network analysis of primate head and neck musculoskeletal system

Network theory is increasingly being used to study morphological modularity and integration. Anatomical network analysis (AnNA) is a framework for quantitatively characterizing the topological organization of anatomical structures and providing an operational way to compare structural integration and modularity. Here we apply AnNA for the first time to study the macroevolution of the musculoskeletal system of the head and neck in primates and their closest living relatives, paying special attention to the evolution of structures associated with facial and vocal communication. We show that well-defined left and right facial modules are plesiomorphic for primates, while anthropoids consistently have asymmetrical facial modules that include structures of both sides, a change likely related to the ability to display more complex, asymmetrical facial expressions. However, no clear trends in network organization were found regarding the evolution of structures related to speech. Remarkably, the increase in the number of head and neck muscles – and thus of musculoskeletal structures – in human evolution led to a decrease in network density and complexity in humans

Protein Sequence Alignment Analysis by Local Covariation: Coevolution Statistics Detect Benchmark Alignment Errors

The use of sequence alignments to understand protein families is ubiquitous in molecular biology. High quality alignments are difficult to build and protein alignment remains one of the largest open problems in computational biology. Misalignments can lead to inferential errors about protein structure, folding, function, phylogeny, and residue importance. Identifying alignment errors is difficult because alignments are built and validated on the same primary criteria: sequence conservation. Local covariation identifies systematic misalignments and is independent of conservation. We demonstrate an alignment curation tool, LoCo, that integrates local covariation scores with the Jalview alignment editor. Using LoCo, we illustrate how local covariation is capable of identifying alignment errors due to the reduction of positional independence in the region of misalignment. We highlight three alignments from the benchmark database, BAliBASE 3, that contain regions of high local covariation, and investigate the causes to illustrate these types of scenarios. Two alignments contain sequential and structural shifts that cause elevated local covariation. Realignment of these misaligned segments reduces local covariation; these alternative alignments are supported with structural evidence. We also show that local covariation identifies active site residues in a validated alignment of paralogous structures. Loco is available at https://sourceforge.net/projects/locoprotein/files

Anatomical Network Comparison of Human Upper and Lower, Newborn and Adult, and Normal and Abnormal Limbs, with Notes on Development, Pathology and Limb Serial Homology vs. Homoplasy

How do the various anatomical parts (modules) of the animal body evolve into very different integrated forms (integration) yet still function properly without decreasing the individual's survival? This long-standing question remains unanswered for multiple reasons, including lack of consensus about conceptual definitions and approaches, as well as a reasonable bias toward the study of hard tissues over soft tissues. A major difficulty concerns the non-trivial technical hurdles of addressing this problem, specifically the lack of quantitative tools to quantify and compare variation across multiple disparate anatomical parts and tissue types. In this paper we apply for the first time a powerful new quantitative tool, Anatomical Network Analysis (AnNA), to examine and compare in detail the musculoskeletal modularity and integration of normal and abnormal human upper and lower limbs. In contrast to other morphological methods, the strength of AnNA is that it allows efficient and direct empirical comparisons among body parts with even vastly different architectures (e.g. upper and lower limbs) and diverse or complex tissue composition (e.g. bones, cartilages and muscles), by quantifying the spatial organization of these parts-their topological patterns relative to each other-using tools borrowed from network theory. Our results reveal similarities between the skeletal networks of the normal newborn/adult upper limb vs. lower limb, with exception to the shoulder vs. pelvis. However, when muscles are included, the overall musculoskeletal network organization of the upper limb is strikingly different from that of the lower limb, particularly that of the more proximal structures of each limb. Importantly, the obtained data provide further evidence to be added to the vast amount of paleontological, gross anatomical, developmental, molecular and embryological data recently obtained that contradicts the long-standing dogma that the upper and lower limbs are serial homologues. In addition, the AnNA of the limbs of a trisomy 18 human fetus strongly supports Pere Alberch's ill-named "logic of monsters" hypothesis, and contradicts the commonly accepted idea that birth defects often lead to lower integration (i.e. more parcellation) of anatomical structures

An Update on MyoD Evolution in Teleosts and a Proposed Consensus Nomenclature to Accommodate the Tetraploidization of Different Vertebrate Genomes

DJM was supported by a Natural Environment Research Council studentship (NERC/S/A/2004/12435).Background: MyoD is a muscle specific transcription factor that is essential for vertebrate myogenesis. In several teleost species, including representatives of the Salmonidae and Acanthopterygii, but not zebrafish, two or more MyoD paralogues are conserved that are thought to have arisen from distinct, possibly lineage-specific duplication events. Additionally, two MyoD paralogues have been characterised in the allotetraploid frog, Xenopus laevis. This has lead to a confusing nomenclature since MyoD paralogues have been named outside of an appropriate phylogenetic framework. Methods and Principal Findings: Here we initially show that directly depicting the evolutionary relationships of teleost MyoD orthologues and paralogues is hindered by the asymmetric evolutionary rate of Acanthopterygian MyoD2 relative to other MyoD proteins. Thus our aim was to confidently position the event from which teleost paralogues arose in different lineages by a comparative investigation of genes neighbouring myod across the vertebrates. To this end, we show that genes on the single myod-containing chromosome of mammals and birds are retained in both zebrafish and Acanthopterygian teleosts in a striking pattern of double conserved synteny. Further, phylogenetic reconstruction of these neighbouring genes using Bayesian and maximum likelihood methods supported a common origin for teleost paralogues following the split of the Actinopterygii and Sarcopterygii. Conclusion: Our results strongly suggest that myod was duplicated during the basal teleost whole genome duplication event, but was subsequently lost in the Ostariophysi ( zebrafish) and Protacanthopterygii lineages. We propose a sensible consensus nomenclature for vertebrate myod genes that accommodates polyploidization events in teleost and tetrapod lineages and is justified from a phylogenetic perspective.Publisher PDFPeer reviewe

Genetic properties of feed efficiency parameters in meat-type chickens

<p>Abstract</p> <p>Background</p> <p>Feed cost constitutes about 70% of the cost of raising broilers, but the efficiency of feed utilization has not kept up the growth potential of today's broilers. Improvement in feed efficiency would reduce the amount of feed required for growth, the production cost and the amount of nitrogenous waste. We studied residual feed intake (RFI) and feed conversion ratio (FCR) over two age periods to delineate their genetic inter-relationships.</p> <p>Methods</p> <p>We used an animal model combined with Gibb sampling to estimate genetic parameters in a pedigreed random mating broiler control population.</p> <p>Results</p> <p>Heritability of RFI and FCR was 0.42-0.45. Thus selection on RFI was expected to improve feed efficiency and subsequently reduce feed intake (FI). Whereas the genetic correlation between RFI and body weight gain (BWG) at days 28-35 was moderately positive, it was negligible at days 35-42. Therefore, the timing of selection for RFI will influence the expected response. Selection for improved RFI at days 28-35 will reduce FI, but also increase growth rate. However, selection for improved RFI at days 35-42 will reduce FI without any significant change in growth rate. The nature of the pleiotropic relationship between RFI and FCR may be dependent on age, and consequently the molecular factors that govern RFI and FCR may also depend on stage of development, or on the nature of resource allocation of FI above maintenance directed towards protein accretion and fat deposition. The insignificant genetic correlation between RFI and BWG at days 35-42 demonstrates the independence of RFI on the level of production, thereby making it possible to study the molecular, physiological and nutrient digestibility mechanisms underlying RFI without the confounding effects of growth. The heritability estimate of FCR was 0.49 and 0.41 for days 28-35 and days 35-42, respectively.</p> <p>Conclusion</p> <p>Selection for FCR will improve efficiency of feed utilization but because of the genetic dependence of FCR and its components, selection based on FCR will reduce FI and increase growth rate. However, the correlated responses in both FI and BWG cannot be predicted accurately because of the inherent problem of FCR being a ratio trait.</p

Identifying and Seeing beyond Multiple Sequence Alignment Errors Using Intra-Molecular Protein Covariation

BACKGROUND: There is currently no way to verify the quality of a multiple sequence alignment that is independent of the assumptions used to build it. Sequence alignments are typically evaluated by a number of established criteria: sequence conservation, the number of aligned residues, the frequency of gaps, and the probable correct gap placement. Covariation analysis is used to find putatively important residue pairs in a sequence alignment. Different alignments of the same protein family give different results demonstrating that covariation depends on the quality of the sequence alignment. We thus hypothesized that current criteria are insufficient to build alignments for use with covariation analyses. METHODOLOGY/PRINCIPAL FINDINGS: We show that current criteria are insufficient to build alignments for use with covariation analyses as systematic sequence alignment errors are present even in hand-curated structure-based alignment datasets like those from the Conserved Domain Database. We show that current non-parametric covariation statistics are sensitive to sequence misalignments and that this sensitivity can be used to identify systematic alignment errors. We demonstrate that removing alignment errors due to 1) improper structure alignment, 2) the presence of paralogous sequences, and 3) partial or otherwise erroneous sequences, improves contact prediction by covariation analysis. Finally we describe two non-parametric covariation statistics that are less sensitive to sequence alignment errors than those described previously in the literature. CONCLUSIONS/SIGNIFICANCE: Protein alignments with errors lead to false positive and false negative conclusions (incorrect assignment of covariation and conservation, respectively). Covariation analysis can provide a verification step, independent of traditional criteria, to identify systematic misalignments in protein alignments. Two non-parametric statistics are shown to be somewhat insensitive to misalignment errors, providing increased confidence in contact prediction when analyzing alignments with erroneous regions because of an emphasis on they emphasize pairwise covariation over group covariation

On the Accessibility of Adaptive Phenotypes of a Bacterial Metabolic Network

The mechanisms by which adaptive phenotypes spread within an evolving population after their emergence are understood fairly well. Much less is known about the factors that influence the evolutionary accessibility of such phenotypes, a pre-requisite for their emergence in a population. Here, we investigate the influence of environmental quality on the accessibility of adaptive phenotypes of Escherichia coli's central metabolic network. We used an established flux-balance model of metabolism as the basis for a genotype-phenotype map (GPM). We quantified the effects of seven qualitatively different environments (corresponding to both carbohydrate and gluconeogenic metabolic substrates) on the structure of this GPM. We found that the GPM has a more rugged structure in qualitatively poorer environments, suggesting that adaptive phenotypes could be intrinsically less accessible in such environments. Nevertheless, on average ∼74% of the genotype can be altered by neutral drift, in the environment where the GPM is most rugged; this could allow evolving populations to circumvent such ruggedness. Furthermore, we found that the normalized mutual information (NMI) of genotype differences relative to phenotype differences, which measures the GPM's capacity to transmit information about phenotype differences, is positively correlated with (simulation-based) estimates of the accessibility of adaptive phenotypes in different environments. These results are consistent with the predictions of a simple analytic theory that makes explicit the relationship between the NMI and the speed of adaptation. The results suggest an intuitive information-theoretic principle for evolutionary adaptation; adaptation could be faster in environments where the GPM has a greater capacity to transmit information about phenotype differences. More generally, our results provide insight into fundamental environment-specific differences in the accessibility of adaptive phenotypes, and they suggest opportunities for research at the interface between information theory and evolutionary biology

Embryonic stem cells in scaffold-free three-dimensional cell culture: osteogenic differentiation and bone generation

Extracorporeal formation of mineralized bone-like tissue is still an unsolved challenge in tissue engineering. Embryonic stem cells may open up new therapeutic options for the future and should be an interesting model for the analysis of fetal organogenesis. Here we describe a technique for culturing embryonic stem cells (ESCs) in the absence of artificial scaffolds which generated mineralized miromasses. Embryonic stem cells were harvested and osteogenic differentiation was stimulated by the addition of dexamethasone, ascorbic acid, and ß-glycerolphosphate (DAG). After three days of cultivation microspheres were formed. These spherical three-dimensional cell units showed a peripheral zone consisting of densely packed cell layers surrounded by minerals that were embedded in the extracellular matrix. Alizarine red staining confirmed evidence of mineralization after 10 days of DAG stimulation in the stimulated but not in the control group. Transmission electron microscopy demonstrated scorching crystallites and collagenous fibrils as early indication of bone formation. These extracellular structures resembled hydroxyl apatite-like crystals as demonstrated by distinct diffraction patterns using electron diffraction analysis. The micromass culture technique is an appropriate model to form three-dimensional bone-like micro-units without the need for an underlying scaffold. Further studies will have to show whether the technique is applicable also to pluripotent stem cells of different origin