On the inference of complex phylogenetic networks by Markov Chain Monte-Carlo

For various species, high quality sequences and complete genomes are nowadays available for many individuals. This makes data analysis challenging, as methods need not only to be accurate, but also time efficient given the tremendous amount of data to process. In this article, we introduce an efficient method to infer the evolutionary history of individuals under the multispecies coalescent model in networks (MSNC). Phylogenetic networks are an extension of phylogenetic trees that can contain reticulate nodes, which allow to model complex biological events such as horizontal gene transfer, hybridization and introgression. We present a novel way to compute the likelihood of biallelic markers sampled along genomes whose evolution involved such events. This likelihood computation is at the heart of a Bayesian network inference method called SnappNet, as it extends the Snapp method inferring evolutionary trees under the multispecies coalescent model, to networks. SnappNet is available as a package of the well-known beast 2 software.Recently, the MCMC_BiMarkers method, implemented in PhyloNet, also extended Snapp to networks. Both methods take biallelic markers as input, rely on the same model of evolution and sample networks in a Bayesian framework, though using different methods for computing priors. However, SnappNet relies on algorithms that are exponentially more time-efficient on non-trivial networks. Using simulations, we compare performances of SnappNet and MCMC_BiMarkers. We show that both methods enjoy similar abilities to recover simple networks, but SnappNet is more accurate than MCMC_BiMarkers on more complex network scenarios. Also, on complex networks, SnappNet is found to be extremely faster than MCMC_BiMarkers in terms of time required for the likelihood computation. We finally illustrate SnappNet performances on a rice data set. SnappNet infers a scenario that is consistent with previous results and provides additional understanding of rice evolution.     

Site-directed Mutagenesis Switching a Dimethylallyl Tryptophan Synthase to a Specific Tyrosine C3-Prenylating Enzyme

The tryptophan prenyltransferases FgaPT2 and 7-DMATS (7-dimethylallyl tryptophan synthase) from Aspergillus fumigatus catalyze C⁴- and C⁷-prenylation of the indole ring, respectively. 7-DMATS was found to accept l-tyrosine as substrate as well and converted it to an O-prenylated derivative. An acceptance of l-tyrosine by FgaPT2 was also observed in this study. Interestingly, isolation and structure elucidation revealed the identification of a C³-prenylated l-tyrosine as enzyme product. Molecular modeling and site-directed mutagenesis led to creation of a mutant FgaPT2_K174F, which showed much higher specificity toward l-tyrosine than l-tryptophan. Its catalytic efficiency toward l-tyrosine was found to be 4.9-fold in comparison with that of non-mutated FgaPT2, whereas the activity toward l-tryptophan was less than 0.4% of that of the wild-type. To the best of our knowledge, this is the first report on an enzymatic C-prenylation of l-tyrosine as free amino acid and altering the substrate preference of a prenyltransferase by mutagenesis.     

Biosynthesis of l-tyrosine O-sulphate from the methyl and ethyl esters of l-tyrosine

1. Rat-liver supernatant preparations are capable of achieving the biological sulphation of l-tyrosine methyl ester, the reaction proceeding maximally at a substrate concentration of 30 mm and at pH 7·0. 2. Two sulphated products are formed, one of which has been identified as l-tyrosine O-sulphate. On the basis of indirect evidence the other product can be assumed to be l-tyrosine O-sulphate methyl ester. 3. An enzyme present in rat-liver supernatant preparations is capable of converting l-tyrosine O-sulphate methyl ester into l-tyrosine O-sulphate. This enzyme is inhibited by l-tyrosine methyl ester. 4. l-Tyrosine ethyl ester also yields two sulphated products when used as an acceptor in the liver sulphating system. One of these has been identified chromatographically as l-tyrosine O-sulphate and the other may be presumed to be l-tyrosine O-sulphate ethyl ester.     

The biological sulphation of l-tyrosyl peptides

1. A rat-liver supernatant preparation can achieve the biological O-sulphation of l-tyrosylglycine and l-tyrosyl-l-alanine at pH7·0. 2. The optimum concentrations of l-tyrosylglycine and l-tyrosyl-l-alanine in this system are 50mm and 60mm respectively. 3. l-Tyrosylglycine yields two sulphated products, whereas l-tyrosyl-l-alanine yields three sulphated products, when used as acceptor for sulphate in the rat-liver system. 4. With both substrates, one of the sulphated products has been identified as the O-sulphate ester of the corresponding parent peptide.     

Specific Duplication and Dorsoventrally Asymmetric Expression Patterns of Cycloidea-Like Genes in Zygomorphic Species of Ranunculaceae

Floral bilateral symmetry (zygomorphy) has evolved several times independently in angiosperms from radially symmetrical (actinomorphic) ancestral states. Homologs of the Antirrhinum majus Cycloidea gene (Cyc) have been shown to control floral symmetry in diverse groups in core eudicots. In the basal eudicot family Ranunculaceae, there is a single evolutionary transition from actinomorphy to zygomorphy in the stem lineage of the tribe Delphinieae. We characterized Cyc homologs in 18 genera of Ranunculaceae, including the four genera of Delphinieae, in a sampling that represents the floral morphological diversity of this tribe, and reconstructed the evolutionary history of this gene family in Ranunculaceae. Within each of the two RanaCyL (Ranunculaceae Cycloidea-like) lineages previously identified, an additional duplication possibly predating the emergence of the Delphinieae was found, resulting in up to four gene copies in zygomorphic species. Expression analyses indicate that the RanaCyL paralogs are expressed early in floral buds and that the duration of their expression varies between species and paralog class. At most one RanaCyL paralog was expressed during the late stages of floral development in the actinomorphic species studied whereas all paralogs from the zygomorphic species were expressed, composing a species-specific identity code for perianth organs. The contrasted asymmetric patterns of expression observed in the two zygomorphic species is discussed in relation to their distinct perianth architecture.     

Enzymatic Biosynthesis of Novel Resveratrol Glucoside and Glycoside Derivatives

A UDP glucosyltransferase from Bacillus licheniformis was overexpressed, purified, and incubated with nucleotide diphosphate (NDP) d- and l-sugars to produce glucose, galactose, 2-deoxyglucose, viosamine, rhamnose, and fucose sugar-conjugated resveratrol glycosides. Significantly higher (90%) bioconversion of resveratrol was achieved with α-d-glucose as the sugar donor to produce four different glucosides of resveratrol: resveratrol 3-O-β-d-glucoside, resveratrol 4′-O-β-d-glucoside, resveratrol 3,5-O-β-d-diglucoside, and resveratrol 3,5,4′-O-β-d-triglucoside. The conversion rates and numbers of products formed were found to vary with the other NDP sugar donors. Resveratrol 3-O-β-d-2-deoxyglucoside and resveratrol 3,5-O-β-d-di-2-deoxyglucoside were found to be produced using TDP-2-deoxyglucose as a donor; however, the monoglycosides resveratrol 4′-O-β-d-galactoside, resveratrol 4′-O-β-d-viosaminoside, resveratrol 3-O-β-l-rhamnoside, and resveratrol 3-O-β-l-fucoside were produced from the respective sugar donors. Altogether, 10 diverse glycoside derivatives of the medically important resveratrol were generated, demonstrating the capacity of YjiC to produce structurally diverse resveratrol glycosides.     

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.     

Simple Y-Autosomal Incompatibilities Cause Hybrid Male Sterility in Reciprocal Crosses Between Drosophila virilis and D. americana

Postzygotic reproductive isolation evolves when hybrid incompatibilities accumulate between diverging populations. Here, I examine the genetic basis of hybrid male sterility between two species of Drosophila, Drosophila virilis and D. americana. From these analyses, I reach several conclusions. First, neither species carries any autosomal dominant hybrid male sterility alleles: reciprocal F₁ hybrid males are perfectly fertile. Second, later generation (backcross and F₂) hybrid male sterility between D. virilis and D. americana is not polygenic. In fact, I identified only three genetically independent incompatibilities that cause hybrid male sterility. Remarkably, each of these incompatibilities involves the Y chromosome. In one direction of the cross, the D. americana Y is incompatible with recessive D. virilis alleles at loci on chromosomes 2 and 5. In the other direction, the D. virilis Y chromosome causes hybrid male sterility in combination with recessive D. americana alleles at a single QTL on chromosome 5. Finally, in contrast with findings from other Drosophila species pairs, the X chromosome has only a modest effect on hybrid male sterility between D. virilis and D. americana.SPECIATION occurs when populations evolve one or more barriers to interbreeding (Dobzhansky 1937; Mayr 1963). One such barrier is intrinsic postzygotic isolation, which typically evolves when diverging populations accumulate different alleles at two or more loci that are incompatible when brought together in hybrid genomes; negative epistasis between these alleles renders hybrids inviable or sterile (Bateson 1909; Dobzhansky 1937; Muller 1942). Classical and recent studies in diverse animal taxa have provided support for two evolutionary patterns that often characterize the genetics of postzygotic isolation (Coyne and Orr 1989a). The first, Haldane''s rule, observes that when there is F₁ hybrid inviability or sterility that affects only one sex, it is almost always the heterogametic sex (Haldane 1922). Over the years, many researchers have tried to account for this pattern, but only two ideas are now thought to provide a general explanation: the “dominance theory,” which posits that incompatibility alleles are generally recessive in hybrids, and the “faster-male theory,” which posits that genes causing hybrid male sterility diverge more rapidly than those causing hybrid female sterility (Muller 1942; Wu and Davis 1993; Turelli and Orr 1995; reviewed in Coyne and Orr 2004). In some cases, however, additional factors might contribute to Haldane''s rule, including meiotic drive, a faster-evolving X chromosome, dosage compensation, and Y chromosome incompatibilities (reviewed in Laurie 1997; Turelli and Orr 2000; Coyne and Orr 2004).The second broad pattern affecting the evolution of postzygotic isolation is the disproportionately large effect of the X chromosome on heterogametic F₁ hybrid sterility (Coyne 1992). This “large X effect” has been documented in genetic analyses of backcross hybrid sterility (e.g., Dobzhansky 1936; Grula and Taylor 1980; Orr 1987; Masly and Presgraves 2007) and inferred from patterns of introgression across natural hybrid zones (e.g., Machado et al. 2002; Saetre et al. 2003; Payseur et al. 2004). However, in only one case has the cause of the large X effect been unambiguously determined: incompatibilities causing hybrid male sterility between Drosophila mauritiana and D. sechellia occur at a higher density on the X than on the autosomes (Masly and Presgraves 2007). Testing the generality of this pattern will require additional high-resolution genetic analyses in diverse taxa (Presgraves 2008). But whatever its causes, there is now general consensus that the X chromosome often plays a special role in the evolution of postzygotic isolation (Coyne and Orr 2004).The contribution of the Y chromosome to animal speciation is less clear. Y chromosomes have far fewer genes than the X or autosomes, and most of these genes are male specific (Lahn and Page 1997; Carvalho et al. 2009). In Drosophila species, the Y chromosome is typically required for male fertility, but not for viability (Voelker and Kojima 1971). How often, then, does the Y chromosome play a role in reproductive isolation? In crosses between Drosophila species, hybrid male sterility is frequently caused by incompatibilities between the X and Y chromosomes (Schafer 1978; Heikkinen and Lumme 1998; Mishra and Singh 2007) or between the Y and heterospecific autosomal alleles (Patterson and Stone 1952; Vigneault and Zouros 1986; Lamnissou et al. 1996). In crosses between D. yakuba and D. santomea, the Y chromosome causes F₁ hybrid male sterility, and accordingly, shows no evidence for recent introgression across a species hybrid zone (Coyne et al. 2004; Llopart et al. 2005). In mammals, reduced introgression of Y-linked loci (relative to autosomal loci) has been shown across natural hybrid zones of mice (Tucker et al. 1992) and rabbits (Geraldes et al. 2008), suggesting that the Y chromosome contributes to reproductive barriers.Here I examine the genetic basis of hybrid male sterility between two species of Drosophila, D. virilis and D. americana. These species show considerable genetic divergence (K_s ∼0.11, Morales-Hojas et al. 2008) and are currently allopatric: D. virilis is a human commensal worldwide with natural populations in Asia, and D. americana is found in riparian habitats throughout much of North America (Throckmorton 1982; McAllister 2002). Nearly 70 years ago, Patterson et al. (1942) showed that incompatibilities between the D. americana Y chromosome and the second and fifth chromosomes from D. virilis cause hybrid male sterility, a result that was confirmed in a more recent study (Lamnissou et al. 1996). Another study suggested that the X chromosome might play the predominant role in causing hybrid male sterility between D. virilis and D. americana (Orr and Coyne 1989). But because previous genetic analyses had to rely on only a few visible markers to map hybrid male sterility, they lacked the resolution to examine the genomic distribution of incompatibility loci.Using the D. virilis genome sequence, I have developed a dense set of molecular markers to investigate the genetic architecture of hybrid male sterility between D. virilis and D. americana. In this study, I perform a comprehensive set of crosses to address several key questions: What is the effect of the X chromosome on hybrid male sterility between D. virilis and D. americana? What is the effect of the Y chromosome? Approximately how many loci contribute to hybrid male sterility between these Drosophila species? Perhaps surprisingly, the answers to these questions differ dramatically from what has been found for other Drosophila species, including the well-studied D. melanogaster group.     

Antagonistic Roles of PP2A-Pab1 and Etd1 in the Control of Cytokinesis in Fission Yeast

In Schizosaccharomyces pombe, Etd1 is a positive regulator of the septation initiation network (SIN), a conserved GTPase-regulated kinase cascade that triggers cytokinesis. Here we show that a mutation in the pab1 gene, which encodes the B-regulatory subunit of the protein phosphatase 2A (PP2A), suppresses mutations in the etd1 gene. Etd1 is required for the function of the GTPase Spg1, a key regulator of SIN signaling. Interestingly, the loss of Pab1 function restored the activity of Spg1 in Etd1-deficient cells. This result suggests that PP2A-Pab1–mediated dephosphorylation inhibits Spg1, thus antagonizing Etd1 function. The loss of pab1 function also rescues the lethality of mutants of other genes in the SIN cascade such as mob1, sid1, and cdc11. Two-hybrid assays indicate that Pab1 physically interacts with Mob1, Sid1, Sid2, and Cdc11, suggesting that the phosphatase 2A B-subunit is a component of the SIN complex. Together, our results indicate that PP2A-Pab1 plays a novel role in cytokinesis, regulating SIN activity at different levels. Pab1 is also required to activate polarized cell growth. Thus, PP2A-Pab1 may be involved in coordinating polar growth and cytokinesis.THE fission yeast Schizosaccharomyces pombe is a leading experimental model for eukaryotic cytokinesis (Bathe and Chang 2009; Pollard and Wu 2010). Fission yeast cells grow in a polarized manner by elongation at the cell ends and divide during cytokinesis by the action of a contractile actomyosin ring assembled in the middle of the cell (Snell and Nurse 1993). At the end of mitosis, when nuclear separation has been completed, actomyosin ring constriction is triggered by the septation initiation network (SIN). This signal transduction cascade is composed of the GTPase Spg1 and three protein kinases—Cdc7, GC-kinase Sid1, and NDR-kinase Sid2 in their presumed order of action—and the associated proteins Cdc14 with Sid1 and Mob1 with Sid2. These proteins are all located at the spindle pole body (SPB) during mitosis on a scaffold composed of the coiled-coil proteins Sid4 and Cdc11 (Krapp et al. 2004). The Sid2-Mob1 protein kinase complex is thought to transmit the division signal from the SPB to the actomyosin ring since it also associates at the division site during septation (Krapp and Simanis 2008). The SIN triggers actomyosin ring contraction coordinated with the synthesis of the primary and secondary septa that will form the new cell wall (Krapp et al. 2004; Wolfe and Gould 2005). The small GTPase Rho1 is known to promote cell-wall formation at the division site by stimulation of Cps1p/Drc1 1,3-β-glucan synthase (Le Goff et al. 1999), but the mechanism remains unclear.SIN activity is tightly regulated during the cell cycle to ensure proper coordination of mitosis and cytokinesis. Mutants that negatively affect SIN function undergo nuclear division in the absence of septation, while increased SIN activity induces septation in interphase cells (Krapp and Simanis 2008). Regulation of the SIN is complex, involving multiple, partially redundant mechanisms, but the nucleotide status of the Ras superfamily small GTPase, Spg1, represents a key step in SIN activity (Lattmann et al. 2009). Cdc16 and Byr4 form a two-component GTPase-activating protein (GAP) for Spg1 that inhibits its activity (Furge et al. 1998; Cerutti and Simanis 1999). Proteins acting as a guanine nucleotide-exchange factor (GEF) for this GTPase have not been identified. In the budding yeast Saccharomyces cerevisiae, the pathway analogous to the SIN is known as the mitotic exit network (MEN) (reviewed in Krapp and Simanis 2008). Contact between the SPB-localized GTPase Tem1 (the Spg1 homolog) with its putative GEF Lte1, which is present only within the bud, has been proposed as a mechanism to ensure that mitotic exit occurs only after the spindle has oriented correctly (Bardin et al. 2000; Pereira et al. 2000). Bfa1-Bub2 (the Cdc16-Byr4 equivalent) are negative regulators of the MEN, acting as a two-component GAP for Tem1 (Geymonat et al. 2002).Etd1 was identified in a genetic screen searching for new regulators of the S. pombe cell division cycle (Jimenez and Oballe 1994). Further characterization indicated that Etd1 acts as a positive regulator of the SIN (Daga et al. 2005). A recent study has established a key role for Etd1 in the timing of cytokinesis via the regulation of Spg1, acting as a potential homolog of budding yeast Lte1 (Garcia-Cortes and McCollum 2009). Loss of Etd1 function can be suppressed by mutations in a number of genes, some of which are involved in morphogenesis (Jimenez and Oballe 1994). Here we show that one of the mutations that bypass the requirement for etd1 in cytokinesis affects the activity of pab1, which encodes the protein phosphatase 2A (PP2A) regulatory subunit B. The characterization of Pab1 and pab1 mutants described in this study reveals a novel role for PP2A-Pab1 in SIN regulation and provides new insight into the mechanism by which Etd1 might regulate SIN signaling. We also show that Pab1 participates in activation of the morphological pathway, suggesting a role for PP2A-Pab1 in the coordination of cytokinesis and morphogenesis.     

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.     

Bilirubin conjugates of human bile. The excretion of bilirubin as the acyl glycosides of aldobiouronic acid,pseudoaldobiouronic acid and hexuronosylhexuronic acid,with a branched-chain hexuronic acid as one of the components of hexuronosylhexuronide

Structure elucidations have been performed on the bilirubin conjugates isolated from human hepatic bile as the phenylazo derivatives. The major bilirubin conjugates are excreted, not as was formerly thought in the form of glucuronides, but as the acyl glycosides of aldobiouronic acid, pseudoaldobiouronic acid and hexuronosylhexuronic acid. The isolated aldobiouronides are proposed to have the structures of an acyl 6-O-hexopyranosyluronic acid-hexopyranoside, an acyl 4-O-hexofuranosyluronic acid-d-glucopyranoside, and an acyl 4-O-β-d-glucofuranosyluronic acid-d-glucopyranoside respectively, with the acyl radicals being those of the phenylazo derivative of bilirubin. The pseudoaldobiouronide is suggested to be the acyl 4-O-α-d-glucofuranosyl-β-d -glucopyranosiduronic acid, with the acyl radical being that of the phenylazo derivative of vinylneoxanthobilirubinic acid. The hexuronosylhexuronide presumably is the acyl 4-O-(3-C-hydroxymethylribofuranosyluronic acid)-β-d-glucopyranosiduronic acid, with the acyl radical being that of the phenylazo derivative of bilirubin. The 3-C-hydroxymethylriburonic acid, isolated as one of the components of the hexuronosylhexuronide, is the first natural branched-chain hexuronic acid to be detected, and the first branched-chain sugar ever detected in humans.     

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.     

Alternative Splicing Modulates Ubx Protein Function in Drosophila melanogaster

Polyploidy is an important aspect of the evolution of flowering plants. The potential of gene copies to diverge and evolve new functions is influenced by meiotic behavior of chromosomes leading to segregation as a single locus or duplicated loci. Switchgrass (Panicum virgatum) linkage maps were constructed using a full-sib population of 238 plants and SSR and STS markers to access the degree of preferential pairing and the structure of the tetraploid genome and as a step toward identification of loci underlying biomass feedstock quality and yield. The male and female framework map lengths were 1645 and 1376 cM with 97% of the genome estimated to be within 10 cM of a mapped marker in both maps. Each map coalesced into 18 linkage groups arranged into nine homeologous pairs. Comparative analysis of each homology group to the diploid sorghum genome identified clear syntenic relationships and collinear tracts. The number of markers with PCR amplicons that mapped across subgenomes was significantly fewer than expected, suggesting substantial subgenome divergence, while both the ratio of coupling to repulsion phase linkages and pattern of marker segregation indicated complete or near complete disomic inheritance. The proportion of transmission ratio distorted markers was relatively low, but the male map was more extensively affected by distorted transmission ratios and multilocus interactions, associated with spurious linkages.POLYPLOIDY is common among plants (Masterson 1994; Levin 2002) and is an important aspect of plant evolution. Widespread paleopolyploidy in flowering plant lineages suggests that ancient polyploidization events have contributed to the radiation of angiosperms (Soltis et al. 2009; Van de Peer et al. 2009a). Whole genome duplications are thought to be the sources of evolutionary novelty (Osborn et al. 2003; Freeling and Thomas 2006; Chen 2007; Hegarty and Hiscock 2008; Flagel and Wendel 2009; Leitch and Leitch 2008). Other attributes of polyploids considered to promote evolutionary success include increased vigor, masking of recessive alleles, and reproductive barriers arising from loss of one of the duplicate genes (Soltis and Soltis 2000; Comai 2005; Otto 2007; Van de Peer et al. 2009b). Among crop species, polyploidy likely contributed to trait improvement under artificial selection (Paterson 2005; Udall and Wendell 2006; Dubcovsky and Dvorak 2007; Hovav et al. 2008).Disomic inheritance in polyploids, in contrast to polysomic inheritance, presents opportunities for duplicated genes to diverge and evolve new functions. The relative age of whole genome duplications and the extent of homology between subgenomes greatly influence chromosomal pairing at meiosis (Soltis and Soltis 1995; Wolfe 2001; Ramsey and Schemske 2002). Polysomic inheritance resulting from random chromosome pairing is associated with doubling of a single set of chromosomes. Disomic inheritance resulting from preferential pairing is often associated with polyploidy arising from combinations of divergent genomes. The evolutionary process of diploidization leads to a shift from random to preferential pairing that is not well understood but is genetically defined in systems such as Ph1 of wheat (Triticum aestivum) and PrBn of Brassica napus (Riley and Chapman 1958; Vega and Feldman 1998; Jenczewski et al. 2003). The degree of preferential pairing also affects allelic diversity and the ability to detect linkage. Accurate information about chromosome pairing and whole or partial genome duplications is thus important for both evolutionary studies and in linkage analysis.Such information is extremely limited in the C4 panicoid species Panicum virgatum (switchgrass), which is now viewed as a promising energy crop in the United States and Europe (Lewandowski et al. 2003; McLaughlin and Kszos 2005) and is planted extensively for forage and soil conservation (Vogel and Jung 2001). Little is known about either its genome structure or inheritance. Much current bioenergy feedstock development is focused on tetraploid cytotypes (2n = 4x = 36) due to their higher yield potentials, and an initial segregation study indicated a high degree of preferential pairing in a single F1 mapping population (Missaoui et al. 2005). A once-dominant component of the tallgrass prairie in North America, switchgrass is largely self-incompatible (Martinez-Reyna and Vogel 2002) with predominantly tetraploid or octoploid cytotypes (Hultquist et al. 1997; Lu et al. 1998). Limited gene flow appears possible between different cytotypes suggested by DNA content variation within collection sites and seed lots (Nielsen 1944; Hultquist et al. 1997; Narasimhamoorthy et al. 2008). True diploids appear to be rare (Nielsen 1944; Young et al. 2010). Multivalents in meiosis have not been observed in tetraploids or F1 hybrids between upland and lowland tetraploids, although rare univalents occurred (Barnett and Carver 1967; Martinez-Reyna et al. 2001). However, polysomic inheritance may occur with random bivalent pairing (Howard and Swaminathan 1953).Sustainable production of switchgrass for bioenergy to meet the goal of reducing greenhouse gas emissions will require advances in feedstock production that include improvements in yield (Carroll and Somerville 2009). Switchgrass has extensive genetic diversity and potential for genetic improvements, but each cycle of phenotypic selection can take several years (McLaughlin and Kszos 2005; Parrish and Fike 2005; Bouton 2007). Detailed understanding of genome structure to enable efficient marker-assisted selection (MAS) can speed this process considerably. Complete linkage maps are therefore required to both understand chromosome pairing and allow MAS.We report the construction of the first complete linkage maps of two switchgrass genotypes. The linkage maps provide genetic evidence for disomic inheritance in lowland, tetraploid switchgrass. Gene-derived markers enabled a comparative analysis to sorghum, revealing syntenic relationships between the diploid sorghum genome and the tetraploid switchgrass subgenomes. Transmission ratio distortion and multilocus interactions were analyzed in detail to document their potential influence on map accuracy and map-based studies in switchgrass.     

Estimating the Parameters of Selection on Nonsynonymous Mutations in Drosophila pseudoobscura and D. miranda

We present the results of surveys of diversity in sets of >40 X-linked and autosomal loci in samples from natural populations of Drosophila miranda and D. pseudoobscura, together with their sequence divergence from D. affinis. Mean silent site diversity in D. miranda is approximately one-quarter of that in D. pseudoobscura; mean X-linked silent diversity is about three-quarters of that for the autosomes in both species. Estimates of the distribution of selection coefficients against heterozygous, deleterious nonsynonymous mutations from two different methods suggest a wide distribution, with coefficients of variation greater than one, and with the average segregating amino acid mutation being subject to only very weak selection. Only a small fraction of new amino acid mutations behave as effectively neutral, however. A large fraction of amino acid differences between D. pseudoobscura and D. affinis appear to have been fixed by positive natural selection, using three different methods of estimation; estimates between D. miranda and D. affinis are more equivocal. Sources of bias in the estimates, especially those arising from selection on synonymous mutations and from the choice of genes, are discussed and corrections for these applied. Overall, the results show that both purifying selection and positive selection on nonsynonymous mutations are pervasive.SURVEYS of DNA sequence diversity and divergence are shedding light on a number of questions in evolutionary genetics (for recent reviews, see Akey 2009; Sella et al. 2009). Two of the most important questions of this kind concern the distribution of selection coefficients against deleterious mutations affecting protein sequences and the proportion of amino acid sequence differences between related species that have been fixed by positive selection. Several different methods have been proposed for studying each of these questions, using different features of data on polymorphism and divergence at nonsynonymous and silent sites.For example, the parameters of the distribution of selection coefficients against deleterious amino acid mutations have been estimated by contrasting the numbers of nonsynonymous and silent within-species polymorphisms and fixed differences between species (Sawyer and Hartl 1992; Bustamante et al. 2002; Piganeau and Eyre-Walker 2003; Sawyer et al. 2007); by fitting the frequency spectra of nonsynonymous and silent variants to models of selection, mutation, and drift (Akashi 1999; Eyre-Walker et al. 2006; Keightley and Eyre-Walker 2007; Kryukov et al. 2007; Boyko et al. 2008; Eyre-Walker and Keightley 2009); or by comparing levels of nonsynonymous and silent diversities between species with different population sizes (Loewe and Charlesworth 2006; Loewe et al. 2006). The results of these different approaches generally agree in suggesting that there is a wide distribution of selection coefficients against nonsynonymous mutations and that the mean selection coefficient against heterozygous carriers of such mutations is very small. The results imply that a typical individual from a human population carries several hundred weakly deleterious mutations (Eyre-Walker et al. 2006; Kryukov et al. 2007; Boyko et al. 2008); for a typical Drosophila population, with its much higher level of variability, the number is probably an order of magnitude greater (Loewe et al. 2006; Keightley and Eyre-Walker 2007).The presence of this large load of slightly deleterious mutations in human and natural populations, most of which are held at low frequencies by natural selection, has many implications. From the point of view of understanding human genetic disease, it means that we have to face the likelihood that susceptibility to a disease can be influenced by variants at many loci, each with small effects (Kryukov et al. 2007). The pervasive presence of deleterious mutations throughout the genome contributes to inbreeding depression (Charlesworth and Willis 2009) and may mean that the effective population size is reduced by background selection effects, even in regions of the genome with normal levels of genetic recombination (Loewe and Charlesworth 2007). Their presence may contribute so strongly to Hill–Robertson effects (Hill and Robertson 1966; Felsenstein 1974) that they cause severely reduced levels of diversity and adaptation in low-recombination regions of the genome (Charlesworth et al. 2010) and create a selective advantage to maintaining nonzero levels of recombination (Keightley and Otto 2006; Charlesworth et al. 2010). In addition, having an estimate of the distribution of selection coefficients against deleterious nonsynonymous mutations allows their contribution to between-species divergence to be predicted, providing a way of estimating the fraction of fixed nonsynonymous differences caused by positive selection (Loewe et al. 2006; Boyko et al. 2008; Eyre-Walker and Keightley 2009).It is thus important to collect data that shed light on the properties of selection against nonsynonymous mutations in a wide range of systems and also to compare the results from different methods of estimation, since they are subject to different sources of difficulty and biases. In a previous study, we proposed the use of a comparison between two related species with different effective population sizes for this purpose (Loewe and Charlesworth 2006; Loewe et al. 2006), using Drosophila miranda and D. pseudoobscura as material. These are well suited for this type of study, as they are closely related, live together in similar habitats, and yet have very different levels of silent nucleotide diversity, indicating different effective population sizes (N_e). This study was hampered by our inability to compare the same set of loci across the two species and by the small number of loci that could be used. We here present the results of a much larger study of DNA variation at X-linked and autosomal loci for these two species, using D. affinis as a basis for estimating divergence. We compare the results, applying the method of Loewe et al. (2006) with that of Eyre-Walker and Keightley (2009) for estimating the distribution of deleterious selection coefficients and with McDonald–Kreitman test-based methods for estimating the proportion of nonsynonymous differences fixed by positive selection. While broadly confirming the conclusions from earlier studies, we note some possible sources of bias and describe methods for minimizing their effects.     

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).