The Genetic Signature of Conditional Expression

Conditionally expressed genes have the property that every individual in a population carries and transmits the gene, but only a fraction, φ, expresses the gene and exposes it to natural selection. We show that a consequence of this pattern of inheritance and expression is a weakening of the strength of natural selection, allowing deleterious mutations to accumulate within and between species and inhibiting the spread of beneficial mutations. We extend previous theory to show that conditional expression in space and time have approximately equivalent effects on relaxing the strength of selection and that the effect holds in a spatially heterogeneous environment even with low migration rates among patches. We support our analytical approximations with computer simulations and delineate the parameter range under which the approximations fail. We model the effects of conditional expression on sequence polymorphism at mutation–selection–drift equilibrium, allowing for neutral sites, and show that sequence variation within and between species is inflated by conditional expression, with the effect being strongest in populations with large effective size. As φ decreases, more sites are recruited into neutrality, leading to pseudogenization and increased drift load. Mutation accumulation diminishes the degree of adaptation of conditionally expressed genes to rare environments, and the mutational cost of phenotypic plasticity, which we quantify as the plasticity load, is greater for more rarely expressed genes. Our theory connects gene-level relative polymorphism and divergence with the spatial and temporal frequency of environments inducing gene expression. Our theory suggests that null hypotheses for levels of standing genetic variation and sequence divergence must be corrected to account for the frequency of expression of the genes under study.IN genetically and ecologically subdivided populations, some individuals will experience a local environment very different from others, making it difficult to evolve a single adaptation adequate for all local conditions. Phenotypic plasticity allows organisms to respond adaptively to spatially and temporally varying environments by developing alternative phenotypes that enhance fitness under local conditions (Scheiner 1993; Via et al. 1995). Examples of alternative phenotypes, i.e., polyphenisms, include the defensive morphologies in Daphnia and algae induced by the presence of predators (e.g., Lively 1986; DeWitt 1998; Harvell 1998; Hazel et al. 2004); the winged and wingless morphs of bean beetles responding to resource variation (e.g., Abouheif and Wray 2002; Roff and Gelinas 2003; Lommen et al. 2005); and bacterial genes involved in traits such as quorum sensing, antibiotic production, biofilm formation, and virulence (Fuqua et al. 1996). The developmental basis of such alternative phenotypes often lies in the inducible expression of some genes in some individuals by environmental variables. That is, all individuals carry and transmit the conditionally expressed genes but only a fraction of individuals, φ, express them when environmental conditions are appropriate.The genes underlying plastic traits should experience relaxed selection due to conditional expression. Wade and co-workers have shown that genes hidden from natural selection in a fraction of individuals in the population by X-linked (Whitlock and Wade 1995; Linksvayer and Wade 2009) or sex-limited expression (Wade 1998; Demuth and Wade 2007) experience relaxed selective constraint. In Drosophila spp., sequence data for genes with maternally limited expression quantitatively support the theoretical predictions both for within-species polymorphism (Barker et al. 2005; Cruickshank and Wade 2008) and for between-species divergence (Barker Et Al 2005; Demuth and Wade 2007; Cruickshank and Wade 2008). Furthermore, male-specific genes in the facultatively sexual pea aphid have been shown to have elevated levels of sequence variation due to relaxed selection (Brisson and Nuzhdin 2008). Genes with spatially restricted expression in a heterogeneous environment should likewise experience relaxed selection. Adaptation to the most common environment in an ecologically subdivided population (Rosenzweig 1987; Holt and Gaines 1992; Holt 1996) allows deleterious mutations to accumulate in traits expressed in rare environments (Kawecki 1994; Whitlock 1996).Here we extend these results by quantifying the consequences of relaxed selection on conditionally expressed genes. Specifically, we show that, with weak selection, spatial and temporal fluctuations in selection intensity generate approximately equivalent effects on mean trait fitness, even with low rates of migration between habitats, resulting in a great simplification of analytical results. Our analytical approximations are supported with deterministic and stochastic simulations, and we note the conditions under which the approximations fail. We then derive general expressions for (1) the expected level of sequence polymorphism within populations under mutation, migration, drift, and purifying selection with conditional gene expression; (2) the rate of sequence divergence among populations, for dominant and recessive mutations; and (3) the reduction in mean population fitness due to accumulation of deleterious mutations at conditionally expressed loci. We find that the rate of accumulation of deleterious mutations for conditionally expressed genes is accelerated and the probability of fixation of beneficial mutations is reduced, causing a reduction in the fitness of conditional traits and an inflation in sequence variation within and between species. Our results suggest that evolutionary null hypotheses must be adjusted to account for the frequency of expression of genes under study, such that signatures of elevated within- or between-species sequence variation are not necessarily evidence of the action of diversifying natural selection. Furthermore, if conditional expression is due to spatial heterogeneity, we show that the level of genetic variation in a sample will often depend on whether or not genotypes were sampled from the selective habitat, the neutral habitat, or both. In the discussion we address the scope and limitations of our theory, as well as its implications for the maintenance of genetic variation, adaptive divergence between species, constraints on phenotypic plasticity, and evolutionary inference from sequence data.     

Inferring Bacterial Genome Flux While Considering Truncated Genes

Bacterial gene content variation during the course of evolution has been widely acknowledged and its pattern has been actively modeled in recent years. Gene truncation or gene pseudogenization also plays an important role in shaping bacterial genome content. Truncated genes could also arise from small-scale lateral gene transfer events. Unfortunately, the information of truncated genes has not been considered in any existing mathematical models on gene content variation. In this study, we developed a model to incorporate truncated genes. Maximum-likelihood estimates (MLEs) of the new model reveal fast rates of gene insertions/deletions on recent branches, suggesting a fast turnover of many recently transferred genes. The estimates also suggest that many truncated genes are in the process of being eliminated from the genome. Furthermore, we demonstrate that the ignorance of truncated genes in the estimation does not lead to a systematic bias but rather has a more complicated effect. Analysis using the new model not only provides more accurate estimates on gene gains/losses (or insertions/deletions), but also reduces any concern of a systematic bias from applying simplified models to bacterial genome evolution. Although not a primary purpose, the model incorporating truncated genes could be potentially used for phylogeny reconstruction using gene family content.GENE content variation as a key feature of bacterial genome evolution has been well recognized (Garcia-Vallvé et al. 2000; Ochman and Jones 2000; Snel et al. 2002; Welch et al. 2002; Kunin and Ouzounis 2003; Fraser-Liggett 2005; Tettelin et al. 2005) and gained increasing attention in recent years. Various methods have been employed to study the variation of gene content in the form of gene insertions/deletions (or gene gains/losses); there are studies of population dynamics (Nielsen and Townsend 2004), birth-and-death evolutionary models (Berg and Kurland 2002; Novozhilov et al. 2005), phylogeny-dependent studies including parsimony methods (Mirkin et al. 2003; Daubin et al. 2003a,b; Hao and Golding 2004), and maximum-likelihood methods (Hao and Golding 2006, 2008b; Cohen et al. 2008; Cohen and Pupko 2010; Spencer and Sangaralingam 2009). The pattern of gene presence/absence also contains phylogenetic signals (Fitz-Gibbon and House 1999; Snel et al. 1999; Tekaia et al. 1999) and has been used for phylogenetic reconstruction (Dutilh et al. 2004; Gu and Zhang 2004; Huson and Steel 2004; Zhang and Gu 2004; Spencer et al. 2007a,b). All these studies make use of the binary information of gene presence or absence and neglect the existence of gene segments or truncated genes.Bacterial genomes are known to harbor pseudogenes. An intracellular species Mycobacterium leprae is an extreme case for both the proportion and the number of pseudogenes: estimated as 40% of the 3.2-Mb genome and 1116 genes (Cole et al. 2001). In free-living bacteria, pseudogenes can make up to 8% of the annotated genes in the genome (Lerat and Ochman 2004). Many pseudogenes result from the degradation of native functional genes (Cole et al. 2001; Mira et al. 2001). Pseudogenes could also result from the degradation of transferred genes and might even be acquired directly via lateral gene transfer. For instance, in plant mitochondrial genomes, which have an α-proteobacterial ancestry, most, if not all, of the laterally transferred genes are pseudogenes (Richardson and Palmer 2007). Furthermore, evidence has been documented that gene transfer could take place at the subgenic level in a wide range of organisms, e.g., among bacteria (Miller et al. 2005; Choi and Kim 2007; Chan et al. 2009), between ancient duplicates in archaea (Archibald and Roger 2002), between different organelles (Hao and Palmer 2009; Hao 2010), and between eukaryotes (Keeling and Palmer 2001). A large fraction of pseudogenes have been shown to arise from failed lateral transfer events (Liu et al. 2004) and most of them are transient in bacterial genomes (Lerat and Ochman 2005). Zhaxybayeva et al. (2007) reported that genomes with truncated homologs might erroneously lead to false inferences of “gene gain” rather than multiple instances of “gene loss.” This raises the question of how a false diagnosis of gene absence affects the estimation of insertion/deletion rates. Recently, we showed that the effect of a false diagnosis of gene absence on estimation of insertion/deletion rates is not systematic, but rather more complicated (Hao and Golding 2008a). To further address the problem, a study incorporating the information of truncated genes is highly desirable. This will not only yield more accurate estimates of the rates of gene insertions/deletions, but also provide a quantitative view of the effect of truncated genes on rate estimation, which has been understudied in bacterial genome evolution.In this study, we developed a model that considers the information of truncated genes and makes use of a parameter-rich time-reversible rate matrix. Rate variation among genes is allowed in the model by incorporating a discrete Γ-distribution. We also allow rates to vary on different parts of the phylogeny (external branches vs. internal branches). Consistent with previous studies, the rates of gene insertions/deletions are comparable to or larger than the rates of nucleotide substitution and the rates of gene insertions/deletions are further inflated in closely related groups and on external branches, suggesting high rates of gene turnover of recently transferred genes. The results from the new model also suggest that many recently truncated genes are in the process of being rapidly deleted from the genome. Some other interesting estimates in the model are also presented and discussed. One implication of the study, though not primary, is that the state of truncated genes could serve as an additional phylogenetic character for phylogenetic reconstruction using gene family content.     

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.     

Utilization of d-Ribitol by Lactobacillus casei BL23 Requires a Mannose-Type Phosphotransferase System and Three Catabolic Enzymes

Lactobacillus casei strains 64H and BL23, but not ATCC 334, are able to ferment d-ribitol (also called d-adonitol). However, a BL23-derived ptsI mutant lacking enzyme I of the phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) was not able to utilize this pentitol, suggesting that strain BL23 transports and phosphorylates d-ribitol via a PTS. We identified an 11-kb region in the genome sequence of L. casei strain BL23 (LCABL_29160 to LCABL_29270) which is absent from strain ATCC 334 and which contains the genes for a GlpR/IolR-like repressor, the four components of a mannose-type PTS, and six metabolic enzymes potentially involved in d-ribitol metabolism. Deletion of the gene encoding the EIIB component of the presumed ribitol PTS indeed prevented d-ribitol fermentation. In addition, we overexpressed the six catabolic genes, purified the encoded enzymes, and determined the activities of four of them. They encode a d-ribitol-5-phosphate (d-ribitol-5-P) 2-dehydrogenase, a d-ribulose-5-P 3-epimerase, a d-ribose-5-P isomerase, and a d-xylulose-5-P phosphoketolase. In the first catabolic step, the protein d-ribitol-5-P 2-dehydrogenase uses NAD⁺ to oxidize d-ribitol-5-P formed during PTS-catalyzed transport to d-ribulose-5-P, which, in turn, is converted to d-xylulose-5-P by the enzyme d-ribulose-5-P 3-epimerase. Finally, the resulting d-xylulose-5-P is split by d-xylulose-5-P phosphoketolase in an inorganic phosphate-requiring reaction into acetylphosphate and the glycolytic intermediate d-glyceraldehyde-3-P. The three remaining enzymes, one of which was identified as d-ribose-5-P-isomerase, probably catalyze an alternative ribitol degradation pathway, which might be functional in L. casei strain 64H but not in BL23, because one of the BL23 genes carries a frameshift mutation.     

Population Structure With Localized Haplotype Clusters

We propose a multilocus version of F_ST and a measure of haplotype diversity using localized haplotype clusters. Specifically, we use haplotype clusters identified with BEAGLE, which is a program implementing a hidden Markov model for localized haplotype clustering and performing several functions including inference of haplotype phase. We apply this methodology to HapMap phase 3 data. With this haplotype-cluster approach, African populations have highest diversity and lowest divergence from the ancestral population, East Asian populations have lowest diversity and highest divergence, and other populations (European, Indian, and Mexican) have intermediate levels of diversity and divergence. These relationships accord with expectation based on other studies and accepted models of human history. In contrast, the population-specific F_ST estimates obtained directly from single-nucleotide polymorphisms (SNPs) do not reflect such expected relationships. We show that ascertainment bias of SNPs has less impact on the proposed haplotype-cluster-based F_ST than on the SNP-based version, which provides a potential explanation for these results. Thus, these new measures of F_ST and haplotype-cluster diversity provide an important new tool for population genetic analysis of high-density SNP data.GENOME-WIDE data sets from worldwide panels of individuals provide an outstanding opportunity to investigate the genetic structure of human populations (Conrad et al. 2006; International Hapmap Consortium 2007; Jakobsson et al. 2008; Auton et al. 2009). Populations around the globe form a continuum rather than discrete units (Serre and Paabo 2004; Weiss and Long 2009). However, notions of discrete populations can be appropriate when, for example, ancestral populations were separated by geographic distance or barriers such that little gene flow occurred.F_ST (Wright 1951; Weir and Cockerham 1984; Holsinger and Weir 2009) is a measure of population divergence. It measures variation between populations vs. within populations. One can calculate a global measure, assuming that all populations are equally diverged from an ancestral population, or one can calculate F_ST for specific populations or for pairs of populations while utilizing data from all populations (Weir and Hill 2002). One use of F_ST is to test for signatures of selection (reviewed in Oleksyk et al. 2010).F_ST may be calculated for single genetic markers. For multiallelic markers, such as microsatellites, this is useful, but single-nucleotide polymorphisms (SNPs) contain much less information when taken one at a time, and thus it is advantageous to calculate averages over windows of markers (Weir et al. 2005) or even over the whole genome. The advantage of windowed F_ST is that it can be used to find regions of the genome that show different patterns of divergence, indicative of selective forces at work during human history.Another measure of human evolutionary history is haplotype diversity. Haplotype diversity may be measured using a count of the number of observed haplotypes in a region or by the expected haplotype heterozygosity based on haplotype frequencies in a region. Application of this regional measure to chromosomal data can be achieved by a haplotype block strategy (Patil et al. 2001) or by windowing (Conrad et al. 2006; Auton et al. 2009).One problem with the analysis of population structure based on genome-wide panels of SNPs is that a large proportion of the SNPs were ascertained in Caucasians, potentially biasing the results of the analyses. Analysis based on haplotypes is less susceptible to such bias (Conrad et al. 2006). This is because haplotypes can be represented by multiple patterns of SNPs; thus lack of ascertainment of a particular SNP does not usually prevent observation of the haplotype. On a chromosome-wide scale, one cannot directly use entire haplotypes, because all the haplotypes in the sample will almost certainly be unique, thus providing no information on population structure. Instead one can use haplotypes on a local basis, either by using windows of adjacent markers or by using localized haplotype clusters, for example those obtained from fastPHASE (Scheet and Stephens 2006) or BEAGLE (Browning 2006; Browning and Browning 2007a).Localized haplotype clusters are a clustering of haplotypes on a localized basis. At the position of each genetic marker, haplotypes are clustered according to their similarity in the vicinity of the position. Both fastPHASE and BEAGLE use hidden Markov modeling to perform the clustering, although the specific models used by the two programs differ.Localized haplotype clusters derived from fastPHASE have been used to investigate haplotype diversity, to create neighbor-joining trees of populations, and to create multidimensional scaling (MDS) plots (Jakobsson et al. 2008). It was found that haplotype clusters showed different patterns of diversity to SNPs, while the neighbor-joining and MDS plots were similar between haplotype clusters and SNPs.In this work, we apply windowed F_ST methods to localized haplotype clusters derived from the BEAGLE program (Browning and Browning 2007a,b, 2009). We consider population-average, population-specific, and pairwise F_ST estimates (Weir and Hill 2002). Population-average F_ST''s either assume that all the populations are equally diverged from a common ancestor, which is not realistic, or represent the average of a set of population-specific values. This can be convenient in that the results are summarized by a single statistic; however, information is lost. A common procedure is to calculate F_ST for each pair of populations, and these values reflect the degree of divergence between the two populations. Different levels of divergence are allowed for each pair of populations but each estimate uses data from only that pair of populations. On the other hand, population-specific F_ST''s allow unequal levels of divergence in a single analysis that makes use of all the data.We compare results from the localized haplotype clusters to those using SNPs directly. The results of applying localized haplotype clusters to population-specific F_ST estimation are very striking, showing better separation of populations and a more realistic pattern of divergence than for population-specific F_ST estimation using SNPs directly. We also use BEAGLE''s haplotype clusters in a haplotype diversity measure and investigate the relationship between this measure of haplotype-cluster diversity and the recombination rate.     

Virulence Gene Regulation by l-Arabinose in Salmonella enterica

Invasion of the intestinal epithelium is a critical step in Salmonella enterica infection and requires functions encoded in the gene cluster known as Salmonella Pathogenicity Island 1 (SPI-1). Expression of SPI-1 genes is repressed by l-arabinose, and not by other pentoses. Transport of l-arabinose is necessary to repress SPI-1; however, repression is independent of l-arabinose metabolism and of the l-arabinose-responsive regulator AraC. SPI-1 repression by l-arabinose is exerted at a single target, HilD, and the mechanism appears to be post-translational. As a consequence of SPI-1 repression, l-arabinose reduces translocation of SPI-1 effectors to epithelial cells and decreases Salmonella invasion in vitro. These observations reveal a hitherto unknown role of l-arabinose in gene expression control and raise the possibility that Salmonella may use L-arabinose as an environmental signal.     

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.     

Accumulation of d-Glucose from Pentoses by Metabolically Engineered Escherichia coli

Escherichia coli that is unable to metabolize d-glucose (with knockouts in ptsG, manZ, and glk) accumulates a small amount of d-glucose (yield of about 0.01 g/g) during growth on the pentoses d-xylose or l-arabinose as a sole carbon source. Additional knockouts in the zwf and pfkA genes, encoding, respectively, d-glucose-6-phosphate 1-dehydrogenase and 6-phosphofructokinase I (E. coli MEC143), increased accumulation to greater than 1 g/liter d-glucose and 100 mg/liter d-mannose from 5 g/liter d-xylose or l-arabinose. Knockouts of other genes associated with interconversions of d-glucose-phosphates demonstrate that d-glucose is formed primarily by the dephosphorylation of d-glucose-6-phosphate. Under controlled batch conditions with 20 g/liter d-xylose, MEC143 generated 4.4 g/liter d-glucose and 0.6 g/liter d-mannose. The results establish a direct link between pentoses and hexoses and provide a novel strategy to increase carbon backbone length from five to six carbons by directing flux through the pentose phosphate pathway.     

Genome Changes After Gene Duplication: Haploidy vs. Diploidy

Since genome size and the number of duplicate genes observed in genomes increase from haploid to diploid organisms, diploidy might provide more evolutionary probabilities through gene duplication. It is still unclear how diploidy promotes genomic evolution in detail. In this study, we explored the evolution of segmental gene duplication in haploid and diploid populations by analytical and simulation approaches. Results show that (1) under the double null recessive (DNR) selective model, given the same recombination rate, the evolutionary trajectories and consequences are very similar between the same-size gene-pool haploid vs. diploid populations; (2) recombination enlarges the probability of preservation of duplicate genes in either haploid or diploid large populations, and haplo-insufficiency reinforces this effect; and (3) the loss of duplicate genes at the ancestor locus is limited under recombination while under complete linkage the loss of duplicate genes is always random at the ancestor and newly duplicated loci. Therefore, we propose a model to explain the advantage of diploidy: diploidy might facilitate the increase of recombination rate, especially under sexual reproduction; more duplicate genes are preserved under more recombination by originalization (by which duplicate genes are preserved intact at a special quasi-mutation-selection balance under the DNR or haplo-insufficient selective model), so genome sizes and the number of duplicate genes in diploid organisms become larger. Additionally, it is suggested that small genomic rearrangements due to the random loss of duplicate genes might be limited under recombination.USUALLY genome size becomes larger from haploid to diploid organisms (Lynch and Conery 2003), and so does the number of duplicate genes observed in genomes (Zhang 2003). It is extensively hypothesized that diploidy might facilitate the preservation and accumulation of duplicate genes, but it is still unclear how diploidy supports the evolution of duplicate genes in detail. The superiority of diploidy is classically attributed to preventing expression of deleterious mutations (Crow and Kimura 1965), but it is also argued that the sheltering of deleterious mutations cannot adequately explain the advantages of diploidy (Perrot et al. 1991).Recombination is a common phenomenon in all three kingdoms of life, Bacteria, Eukarya, and Archaea. It has been reported that recombination influences the loss of duplicate genes (Zhang and Kishino 2004; Xue et al. 2010). In diploid organisms, if recombination between the ancestor locus and the newly duplicated locus is free, the rate of recombination is maximally 0.5, which is commonly observed especially when the two loci are located on different chromosomes. Although recombination should not be regarded as an exception in haploid organisms (Fraser et al. 2007), recombination events usually occur more frequently in diploid populations than they do in haploid populations. In other words, diploidy might facilitate the occurrence of recombination. The difference of recombination behaviors between haploid and diploid organisms is an obvious and important feature during genomic evolution.In our recent studies of genomic duplication, we proposed a new possible way of preserving and accumulating duplicate genes in genomes—originalization (Xue and Fu 2009a). As is well known, for a locus in an infinite diploid population, the frequencies of wild-type and degenerative alleles will move to an equilibrium under purifying selection and mutation, which is known as the mutation–selection balance. After genomic duplication, under two simple selective models, double null recessive (DNR, under which valid individuals require at least one active wild-type allele on the ancestor and newly duplicated loci) and haplo-insufficient (HI or partial dominant, under which valid individuals require at least two active wild-type alleles on both loci) models, a special equilibrium of allele frequencies at the ancestor and newly duplicated loci will be reached under recombination, in which the frequency of wild-type allele is kept high at both loci. Under the HI selective model this balance becomes so stable and flexible that the fixation of a degenerative allele at one of these two loci (or the balance being broken) becomes very difficult even in a modest population (Xue and Fu 2009a,b). However, if the two loci are tightly linked (recombination rate r = 0), this balance of allele frequencies does not appear. As r increases, the balance becomes more stable and the frequency of the wild-type allele at two loci becomes higher. High frequency of the wild-type allele at both loci means that duplicate genes are preserved intact in genomes, so this phenomenon was named originalization.Although many duplicate genes originated from genomic duplications in some species, such as yeast, maize, and fish (Li et al. 2005), those from segmental duplications are also very popular (Zhang et al. 2000; Leister 2004). In haploid populations, most duplication events are small segmental duplications. Therefore, to understand genomic evolution comprehensively, it is necessary to explore the evolution of segmental genomic duplication.Lynch et al. (2001) and Tanaka et al. (2009) have studied the evolution of segmental gene duplication in diploid populations theoretically. However, in this study, we further compared the evolution of segmental gene duplication in haploid vs. diploid populations by numerical and simulation approaches under the DNR and HI selective models. We observed that haploid and diploid populations with the same-size gene pool are very similar under the DNR model and the same recombination rate. Recombination enlarges the probability of preservation of duplicate genes in either haploid or diploid populations via originalization, and haplo-insufficiency reinforces this effect. The loss of duplicate genes at the ancestor locus might be limited under recombination, while under complete linkage, the loss of duplicate genes is random at the ancestor and newly duplicated loci. According to these results, we propose a model with which to explain the revolutionary genomic transition from haploidy to diploidy.     

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.     

Identification of an l-Arabinose Reductase Gene in Aspergillus niger and Its Role in l-Arabinose Catabolism

The first enzyme in the pathway for l-arabinose catabolism in eukaryotic microorganisms is a reductase, reducing l-arabinose to l-arabitol. The enzymes catalyzing this reduction are in general nonspecific and would also reduce d-xylose to xylitol, the first step in eukaryotic d-xylose catabolism. It is not clear whether microorganisms use different enzymes depending on the carbon source. Here we show that Aspergillus niger makes use of two different enzymes. We identified, cloned, and characterized an l-arabinose reductase, larA, that is different from the d-xylose reductase, xyrA. The larA is up-regulated on l-arabinose, while the xyrA is up-regulated on d-xylose. There is however an initial up-regulation of larA also on d-xylose but that fades away after about 4 h. The deletion of the larA gene in A. niger results in a slow growth phenotype on l-arabinose, whereas the growth on d-xylose is unaffected. The l-arabinose reductase can convert l-arabinose and d-xylose to their corresponding sugar alcohols but has a higher affinity for l-arabinose. The K_m for l-arabinose is 54 ± 6 mm and for d-xylose 155 ± 15 mm.     

Association Mapping of Quantitative Disease Resistance in a Natural Population of Loblolly Pine (Pinus taeda L.)

Genetic resistance to disease incited by necrotrophic pathogens is not well understood in plants. Whereas resistance is often quantitative, there is limited information on the genes that underpin quantitative variation in disease resistance. We used a population genomic approach to identify genes in loblolly pine (Pinus taeda) that are associated with resistance to pitch canker, a disease incited by the necrotrophic pathogen Fusarium circinatum. A set of 498 largely unrelated, clonally propagated genotypes were inoculated with F. circinatum microconidia and lesion length, a measure of disease resistance, data were collected 4, 8, and 12 weeks after inoculation. Best linear unbiased prediction was used to adjust for imbalance in number of observations and to identify highly susceptible and highly resistant genotypes (“tails”). The tails were reinoculated to validate the results of the full population screen. Significant associations were detected in 10 single nucleotide polymorphisms (SNPs) (out of 3938 tested). As hypothesized for genes involved in quantitative resistance, the 10 SNPs had small effects and proposed roles in basal resistance, direct defense, and signal transduction. We also discovered associated genes with unknown function, which would have remained undetected in a candidate gene approach constrained by annotation for disease resistance or stress response.GENETIC interactions between host and pathogen populations result in abundant natural variation in the genes involved in host disease resistance. Most of the studies leading to identification and cloning of disease resistance genes are focused on major gene disease resistance (Johal and Briggs 1992; Dangl and Jones 2001; Jones and Dangl 2006). In cases where resistance is associated with single genes, genetic effects are large in magnitude and detection is straightforward. In contrast, quantitative disease resistance is typically conditioned by many genes with relatively small effects. Quantitative resistance is generally considered to be more durable but also more difficult to investigate relative to major gene resistance, since the effects of individual genes are small and phenotyping experiments must be performed with high levels of precision. As a consequence, the genes and mechanisms of quantitative disease resistance are poorly understood, in part due to the smaller effect of individual genes on the resistance phenotype. Interactions between plants and necrotrophic pathogens often exhibit quantitative resistance (Balint-Kurti et al. 2008; Poland et al. 2009).Pitch canker disease of loblolly pine and other pine species is incited by the necrotrophic pathogen Fusarium circinatum and is manifest as resinous lesions in stems and branches (Dwinell et al. 1985; Enebak and Stanosz 2003; Carey et al. 2005; Sakamoto and Gordon 2006). There is evidence for heritable resistance to pitch canker in loblolly pine (Kayihan et al. 2005) as well as other pine species (Hodge and Dvorak 2000, 2007). In this article we report the first population-wide phenotypic screen of a clonally propagated population of loblolly pine for association testing (Eckert et al. 2010). Clonal propagation of this population enabled precise phenotyping, which was required to obtain the resolution needed to identify candidates for quantitative disease resistance loci.Pine species in general exhibit high levels of nucleotide variation and low linkage disequilibrium (LD) (Brown et al. 2004). An association genetic approach relies on the premise that historical, unrecorded recombination events over many generations have reduced LD between markers and quantitative trait loci such that only those marker-trait pairs that are tightly linked remain detectable; this may enable “fine mapping” to identify genes underlying quantitative variation (Flint-Garcia et al. 2003; Neale and Savolainen 2004). Association-based approaches have been used to identify candidate genes underlying traits in plants (Zhao et al. 2007; Stich et al. 2008; Wang et al. 2008; Yahiaoui et al. 2008; Inostroza et al. 2009; Stracke et al. 2009), based in part on applications in humans (D''alfonso et al. 2002; McGuffin et al. 2003; Easton et al. 2007; Lee et al. 2007), livestock (Martinez et al. 2006; Charlier et al. 2008; Goddard and Hayes 2009), and Drosophila (Kennington et al. 2007; Norry et al. 2007; Jiang et al. 2009). Recent association studies in tree species have evaluated single candidate genes or a modest number of candidate genes for association (Thumma et al. 2005; Gonzalez-Martinez et al. 2007, 2008; Ingvarsson et al. 2008; Eckert et al. 2009a). Association mapping has been used to identify disease resistance genes in several crop species including sugarcane, maize, barley, and potato (Flint-Garcia et al. 2005; Wei et al. 2006; Yu and Buckler 2006; Malosetti et al. 2007; Stich et al. 2008; Inostroza et al. 2009; Murray et al. 2009). The population analyzed in this study was genotyped at 3938 SNP loci that were selected without regard to the functional annotation of ESTs from which they were derived. Thus, we reasoned that the status of any particular marker as a candidate disease resistance gene would be determined by association testing, as opposed to previous studies in which markers were typically evaluated on the basis of their presumed roles in disease resistance in other species.Several different, but not mutually exclusive hypotheses have been proposed regarding the genetic origins of quantitative resistance (Poland et al. 2009), providing a useful framework for understanding evolution of resistance to necrotrophic pathogens. These six hypotheses proposed by Poland et al. (2009) predict that quantitative disease resistance is conditioned by: (1) genes regulating morphological and developmental phenotypes; (2) mutations in genes involved in basal defense causing small, incremental levels of resistance; (3) components of chemical warfare, through the action of genes producing antibiotic or antifungal compounds; (4) genes involved in defense signal transduction pathways; (5) weak forms of defeated R genes; and/or (6) genes not yet known to be involved in disease resistance.In this study, our main objective was to evaluate the genetic architecture of pitch canker disease resistance: to quantify the extent to which genes contribute to variation in the disease phenotype, to evaluate the hypothesis that disease resistance was quantitative, and to identify candidate genes for resistance as well as quantify their magnitude of effect. In the process of identifying candidate genes for resistance we were also able to evaluate support for hypotheses recently put forth by Poland et al. (2009) regarding the biological roles and origins of quantitative resistance genes.     

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?     

Association Studies Identify Natural Variation at PHYC Linked to Flowering Time and Morphological Variation in Pearl Millet

The identification of genes selected during and after plant domestication is an important research topic to enhance knowledge on adaptative evolution. Adaptation to different climates was a key factor in the spread of domesticated crops. We conducted a study to identify genes responsible for these adaptations in pearl millet and developed an association framework to identify genetic variations associated with the phenotype in this species. A set of 90 inbred lines genotyped using microsatellite loci and AFLP markers was used. The population structure was assessed using two different Bayesian approaches that allow inbreeding or not. Association studies were performed using a linear mixed model considering both the population structure and familial relationships between inbred lines. We assessed the ability of the method to limit the number of false positive associations on the basis of the two different Bayesian methods, the number of populations considered and different morphological traits while also assessing the power of the methodology to detect given additive effects. Finally, we applied this methodology to a set of eight pearl millet genes homologous to cereal flowering pathway genes. We found significant associations between several polymorphisms of the pearl millet PHYC gene and flowering time, spike length, and stem diameter in the inbred line panel. To validate this association, we performed a second association analysis in a different set of pearl millet individuals from Niger. We confirmed a significant association between genetic variation in this gene and these characters.DOMESTICATION and dispersion of cultivated plants were associated with their adaptation to the agricultural environment. These adaptations led to genetic changes shared by all individuals of a cultivated species (domestication genes) or to variations between varieties within a cultivated species (genes controlling varietal differences). Domestication genes like tb1 (Doebley et al. 1997; Wang et al. 1999) in maize (Zea mays) were selected very early by human populations (Jaenicke-Després et al. 2003). After the first early selection, adaptation of the flowering phenotype to different climatic conditions was certainly a key innovation that enabled colonization of new environments. One of the most well-known examples was the adaptation of maize—a tropical plant—to northern climates. Maize cultivation spread late to northeastern America. By 1000 YBP, only maize was an established staple crop (Fritz 1995). A genetic variant of the Dwarf8 gene led to an earlier flowering phenotype (Thornsberry et al. 2001). This early allele was present at a high frequency in North America and was certainly selected after the domestication of maize under northern climatic conditions (Camus-Kulandaivelu et al. 2006).Pearl millet (Pennisetum glaucum [(L.) R. Br.]), one of the most important West African cereals, was most likely domesticated once in the Sahelian zone of West Africa (Oumar et al. 2008). By 3500 YBP, it was already being cultivated throughout Sahelian and tropical West African countries (D''Andrea et al. 2001; D''Andrea and Casey 2002). The adaptation of pearl millet in West Africa was also associated with an environmental gradient (Haussmann et al. 2006). Pearl millet varieties from tropical coastal West Africa flower very late (up to 160 days from planting to female flowering) as compared to varieties from Sahelian West Africa, which may have a flowering time as short as 45 days (Haussmann et al. 2006). The genetic factors underlying the differences between these varieties are still unknown.Association studies offer new opportunities for assessing the role of a particular gene on a phenotype. Contrary to QTL analysis, association studies have the challenging task of taking an unknown evolutionary history of studied individuals into account. For example, population structure is a common confounding effect in association studies (Pritchard et al. 2000a). Allele frequencies evolve between divergent structured populations via drift, mutation, and selection. Differences in allele frequencies may be correlated with any morphological traits that differentiate two populations. Then a statistical correlation between a gene and a trait is not necessarily associated with a “causative” relationship between the gene and the morphology, which can lead to a high number of false positives. The use of population structure to correct the number of false positives was a significant breakthrough in plant studies (Thornsberry et al. 2001). This approach was recently further refined by also using a matrix of kinship coefficients, which proves efficient when there is a complex structure and familial relationship between individuals (Yu et al. 2006; Kang et al. 2008; Stich et al. 2008). Complex structures and familial relationships are common in inbred cultivated crop material. In the current association study framework (Thornsberry et al. 2001; Yu et al. 2006; Casa et al. 2008; Kang et al. 2008; Stich et al. 2008), population structure was assessed using STRUCTURE software (Pritchard et al. 2000b). This tool is not implemented to deal with selfed inbred materials or inbred species (Pritchard et al. 2000b). Through new methodological developments, population structure analysis can now be performed using Bayesian methods in these particular cases (Gao et al. 2007). The extent to which the power of association studies will differ when dealing with inbred material or selfing species using either Bayesian method has yet to be evaluated.In this study, we developed an association framework for pearl millet to assess the role of flowering pathway genes. We assessed the ability of the method to control the number of false positives, while taking different methodological inferences of population structure that allow inbreeding or not into account. We also assessed the power of the association framework to detect given additive genetic effects. Finally, we applied this method to a set of eight flowering time gene homologs sequenced in pearl millet. We assessed sequence variation in light perception genes (PHYA, PHYB, PHYC, and CRY2) and downstream regulators of flowering (GI, Hd6, Hd1, and FLORICAULA). Variation was detected in the PHYC gene associated with variations in flowering time and morphological traits. This association was noted in two different data sets.     

Patterns and Processes of Genome-Wide Divergence Between North American and African Drosophila melanogaster

Genomic tools and analyses are now being widely used to understand genome-wide patterns and processes associated with speciation and adaptation. In this article, we apply a genomics approach to the model organism Drosophila melanogaster. This species originated in Africa and subsequently spread and adapted to temperate environments of Eurasia and the New World, leading some populations to evolve reproductive isolation, especially between cosmopolitan and Zimbabwean populations. We used tiling arrays to identify highly differentiated regions within and between North America (the United States and Caribbean) and Africa (Cameroon and Zimbabwe) across 63% of the D. melanogaster genome and then sequenced representative fragments to study their genetic divergence. Consistent with previous findings, our results showed that most differentiation was between populations living in Africa vs. outside of Africa (i.e., “out-of-Africa” divergence), with all other geographic differences being less substantial (e.g., between cosmopolitan and Zimbabwean races). The X chromosome was much more strongly differentiated than the autosomes between North American and African populations (i.e., greater X divergence). Overall differentiation was positively associated with recombination rates across chromosomes, with a sharp reduction in regions near centromeres. Fragments surrounding these high F_ST sites showed reduced haplotype diversity and increased frequency of rare and derived alleles in North American populations compared to African populations. Nevertheless, despite sharp deviation from neutrality in North American strains, a small set of bottleneck/expansion demographic models was consistent with patterns of variation at the majority of our high F_ST fragments. Although North American populations were more genetically variable compared to Europe, our simulation results were generally consistent with those previously based on European samples. These findings support the hypothesis that most differentiation between North America and Africa was likely driven by the sorting of African standing genetic variation into the New World via Europe. Finally, a few exceptional loci were identified, highlighting the need to use an appropriate demographic null model to identify possible cases of selective sweeps in species with complex demographic histories.THE study of genetic differentiation between populations and species has recently been empowered by the use of genomic techniques and analysis (e.g., Noor and Feder 2006; Stinchcombe and Hoekstra 2008). In the past decade, genetic studies of adaptation and speciation have taken advantage of emerging molecular techniques to scan the genomes of diverging populations for highly differentiated genetic regions (e.g., Wilding et al. 2001; Emelianov et al. 2003; Beaumont and Balding 2004; Campbell and Bernatchez 2004; Scotti-Saintagne et al. 2004; Achere et al. 2005; Turner et al. 2005; Vasemagi et al. 2005; Bonin et al. 2006, 2007; Murray and Hare 2006; Savolainen et al. 2006; Yatabe et al. 2007; Nosil et al. 2008, 2009; Turner et al. 2008a,b; Kulathinal et al. 2009). As a result, genome scans can identify candidate regions that may be associated with adaptive evolution between diverging populations and, more broadly, are able to describe genome-wide patterns and processes of population differentiation (Begun et al. 2007; Stinchcombe and Hoekstra 2008).Genome scans in well-studied genetic model species such as Drosophila melanogaster gain particular power because differentiated loci are mapped to a well-annotated genome. Moreover, the evolutionary history of D. melanogaster is rich with adaptive and demographic events with many parallels to human evolution. Most notable is the historical out-of-Africa migration and subsequent adaptation to temperate ecological environments of Europe, Asia, North America, and Australia. This has resulted in widespread genetic and phenotypic divergence between African and non-African populations (e.g., David and Capy 1988; Begun and Aquadro 1993; Capy et al. 1994; Colegrave et al. 2000; Rouault et al. 2001; Takahashi et al. 2001; Caracristi and Schlötterer 2003; Baudry et al. 2004; Pool and Aquadro 2006; Schmidt et al. 2008; Yukilevich and True 2008a,b). Further, certain populations in Africa and in the Caribbean vary in their degree of reproductive isolation from populations in more temperate regions (Wu et al. 1995; Hollocher et al. 1997; Yukilevich and True 2008a,b). In particular, the Zimbabwe and nearby populations of southern Africa are strongly sexually isolated from all other populations, designating them as a distinct behavioral race (Wu et al. 1995).D. melanogaster has received a great deal of attention from the population geneticists in studying patterns of sequence variation across African and non-African populations. Many snapshots have been taken of random microsatellite and SNP variants spread across X and autosomes, and these have generated several important conclusions. Polymorphism patterns in European populations are characterized by reduced levels of nucleotide and haplotype diversity, an excess of high frequency-derived polymorphisms, and elevated levels of linkage disequilibrium relative to African populations (e.g., Begun and Aquadro 1993; Andolfatto 2001; Glinka et al. 2003; Haddrill et al. 2005; Ometto et al. 2005; Thornton and Andolfatto 2006; Hutter et al. 2007; Singh et al. 2007). These results have been generally interpreted as compatible with population size reduction/bottlenecks followed by recent population expansions. On the other hand, African populations are generally assumed either to have been relatively constant in size over time or to have experienced population size expansions. They generally show higher levels of nucleotide and haplotype diversity, an excess of rare variants, and a deficit of high frequency-derived alleles (Glinka et al. 2003; Ometto et al. 2005; Pool and Aquadro 2006; Hutter et al. 2007; but see Haddrill et al. 2005 for evidence of bottlenecks in Africa).Previous work also shows that the ratio of X-linked to autosomal polymorphism deviates from neutral expectations in opposite directions in African and European populations with more variation on the X than expected in Africa and less variation on the X than expected in Europe (Andolfatto 2001; Kauer et al. 2002; Hutter et al. 2007; Singh et al. 2007). The deviation from neutrality in the ratio of X-autosome polymorphism may be explained by positive selection being more prevalent on the X in Europe and/or by a combination of bottlenecks and male-biased sex ratios in Europe and female-biased sex ratios in Africa (Charlesworth 2001; Hutter et al. 2007; Singh et al. 2007). The selective explanation stems from the argument that, under the hitchhiking selection model, X-linked loci are likely to be more affected by selective sweeps than autosomal loci (Maynard Smith and Haigh 1974; Charlesworth et al. 1987; Vicoso and Charlesworth 2006, 2009).The relative contribution of selective and demographic processes in shaping patterns of genomic variation and differentiation is highly debated (Wall et al. 2002; Glinka et al. 2003; Haddrill et al. 2005; Ometto et al. 2005; Schöfl and Schlötterer 2004; Thornton and Andolfatto 2006; Hutter et al. 2007; Singh et al. 2007; Shapiro et al. 2007; Stephan and Li 2007; Hahn 2008; Macpherson et al. 2008; Noor and Bennett 2009; Sella et al. 2009). This is especially the case in D. melanogaster because derived non-African populations have likely experienced a complex set of demographic events during their migration out of Africa (e.g., Thornton and Andolfatto 2006; Singh et al. 2007; Stephan and Li 2007), making population genetics signatures of demography and selection difficult to tease apart (e.g., Macpherson et al. 2008). Thus it is still unclear what role selection has played in shaping overall patterns of genomic variation and differentiation relative to demographic processes in this species.While there is a long tradition in studying arbitrarily or opportunistically chosen sequences in D. melanogaster, genomic scans that focus particularly on highly differentiated sites across the genome have received much less attention. Such sites are arguably the best candidates to resolve the debate on which processes have shaped genomic differentiation within species (e.g., Przeworski 2002). Recently, a genome-wide scan of cosmopolitan populations in the United States and in Australia was performed to investigate clinal genomic differentiation on the two continents (Turner et al. 2008a). Many single feature polymorphisms differentiating Northern and Southern Hemisphere populations were identified. Among the most differentiated loci in common between continents, 80% were differentiated in the same orientation relative to the Equator, implicating selection as the likely explanation (Turner et al. 2008a). Larger regions of genomic differentiation within and between African and non-African populations have also been discovered, some of them possibly being driven by divergent selection (e.g., Dopman and Hartl 2007; Emerson et al. 2008; Turner et al. 2008a, Aguade 2009). Despite this recent progress, we still know relatively little about large-scale patterns of genomic differentiation in this species, especially between African and non-African populations, and whether most of this differentiation is consistent with demographic processes alone or if it requires selective explanations.In this work, we explicitly focus on identifying differentiated sites across the genome between U.S., Caribbean, West African, and Zimbabwean populations. This allows us to address several fundamental questions related to genomic evolution in D. melanogaster, such as the following: (1) Do genome-wide patterns of differentiation reflect patterns of reproductive isolation? (2) Is genomic differentiation random across and within chromosomes or are some regions overrepresented? (3) What are the population genetics properties of differentiated sites and their surrounding sequences? (4) Can demographic historical processes alone explain most of the observed differentiation on a genome-wide level or is it necessary to involve selection in their explanation?In general, our findings revealed that most genomic differentiation within D. melanogaster shows an out-of-Africa genetic signature. These results are inconsistent with the notion that most genomic differentiation occurs between cosmopolitan and Zimbabwean reproductively isolated races. Further, we found that the X is more differentiated between North American and African populations and more strongly deviates from pure neutrality in North American populations relative to autosomes. Nevertheless, our article shows that much of this deviation from neutrality is broadly consistent with several demographic null models, with a few notable exceptions. Athough this does not exclude selection as a possible alternative mechanism for the observed patterns, it supports the idea that most differentiation in D. melanogaster was likely driven by the sorting of African standing genetic variation into the New World.     

High-Level Heterologous Production of d-Cycloserine by Escherichia coli

Previously, we successfully cloned a d-cycloserine (d-CS) biosynthetic gene cluster consisting of 10 open reading frames (designated dcsA to dcsJ) from d-CS-producing Streptomyces lavendulae ATCC 11924. In this study, we put four d-CS biosynthetic genes (dcsC, dcsD, dcsE, and dcsG) in tandem under the control of the T7 promoter in an Escherichia coli host. SDS-PAGE analysis demonstrated that the 4 gene products were simultaneously expressed in host cells. When l-serine and hydroxyurea (HU), the precursors of d-CS, were incubated together with the E. coli resting cell suspension, the cells produced significant amounts of d-CS (350 ± 20 μM). To increase the productivity of d-CS, the dcsJ gene, which might be responsible for the d-CS excretion, was connected downstream of the four genes. The E. coli resting cells harboring the five genes produced d-CS at 660 ± 31 μM. The dcsD gene product, DcsD, forms O-ureido-l-serine from O-acetyl-l-serine (OAS) and HU, which are intermediates in d-CS biosynthesis. DcsD also catalyzes the formation of l-cysteine from OAS and H₂S. To repress the side catalytic activity of DcsD, the E. coli chromosomal cysJ and cysK genes, encoding the sulfite reductase α subunit and OAS sulfhydrylase, respectively, were disrupted. When resting cells of the double-knockout mutant harboring the four d-CS biosynthetic genes, together with dcsJ, were incubated with l-serine and HU, the d-CS production was 980 ± 57 μM, which is comparable to that of d-CS-producing S. lavendulae ATCC 11924 (930 ± 36 μM).