Safe and complete contig assembly via omnitigs

Contig assembly is the first stage that most assemblers solve when reconstructing a genome from a set of reads. Its output consists of contigs -- a set of strings that are promised to appear in any genome that could have generated the reads. From the introduction of contigs 20 years ago, assemblers have tried to obtain longer and longer contigs, but the following question was never solved: given a genome graph $G$ (e.g. a de Bruijn, or a string graph), what are all the strings that can be safely reported from $G$ as contigs? In this paper we finally answer this question, and also give a polynomial time algorithm to find them. Our experiments show that these strings, which we call omnitigs, are 66% to 82% longer on average than the popular unitigs, and 29% of dbSNP locations have more neighbors in omnitigs than in unitigs.Comment: Full version of the paper in the proceedings of RECOMB 201

Body odor quality predicts behavioral attractiveness in humans

Growing effort is being made to understand how different attractive physical traits co-vary within individuals, partly because this might indicate an underlying index of genetic quality. In humans, attention has focused on potential markers of quality such as facial attractiveness, axillary odor quality, the second-to-fourth digit (2D:4D) ratio and body mass index (BMI). Here we extend this approach to include visually-assessed kinesic cues (nonverbal behavior linked to movement) which are statistically independent of structural physical traits. The utility of such kinesic cues in mate assessment is controversial, particularly during everyday conversational contexts, as they could be unreliable and susceptible to deception. However, we show here that the attractiveness of nonverbal behavior, in 20 male participants, is predicted by perceived quality of their axillary body odor. This finding indicates covariation between two desirable traits in different sensory modalities. Depending on two different rating contexts (either a simple attractiveness rating or a rating for long-term partners by 10 female raters not using hormonal contraception), we also found significant relationships between perceived attractiveness of nonverbal behavior and BMI, and between axillary odor ratings and 2D:4D ratio. Axillary odor pleasantness was the single attribute that consistently predicted attractiveness of nonverbal behavior. Our results demonstrate that nonverbal kinesic cues could reliably reveal mate quality, at least in males, and could corroborate and contribute to mate assessment based on other physical traits

Complex trait subtypes identification using transcriptome profiling reveals an interaction between two QTL affecting adiposity in chicken

<p>Abstract</p> <p>Background</p> <p>Integrative genomics approaches that combine genotyping and transcriptome profiling in segregating populations have been developed to dissect complex traits. The most common approach is to identify genes whose eQTL colocalize with QTL of interest, providing new functional hypothesis about the causative mutation. Another approach includes defining subtypes for a complex trait using transcriptome profiles and then performing QTL mapping using some of these subtypes. This approach can refine some QTL and reveal new ones.</p> <p>In this paper we introduce Factor Analysis for Multiple Testing (FAMT) to define subtypes more accurately and reveal interaction between QTL affecting the same trait. The data used concern hepatic transcriptome profiles for 45 half sib male chicken of a sire known to be heterozygous for a QTL affecting abdominal fatness (AF) on chromosome 5 distal region around 168 cM.</p> <p>Results</p> <p>Using this methodology which accounts for hidden dependence structure among phenotypes, we identified 688 genes that are significantly correlated to the AF trait and we distinguished 5 subtypes for AF trait, which are not observed with gene lists obtained by classical approaches. After exclusion of one of the two lean bird subtypes, linkage analysis revealed a previously undetected QTL on chromosome 5 around 100 cM. Interestingly, the animals of this subtype presented the same q paternal haplotype at the 168 cM QTL. This result strongly suggests that the two QTL are in interaction. In other words, the "q configuration" at the 168 cM QTL could hide the QTL existence in the proximal region at 100 cM. We further show that the proximal QTL interacts with the previous one detected on the chromosome 5 distal region.</p> <p>Conclusion</p> <p>Our results demonstrate that stratifying genetic population by molecular phenotypes followed by QTL analysis on various subtypes can lead to identification of novel and interacting QTL.</p

Human-animal chimeras for vaccine development: an endangered species or opportunity for the developing world?

<p>Abstract</p> <p>Background</p> <p>In recent years, the field of vaccines for diseases such as Human Immunodeficiency Virus (HIV) which take a heavy toll in developing countries has faced major failures. This has led to a call for more basic science research, and development as well as evaluation of new vaccine candidates. Human-animal chimeras, developed with a 'humanized' immune system could be useful to study infectious diseases, including many neglected diseases. These would also serve as an important tool for the efficient testing of new vaccine candidates to streamline promising candidates for further trials in humans. However, developing human-animal chimeras has proved to be controversial.</p> <p>Discussion</p> <p>Development of human-animal chimeras for vaccine development has been slowed down because of opposition by some philosophers, ethicists and policy makers in the west-they question the moral status of such animals, and also express discomfort about transgression of species barriers. Such opposition often uses a contemporary western world view as a reference point. Human-animal chimeras are often being created for diseases which cause significantly higher morbidity and mortality in the developing world as compared to the developed world. We argue in our commentary that given this high disease burden, we should look at socio-cultural perspectives on human-animal chimera like beings in the developing world. On examination, it's clear that such beings have been part of mythology and cultural descriptions in many countries in the developing world.</p> <p>Summary</p> <p>To ensure that important research on diseases afflicting millions like malaria, HIV, Hepatitis-C and dengue continues to progress, we recommend supporting human-animal chimera research for vaccine development in developing countries (especially China and India which have growing technical expertise in the area). The negative perceptions in some parts of the west about human-animal chimeras can be used as an opportunity for nurturing important vaccine development research in the developing world.</p

Genetic variation in the pleiotropic association between physical activity and body weight in mice

<p>Abstract</p> <p>Background</p> <p>A sedentary lifestyle is often assumed to lead to increases in body weight and potentially obesity and related diseases but in fact little is known about the genetic association between physical activity and body weight. We tested for such an association between body weight and the distance, duration, and speed voluntarily run by 310 mice from the F₂generation produced from an intercross of two inbred lines that differed dramatically in their physical activity levels.</p> <p>Methods</p> <p>We used a conventional interval mapping approach with SNP markers to search for QTLs that affected both body weight and activity traits. We also conducted a genome scan to search for relationship QTLs (<it>rel</it>QTLs), or chromosomal regions that affected an activity trait variably depending on the phenotypic value of body weight.</p> <p>Results</p> <p>We uncovered seven quantitative trait loci (QTLs) affecting body weight, but only one co-localized with another QTL previously found for activity traits. We discovered 19 <it>rel</it>QTLs that provided evidence for a genetic (pleiotropic) association of physical activity and body weight. The three genotypes at each of these loci typically exhibited a combination of negative, zero, and positive regressions of the activity traits on body weight, the net effect of which was to produce overall independence of body weight from physical activity. We also demonstrated that the <it>rel</it>QTLs produced these varying associations through differential epistatic interactions with a number of other epistatic QTLs throughout the genome.</p> <p>Conclusion</p> <p>It was concluded that individuals with specific combinations of genotypes at the <it>rel</it>QTLs and <it>epi</it>QTLs might account for some of the variation typically seen in plots of the association of physical activity with body weight.</p

Homozygosity Mapping on Homozygosity Haplotype Analysis to Detect Recessive Disease-Causing Genes from a Small Number of Unrelated, Outbred Patients

Genes involved in disease that are not common are often difficult to identify; a method that pinpoints them from a small number of unrelated patients will be of great help. In order to establish such a method that detects recessive genes identical-by-descent, we modified homozygosity mapping (HM) so that it is constructed on the basis of homozygosity haplotype (HM on HH) analysis. An analysis using 6 unrelated patients with Siiyama-type α1-antitrypsin deficiency, a disease caused by a founder gene, the correct gene locus was pinpointed from data of any 2 patients (length: 1.2–21.8 centimorgans, median: 1.6 centimorgans). For a test population in which these 6 patients and 54 healthy subjects were scrambled, the approach accurately identified these 6 patients and pinpointed the locus to a 1.4-centimorgan fragment. Analyses using synthetic data revealed that the analysis works well for IBD fragment derived from a most recent common ancestor (MRCA) who existed less than 60 generations ago. The analysis is unsuitable for the genes with a frequency in general population more than 0.1. Thus, HM on HH analysis is a powerful technique, applicable to a small number of patients not known to be related, and will accelerate the identification of disease-causing genes for recessive conditions

Quantitative trait locus analysis of hybrid pedigrees: variance-components model, inbreeding parameter, and power

<p>Abstract</p> <p>Background</p> <p>For the last years reliable mapping of quantitative trait loci (QTLs) has become feasible through linkage analysis based on the variance-components method. There are now many approaches to the QTL analysis of various types of crosses within one population (breed) as well as crosses between divergent populations (breeds). However, to analyse a complex pedigree with dominance and inbreeding, when the pedigree's founders have an inter-population (hybrid) origin, it is necessary to develop a high-powered method taking into account these features of the pedigree.</p> <p>Results</p> <p>We offer a universal approach to QTL analysis of complex pedigrees descended from crosses between outbred parental lines with different QTL allele frequencies. This approach improves the established variance-components method due to the consideration of the genetic effect conditioned by inter-population origin and inbreeding of individuals. To estimate model parameters, namely additive and dominant effects, and the allelic frequencies of the QTL analysed, and also to define the QTL positions on a chromosome with respect to genotyped markers, we used the maximum-likelihood method. To detect linkage between the QTL and the markers we propose statistics with a non-central χ²-distribution that provides the possibility to deduce analytical expressions for the power of the method and therefore, to estimate the pedigree's size required for 80% power. The method works for arbitrarily structured pedigrees with dominance and inbreeding.</p> <p>Conclusion</p> <p>Our method uses the phenotypic values and the marker information for each individual of the pedigree under observation as initial data and can be valuable for fine mapping purposes. The power of the method is increased if the QTL effects conditioned by inter-population origin and inbreeding are enhanced. Several improvements can be developed to take into account fixed factors affecting trait formation, such as age and sex.</p

A Bayesian Framework to Account for Complex Non-Genetic Factors in Gene Expression Levels Greatly Increases Power in eQTL Studies

Gene expression measurements are influenced by a wide range of factors, such as the state of the cell, experimental conditions and variants in the sequence of regulatory regions. To understand the effect of a variable of interest, such as the genotype of a locus, it is important to account for variation that is due to confounding causes. Here, we present VBQTL, a probabilistic approach for mapping expression quantitative trait loci (eQTLs) that jointly models contributions from genotype as well as known and hidden confounding factors. VBQTL is implemented within an efficient and flexible inference framework, making it fast and tractable on large-scale problems. We compare the performance of VBQTL with alternative methods for dealing with confounding variability on eQTL mapping datasets from simulations, yeast, mouse, and human. Employing Bayesian complexity control and joint modelling is shown to result in more precise estimates of the contribution of different confounding factors resulting in additional associations to measured transcript levels compared to alternative approaches. We present a threefold larger collection of cis eQTLs than previously found in a whole-genome eQTL scan of an outbred human population. Altogether, 27% of the tested probes show a significant genetic association in cis, and we validate that the additional eQTLs are likely to be real by replicating them in different sets of individuals. Our method is the next step in the analysis of high-dimensional phenotype data, and its application has revealed insights into genetic regulation of gene expression by demonstrating more abundant cis-acting eQTLs in human than previously shown. Our software is freely available online at http://www.sanger.ac.uk/resources/software/peer/