The long-awaited public stock offer from third-generation sequencing technology company Pacific Biosciences has finally arrived (here’s the SEC filing, and coverage from Matthew Herper and GenomeWeb).
PacBio has already raised almost US$400 million in venture capital, and aims to increase this by up to US$200 million from its share offering. The sheer scale of these figures gives you a sense of just how much money is being thrown around in the race towards ever cheaper, faster and more accurate DNA sequencing machines. Of course, whether investors will decide to throw more money in PacBio’s direction remains to be seen.
The next few months for the company will be crucial: given the sheer amount of money now sloshing around behind PacBio, it needs to sell a
lot of machines to keep its investors happy.
The sequencing arms race has been heating up in other areas too: two weeks ago, sequencing service provider Complete Genomics also filed for a public offering – and was promptly sued by established competitor Illumina for alleged infringement of two patents.
Why did Illumina wait until after Complete’s announcement to launch the law-suit? Obviously we don’t know for sure, but the timing sure looks like a seriously impressive piece of evil genius: unsettled litigation is liable to make investors nervous, but restrictions on the public statements companies can make after announcing a public offering mean Complete is basically unable to defend itself in the media.
Interesting times indeed…
Read the comments on this post…
Forest trees are ideally suited to association mapping due to their high levels of diversity and low genomic linkage disequilibrium. Using an association mapping approach, single-nucleotide polymorphism (SNP) markers influencing quantitative variation in wood quality were identified in a natural population of Pinus radiata. Of 149 sites examined, 10 demonstrated significant associations (P < 0.05, q < 0.1) with one or more traits after accounting for population structure and experimentwise error. Without accounting for marker interactions, phenotypic variation attributed to individual SNPs ranged from 2 to 6.5%. Undesirable negative correlations between wood quality and growth were not observed, indicating potential to break negative correlations by selecting for individual SNPs in breeding programs. Markers that yielded significant associations were reexamined in an Australian land race. SNPs from three genes (PAL1, PCBER, and SUSY) yielded significant associations. Importantly, associations with two of these genes validated associations with density previously observed in the discovery population. In both cases, decreased wood density was associated with the minor allele, suggesting that these SNPs may be under weak negative purifying selection for density in the natural populations. These results demonstrate the utility of LD mapping to detect associations, even when the power to detect SNPs with small effect is anticipated to be low.
Efficient genomic selection in animals or crops requires the accurate prediction of the agronomic performance of individuals from their high-density molecular marker profiles. Using a training data set that contains the genotypic and phenotypic information of a large number of individuals, each marker or marker allele is associated with an estimated effect on the trait under study. These estimated marker effects are subsequently used for making predictions on individuals for which no phenotypic records are available. As most plant and animal breeding programs are currently still phenotype driven, the continuously expanding collection of phenotypic records can only be used to construct a genomic prediction model if a dense molecular marker fingerprint is available for each phenotyped individual. However, as the genotyping budget is generally limited, the genomic prediction model can only be constructed using a subset of the tested individuals and possibly a genome-covering subset of the molecular markers. In this article, we demonstrate how an optimal selection of individuals can be made with respect to the quality of their available phenotypic data. We also demonstrate how the total number of molecular markers can be reduced while a maximum genome coverage is ensured. The third selection problem we tackle is specific to the construction of a genomic prediction model for a hybrid breeding program where only molecular marker fingerprints of the homozygous parents are available. We show how to identify the set of parental inbred lines of a predefined size that has produced the highest number of progeny. These three selection approaches are put into practice in a simulation study where we demonstrate how the trade-off between sample size and sample quality affects the prediction accuracy of genomic prediction models for hybrid maize.
The genomics revolution has spurred the undertaking of HapMap studies of numerous species, allowing for population genomics to increase the understanding of how selection has created genetic differences between subspecies populations. The objectives of this study were to (1) develop an approach to detect signatures of selection in subsets of phenotypically similar breeds of livestock by comparing single nucleotide polymorphism (SNP) diversity between the subset and a larger population, (2) verify this method in breeds selected for simply inherited traits, and (3) apply this method to the dairy breeds in the International Bovine HapMap (IBHM) study. The data consisted of genotypes for 32,689 SNPs of 497 animals from 19 breeds. For a given subset of breeds, the test statistic was the parametric composite log likelihood (CLL) of the differences in allelic frequencies between the subset and the IBHM for a sliding window of SNPs. The null distribution was obtained by calculating CLL for 50,000 random subsets (per chromosome) of individuals. The validity of this approach was confirmed by obtaining extremely large CLLs at the sites of causative variation for polled (BTA1) and black-coat-color (BTA18) phenotypes. Across the 30 bovine chromosomes, 699 putative selection signatures were detected. The largest CLL was on BTA6 and corresponded to KIT, which is responsible for the piebald phenotype present in four of the five dairy breeds. Potassium channel-related genes were at the site of the largest CLL on three chromosomes (BTA14, -16, and -25) whereas integrins (BTA18 and -19) and serine/arginine rich splicing factors (BTA20 and -23) each had the largest CLL on two chromosomes. On the basis of the results of this study, the application of population genomics to farm animals seems quite promising. Comparisons between breed groups have the potential to identify genomic regions influencing complex traits with no need for complex equipment and the collection of extensive phenotypic records and can contribute to the identification of candidate genes and to the understanding of the biological mechanisms controlling complex traits.
A novel method, called linkage disequilibrium multilocus iterative peeling (LDMIP), for the imputation of phase and missing genotypes is developed. LDMIP performs an iterative peeling step for every locus, which accounts for the family data, and uses a forward–backward algorithm to accumulate information across loci. Marker similarity between haplotype pairs is used to impute possible missing genotypes and phases, which relies on the linkage disequilibrium between closely linked markers. After this imputation step, the combined iterative peeling/forward–backward algorithm is applied again, until convergence. The calculations per iteration scale linearly with number of markers and number of individuals in the pedigree, which makes LDMIP well suited to large numbers of markers and/or large numbers of individuals. Per iteration calculations scale quadratically with the number of alleles, which implies biallelic markers are preferred. In a situation with up to 15% randomly missing genotypes, the error rate of the imputed genotypes was <1% and ~99% of the missing genotypes were imputed. In another example, LDMIP was used to impute whole-genome sequence data consisting of 17,321 SNPs on a chromosome. Imputation of the sequence was based on the information of 20 (re)sequenced founder individuals and genotyping their descendants for a panel of 3000 SNPs. The error rate of the imputed SNP genotypes was 10%. However, if the parents of these 20 founders are also sequenced, >99% of missing genotypes are imputed correctly.
S1 is the most important locus acting as a reproductive barrier between Oryza sativa and O. glaberrima. It is a complex locus, with factors that may affect male and female fertility separately. Recently, the component causing the allelic elimination of pollen was fine mapped. However, the position and nature of the component causing female sterility remains unknown. To fine map the factor of the S1 locus affecting female fertility, we developed a mapping approach based on the evaluation of the degree of female transmission ratio distortion (fTRD) of markers. Through implementing this methodology in four O. sativa x O. glaberrima crosses, the female component of the S1 locus was mapped into a 27.8-kb (O. sativa) and 50.3-kb (O. glaberrima) region included within the interval bearing the male component of the locus. Moreover, evidence of additional factors interacting with S1 was also found. In light of the available data, a model where incompatibilities in epistatic interactions between S1 and the additional factors are the cause of the female sterility barrier between O. sativa and O. glaberrima was developed to explain the female sterility and the TRD mediated by S1. According to our model, the recombination ratio and allelic combinations between these factors would determine the final allelic frequencies observed for a given cross.
Loci involved in local adaptation can potentially be identified by an unusual correlation between allele frequencies and important ecological variables or by extreme allele frequency differences between geographic regions. However, such comparisons are complicated by differences in sample sizes and the neutral correlation of allele frequencies across populations due to shared history and gene flow. To overcome these difficulties, we have developed a Bayesian method that estimates the empirical pattern of covariance in allele frequencies between populations from a set of markers and then uses this as a null model for a test at individual SNPs. In our model the sample frequencies of an allele across populations are drawn from a set of underlying population frequencies; a transform of these population frequencies is assumed to follow a multivariate normal distribution. We first estimate the covariance matrix of this multivariate normal across loci using a Monte Carlo Markov chain. At each SNP, we then provide a measure of the support, a Bayes factor, for a model where an environmental variable has a linear effect on the transformed allele frequencies compared to a model given by the covariance matrix alone. This test is shown through power simulations to outperform existing correlation tests. We also demonstrate that our method can be used to identify SNPs with unusually large allele frequency differentiation and offers a powerful alternative to tests based on pairwise or global FST. Software is available at http://www.eve.ucdavis.edu/gmcoop/.
Mutation-accumulation experiments are widely used to estimate parameters of spontaneous mutations affecting fitness. In many experiments only one component of fitness is measured. In a previous study involving the diploid yeast Saccharomyces cerevisiae, we measured the growth rate of 151 mutation-accumulation lines to estimate parameters of mutation. We found that an unexpectedly high frequency of fitness-altering mutations was beneficial. Here, we build upon our previous work by examining sporulation efficiency, spore viability, and haploid growth rate and find that these components of fitness also show a high frequency of beneficial mutations. We also examine whether mutation-acycumulation (MA) lines show any evidence of pleiotropy among accumulated mutations and find that, for most, there is none. However, MA lines that have zero fitness (i.e., lethality) for any one fitness component do show evidence for pleiotropy among accumulated mutations. We also report estimates of other parameters of mutation based on each component of fitness.
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.
Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have “fast-forwarded” a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.
The general coalescent tree framework is a family of models for determining ancestries among random samples of DNA sequences at a nonrecombining locus. The ancestral models included in this framework can be derived under various evolutionary scenarios. Here, a computationally tractable full-likelihood-based inference method for neutral polymorphisms is presented, using the general coalescent tree framework and the infinite-sites model for mutations in DNA sequences. First, an exact sampling scheme is developed to determine the topologies of conditional ancestral trees. However, this scheme has some computational limitations and to overcome these limitations a second scheme based on importance sampling is provided. Next, these schemes are combined with Monte Carlo integrations to estimate the likelihood of full polymorphism data, the ages of mutations in the sample, and the time of the most recent common ancestor. In addition, this article shows how to apply this method for estimating the likelihood of neutral polymorphism data in a sample of DNA sequences completely linked to a mutant allele of interest. This method is illustrated using the data in a sample of DNA sequences at the APOE gene locus.
Much of population genetics is based on the diffusion limit of the Wright–Fisher model, which assumes a fixed population size. This assumption is violated in most natural populations, particularly for microbes. Here we study a more realistic model that decouples birth and death events and allows for a stochastically varying population size. Under this model, classical quantities such as the probability of and time before fixation of a mutant allele can differ dramatically from their Wright–Fisher expectations. Moreover, inferences about natural selection based on Wright–Fisher assumptions can yield erroneous and even contradictory conclusions: at small population densities one allele will appear superior, whereas at large densities the other allele will dominate. Consequently, competition assays in laboratory conditions may not reflect the outcome of long-term evolution in the field. These results highlight the importance of incorporating demographic stochasticity into basic models of population genetics.
We propose a multilocus version of FST 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 FST 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 FST than on the SNP-based version, which provides a potential explanation for these results. Thus, these new measures of FST and haplotype-cluster diversity provide an important new tool for population genetic analysis of high-density SNP data.
Juvenile hormone (JH) is critical for multiple aspects of insect development and physiology. Although roles for the hormone have received considerable study, an understanding of the molecules necessary for JH action in insects has been frustratingly slow to evolve. Methoprene-tolerant (Met) in Drosophila melanogaster fulfills many of the requirements for a hormone receptor gene. A paralogous gene, germ-cell expressed (gce), possesses homology and is a candidate as a Met partner in JH action. Expression of gce was found to occur at multiple times and in multiple tissues during development, similar to that previously found for Met. To probe roles of this gene in JH action, we carried out in vivo gce over- and underexpression studies. We show by overexpression studies that gce can substitute in vivo for Met, alleviating preadult but not adult phenotypic characters. We also demonstrate that RNA interference-driven knockdown of gce expression in transgenic flies results in preadult lethality in the absence of MET. These results show that (1) unlike Met, gce is a vital gene and shows functional flexibility and (2) both gene products appear to promote JH action in preadult but not adult development.
Regulation of cytoskeletal structure and dynamics is essential for multiple aspects of cellular behavior, yet there is much to learn about the molecular machinery underlying the coordination between the cytoskeleton and its effector systems. One group of proteins that regulate microtubule behavior and its interaction with other cellular components, such as actin-regulatory proteins and transport machinery, is the plus-end tracking proteins (MT+TIPs). In particular, evidence suggests that the MT+TIP, CLASP, may play a pivotal role in the coordination of microtubules with other cellular structures in multiple contexts, although the molecular mechanism by which it functions is still largely unknown. To gain deeper insight into the functional partners of CLASP, we conducted parallel genetic and proteome-wide screens for CLASP interactors in Drosophila melanogaster. We identified 36 genetic modifiers and 179 candidate physical interactors, including 13 that were identified in both data sets. Grouping interactors according to functional classifications revealed several categories, including cytoskeletal components, signaling proteins, and translation/RNA regulators. We focused our initial investigation on the MT+TIP Minispindles (Msps), identified among the cytoskeletal effectors in both genetic and proteomic screens. Here, we report that Msps is a strong modifier of CLASP and Abl in the retina. Moreover, we show that Msps functions during axon guidance and antagonizes both CLASP and Abl activity. Our data suggest a model in which CLASP and Msps converge in an antagonistic balance in the Abl signaling pathway.
Sensory communication depends on the precise matching between the emission and the perception of sex- and species-specific signals; understanding both the coevolutionary process and the genes involved in both production and detection is a major challenge. desat1 determines both aspects of communication—a mutation in desat1 simultaneously alters both sex pheromone emission and perception in Drosophila melanogaster flies. We investigated whether the alteration of pheromonal perception is a consequence of the altered production of pheromones or if the two phenotypes are independently controlled by the same locus. Using several genetic tools, we were able to separately manipulate the two pheromonal phenotypes, implying that desat1 is the sole gene responsible, exerting a pleiotropic effect on both transmission and detection. The levels of the five desat1 trancripts, measured in the head and body of manipulated flies, were related to variation in pheromone production. This suggests that the pleiotropic action of desat1 on pheromonal communication depends on the fine regulation of its transcriptional activity.
The cell surface receptor kinase BRASSINOSTEROID-INSENSITIVE-1 (BRI1) is the major receptor for steroid hormones in Arabidopsis. Plants homozygous for loss-of-function mutations in BRI1 display a reduction in the size of vegetative organs, resulting in dwarfism. The recessive bri1-5 mutation produces receptors that do not accumulate to wild-type levels and are retained mainly in the endoplasmic reticulum. We have isolated a dominant suppressor of the dwarf phenotype of bri1-5 plants. We show that this suppression is caused by a second-site mutation in BRI1, bri1-5R1. The bri1-5R1 mutation partially rescues the phenotypes of bri1-5 in many tissues and enhances bri1-5 phenotypes above wild-type levels in several other tissues. We demonstrate that the phenotypes of bri1-5R1 plants are due to both increased cell expansion and increased cell division. To test the mechanism of bri1-5 suppression, we assessed whether the phenotypic suppression in bri1-5R1 was dependent on ligand availability and the integrity of the signaling pathway. Our results indicate that the suppression of the dwarf phenotypes associated with bri1-5R1 requires both BR biosynthesis and the receptor kinase BRI1-ASSOCIATED KINASE-1 (BAK1). Finally, we show that bri1-5R1 partially restores the accumulation and plasma membrane localization of BRI1. Collectively, our results point toward a model in which bri1-R1 compensates for the protein-folding abnormalities caused by bri1-5, restoring accumulation of the receptor and its delivery to the cell surface.
Meiosis is a highly regulated process in eukaryotic species. The filamentous fungus Neurospora crassa has been shown to be missing homologs of a number of meiotic initiation genes conserved in Saccharomyces cerevisiae, but has three homologs of the well-characterized middle meiotic transcriptional regulator NDT80. In this study, we evaluated the role of all three NDT80 homologs in the formation of female reproductive structures, sexual development, and meiosis. We found that none of the NDT80 homologs were required for meiosis and that even the triple mutant was unaffected. However, strains containing mutations in NCU09915 (fsd-1) were defective in female sexual development and ascospore maturation. vib-1 was a major regulator of protoperithecial development in N. crassa, and double mutants carrying deletions of both vib-1 (NCU03725) and fsd-1 exhibited a synergistic effect on the timing of female reproductive structure (protoperithecia) formation. We further evaluated the role of the N. crassa homolog of IME2, a kinase involved in initiation of meiosis in S. cerevisiae. Strains containing mutations in ime-2 showed unregulated development of protoperithecia. Genetic analysis indicated that mutations in vib-1 were epistatic to ime-2, suggesting that IME-2 may negatively regulate VIB-1 activity. Our data indicate that the IME2/NDT80 pathway is not involved in meiosis in N. crassa, but rather regulates the formation of female reproductive structures.
While mitochondria are renowned for their role in energy production, they also perform several other integral functions within the cell. Thus, it is not surprising that mitochondrial dysfunction can negatively impact cell viability. Although mitochondria have received an increasing amount of attention in recent years, there is still relatively little information about how proper maintenance of mitochondria and its genomes is achieved. The Neurospora crassa mus-10 mutant was first identified through its increased sensitivity to methyl methanesulfonate (MMS) and was thus believed to be defective in some aspect of DNA repair. Here, we report that mus-10 harbors fragmented mitochondria and that it accumulates deletions in its mitochondrial DNA (mtDNA), suggesting that the mus-10 gene product is involved in mitochondrial maintenance. Interestingly, mus-10 begins to senesce shortly after deletions are visualized in its mtDNA. To uncover the function of MUS-10, we used a gene rescue approach to clone the mus-10 gene and discovered that it encodes a novel F-box protein. We show that MUS-10 interacts with a core component of the Skp, Cullin, F-box containing (SCF) complex, SCON-3, and that its F-box domain is essential for its function in vivo. Thus, we provide evidence that MUS-10 is part of an E3 ubiquitin ligase complex involved in maintaining the integrity of mitochondria and may function to prevent cellular senescence.
Sister chromatid cohesion refers to the process by which sister chromatids are tethered together until the metaphase-to-anaphase transition. The evolutionarily conserved cohesin complex mediates sister chromatid cohesion. Cohesin not only ensures proper chromosome segregation, but also promotes high-fidelity DNA repair and transcriptional regulation. Two subunits of cohesin (Smc1p, Smc3p) are members of the structural maintenance of chromosomes (SMC) family. The SMC family is recognized by their large coiled-coil arms and conserved ATP-binding cassette-like ATPase domain. While both Smc1p and Smc3p ATP binding and hydrolysis are essential for cohesin function in vivo, little is known about how this core enzymatic activity is regulated to facilitate sister chromatid cohesion. Here we use SMC mutant proteins to block specific steps in cohesin’s ATPase cycle in Saccharomyces cerevisiae. We show that blocking Smc3p-mediated ATP binding or Smc3p ATP hydrolysis traps unique functional states in cohesion. Finally, we provide evidence that Smc3p acetylation, which has an essential role in cohesion establishment, modulates the Smc3p ATP-bound state.
