Allelic penetrance approach as a tool to model two-locus interaction in complex binary traits

Many binary phenotypes do not follow a classical Mendelian inheritance pattern. Interaction between genetic and environmental factors is thought to contribute to the incomplete penetrance phenomena often observed in these complex binary traits. Several two-locus models for penetrance have been proposed to aid the genetic dissection of binary traits. Such models assume linear genetic effects of both loci in different mathematical scales of penetrance, resembling the analytical framework of quantitative traits. However, changes in phenotypic scale are difficult to envisage in binary traits and limited genetic interpretation is extractable from current modeling of penetrance. To overcome this limitation, we derived an allelic penetrance approach that attributes incomplete penetrance to the stochastic expression of the alleles controlling the phenotype, the genetic background and environmental factors. We applied this approach to formulate dominance and recessiveness in a single diallelic locus and to model different genetic mechanisms for the joint action of two diallelic loci. We fit the models to data on the genetic susceptibility of mice following infections with Listeria monocytogenes and Plasmodium berghei. These models gain in genetic interpretation, because they specify the alleles that are responsible for the genetic (inter)action and their genetic nature (dominant or recessive), and predict genotypic combinations determining the phenotype. Further, we show via computer simulations that the proposed models produce penetrance patterns not captured by traditional two-locus models. This approach provides a new analysis framework for dissecting mechanisms of interlocus joint action in binary traits using genetic crosses.     

From complex traits to complex alleles

Several recent studies using natural populations of Drosophila show that one must be very careful when sorting among the large number of molecular polymorphisms found at most loci to identify the nucleotide changes responsible for phenotypic variation in complex traits. Indeed, several mutations within a single allele can interact to generate the overall observed effect. The results are instructive both for those interested in the genetics of evolutionary change and for those attempting to ferret out the genetic basis of complex human diseases.     

ForSim: a tool for exploring the genetic architecture of complex traits with controlled truth

Many important problems in biology involve complex traits affected by multiple coding or regulatory parts of the genome. How well the underlying genetic architecture can be inferred by statistical methods such as mapping and association studies are active research areas. ForSim is a flexible forward evolutionary simulation tool for exploring the consequences of evolution by phenotype, whereby demographic, chance, behavioral and selective effects mold genetic architecture. Simulation is useful for exploring these issues as well as the choice of study design inferential methods. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

A tooth,a toe,and a vertebra: The genetic dimensions of complex morphological traits

Most readers probably share with me the profoundly affecting experience of wandering through a natural history museum and being surrounded by the skeletal remains of our vertebrate predecessors. Who does not stand in awe before these scaffolds of the great dinosaurs or ponder the groping stages through which our primate ancestors explored the skeletal and dental opportunities of forest life (Fig. 1)? How did this proliferation of limbs, teeth, and vertebrae, not to mention the complexity we can infer about the overlying anatomy and physiology, arise from the simple early forms of animal life? How did DNA evolve to contain the program for such complexity?     

The genetic dissection of complex traits in a founder population

We estimated broad heritabilities (H(2)) and narrow heritabilities (h(2)) and conducted genomewide screens, using a novel association-based mapping approach for 20 quantitative trait loci (QTLs) among the Hutterites, a founder population that practices a communal lifestyle. Heritability estimates ranged from.21 for diastolic blood pressure (DBP) to.99 for whole-blood serotonin levels. Using a multipoint method to detect association under a recessive model we found evidence of major QTLs for six traits: low-density lipoprotein (LDL), triglycerides, lipoprotein (a) (Lp[a]), systolic blood pressure (SBP), serum cortisol, and whole-blood serotonin. Second major QTLs for Lp(a) and for cortisol were identified using a single-point method to detect association under a general two-allele model. The heritabilities for these six traits ranged from.37 for triglycerides to.99 for serotonin, and three traits (LDL, SBP, and serotonin) had significant dominance variances (i.e., H(2) > h(2)). Surprisingly, there was little correlation between measures of heritability and the strength of association on a genomewide screen (P>.50), suggesting that heritability estimates per se do not identify phenotypes that are influenced by genes with major effects. The present study demonstrates the feasibility of genomewide association studies for QTL mapping. However, even in this young founder population that has extensive linkage disequilibrium, map densities <5 cM may be required to detect all major QTLs.     

Epistasis,complex traits,and mapping genes

Using a three-locus model wherein two loci regulate a third, candidate locus, I examine physiological epistasis from the gene's eye view of the regulated locus. I show that, depending upon genetic background at the regulatory loci, an allele at the candidate locus can be dominant, additive, recessive, neutral, over-dominant, or under-dominant in its effects on fitness. This kind of variation in allelic effect caused by variation in genetic background from population to population, from time to time in the same population, or sample to sample makes finding and mapping the genes underlying a complex phenotype difficult. The rate of evolution of such genes can also be slowed, especially in genetically subdivided metapopulations with migration. Nevertheless, understanding how variation in genetic background causes variation in allelic effects permits the genetic architecture of such complex traits to be dissected into the interacting component genes. While some backgrounds diminish allelic effects and make finding and mapping genes difficult, other backgrounds enhance allelic effects and facilitate gene mapping.     

Mining the HapMap to dissect complex traits

A report on the Keystone Symposium 'Genome Sequence Variation and the Inherited Basis of Common Disease and Complex Traits', Big Sky, USA, 8-13 January 2006.     

Approaches to mapping genetically correlated complex traits

Our Markov chain Monte Carlo (MCMC) methods were used in linkage analyses of the Framingham Heart Study data using all available pedigrees. Our goal was to detect and map loci associated with covariate-adjusted traits log triglyceride (lnTG) and high-density lipoprotein cholesterol (HDL) using multipoint LOD score analysis, Bayesian oligogenic linkage analysis and identity-by-descent (IBD) scoring methods. Each method used all marker data for all markers on a chromosome. Bayesian linkage analysis detected a linkage signal on chromosome 7 for lnTG and HDL, corroborating previously published results. However, these results were not replicated in a classical linkage analysis of the data or by using IBD scoring methods.We conclude that Bayesian linkage analysis provides a powerful paradigm for mapping trait loci but interpretation of the Bayesian linkage signals is subjective. In the absence of a LOD score method accommodating genetically complex traits and linkage heterogeneity, validation of these signals remains elusive.     

The genome architecture of the Collaborative Cross mouse genetic reference population

The Collaborative Cross Consortium reports here on the development of a unique genetic resource population. The Collaborative Cross (CC) is a multiparental recombinant inbred panel derived from eight laboratory mouse inbred strains. Breeding of the CC lines was initiated at multiple international sites using mice from The Jackson Laboratory. Currently, this innovative project is breeding independent CC lines at the University of North Carolina (UNC), at Tel Aviv University (TAU), and at Geniad in Western Australia (GND). These institutions aim to make publicly available the completed CC lines and their genotypes and sequence information. We genotyped, and report here, results from 458 extant lines from UNC, TAU, and GND using a custom genotyping array with 7500 SNPs designed to be maximally informative in the CC and used a novel algorithm to infer inherited haplotypes directly from hybridization intensity patterns. We identified lines with breeding errors and cousin lines generated by splitting incipient lines into two or more cousin lines at early generations of inbreeding. We then characterized the genome architecture of 350 genetically independent CC lines. Results showed that founder haplotypes are inherited at the expected frequency, although we also consistently observed highly significant transmission ratio distortion at specific loci across all three populations. On chromosome 2, there is significant overrepresentation of WSB/EiJ alleles, and on chromosome X, there is a large deficit of CC lines with CAST/EiJ alleles. Linkage disequilibrium decays as expected and we saw no evidence of gametic disequilibrium in the CC population as a whole or in random subsets of the population. Gametic equilibrium in the CC population is in marked contrast to the gametic disequilibrium present in a large panel of classical inbred strains. Finally, we discuss access to the CC population and to the associated raw data describing the genetic structure of individual lines. Integration of rich phenotypic and genomic data over time and across a wide variety of fields will be vital to delivering on one of the key attributes of the CC, a common genetic reference platform for identifying causative variants and genetic networks determining traits in mammals.     

The Collaborative Cross,developing a resource for mammalian systems genetics: A status report of the Wellcome Trust cohort

We report on the progress of a project funded by the Wellcome Trust to produce over 100 recombinant inbred mouse lines as part of the Collaborative Cross (CC) genetic reference panel. These new strains of mice are being derived from a set of eight genetically diverse founders. The genomes of the finished strains will be mosaics of the founder strains’ genomes with a high density of independent recombination breakpoints. The CC mice will be available for distribution free of any intellectual property constraints to serve as a community resource for systems genetics studies.     

Epigrass: a tool to study disease spread in complex networks

Background  

The construction of complex spatial simulation models such as those used in network epidemiology, is a daunting task due to the large amount of data involved in their parameterization. Such data, which frequently resides on large geo-referenced databases, has to be processed and assigned to the various components of the model. All this just to construct the model, then it still has to be simulated and analyzed under different epidemiological scenarios. This workflow can only be achieved efficiently by computational tools that can automate most, if not all, these time-consuming tasks. In this paper, we present a simulation software, Epigrass, aimed to help designing and simulating network-epidemic models with any kind of node behavior.     

Engineering and control of biological systems: A new way to tackle complex diseases

The ongoing merge between engineering and biology has contributed to the emerging field of synthetic biology. The defining features of this new discipline are abstraction and standardisation of biological parts, decoupling between parts to prevent undesired cross-talking, and the application of quantitative modelling of synthetic genetic circuits in order to guide their design. Most of the efforts in the field of synthetic biology in the last decade have been devoted to the design and development of functional gene circuits in prokaryotes and unicellular eukaryotes. Researchers have used synthetic biology not only to engineer new functions in the cell, but also to build simpler models of endogenous gene regulatory networks to gain knowledge of the rules governing their wiring diagram. However, the need for innovative approaches to study and modify complex signalling and regulatory networks in mammalian cells and multicellular organisms has prompted advances of synthetic biology also in these species, thus contributing to develop innovative ways to tackle human diseases. In this work, we will review the latest progress in synthetic biology and the most significant developments achieved so far, both in unicellular and multicellular organisms, with emphasis on human health.     

The evolution of quantitative traits in complex environments

Species inhabit complex environments and respond to selection imposed by numerous abiotic and biotic conditions that vary in both space and time. Environmental heterogeneity strongly influences trait evolution and patterns of adaptive population differentiation. For example, heterogeneity can favor local adaptation, or can promote the evolution of plastic genotypes that alter their phenotypes based on the conditions they encounter. Different abiotic and biotic agents of selection can act synergistically to either accelerate or constrain trait evolution. The environmental context has profound effects on quantitative genetic parameters. For instance, heritabilities measured in controlled conditions often exceed those measured in the field; thus, laboratory experiments could overestimate the potential for a population to respond to selection. Nevertheless, most studies of the genetic basis of ecologically relevant traits are conducted in simplified laboratory environments, which do not reflect the complexity of nature. Here, we advocate for manipulative field experiments in the native ranges of plant species that differ in mating system, life-history strategy and growth form. Field studies are vital to evaluate the roles of disparate agents of selection, to elucidate the targets of selection and to develop a nuanced perspective on the evolution of quantitative traits. Quantitative genetics field studies will also shed light on the potential for natural populations to adapt to novel climates in highly fragmented landscapes. Drawing from our experience with the ecological model system Boechera (Brassicaceae), we discuss advancements possible through dedicated field studies, highlight future research directions and examine the challenges associated with field studies.     

The Collaborative Cross as a Resource for Modeling Human Disease: CC011/Unc,a New Mouse Model for Spontaneous Colitis

Inflammatory bowel disease (IBD) is an immune-mediated condition driven by improper responses to intestinal microflora in the context of environmental and genetic background. GWAS in humans have identified many loci associated with IBD, but animal models are valuable for dissecting the underlying molecular mechanisms, characterizing environmental and genetic contributions and developing treatments. Mouse models rely on interventions such as chemical treatment or introduction of an infectious agent to induce disease. Here, we describe a new model for IBD in which the disease develops spontaneously in 20-week-old mice in the absence of known murine pathogens. The model is part of the Collaborative Cross and came to our attention due to a high incidence of rectal prolapse in an incompletely inbred line. Necropsies revealed a profound proliferative colitis with variable degrees of ulceration and vasculitis, splenomegaly and enlarged mesenteric lymph nodes with no discernible anomalies of other organ systems. Phenotypic characterization of the CC011/Unc mice with homozygosity ranging from 94.1 to 99.8 % suggested that the trait was fixed and acted recessively in crosses to the colitis-resistant C57BL/6J inbred strain. Using a QTL approach, we identified four loci, Ccc1, Ccc2, Ccc3 and Ccc4 on chromosomes 12, 14, 1 and 8 that collectively explain 27.7 % of the phenotypic variation. Surprisingly, we also found that minute levels of residual heterozygosity in CC011/Unc have significant impact on the phenotype. This work demonstrates the utility of the CC as a source of models of human disease that arises through new combinations of alleles at susceptibility loci.     

Response of a complex foraging phenotype to artificial selection on its component traits

In nature, where predators must often track dynamic and dispersed prey populations, predator consumption rate, conversion efficiency, dispersal, and prey finding are likely to be important links between foraging and predator–prey population dynamics. Small differences in predator foraging caused by variation in any of the abovementioned traits might lead to significant differences in predator success as well as population dynamics. We used artificial selection to create lines of the predatory mite, Phytoseiulus persimilis in order to determine the potential for or constraints on the evolution of predator foraging behaviors. All four foraging traits demonstrated considerable phenotypic variation. They also exhibited significant realized heritabilities after artificial selection, except that prey finding did not respond to downward selection. Lines that responded to selection did so rapidly, and high-consumption, high-conversion efficiency, and high- and low-dispersal were stable for at least four generations after artificial selection was relaxed. There were some indirect responses to selection among the foraging traits. For example, there was positive correlation between consumption and dispersal. However, none of the correlated responses were of the magnitude of the direct responses we measured on the same trait. We also observed some correlations between foraging traits and life-history traits such as low-consumption and development time (negative), high-consumption and fecundity (positive), and high-conversion efficiency and fecundity (positive), but these were more likely to represent non-genetic constraints. Intrinsic rates of increase in low-consumption and low-conversion efficiency lines were lower than in their respective high lines and the unselected control, whereas rates of increase in dispersal and olfactory response lines did not differ from the unselected control. Thus, traits that make up foraging share partially overlapping genetic architectures with highly heritable phenotypic components, suggesting that each foraging trait will be able to respond rapidly to changes in the density and distribution of resources.     

DNA barcoding for community ecology--how to tackle a hyperdiverse, mostly undescribed Melanesian fauna

Background

Trigonopterus weevils are widely distributed throughout Melanesia and hyperdiverse in New Guinea. They are a dominant feature in natural forests, with narrow altitudinal zonation. Their use in community ecology has been precluded by the “taxonomic impediment”.

Methodology/Principal Findings

We sampled >6,500 specimens from seven areas across New Guinea; 1,002 specimens assigned to 270 morphospecies were DNA sequenced. Objective clustering of a refined dataset (excluding nine cryptic species) at 3% threshold revealed 324 genetic clusters (DNA group count relative to number of morphospecies = 20.0% overestimation of species diversity, or 120.0% agreement) and 85.6% taxonomic accuracy (the proportion of DNA groups that “perfectly” agree with morphology-based species hypotheses). Agreement and accuracy were best at an 8% threshold. GMYC analysis revealed 328 entities (21.5% overestimation) with 227 perfect GMYC entities (84.1% taxonomic accuracy). Both methods outperform the parataxonomist (19% underestimation; 31.6% taxonomic accuracy). The number of species found in more than one sampling area was highest in the Eastern Highlands and Huon (Sørensen similarity index 0.07, 4 shared species); ⅓ of all areas had no species overlap. Success rates of DNA barcoding methods were lowest when species showed a pronounced geographical structure. In general, Trigonopterus show high α and β-diversity across New Guinea.

Conclusions/Significance

DNA barcoding is an excellent tool for biodiversity surveys but success rates might drop when closer localities are included. Hyperdiverse Trigonopterus are a useful taxon for evaluating forest remnants in Melanesia, allowing finer-grained analyses than would be possible with vertebrate taxa commonly used to date. Our protocol should help establish other groups of hyperdiverse fauna as target taxa for community ecology. Sequencing delivers objective data on taxa of incredible diversity but mostly without a solid taxonomic foundation and should help pave the road for the eventual formal naming of new species.     

Bay-0 x Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis

Natural genetic variation in Arabidopsis is considerable, but has not yet been used extensively as a source of variants to identify new genes of interest. From the cross between two genetically distant ecotypes, Bay-0 and Shahdara, we generated a Recombinant Inbred Line (RIL) population dedicated to Quantitative Trait Locus (QTL) mapping. A set of 38 physically anchored microsatellite markers was created to construct a robust genetic map from the 420 F6 lines. These markers, evenly distributed throughout the five chromosomes, revealed a remarkable equilibrium in the segregation of parental alleles in the genome. As a model character, we have analysed the genetic basis of variation in flowering time in two different environments. The simultaneous mapping of both large- and small-effect QTLs responsible for this variation explained 90% of the total genotypic variance. Two of the detected QTLs colocalize very precisely with FRIGIDA and FLOWERING LOCUS C genes; we provide information on the polymorphism of genes confirming this hypothesis. Another QTL maps in a region where no QTL had been found previously for this trait. This confirms the accuracy of QTL detection using the Bay-0 x Shahdara RIL population, which constitutes the largest in size available so far in Arabidopsis. As an alternative to mutant analysis, this population represents a powerful tool which is currently being used to undertake the genetic dissection of complex metabolic pathways.