A Hidden Markov Model Combining Linkage and Linkage Disequilibrium Information for Haplotype Reconstruction and Quantitative Trait Locus Fine Mapping

Faithful reconstruction of haplotypes from diploid marker data (phasing) is important for many kinds of genetic analyses, including mapping of trait loci, prediction of genomic breeding values, and identification of signatures of selection. In human genetics, phasing most often exploits population information (linkage disequilibrium), while in animal genetics the primary source of information is familial (Mendelian segregation and linkage). We herein develop and evaluate a method that simultaneously exploits both sources of information. It builds on hidden Markov models that were initially developed to exploit population information only. We demonstrate that the approach improves the accuracy of allele phasing as well as imputation of missing genotypes. Reconstructed haplotypes are assigned to hidden states that are shown to correspond to clusters of genealogically related chromosomes. We show that these cluster states can directly be used to fine map QTL. The method is computationally effective at handling large data sets based on high-density SNP panels.ARRAY technology now allows genotyping of large cohorts for thousands to millions of single nucleotide polymorphisms (SNPs), which are becoming available for a growing list of organisms including human and domestic animals. Among other applications, these advances permit systematic scanning of the genome to map trait loci by association (e.g., Wellcome Trust Case Control Consortium 2007; Charlier et al. 2008), to predict genomic breeding values for complex traits (Meuwissen et al. 2001; Goddard and Hayes 2009), or to identify signatures of selection (e.g., Voight et al. 2006).Present-day genotyping platforms do not directly provide information about linkage phase; i.e., co-inherited alleles at adjacent heterozygous markers (haplotypes) are not identified as such. As haplotype information may considerably empower genetic analyses, indirect phasing strategies have been devised: haplotypes can be reconstructed from unphased genotypes using either familial information (Mendelian segregation and linkage) and/or population information (linkage disequilibrium, LD, and surrogate parents) (e.g., Windig and Meuwissen 2004; Scheet and Stephens 2006; Kong et al. 2008).Haplotype-based approaches are routinely applied in animal genetics for combined linkage and LD mapping of QTL (e.g., Meuwissen and Goddard 2000; Blott et al. 2003). In these studies, phasing has so far relied on familial information provided by the extended pedigrees typical of livestock (e.g., Windig and Meuwissen 2004). This approach, however, leaves a nonnegligible proportion of genotypes unphased, especially for the less connected individuals. After phasing, identity-by-descent (IBD) probabilities conditional on haplotype data—needed for QTL mapping—are computed for all chromosome pairs, using familial as well as population information (hence combined linkage and LD mapping – L + LD) (e.g., Meuwissen and Goddard 2001). However, the use of high-density SNP chips and the analysis of ever larger cohorts render the computation of pairwise IBD probabilities a bottleneck.We herein propose a more efficient, heuristic approach based on hidden Markov models (HMM). It simultaneously phases and sorts haplotypes in clusters that can be used directly for mapping or other purposes. The proposed method exploits familial as well as population information, and imputes missing genotypes. We herein describe the accuracy of the proposed method and its use for L + LD mapping of QTL.     

Bayesian Inference of Genetic Parameters Based on Conditional Decompositions of Multivariate Normal Distributions

It is widely recognized that the mixed linear model is an important tool for parameter estimation in the analysis of complex pedigrees, which includes both pedigree and genomic information, and where mutually dependent genetic factors are often assumed to follow multivariate normal distributions of high dimension. We have developed a Bayesian statistical method based on the decomposition of the multivariate normal prior distribution into products of conditional univariate distributions. This procedure permits computationally demanding genetic evaluations of complex pedigrees, within the user-friendly computer package WinBUGS. To demonstrate and evaluate the flexibility of the method, we analyzed two example pedigrees: a large noninbred pedigree of Scots pine (Pinus sylvestris L.) that includes additive and dominance polygenic relationships and a simulated pedigree where genomic relationships have been calculated on the basis of a dense marker map. The analysis showed that our method was fast and provided accurate estimates and that it should therefore be a helpful tool for estimating genetic parameters of complex pedigrees quickly and reliably.MUCH effort in genetics has been devoted to revealing the underlying genetic architecture of quantitative or complex traits. Traditionally, the polygenic model has been used extensively to estimate genetic variances and breeding values of natural and breeding populations, where an infinite number of genes is assumed to code for the trait of interest (Bulmer 1971; Falconer and Mackay 1996). The genetic variance of a quantitative trait can be decomposed into an additive part that corresponds to the effects of individual alleles and a part that is nonadditive because of interactions between alleles. Attention has generally been focused on the estimation of additive genetic variance (and heritability), since additive variation is directly proportional to the response of selection via the breeder''s equation (Falconer and Mackay 1996, Chap. 11). However, to estimate additive genetic variation and heritability accurately, it can be important to identify potential nonadditive sources in genetic evaluations (Misztal 1997; Ovaskainen et al. 2008; Waldmann et al. 2008), especially if the pedigree being analyzed contains a large proportion of full-sibs and clones, as these in particular give rise to nonadditive genetic relationships (Lynch and Walsh 1998, pp. 145). The polygenic model using pedigree and phenotypic information, i.e., the animal model (Henderson 1984), has been the model of choice for estimating genetic parameters in breeding and natural populations (Abney et al. 2000; Sorensen and Gianola 2002; O′Hara et al. 2008).Recent breakthroughs in molecular techniques have made it possible to create genome-wide, single nucleotide polymorphism (SNP) maps. These maps have helped to uncover a vast amount of new loci responsible for trait expression and have provided general insights into the genetic architecture of quantitative traits (e.g., Valdar et al. 2006; Visscher 2008; Flint and Mackay 2009). These insights can help when calculating disease risks in humans, when attempting to increase the yield from breeding programs, and when estimating relatedness in conservation programs. High-density SNPs of many species of importance to science and agriculture can now be scored quickly and relatively cheaply, for example, in mice (Valdar et al. 2006), chickens (Muir et al. 2008), and dairy cattle (VanRaden et al. 2009).In the analysis of populations of breeding stock, the inclusion of dense marker data has improved the predictive ability (i.e., reliability) of genetic evaluations compared to the traditional phenotype model, both in simulations (Meuwissen et al. 2001; Calus et al. 2008; Hayes et al. 2009) and when using real data (Legarra et al. 2008; VanRaden et al. 2009; González-Recio et al. 2009). Meuwissen et al. (2001) suggested that the effect of all markers should first be estimated, and then summed, to obtain genomic estimated breeding values (GEBVs). An alternative procedure, where all markers are used to compute the genomic relationship matrix (in place of the additive polygenic relationship matrix) has also been suggested (e.g., Villanueva et al. 2005; VanRaden 2008; Hayes et al. 2009); this matrix is then incorporated into the statistical analysis to estimate GEBVs. A comparison of both procedures (VanRaden 2008) yielded similar estimates of GEBVs in cases where the effect of an individual allele was small. In addition, if not all pedigree members have marker information, a combined relationship matrix derived from both genotyped and ungenotyped individuals could be computed; this has been shown to increase the accuracy of GEBVs (Legarra et al. 2009; Misztal et al. 2009). Another plausible option to incorporate marker information is to use low-density SNP panels within families and to trace the effect of SNPs from high-density genotyped ancestors, as suggested by Habier et al. (2009) and Weigel et al. (2009). However, fast and powerful computer algorithms, which can use the marker information as efficiently as possible in the analysis of quantitative traits, are needed to obtain accurate GEBVs from genome-wide marker data.This study describes the development of an efficient Bayesian method for incorporating general relationships into the genetic evaluation procedure. The method is based on expressing the multivariate normal prior distribution as a product of one-dimensional normal distributions, each conditioned on the descending variables. When evaluating the genetic parameters of natural and breeding populations, high-dimensional distributions are often used as prior distributions of various genetic effects, such as the additive polygenic effect (Wang et al. 1993), multivariate additive polygenic effects (Van Tassell and Van Vleck 1996), and quantitative trait loci (QTL) effects via the identical-by-decent matrix (Yi and Xu 2000). A Bayesian framework is adopted to obtain posterior distributions of all unknown parameters, estimated by using Markov chain Monte Carlo (MCMC) sampling algorithms in the software package WinBUGS (Lunn et al. 2000, 2009). By performing prior calculations in the form of the factorized product of simple univariate conditional distributions, the computational time of the MCMC estimation procedure is reduced considerably. This feature permits rapid inference for both the polygenic model and the genomic relationship model. Moreover, the decomposition allows for inbreeding of varying degree, since the correct genetic covariance structure can be inferred into the analysis. In this article, we test the method on two previously published pedigree data sets: phenotype data from a large pedigree of Scots pine, incorporation of information on both additive and dominance genetic relationships (Waldmann et al. 2008); and genomic information obtained from a genome-wide scan of a simulated animal population (Lund et al. 2009).     

A Two-Stage Approximation for Analysis of Mixture Genetic Models in Large Pedigrees

Information from cosegregation of marker and QTL alleles, in addition to linkage disequilibrium (LD), can improve genomic selection. Variance components linear models have been proposed for this purpose, but accommodating dominance and epistasis is not straightforward with them. A full-Bayesian analysis of a mixture genetic model is favorable in this respect, but is computationally infeasible for whole-genome analyses. Thus, we propose an approximate two-step approach that neglects information from trait phenotypes in inferring ordered genotypes and segregation indicators of markers. Quantitative trait loci (QTL) fine-mapping scenarios, using high-density markers and pedigrees of five generations without genotyped females, were simulated to test this strategy against an exact full-Bayesian approach. The latter performed better in estimating QTL genotypes, but precision of QTL location and accuracy of genomic breeding values (GEBVs) did not differ for the two methods at realistically low LD. If, however, LD was higher, the exact approach resulted in a slightly higher accuracy of GEBVs. In conclusion, the two-step approach makes mixture genetic models computationally feasible for high-density markers and large pedigrees. Furthermore, markers need to be sampled only once and results can be used for the analysis of all traits. Further research is needed to evaluate the two-step approach for complex pedigrees and to analyze alternative strategies for modeling LD between QTL and markers.DUE to advances in molecular genetics, high-density single-nucleotide polymorphisms (SNPs) are becoming available in animal and plant breeding. These can be used for whole-genome analyses such as prediction of genomic breeding values (GEBVs) and fine mapping of quantitative trait loci (QTL). Genomic selection (GS) (Meuwissen et al. 2001) is promising to improve response to selection by exploiting linkage disequilibrium (LD) between SNPs and QTL (Hayes et al. 2009; Vanraden et al. 2009), but the accuracy of GEBVs depends on additive-genetic relationships between the individuals used to estimate SNP effects and selection candidates (Habier et al. 2007, 2010). Use of cosegregation information, in addition to LD, may reduce this dependency and improve GS. Calus et al. (2008) used a variance components linear model for this purpose in which random QTL effects are modeled conditional on marker haplotypes. The covariance between founder haplotypes allows one to include LD (Meuwissen and Goddard 2000), and the covariance between nonfounder haplotypes computed as in Fernando and Grossman (1989) allows one to include cosegregation. The resulting covariance matrices, however, can be nonpositive definite, which necessitates bending with the effect that information can be lost (Legarra and Fernando 2009). Furthermore, accommodating dominance and epistasis is not straightforward with linear models, especially for crossbred data. In contrast with mixture genetic models, genetic covariance matrices do not enter into the analysis, and accommodating dominance and epistasis is more straightforward (Goddard 1998; Pong-Wong et al. 1998; Stricker and Fernando 1998; Du et al. 1999; Du and Hoeschele 2000; Hoeschele 2001; Yi and Xu 2002; Pérez-Enciso 2003; Yi et al. 2003, 2005).Mixture model analyses, however, are more computationally demanding because the unknowns to be estimated in these analyses include the effects of unobservable QTL genotypes. In linear model analyses, in contrast, it is effects of observable marker genotypes that are estimated. The mixture model analysis can be thought of as a weighted sum of linear model analyses corresponding to each possible state for the unobservable QTL genotypes, where the weights are the probabilities of the QTL genotype states conditional on the observed marker genotypes and trait phenotypes. In practice, the analysis needs to consider all possible haplotypes at the markers also because even when all marker genotypes are observed, some of the marker haplotypes may not be known. As a result, the computational burden of these analyses stems from the number of unknown genotype and haplotype states that need to be summed over being exponentially related to the number of individuals in the pedigree and the number of loci.It can be shown that conditional on the genotypes of their parents, genotypes of offspring are independent of the genotypes of all their ancestors. This conditional independence can be exploited to efficiently compute the weighted summation in the mixture model analysis, provided the pedigree is not too complex (Lauritzen and Sheehan 2003). In genetics, this strategy is called peeling (Elston and Stewart 1971; Cannings et al. 1978) and is equivalent to variable elimination in graphical models (Lauritzen and Sheehan 2003). This approach, however, becomes infeasible when the pedigree is complex and the number of loci is large. Another strategy for analysis of mixture models is based on using Markov chain Monte Carlo (MCMC) theory to draw samples of QTL genotypes and marker haplotypes conditional on the observed marker genotypes and trait phenotypes. Pérez-Enciso (2003) developed an MCMC-based Bayesian analysis for a mixture genetic model that uses information from both LD and cosegregation to fine map a single QTL, but this approach becomes computationally infeasible for whole-genome analyses without approximations.In this article, we investigate a two-stage, approximate analysis that uses information from both LD and cosegregation. In the first stage, ordered genotypes of markers are sampled conditional only on the observed, unordered marker genotypes, ignoring information from the trait phenotypes. These samples are drawn using a Gibbs sampler with overlapping blocks (Thomas et al. 2000; Abraham et al. 2007) in which peeling is performed within a block while conditioning on variables outside the block. From these samples, founder haplotype probabilities and segregation probabilities for the QTL, also called probabilities of descent of QTL (PDQs) alleles, are calculated. In the second stage, these probabilities are used to sample QTL genotypes conditional on the trait phenotypes. In this analysis, information from LD is incorporated by allowing the QTL allele frequencies in founders to be dependent on the marker haplotypes, and information from cosegregation is incorporated by using the PDQs from the first stage to sample QTL alleles in nonfounders. The approximation comes from ignoring trait phenotypes in sampling ordered marker genotypes. A major advantage of the two-step approach is that markers have to be sampled only once and can then be used to analyze all quantitative traits with a mixture model.The objective of this study is to test the hypothesis that this approximation is negligible given high-density SNPs. To test this hypothesis, results from the two-stage, approximate analysis are compared to a full-Bayesian analysis that does not ignore the information from the trait phenotypes in sampling the ordered marker genotypes. The full-Bayesian approach was selected, because it is considered to be the ideal statistical model as it accounts for all uncertainties (Hoeschele 2001). Because the full-Bayesian approach is computationally too demanding for application to GS, the approximate and full-Bayesian analyses are used to fine map within a simulated chromosomal region that is known to contain a QTL to make the comparison computationally feasible. If the consequences of ignoring trait phenotypes to sample ordered marker genotypes are negligible, further research to apply mixture genetic models to GS and comparisons with linear models are justifiable.     

Phenotypic Consequences of Aneuploidy in Arabidopsis thaliana

Aneuploid cells are characterized by incomplete chromosome sets. The resulting imbalance in gene dosage has phenotypic consequences that are specific to each karyotype. Even in the case of Down syndrome, the most viable and studied form of human aneuploidy, the mechanisms underlying the connected phenotypes remain mostly unclear. Because of their tolerance to aneuploidy, plants provide a powerful system for a genome-wide investigation of aneuploid syndromes, an approach that is not feasible in animal systems. Indeed, in many plant species, populations of aneuploid individuals can be easily obtained from triploid individuals. We phenotyped a population of Arabidopsis thaliana aneuploid individuals containing 25 different karyotypes. Even in this highly heterogeneous population, we demonstrate that certain traits are strongly associated with the dosage of specific chromosome types and that chromosomal effects can be additive. Further, we identified subtle developmental phenotypes expressed in the diploid progeny of aneuploid parent(s) but not in euploid controls from diploid lineages. These results indicate long-term phenotypic consequences of aneuploidy that can persist after chromosomal balance has been restored. We verified the diploid nature of these individuals by whole-genome sequencing and discuss the possibility that trans-generational phenotypic effects stem from epigenetic modifications passed from aneuploid parents to their diploid progeny.THE genome of aneuploid individuals contains incomplete chromosome sets. The balance between chromosome types, and the genes they encode, is compromised, resulting in altered expression of many genes, including genes with dosage-sensitive effects on phenotypes. In humans, only a few types of aneuploid karyotypes are viable (Hassold and Hunt 2001), highlighting the deleterious effect of chromosome imbalance. The most commonly known viable form of aneuploidy in humans is Down syndrome, which results from a trisomy of chromosome 21 in an otherwise diploid background. Down syndrome patients exhibit many specific phenotypes, sometimes visible only in a subset of patients (Antonarakis et al. 2004). For phenotypes found in all Down syndrome patients, the penetrance of each phenotype varies between patients (Antonarakis et al. 2004). Despite the increasing amount of information available about the human genome and the availability of a mouse model for Down syndrome (O''Doherty et al. 2005), the genes responsible for most of the phenotypes associated with Down syndrome are still unknown (Patterson 2007; Korbel et al. 2009; Patterson 2009). Recently, detailed phenotypic analyses of as many as 30 aneuploid patients have allowed the identification of susceptibility regions for several specific phenotypes (Patterson 2007, 2009; Korbel et al. 2009; Lyle et al. 2009), but the specific genes remain to be identified. Understanding the physiology of aneuploidy is not only relevant to those individuals with aneuploid genomes but also to understanding cancer since most cancerous cells are aneuploid (Matzke et al. 2003; Pihan and Doxsey 2003; Storchova and Pellman 2004; Holland and Cleveland 2009; Williams and Amon 2009) or the consequences of copy number variation and dosage sensitivity (Dear 2009; Henrichsen et al. 2009).Plants are more tolerant of aneuploidy than animals (Matzke et al. 2003) for reasons that remain unclear. Since the discovery of the Datura trisomic “chromosome mutants” by Blakeslee (1921, 1922), viable trisomics of each chromosome type have been described in numerous species. Trisomics exhibit phenotypes specific to the identity of the triplicated chromosome (Blakeslee 1922; Khush 1973; Koornneef and Van der Veen 1983; Singh 2003). More complex aneuploids, i.e., individuals carrying more than one additional chromosome, can be viable as well and have been observed in many plants species, especially among the progeny of triploid individuals (McClintock 1929; Levan 1942; Johnsson 1945; Khush 1973). Some species appear to be more tolerant of complex aneuploidies than others, suggesting a genetic basis for aneuploidy tolerance (Satina and Blakeslee 1938; Khush 1973; Ramsey and Schemske 2002; Henry et al. 2009). Aneuploid individuals frequently appear spontaneously within polyploid plant populations, presumably due to a failure to equally partition the multiple chromosome sets at meiosis (Randolph 1935; Doyle 1986). These aneuploids exhibit few or subtle phenotypic abnormalities and can often compete with their euploid progenitors (Ramsey and Schemske 1998). Plants therefore provide an excellent opportunity for a genome-wide investigation of aneuploid syndromes: sample size is not limited, phenotypes can be described and assessed in detail, and plant aneuploid populations provide a complex mixture of viable karyotypes.In this article, we report our investigation of the relationship between phenotype and karyotype in populations of aneuploid Arabidopsis thaliana plants. All simple trisomics of A. thaliana have been previously isolated and phenotypically characterized (Steinitz-Sears 1962; Lee-Chen and Steinitz-Sears 1967; Steinitz-Sears and Lee-Chen 1970; Koornneef and Van der Veen 1983), demonstrating that they are tolerated in A. thaliana. We previously reported that aneuploid swarms—populations of aneuploid individuals of varying aneuploid karyotypes—could be obtained from the progeny of triploid A. thaliana individuals (Henry et al. 2005, 2009). Using a combination of a quantitative PCR-based method and flow cytometry, we were able to derive the full aneuploid karyotype of each of these individuals (Henry et al. 2006). We further crossed triploid A. thaliana to diploid or tetraploid individuals and demonstrated that at least 44 of the 60 possible aneuploid karyotypes that could result from these crosses (aneuploid individuals carrying between 11 and 19 chromosomes) were viable and successfully produced adult plants. Taken together, these populations and methods make it possible to explore the basis of aneuploid syndromes in A. thaliana. In this study, we were able to phenotypically characterize at least one individual from 25 different aneuploid karyotypes falling between diploidy and tetraploidy. We demonstrated that specific phenotypes are affected by the dosage of specific chromosome types. The effect of the dosage of specific chromosome types on traits was additive and could be used to predict the observed phenotype. The availability of multiple generations of aneuploid and euploid individuals allowed us to investigate potential long-term effects of aneuploidy as well as parent-of-origin effects on aneuploid phenotypes.     

The Effect of Recent Admixture on Inference of Ancient Human Population History

Despite the widespread study of genetic variation in admixed human populations, such as African-Americans, there has not been an evaluation of the effects of recent admixture on patterns of polymorphism or inferences about population demography. These issues are particularly relevant because estimates of the timing and magnitude of population growth in Africa have differed among previous studies, some of which examined African-American individuals. Here we use simulations and single-nucleotide polymorphism (SNP) data collected through direct resequencing and genotyping to investigate these issues. We find that when estimating the current population size and magnitude of recent growth in an ancestral population using the site frequency spectrum (SFS), it is possible to obtain reasonably accurate estimates of the parameters when using samples drawn from the admixed population under certain conditions. We also show that methods for demographic inference that use haplotype patterns are more sensitive to recent admixture than are methods based on the SFS. The analysis of human genetic variation data from the Yoruba people of Ibadan, Nigeria and African-Americans supports the predictions from the simulations. Our results have important implications for the evaluation of previous population genetic studies that have considered African-American individuals as a proxy for individuals from West Africa as well as for future population genetic studies of additional admixed populations.STUDIES of archeological and genetic data show that anatomically modern humans originated in Africa and more recently left Africa to populate the rest of the world (Tishkoff and Williams 2002; Barbujani and Goldstein 2004; Garrigan and Hammer 2006; Reed and Tishkoff 2006; Campbell and Tishkoff 2008; Jakobsson et al. 2008; Li et al. 2008). Given the central role Africa has played in the origin of diverse human populations, understanding patterns of genetic variation and the demographic history of populations within Africa is important for understanding the demographic history of global human populations. The availability of large-scale single-nucleotide polymorphism (SNP) data sets coupled with recent advances in statistical methodology for inferring parameters in population genetic models provides a powerful means of accomplishing these goals (Keinan et al. 2007; Boyko et al. 2008; Lohmueller et al. 2009; Nielsen et al. 2009).It is important to realize that studies of African demographic history using genetic data have come to qualitatively different conclusions regarding important parameters. Some recent studies have found evidence for ancient (>100,000 years ago) two- to fourfold growth in African populations (Adams and Hudson 2004; Marth et al. 2004; Keinan et al. 2007; Boyko et al. 2008). Other studies have found evidence of very recent growth (Pluzhnikov et al. 2002; Akey et al. 2004; Voight et al. 2005; Cox et al. 2009; Wall et al. 2009) or could not reject a model with a constant population size (Pluzhnikov et al. 2002; Voight et al. 2005). It is unclear why studies found such different parameter estimates. However, these studies all differ from each other in the amount of data considered, the types of data used (e.g., SNP genotypes vs. full resequencing), the genomic regions studied (e.g., noncoding vs. coding SNPs), and the types of demographic models considered (e.g., including migration vs. not including migration postseparation of African and non-African populations).Another important way in which studies of African demographic history differ from each other is in the populations sampled. Some studies have focused on genetic data from individuals sampled from within Africa (Pluzhnikov et al. 2002; Adams and Hudson 2004; Voight et al. 2005; Keinan et al. 2007; Cox et al. 2009; Wall et al. 2009), while other studies included American individuals with African ancestry (Adams and Hudson 2004; Akey et al. 2004; Marth et al. 2004; Boyko et al. 2008). While there is no clear correspondence between those studies which sampled native African individuals (as opposed to African-Americans) and particular growth scenarios, it is clear from previous studies that African-American populations do differ from African populations in their recent demographic history. In particular, genetic studies suggest that there is wide variation in the degree of European admixture in most African-American individuals in the United States and that they have, on average, ∼80% African ancestry and 20% European ancestry (Parra et al. 1998; Pfaff et al. 2001; Falush et al. 2003; Patterson et al. 2004; Tian et al. 2006; Lind et al. 2007; Reiner et al. 2007; Price et al. 2009; Bryc et al. 2010). Furthermore, both historical records and genetic evidence suggest that the admixture process began quite recently, within the last 20 generations (Pfaff et al. 2001; Patterson et al. 2004; Seldin et al. 2004; Tian et al. 2006). Recent population admixture can alter patterns of genetic variation in a discernible and predictable way. For example, recently admixed populations will exhibit correlation in allele frequencies (i.e., linkage disequilibrium) among markers that differ in frequency between the parental populations. This so-called admixture linkage disequilibrium (LD) (Chakraborty and Weiss 1988) can extend over long physical distances (Lautenberger et al. 2000) and decays exponentially with time the since the admixture process began (i.e., recently admixed populations typically exhibit LD over a longer physical distance than anciently admixed populations).While it is clear that African-American populations have a different recent demographic history than do African populations from within Africa and that admixture tracts can be identified in admixed individuals (Falush et al. 2003; Patterson et al. 2004; Tang et al. 2006; Sankararaman et al. 2008a,b; Price et al. 2009; Bryc et al. 2010), the effect that admixture has on other patterns of genetic variation remains unclear. For example, Xu et al. (2007) found similar LD decay patterns when comparing African-American and African populations. It is also unclear whether the recent admixture affects our ability to reconstruct ancient demographic events (such as expansions that predate the spread of humans out of Africa) from whole-genome SNP data. Most studies of demographic history have summarized the genome-wide SNP data by allele frequency or haplotype summary statistics. If these summary statistics are not sensitive to the recent European admixture, then the African-American samples may yield estimates of demographic parameters that are close to the true demographic parameters for the ancestral, unsampled, African populations. This would suggest that the differences in growth parameter estimates obtained from African populations cannot be explained by certain studies sampling African-American individuals and others sampling African individuals from within Africa. However, if these statistics are sensitive to recent admixture, then they may give biased estimates of growth parameters.Here, we examine the effect of recent admixture on the estimation of population demography. In particular, we estimate growth parameters from simulated data sets using SNP frequencies as well as a recently developed haplotype summary statistic (Lohmueller et al. 2009). We compare the demographic parameter estimates made from the admixed and nonadmixed populations and find that some parameter estimates are qualitatively similar between the two populations when inferred using allele frequencies. Inferences of growth using haplotype-based approaches appear to be more sensitive to recent admixture than inferences based on SNP frequencies. We discuss implications that our results have for interpreting studies of demography in admixed populations.     

The Conserved miR-51 microRNA Family Is Redundantly Required for Embryonic Development and Pharynx Attachment in Caenorhabditis elegans

A major question about cytokinesis concerns the role of the septin proteins, which localize to the division site in all animal and fungal cells but are essential for cytokinesis only in some cell types. For example, in Schizosaccharomyces pombe, four septins localize to the division site, but deletion of the four genes produces only a modest delay in cell separation. To ask if the S. pombe septins function redundantly in cytokinesis, we conducted a synthetic-lethal screen in a septin-deficient strain and identified seven mutations. One mutation affects Cdc4, a myosin light chain that is an essential component of the cytokinetic actomyosin ring. Five others cause frequent cell lysis during cell separation and map to two loci. These mutations and their dosage suppressors define a signaling pathway (including Rho1 and a novel arrestin) for repairing cell-wall damage. The seventh mutation affects the poorly understood RNA-binding protein Scw1 and severely delays cell separation when combined either with a septin mutation or with a mutation affecting the septin-interacting, anillin-like protein Mid2, suggesting that Scw1 functions in a pathway parallel to that of the septins. Taken together, our results suggest that the S. pombe septins participate redundantly in one or more pathways that cooperate with the actomyosin ring during cytokinesis and that a septin defect causes septum defects that can be repaired effectively only when the cell-integrity pathway is intact.THE fission yeast Schizosaccharomyces pombe provides an outstanding model system for studies of cytokinesis (McCollum and Gould 2001; Balasubramanian et al. 2004; Pollard and Wu 2010). As in most animal cells, successful cytokinesis in S. pombe requires an actomyosin ring (AMR). The AMR begins to assemble at the G₂/M transition and involves the type II myosin heavy chains Myo2 and Myp2 and the light chains Cdc4 and Rlc1 (Wu et al. 2003). Myo2 and Cdc4 are essential for cytokinesis under all known conditions, Rlc1 is important at all temperatures but essential only at low temperatures, and Myp2 is essential only under stress conditions. As the AMR constricts, a septum of cell wall is formed between the daughter cells. The primary septum is sandwiched by secondary septa and subsequently digested to allow cell separation (Humbel et al. 2001; Sipiczki 2007). Because of the internal turgor pressure of the cells, the proper assembly and structural integrity of the septal layers are essential for cell survival.Septum formation involves the β-glucan synthases Bgs1/Cps1/Drc1, Bgs3, and Bgs4 (Ishiguro et al. 1997; Le Goff et al. 1999; Liu et al. 1999, 2002; Martín et al. 2003; Cortés et al. 2005) and the α-glucan synthase Ags1/Mok1 (Hochstenbach et al. 1998; Katayama et al. 1999). These synthases are regulated by the Rho GTPases Rho1 and Rho2 and the protein kinase C isoforms Pck1 and Pck2 (Arellano et al. 1996, 1997, 1999; Nakano et al. 1997; Hirata et al. 1998; Calonge et al. 2000; Sayers et al. 2000; Ma et al. 2006; Barba et al. 2008; García et al. 2009b). The Rho GTPases themselves appear to be regulated by both GTPase-activating proteins (GAPs) and guanine-nucleotide-exchange factors (GEFs) (Nakano et al. 2001; Calonge et al. 2003; Iwaki et al. 2003; Tajadura et al. 2004; Morrell-Falvey et al. 2005; Mutoh et al. 2005; García et al. 2006, 2009a,b). In addition, septum formation and AMR function appear to be interdependent. In the absence of a normal AMR, cells form aberrant septa and/or deposit septal materials at random locations, whereas a mutant defective in septum formation (bgs1) is also defective in AMR constriction (Gould and Simanis 1997; Le Goff et al. 1999; Liu et al. 1999, 2000). Both AMR constriction and septum formation also depend on the septation initiation network involving the small GTPase Spg1 (McCollum and Gould 2001; Krapp and Simanis 2008). Despite this considerable progress, many questions remain about the mechanisms and regulation of septum formation and its relationships to the function of the AMR.One major question concerns the role(s) of the septins. Proteins of this family are ubiquitous in fungal and animal cells and typically localize to the cell cortex, where they appear to serve as scaffolds and diffusion barriers for other proteins that participate in a wide variety of cellular processes (Longtine et al. 1996; Gladfelter et al. 2001; Hall et al. 2008; Caudron and Barral 2009). Despite the recent progress in elucidating the mechanisms of septin assembly (John et al. 2007; Sirajuddin et al. 2007; Bertin et al. 2008; McMurray and Thorner 2008), the details of septin function remain obscure. However, one prominent role of the septins and associated proteins is in cytokinesis. Septins concentrate at the division site in every cell type that has been examined, and in Saccharomyces cerevisiae (Hartwell 1971; Longtine et al. 1996; Lippincott et al. 2001; Dobbelaere and Barral 2004) and at least some Drosophila (Neufeld and Rubin 1994; Adam et al. 2000) and mammalian (Kinoshita et al. 1997; Surka et al. 2002) cell types, the septins are essential for cytokinesis. In S. cerevisiae, the septins are required for formation of the AMR (Bi et al. 1998; Lippincott and Li 1998). However, this cannot be their only role, because the AMR itself is not essential for cytokinesis in this organism (Bi et al. 1998; Korinek et al. 2000; Schmidt et al. 2002). Moreover, there is no evidence that the septins are necessary for AMR formation or function in any other organism. A further complication is that in some cell types, including most Caenorhabditis elegans cells (Nguyen et al. 2000; Maddox et al. 2007) and some Drosophila cells (Adam et al. 2000; Field et al. 2008), the septins do not appear to be essential for cytokinesis even though they localize to the division site.S. pombe has seven septins, four of which (Spn1, Spn2, Spn3, and Spn4) are expressed in vegetative cells and localize to the division site shortly before AMR constriction and septum formation (Longtine et al. 1996; Berlin et al. 2003; Tasto et al. 2003; Wu et al. 2003; An et al. 2004; Petit et al. 2005; Pan et al. 2007; Onishi et al. 2010). Spn1 and Spn4 appear to be the core members of the septin complex (An et al. 2004; McMurray and Thorner 2008), and mutants lacking either of these proteins do not assemble the others at the division site. Assembly of a normal septin ring also depends on the anillin-like protein Mid2, which colocalizes with the septins (Berlin et al. 2003; Tasto et al. 2003). Surprisingly, mutants lacking the septins are viable and form seemingly complete septa with approximately normal timing. These mutants do, however, display a variable delay in separation of the daughter cells, suggesting that the septins play some role(s) in the proper completion of the septum or in subsequent processes necessary for cell separation (Longtine et al. 1996; An et al. 2004; Martín-Cuadrado et al. 2005).It is possible that the septins localize to the division site and yet are nonessential for division in some cell types because their role is redundant with that of some other protein(s) or pathway(s). To explore this possibility in S. pombe, we screened for mutations that were lethal in combination with a lack of septins. The results suggest that the septins cooperate with the AMR during cytokinesis and that, in the absence of septin function, the septum is not formed properly, so that an intact system for recognizing and repairing cell-wall damage becomes critical for cell survival.     

Polymorphic Genes of Major Effect: Consequences for Variation,Selection and Evolution in Arabidopsis thaliana

The importance of genes of major effect for evolutionary trajectories within and among natural populations has long been the subject of intense debate. For example, if allelic variation at a major-effect locus fundamentally alters the structure of quantitative trait variation, then fixation of a single locus can have rapid and profound effects on the rate or direction of subsequent evolutionary change. Using an Arabidopsis thaliana RIL mapping population, we compare G-matrix structure between lines possessing different alleles at ERECTA, a locus known to affect ecologically relevant variation in plant architecture. We find that the allele present at ERECTA significantly alters G-matrix structure—in particular the genetic correlations between branch number and flowering time traits—and may also modulate the strength of natural selection on these traits. Despite these differences, however, when we extend our analysis to determine how evolution might differ depending on the ERECTA allele, we find that predicted responses to selection are similar. To compare responses to selection between allele classes, we developed a resampling strategy that incorporates uncertainty in estimates of selection that can also be used for statistical comparisons of G matrices.THE structure of the genetic variation that underlies phenotypic traits has important consequences for understanding the evolution of quantitative traits (Fisher 1930; Lande 1979; Bulmer 1980; Kimura 1983; Orr 1998; Agrawal et al. 2001). Despite the infinitesimal model''s allure and theoretical tractability (see Orr and Coyne 1992; Orr 1998, 2005a,b for reviews of its influence), evidence has accumulated from several sources (artificial selection experiments, experimental evolution, and QTL mapping) to suggest that genes of major effect often contribute to quantitative traits. Thus, the frequency and role of genes of major effect in evolutionary quantitative genetics have been a subject of intense debate and investigation for close to 80 years (Fisher 1930; Kimura 1983; Orr 1998, 2005a,b). Beyond the conceptual implications, the prevalence of major-effect loci also affects our ability to determine the genetic basis of adaptations and species differences (e.g., Bradshaw et al. 1995, 1998).Although the existence of genes of major effect is no longer in doubt, we still lack basic empirical data on how segregating variation at such genes affects key components of evolutionary process (but see Carrière and Roff 1995). In other words, How does polymorphism at genes of major effect alter patterns of genetic variation and covariation, natural selection, and the likely response to selection? The lack of data stems, in part, from the methods used to detect genes of major effect: experimental evolution (e.g., Bull et al. 1997; Zeyl 2005) and QTL analysis (see Erickson et al. 2004 for a review) often detect such genes retrospectively after they have become fixed in experimental populations or the species pairs used to generate the mapping population. The consequences of polymorphism at these genes on patterns of variation, covariation, selection, and the response to selection—which can be transient (Agrawal et al. 2001)—are thus often unobserved.A partial exception to the absence of data on the effects of major genes comes from artificial selection experiments, in which a substantial evolutionary response to selection in the phenotype after a plateau is often interpreted as evidence for the fixation of a major-effect locus (Frankham et al. 1968; Yoo 1980a,b; Frankham 1980; Shrimpton and Robertson 1988a,b; Caballero et al. 1991; Keightley 1998; see Mackay 1990 and Hill and Caballero 1992 for reviews). However, many of these experiments report only data on the selected phenotype (e.g., bristle number) or, alternatively, the selected phenotype and some measure of fitness (e.g., Frankham et al. 1968, Yoo 1980b; Caballero et al. 1991; Mackay et al. 1994; Fry et al. 1995; Nuzhdin et al. 1995; Zur Lage et al. 1997), making it difficult to infer how a mutation will affect variation, covariation, selection, and evolutionary responses for a suite of traits that might affect fitness themselves. One approach is to document how variation at individual genes of major effect affects the genetic variance–covariance matrix (“G matrix”; Lande 1979), which represents the additive genetic variance and covariance between traits.Although direct evidence for variation at major-effect genes altering patterns of genetic variation, covariation, and selection is rare, there is abundant evidence for the genetic mechanisms that could produce these dynamics. A gene of major effect could have these consequences due to any of at least three genetic mechanisms: (1) pleiotropy, where a gene of major effect influences several traits, including potentially fitness, simultaneously, (2) physical linkage or linkage disequilibrium (LD), in which a gene of major effect is either physically linked or in LD with other genes that influence other traits under selection, and (3) epistasis, in which the allele present at a major-effect gene alters the phenotypic effect of other loci and potentially phenotypes under selection. Evidence for these three evolutionary genetic mechanisms leading to changes in suites of traits comes from a variety of sources, including mutation accumulation experiments (Clark et al. 1995; Fernandez and Lopez-Fanjul 1996), mutation induction experiments (Keightley and Ohnishi 1998), artificial selection experiments (Long et al. 1995), and transposable element insertions (Rollmann et al. 2006). For pleiotropy in particular, major-effect genes that have consequences on several phenotypic traits are well known from the domestication and livestock breeding literature [e.g., myostatin mutations in Belgian blue cattle and whippets (Arthur 1995; Grobet et al. 1997; Mosher et al. 2007), halothane genes in pigs (Christian and Rothschild 1991; Fujii et al. 1991), and Booroola and Inverdale genes in sheep (Amer et al. 1999; Visscher et al. 2000)]. While these data suggest that variation at major-effect genes could—and probably does—influence variation, covariation, and selection on quantitative traits, data on the magnitude of these consequences remain lacking.Recombinant inbred line (RIL) populations are a promising tool for investigating the influence of major-effect loci. During advancement of the lines from F₂''s to RILs, alternate alleles at major-effect genes (and most of the rest of the genome) will be made homozygous, simplifying comparisons among genotypic classes. Because of the high homozygosity, individuals within RILs are nearly genetically identical, facilitating phenotyping of many genotypes under a range of environments. In addition, because of recombination, alternative alleles are randomized across genetic backgrounds—facilitating robust comparisons between sets of lines differing at a major-effect locus.Here we investigate how polymorphism at an artificially induced mutation, the erecta locus in Arabidopsis thaliana, affects the magnitude of these important evolutionary genetic parameters under ecologically realistic field conditions. We use the Landsberg erecta (Ler) × Columbia (Col) RIL population of A. thaliana to examine how variation at a gene of major effect influences genetic variation, covariation, and selection on quantitative traits in a field setting. The Ler × Col RIL population is particularly suitable, because it segregates for an artificially induced mutation at the erecta locus, which has been shown to influence a wide variety of plant traits. The Ler × Col population thus allows a powerful test of the effects of segregating variation at a gene—chosen a priori—with numerous pleiotropic effects. The ERECTA gene is a leucine-rich receptor-like kinase (LRR-RLK) (Torii et al. 1996) and has been shown to affect plant growth rates (El-Lithy et al. 2004), stomatal patterning and transpiration efficiency (Masle et al. 2005; Shpak et al. 2005), bacterial pathogen resistance (Godiard et al. 2003), inflorescence and floral organ size and shape (Douglas et al. 2002; Shpak et al. 2003, 2004), and leaf polarity (Xu et al. 2003; Qi et al. 2004).Specifically, we sought to answer the following questions: (1) Is variation at erecta significantly associated with changes to the G matrix? (2) Is variation at erecta associated with changes in natural selection on genetically variable traits? And (3) is variation at erecta associated with significantly different projected evolutionary responses to selection?     

Prediction of Genetic Values of Quantitative Traits in Plant Breeding Using Pedigree and Molecular Markers

The availability of dense molecular markers has made possible the use of genomic selection (GS) for plant breeding. However, the evaluation of models for GS in real plant populations is very limited. This article evaluates the performance of parametric and semiparametric models for GS using wheat (Triticum aestivum L.) and maize (Zea mays) data in which different traits were measured in several environmental conditions. The findings, based on extensive cross-validations, indicate that models including marker information had higher predictive ability than pedigree-based models. In the wheat data set, and relative to a pedigree model, gains in predictive ability due to inclusion of markers ranged from 7.7 to 35.7%. Correlation between observed and predictive values in the maize data set achieved values up to 0.79. Estimates of marker effects were different across environmental conditions, indicating that genotype × environment interaction is an important component of genetic variability. These results indicate that GS in plant breeding can be an effective strategy for selecting among lines whose phenotypes have yet to be observed.PEDIGREE-BASED prediction of genetic values based on the additive infinitesimal model (Fisher 1918) has played a central role in genetic improvement of complex traits in plants and animals. Animal breeders have used this model for predicting breeding values either in a mixed model (best linear unbiased prediction, BLUP) (Henderson 1984) or in a Bayesian framework (Gianola and Fernando 1986). More recently, plant breeders have incorporated pedigree information into linear mixed models for predicting breeding values (Crossa et al. 2006, 2007; Oakey et al. 2006; Burgueño et al. 2007; Piepho et al. 2007).The availability of thousands of genome-wide molecular markers has made possible the use of genomic selection (GS) for prediction of genetic values (Meuwissen et al. 2001) in plants (e.g., Bernardo and Yu 2007; Piepho 2009; Jannink et al. 2010) and animals (Gonzalez-Recio et al. 2008; VanRaden et al. 2008; Hayes et al. 2009; de los Campos et al. 2009a). Implementing GS poses several statistical and computational challenges, such as how models can cope with the curse of dimensionality, colinearity between markers, or the complexity of quantitative traits. Parametric (e.g., Meuwissen et al. 2001) and semiparametric (e.g., Gianola et al. 2006; Gianola and van Kaam 2008) methods address these problems differently.In standard genetic models, phenotypic outcomes, , are viewed as the sum of a genetic value, , and a model residual, ; that is, . In parametric models for GS, is described as a regression on marker covariates (j = 1,  …  , p molecular markers) of the form , such that(or , in matrix notation), where is the regression of on the jth marker covariate .Estimation of via multiple regression by ordinary least squares (OLS) is not feasible when p > n. A commonly used alternative is to estimate marker effects jointly using penalized methods such as ridge regression (Hoerl and Kennard 1970) or the Least Absolute Shrinkage and Selection Operator (LASSO) (Tibshirani 1996) or their Bayesian counterpart. This approach yields greater accuracy of estimated genetic values and can be coupled with geostatistical techniques commonly used in plant breeding to model multienvironments trials (Piepho 2009).In ridge regression (or its Bayesian counterpart) the extent of shrinkage is homogeneous across markers, which may not be appropriate if some markers are located in regions that are not associated with genetic variance, while markers in other regions may be linked to QTL (Goddard and Hayes 2007). To overcome this limitation, many authors have proposed methods that use marker-specific shrinkage. In a Bayesian setting, this can be implemented using priors of marker effects that are mixtures of scaled-normal densities. Examples of this are methods Bayes A and Bayes B of Meuwissen et al. (2001) and the Bayesian LASSO of Park and Casella (2008).An alternative to parametric regressions is to use semiparametric methods such as reproducing kernel Hilbert spaces (RKHS) regression (Gianola and van Kaam 2008). The Bayesian RKHS regression regards genetic values as random variables coming from a Gaussian process centered at zero and with a (co)variance structure that is proportional to a kernel matrix K (de los Campos et al. 2009b); that is, , where , are vectors of marker genotypes for the ith and jth individuals, respectively, and is a positive definite function evaluated in marker genotypes. In a finite-dimensional setting this amounts to modeling the vector of genetic values, , as multivariate normal; that is, where is a variance parameter. One of the most attractive features of RKHS regression is that the methodology can be used with almost any information set (e.g., covariates, strings, images, graphs). A second advantage is that with RKHS the model is represented in terms of n unknowns, which gives RKHS a great computational advantage relative to some parametric methods, especially when p ≫ n.This study presents an evaluation of several methods for GS, using two extensive data sets. One contains phenotypic records of a series of wheat trials and recently generated genomic data. The other data set pertains to international maize trials in which different traits were measured in maize lines evaluated under severe drought and well-watered conditions.     

Age Structure,Changing Demography and Effective Population Size in Atlantic Salmon (Salmo salar)

Effective population size (N_e) is a central evolutionary concept, but its genetic estimation can be significantly complicated by age structure. Here we investigate N_e in Atlantic salmon (Salmo salar) populations that have undergone changes in demography and population dynamics, applying four different genetic estimators. For this purpose we use genetic data (14 microsatellite markers) from archived scale samples collected between 1951 and 2004. Through life table simulations we assess the genetic consequences of life history variation on N_e. Although variation in reproductive contribution by mature parr affects age structure, we find that its effect on N_e estimation may be relatively minor. A comparison of estimator models suggests that even low iteroparity may upwardly bias N_e estimates when ignored (semelparity assumed) and should thus empirically be accounted for. Our results indicate that N_e may have changed over time in relatively small populations, but otherwise remained stable. Our ability to detect changes in N_e in larger populations was, however, likely hindered by sampling limitations. An evaluation of N_e estimates in a demographic context suggests that life history diversity, density-dependent factors, and metapopulation dynamics may all affect the genetic stability of these populations.THE effective size of a population (N_e) is an evolutionary parameter that can be informative on the strength of stochastic evolutionary processes, the relevance of which relative to deterministic forces has been debated for decades (e.g., Lande 1988). Stochastic forces include environmental, demographic, and genetic components, the latter two of which are thought to be more prominent at reduced population size, with potentially detrimental consequences for average individual fitness and population persistence (Newman and Pilson 1997; Saccheri et al. 1998; Frankham 2005). The quantification of N_e in conservation programs is thus frequently advocated (e.g., Luikart and Cornuet 1998; Schwartz et al. 2007), although gene flow deserves equal consideration given its countering effects on genetic stochasticity (Frankham et al. 2003; Palstra and Ruzzante 2008).Effective population size is determined mainly by the lifetime reproductive success of individuals in a population (Wright 1938; Felsenstein 1971). Variance in reproductive success, sex ratio, and population size fluctuations can reduce N_e below census population size (Frankham 1995). Given the difficulty in directly estimating N_e through quantification of these demographic factors (reviewed by Caballero 1994), efforts have been directed at inferring N_e indirectly through measurement of its genetic consequences (see Leberg 2005, Wang 2005, and Palstra and Ruzzante 2008 for reviews). Studies employing this approach have quantified historical levels of genetic diversity and genetic threats to population persistence (e.g., Nielsen et al. 1999b; Miller and Waits 2003; Johnson et al. 2004). N_e has been extensively studied in (commercially important) fish species, due to the common availability of collections of archived samples that facilitate genetic estimation using the temporal method (e.g., Hauser et al. 2002; Shrimpton and Heath 2003; Gomez-Uchida and Banks 2006; Saillant and Gold 2006).Most models relating N_e to a population''s genetic behavior make simplifying assumptions regarding population dynamics. Chiefly among these is the assumption of discrete generations, frequently violated in practice given that most natural populations are age structured with overlapping generations. Here, theoretical predictions still apply, provided that population size and age structure are constant (Felsenstein 1971; Hill 1972). Ignored age structure can introduce bias into temporal genetic methods for the estimation of N_e, especially for samples separated by time spans that are short relative to generation interval (Jorde and Ryman 1995; Waples and Yokota 2007; Palstra and Ruzzante 2008). Moreover, estimation methods that do account for age structure (e.g., Jorde and Ryman 1995) still assume this structure to be constant. Population dynamics will, however, likely be altered as population size changes, thus making precise quantifications of the genetic consequences of acute population declines difficult (Nunney 1993; Engen et al. 2005; Waples and Yokota 2007). This problem may be particularly relevant when declines are driven by anthropogenic impacts, such as selective harvesting regimes, that can affect age structure and N_e simultaneously (Ryman et al. 1981; Allendorf et al. 2008). Demographic changes thus have broad conservation implications, as they can affect a population''s sensitivity to environmental stochasticity and years of poor recruitment (Warner and Chesson 1985; Ellner and Hairston 1994; Gaggiotti and Vetter 1999). Consequently, although there is an urgent need to elucidate the genetic consequences of population declines, relatively little is understood about the behavior of N_e when population dynamics change (but see Engen et al. 2005, 2007).Here we focus on age structure and N_e in Atlantic salmon (Salmo salar) river populations in Newfoundland and Labrador. The freshwater habitat in this part of the species'' distribution range is relatively pristine (Parrish et al. 1998), yet Atlantic salmon in this area have experienced demographic declines, associated with a commercial marine fishery, characterized by high exploitation rates (40–80% of anadromous runs; Dempson et al. 2001). A fishery moratorium was declared in 1992, with rivers displaying differential recovery patterns since then (Dempson et al. 2004b), suggesting a geographically variable impact of deterministic and stochastic factors, possibly including genetics. An evaluation of those genetic consequences thus requires accounting for potential changes in population dynamics as well as in life history. Life history in Atlantic salmon can be highly versatile (Fleming 1996; Hutchings and Jones 1998; Fleming and Reynolds 2004), as exemplified by the high variation in age-at-maturity displayed among and within populations (Hutchings and Jones 1998), partly reflecting high phenotypic plasticity (Hutchings 2004). This diversity is particularly evident in the reproductive biology of males, which can mature as parr during juvenile freshwater stages (Jones and King 1952; Fleming and Reynolds 2004) and/or at various ages as anadromous individuals, when returning to spawn in freshwater from ocean migration. Variability in life history strategies is further augmented by iteroparity, which can be viewed as a bet-hedging strategy to deal with environmental uncertainty (e.g., Orzack and Tuljapurkar 1989; Fleming and Reynolds 2004). Life history diversity and plasticity may allow salmonid fish populations to alter and optimize their life history under changing demography and population dynamics, potentially acting to stabilize N_e. Reduced variance in individual reproductive success at low breeder abundance (genetic compensation) will achieve similar effects and might be a realistic aspect of salmonid breeding systems (Ardren and Kapuscinski 2003; Fraser et al. 2007b). Little is currently known about the relationships between life history plasticity, demographic change and N_e, partly due to scarcity of the multivariate data required for these analyses.Our objective in this article is twofold. First, we use demographic data for rivers in Newfoundland to quantify how life history variation influences age structure in Atlantic salmon and hence N_e and its empirical estimation from genetic data. We find that variation in reproductive contribution by mature parr has a much smaller effect on the estimation of N_e than is often assumed. Second, we use temporal genetic data to estimate N_e and quantify the genetic consequences of demographic changes. We attempt to account for potential sources of bias, associated with (changes in) age structure and life history, by using four different analytical models to estimate N_e: a single-sample estimator using the linkage disequilibrium method (Hill 1981), the temporal model assuming discrete generations (Nei and Tajima 1981; Waples 1989), and two temporal models for species with overlapping generations (Waples 1990a,b; Jorde and Ryman 1995) that differ principally in assumptions regarding iteroparity. A comparison of results from these different estimators suggests that iteroparity may often warrant analytical consideration, even when it is presumably low. Although sometimes limited by statistical power, a quantification and comparison of temporal changes in N_e among river populations suggests a more prominent impact of demographic changes on N_e in relatively small river populations.     

Deletion of a Novel F-Box Protein,MUS-10, in Neurospora crassa Leads to Altered Mitochondrial Morphology,Instability of mtDNA and Senescence

While mitochondria are renowned for their role in energy production, they also perform several other integral functions within the cell. Thus, it is not surprising that mitochondrial dysfunction can negatively impact cell viability. Although mitochondria have received an increasing amount of attention in recent years, there is still relatively little information about how proper maintenance of mitochondria and its genomes is achieved. The Neurospora crassa mus-10 mutant was first identified through its increased sensitivity to methyl methanesulfonate (MMS) and was thus believed to be defective in some aspect of DNA repair. Here, we report that mus-10 harbors fragmented mitochondria and that it accumulates deletions in its mitochondrial DNA (mtDNA), suggesting that the mus-10 gene product is involved in mitochondrial maintenance. Interestingly, mus-10 begins to senesce shortly after deletions are visualized in its mtDNA. To uncover the function of MUS-10, we used a gene rescue approach to clone the mus-10 gene and discovered that it encodes a novel F-box protein. We show that MUS-10 interacts with a core component of the Skp, Cullin, F-box containing (SCF) complex, SCON-3, and that its F-box domain is essential for its function in vivo. Thus, we provide evidence that MUS-10 is part of an E3 ubiquitin ligase complex involved in maintaining the integrity of mitochondria and may function to prevent cellular senescence.THE mus-10 mutant was isolated from a screen aimed at identifying Neurospora crassa strains that were sensitive to MMS and therefore likely to lack proper DNA repair mechanisms (Kafer and Perlmutter 1980). Epistasis analyses involving mus-10 suggested that it belonged to the uvs-6 epistasis group, which functions in recombination repair (Kafer and Perlmutter 1980; Kafer 1983). However, mus-10 did not display several phenotypes common to other members of the uvs-6 epistasis group: chromosomal instability, a high sensitivity to histidine, and the inability to produce viable ascospores in homozygous crosses (Newmeyer et al. 1978; Newmeyer and Galeazzi 1978; Kafer and Perlmutter 1980; Kafer 1981; Schroeder 1986; Watanabe et al. 1997; Handa et al. 2000; Sakuraba et al. 2000). Furthermore, the frequencies of spontaneous and radiation-induced mutation observed in mus-10 were similar to those of a wild-type strain (Kafer 1981). Past efforts to uncover the nature of these discrepancies or the function of the mus-10 gene product have been uninformative.The majority of cellular ATP is produced in mitochondria through aerobic respiration, which couples electron flow through respiratory complexes within the mitochondrial inner membrane with oxidative phosphorylation. Besides their role in ATP synthesis, mitochondria are also involved in many other cellular processes including beta-oxidation (Bartlett and Eaton 2004), calcium homeostasis (Gunter et al. 2004; Rimessi et al. 2008), production of iron-sulfur clusters (Zheng et al. 1998; Gerber and Lill 2002; Lill and Muhlenhoff 2005; Rouault and Tong 2005), and apoptosis (Green 2005; Antignani and Youle 2006; Xu and Shi 2007). Although virtually all mitochondrial proteins are encoded within the nucleus, a small number of proteins are encoded by mitochondrial DNA (mtDNA). The integrity of the mitochondrial genome may affect cell survival as mutations in mtDNA accumulate in patients suffering from severe neurological diseases including Alzheimer''s, Huntington''s and Parkinson''s, as well as several types of cancer (Chatterjee et al. 2006; Higuchi 2007; Krishnan et al. 2007; Reeve et al. 2008). The number of mtDNA mutations also increases with age, suggesting a link between mitochondrial dysfunction and ageing (Cortopassi and Arnheim 1990; Corral-Debrinski et al. 1992; Cortopassi et al. 1992; Simonetti et al. 1992; Reeve et al. 2008). Contrary to the single genome in the nucleus, there are several copies of mtDNA in each mitochondrion. Thus, defects in a few mitochondrial genomes do not necessarily lead to mitochondrial dysfunction. Many patients suffering from mitochondrial diseases exhibit heteroplasmy, a phenomenon in which a mixture of wild-type and mutant mtDNAs exist in a single cell. The ratio of wild-type to mutant mtDNAs is critical in determining the penetrance of the genetic defect, where mutant loads >60% are required to cause respiratory chain dysfunction within an individual cell (Boulet et al. 1992; Chomyn et al. 1992; Sciacco et al. 1994).Even though N. crassa strains are generally deemed immortal if they can be subcultured ∼50 times, a wild-type strain was recently reported to senesce after 12,000 hr of growth, implying that this fungus undergoes natural or programmed ageing (Maheshwari and Navaraj 2008; Kothe et al. 2010). However, replicative life span is also influenced by genetic background as certain mutations can cause progressive deterioration of growth, ultimately leading to death. One such example is the nuclear-encoded natural death (nd), which when mutant causes a senescence phenotype correlating with the accumulation of multiple mtDNA deletions (Sheng 1951; Seidel-Rogol et al. 1989). The deletions of mtDNA in nd occurred between two 70- to 701-bp direct repeats, suggesting that the nd gene product regulates recombination, repair, or replication of mtDNA (Bertrand et al. 1993). Another nuclear mutation, senescence (sen), was isolated from N. intermedia and introgressed into N. crassa (Navaraj et al. 2000). Deletions were also observed in the mtDNA of sen mutants, but unlike those occurring in nd were flanked by 6- to 10-bp repeats typically associated with GC-rich palindromic sequences (D''Souza et al. 2005). The nature of the sequences that flanked the mtDNA deletions in these two mutants supported the existence of two distinct systems of mtDNA recombination in N. crassa: a general system of homologous recombination (system I) and a site-specific mechanism (system II), mediated in part by nd and sen, respectively (Bertrand et al. 1993; D''Souza et al. 2005). The nd and sen mutations have been mapped to linkage groups I and V, respectively, but neither gene has been cloned and the precise function of their gene products remains unclear. Two ultraviolet (UV)-sensitive mutants, uvs-4 and uvs-5, are thought to undergo senescence, but unfortunately, these strains have not been studied in great detail (Schroeder 1970; Perkins et al. 1993; Hausner et al. 2006). Premature senescence has also been observed in cytoplasmic mutants of N. crassa including the E35 and ER-3 stopper mutants that harbor large mtDNA deletions, as well as strains that accumulate mitochondrial plasmids capable of inserting into mtDNA through homologous recombination (de Vries et al. 1986; Akins et al. 1989; Myers et al. 1989; Niagro and Mishra 1989; Court et al. 1991; Alves and Videira 1998).While trying to establish the role of MUS-10 in DNA repair, we discovered that the mus-10 mutant exhibited a shortened life span, an abnormal mitochondrial morphology and mtDNA instability. We cloned the mus-10 gene through its ability to complement the MMS sensitivity of the mus-10 mutant and revealed that it encoded a novel F-box protein. This suggested that MUS-10 is part of an Skp, Cullin, F-box containing (SCF) E3 ubiquitin ligase complex that targets proteins for degradation by the 26S proteasome. The data we present in this article offer proof that an SCF complex can regulate both mitochondrial maintenance and cellular senescence.     

Binary Trait Mapping in Experimental Crosses With Selective Genotyping

Selective genotyping is an efficient strategy for mapping quantitative trait loci. For binary traits, where there are only two distinct phenotypic values (e.g., affected/unaffected or present/absent), one may consider selective genotyping of affected individuals, while genotyping none or only some of the unaffecteds. If selective genotyping of this sort is employed, the usual method for binary trait mapping, which considers phenotypes conditional on genotypes, cannot be used. We present an alternative approach, instead considering genotypes conditional on phenotypes, and compare this to the more standard method of analysis, both analytically and by example. For studies of rare binary phenotypes, we recommend performing an initial genome scan with all affected individuals and an equal number of unaffecteds, followed by genotyping the full cross in genomic regions of interest to confirm results from the initial screen.WE consider the problem of mapping genetic loci contributing to a binary trait in an experimental cross with selective genotyping. There are two clear approaches for linkage analysis with a binary trait. Typically, we compare the proportion of affected individuals across genotype groups (Xu and Atchley 1996). Alternatively, we can compare genotype frequencies between affected and unaffected individuals, similar to Henshall and Goddard (1999). Beyond these two basic approaches, binary trait mapping has seen fundamental advances in regression models (McIntyre et al. 2001; Deng et al. 2006), extensions to multiple-QTL mapping (Coffman et al. 2005; Chen and Liu 2009), and the development of Bayesian algorithms (Yi and Xu 2000; Huang et al. 2007). However, the original data structure and approach have remained intact. Existing methods for binary trait mapping largely require the availability of genotype and phenotype data for a representative sample of both affected and unaffected individuals, and we have not yet seen a well-developed framework for binary trait mapping in the presence of selective genotyping.It is not uncommon to see genotype data on affected individuals only, in which case the above methods cannot be used. Instead, we can compare observed genotype frequencies to the expected segregation ratios given the cross type, in a test for segregation distortion (see Faris et al. 1998; Lambrides et al. 2004). For example, the expected segregation proportions for an intercross are 1:2:1. The observed genotypes can then be described by a multinomial model, and statistically significant deviation from the expected segregation ratios among the genotyped affected individuals would suggest genotype–phenotype association. Gene mapping approaches that model genotypes rather than phenotypes have been developed extensively in the analysis of affected human relative pairs (see, for example, Risch 1990; Holmans 1993; Hauser and Boehnke 1998). In the analysis of experimental crosses, however, this type of approach has been developed primarily for the identification of monogenic mutants (Moran et al. 2006).Once all affected individuals are genotyped, an investigator may go on to genotype unaffected individuals. With this genotyping strategy in mind, we present several potential methods of analysis that might be applied in this context. First, we consider a standard analysis of the genotyped individuals, with disease proportions compared across genotype groups (Xu and Atchley 1996). Having omitted ungenotyped individuals, this method of analysis appears invalid because the estimated disease proportions are biased upward, reflecting an overrepresentation of affecteds in the set of genotyped individuals under consideration. As an alternative, we develop a reverse approach with genotype frequencies compared across phenotype groups. Because selective genotyping does provide a representative sample of genotypes for each phenotype group, this reverse approach does not face the bias in parameter estimation seen with the standard approach. We further extend the reverse approach to incorporate a segregation assumption, as is necessary for an affecteds only analysis. Finally, we present a full-likelihood analysis accounting for selective genotyping, similar to that suggested by Lander and Botstein (1989) for quantitative traits. We develop the full-likelihood approach both with and without incorporating an assumption on the genotype segregation proportions.Having put forth each of these methods, we derive analytic relationships among them. These relationships provide important insight regarding application of the presented methods under selective genotyping. Most notably, we find that making a segregation assumption can lead to spurious evidence of a QTL, but is necessary to treat the case of affecteds only genotyping. We demonstrate properties of the methods in an analysis of recovery from infection by Listeria monocytogenes in intercross mice and further compare power of the methods through computer simulations. Finally, we synthesize our analytical and simulation results to offer more general suggestions for the analysis of binary trait data with selective genotyping.     

Mating-System Variation,Demographic History and Patterns of Nucleotide Diversity in the Tristylous Plant Eichhornia paniculata

Inbreeding in highly selfing populations reduces effective size and, combined with demographic conditions associated with selfing, this can erode genetic diversity and increase population differentiation. Here we investigate the role that variation in mating patterns and demographic history play in shaping the distribution of nucleotide variation within and among populations of the annual neotropical colonizing plant Eichhornia paniculata, a species with wide variation in selfing rates. We sequenced 10 EST-derived nuclear loci in 225 individuals from 25 populations sampled from much of the geographic range and used coalescent simulations to investigate demographic history. Highly selfing populations exhibited moderate reductions in diversity but there was no significant difference in variation between outcrossing and mixed mating populations. Population size interacted strongly with mating system and explained more of the variation in diversity within populations. Bayesian structure analysis revealed strong regional clustering and selfing populations were highly differentiated on the basis of an analysis of F_st. There was no evidence for a significant loss of within-locus linkage disequilibrium within populations, but regional samples revealed greater breakdown in Brazil than in selfing populations from the Caribbean. Coalescent simulations indicate a moderate bottleneck associated with colonization of the Caribbean from Brazil ∼125,000 years before the present. Our results suggest that the recent multiple origins of selfing in E. paniculata from diverse outcrossing populations result in higher diversity than expected under long-term equilibrium.THE rate of self-fertilization in hermaphrodite organisms is expected to affect a number of important features of population genetic structure and diversity. Most directly, homozygosity increases as a function of the selfing rate and thus reduces the effective population size (N_e), up to twofold with complete selfing (Pollak 1987; Charlesworth et al. 1993; Nordborg 2000). Further, because of increased homozygosity, crossing over rarely occurs between heterozygous sites, thus increasing linkage disequilibrium (LD). Higher LD causes stronger hitchhiking effects such as selective sweeps, background selection, and Hill–Robertson interference, all of which are expected to further reduce the amount of neutral genetic variation within populations (reviewed in Charlesworth and Wright 2001).Population genetic processes resulting from inbreeding may be further augmented by demographic and life-history characteristics associated with the selfing habit. In particular, selfing populations can be founded by single individuals, resulting in striking reductions in diversity as a result of genetic bottlenecks and reproductive isolation. The capacity for uniparental reproduction gives many selfers prolific colonizing ability and the capacity to establish after long-distance dispersal, especially in comparison with obligate outcrossers (Baker 1955; Pannell and Barrett 1998). The colonization–extinction dynamics typical of many selfing species and limited pollen-mediated gene flow also increase differentiation among populations, resulting in considerable population subdivision (Hamrick and Godt 1990, 1996; Schoen and Brown 1991). Although the total amounts of among-population variation may be less affected by these processes (Pannell and Charlesworth 1999; Ingvarsson 2002), the demographic and life-history characteristics of many selfing species are likely to result in nonequilibrium conditions occurring in selfing populations.In many taxa where selfing has evolved it may be of relatively recent origin (Schoen et al. 1997; Takebayashi and Morrell 2001; Foxe et al. 2009; Guo et al. 2009). Where selfing has recently established, demographic forces associated with colonization may be as important as the mating system per se in structuring patterns of diversity. For example, if selfing originates through the establishment of a small number of founders, we would expect a sharp reduction in diversity relative to the outcrossing progenitor and a strong signature of a genetic bottleneck. In contrast, if selfing has evolved recently through the spread of genetic modifiers of small effect, newly established populations may retain significant amounts of ancestral polymorphism from their outcrossing progenitors. In this latter case populations may retain considerably more variation than expected under long-term equilibrium predictions.Molecular evidence for reduced nucleotide diversity and greater differentiation among populations of selfing taxa compared to populations of related outcrossing taxa has been reported from Leavenworthia (Liu et al. 1998, 1999), Arabidopsis (Savolainen et al. 2000; Wright et al. 2002), Solanum (Baudry et al. 2001), Mimulus (Sweigart and Willis 2003), Amsinckia (Perusse and Schoen 2004), and Caenorhabditis (Graustein et al. 2002; Cutter et al. 2006; Cutter 2008). In each case the reduction in diversity was more severe than the twofold reduction predicted for selfing populations at equilibrium. This indicates that factors in addition to the mating system are reducing diversity, but it has been difficult to uncouple the relative importance of genetic hitchhiking from the ecology and demographic history of selfing taxa. This challenge parallels similar difficulties in efforts to distinguish selective from demographic explanations in population genetic studies of Drosophila (Haddrill et al. 2005; Ometto et al. 2005; Thornton and Andolfatto 2006; Jensen et al. 2008). However, in many plant populations, especially those with annual life histories and small structured populations, demographic processes may play a more prominent role in causing reduced diversity than increased hitchhiking associated with selfing.Molecular population genetic studies of selfing in plants have generally focused on either small samples from a large number of populations (e.g., Sweigart and Willis 2003; Nordborg et al. 2005) or relatively large within-population samples from a small number of populations (e.g., Baudry et al. 2001). Ideally, a deeper sampling both within and among populations combined with independent ecological and historical information is required to improve understanding of the interplay of demographic and selective factors. Here we address these issues by examining patterns of nucleotide diversity within a large sample of populations of Eichhornia paniculata (Pontederiaceae), an annual species for which there is considerable ecological and demographic information (reviewed in Barrett and Husband 1997).E. paniculata occurs primarily in northeastern (N.E.) Brazil and the Caribbean islands of Cuba and Jamaica. Various lines of evidence suggest that Brazil is the original source region for Caribbean populations (reviewed in Barrett et al. 2009). Populations of E. paniculata exhibit striking mating-system diversity, ranging from predominantly outcrossing to those that are highly selfing (outcrossing rate, t = 0.002–0.96; n = 54 populations) (Barrett and Husband 1990; Barrett et al. 1992). Variation in mating system is associated with the evolutionary breakdown of the species'' tristylous genetic polymorphism and the spread and fixation of selfing variants capable of autonomous self-pollination (Barrett et al. 1989). Populations of E. paniculata are characterized by three morph structures: trimorphic with long-, mid-, and short-styled morphs (hereafter L-, M-, and S-morphs); dimorphic, with two floral morphs, most commonly the L- and M-morphs; and monomorphic, primarily composed of selfing variants of the M-morph. The morph structure and presence of selfing variants within populations explain ∼60% of the variation in outcrossing rates among populations (Barrett and Husband 1990). Trimorphic populations are largely outcrossing, dimorphic populations display mixed mating, and monomorphic populations are highly selfing. Patterns of allozyme variation indicate a reduction in diversity with increased selfing rates and greater among-population differentiation (Glover and Barrett 1987; Barrett and Husband 1990; Husband and Barrett 1993). Finally, studies of the inheritance of mating-system modifiers (Fenster and Barrett 1994; Vallejo-Marín and Barrett 2009) in combination with allozyme (Husband and Barrett 1993) and molecular evidence (Barrett et al. 2009) indicate that the transition from outcrossing to selfing in E. paniculata has occurred on multiple occasions.The goal of our study was to investigate the relation between mating-system variation and neutral molecular diversity for a large sample of E. paniculata populations encompassing most of the geographical range. This was accomplished by collecting multilocus nucleotide sequence data from 225 individuals sampled from 25 populations including trimorphic, dimorphic, and monomorphic populations. Because it has been previously demonstrated that this sequence of morph structures is strongly associated with increasing rates of self-fertilization (see Barrett and Husband 1990), we predicted a decrease in neutral diversity and increases in F_st and linkage disequilibrium from floral trimorphism to monomorphism. This extensive population-level sampling across a wide range of selfing rates allowed us to investigate the relative importance of mating system, geography, and current population size in structuring genetic variation. We also applied the approaches of Bayesian clustering (Pritchard et al. 2000; Falush et al. 2003; Gao et al. 2007) and divergence population genetics (Wakeley and Hey 1997; Hey and Nielsen 2004; Becquet and Przeworski 2007) to investigate the demographic history of E. paniculata and to provide a framework for understanding island colonization and the transition from outcrossing to selfing.