GeneticMap

Morgan's three students (Muller, Sturtevant, and Bridges) introduced reductionist empirical methods to the study of the chromosomal theory of heredity. Herman J. Muller concentrated on mutations, namely changes in the heterocatalytic properties of genes, without losing their autocatalytic (self-replication) properties. Experimental induction of mutations allowed quantitative analyses of genes' parameters, but hopes to deduce their chemicophysical character were never fulfilled. Once the model for DNA structure was proposed, the reductionist notions of mutation analysis were successfully applied to the molecular genes. However, it was soon realized that the concept of the particulate gene was inadequate. The more the molecular analysis of the genome advanced, the clearer it became that the entities of heredity must be conceived within systems' perspectives, for which special tools for handling large number of variables were developed. Analytic mutagenesis, however, continues to be a major strategy for the study of the cellular and chromosomal mechanisms that control mutation inductions.

The variation of expression pattern exhibited by a transgene as a result of random integration, known as position effect, is, among other mechanisms, a particular challenge to reverse genetics. We present a strategy to counteract position effect in Arabidopsis thaliana by flanking the transgenes with the gypsy insulator from Drosophila melanogaster. In addition, Suppressor of Hairy-wing [Su(Hw)], the binding protein of the gypsy insulator, was coexpressed. Results indicated that the gypsy insulators could efficiently improve the expression levels of reporter genes driven by various kinds of promoters by 8- to 13-fold. Coexpression of the Su(Hw) protein led to a more uniform expression level of transgenes, as the coefficient of variation of expression levels was reduced further. The gypsy-Su(Hw) system enhanced expression levels, but did not alter the specificity of promoter activities, as experimentally evidenced by the promoters of the PIN and the AFB gene families. Interestingly, the gypsy insulator was also able to improve the expression of a selectable marker gene outside the insulated region, which facilitated the screen of transformants. Our system will likely decrease the number of lines that experimenters need to create and examine for a given transgene by contributing to relatively high and precise expression of transgenes in plants. Certain features of the gypsy insulator in Arabidopsis also provide new perspectives on the insulator field.

DNA replication and the correct packaging of DNA into different states of chromatin are both essential processes in all eukaryotic cells. High-fidelity replication of DNA is essential for the transmission of genetic material to cells. Likewise the maintenance of the epigenetic chromatin states is essential to the faithful reproduction of the transcriptional state of the cell. It is becoming more apparent that these two processes are linked through interactions between DNA replication proteins and chromatin-associated proteins. In addition, more proteins are being discovered that have dual roles in both DNA replication and the maintenance of epigenetic states. We present an analysis of two Drosophila mutants in the conserved DNA replication protein Mcm10. A hypomorphic mutant demonstrates that Mcm10 has a role in heterochromatic silencing and chromosome condensation, while the analysis of a novel C-terminal truncation allele of Mcm10 suggests that an interaction with Mcm2 is not required for chromosome condensation and heterochromatic silencing but is important for DNA replication.

Two mechanisms for chloroplast DNA replication have been revealed through the study of an unusual heteroplasmic strain of the green alga Chlamydomonas reinhardtii. Heteroplasmy is a state in which more than one genome type occurs in a mitochondrion or chloroplast. The Chlamydomonas strain spa19 bears two distinct chloroplast genomes, termed PS+ and PS–. PS+ genomes predominate and are stably maintained in vegetative cells, despite their lack of known replication origins. In sexual crosses with spa19 as the mating type plus parent, however, PS+ genomes are transmitted in only ~25% of tetrads, whereas the PS– genomes are faithfully inherited in all progeny. In this research, we have explored the mechanism underlying this biased uniparental inheritance. We show that the relative reduction and dilution of PS+ vs. PS– genomes takes place during gametogenesis. Bromodeoxyuridine labeling, followed by immunoprecipitation and PCR, was used to compare replication activities of PS+ and PS– genomes. We found that the replication of PS+ genomes is specifically suppressed during gametogenesis and germination of zygospores, a phenomenon that also was observed when spa19 cells were treated with rifampicin, an inhibitor of the chloroplast RNA polymerase. Furthermore, when bromodeoxyuridine incorporation was compared at 11 sites within the chloroplast genome between vegetative cells, gametes, and rifampicin-treated cells by quantitative PCR, we found that incorporation was often reduced at the same sites in gametes that were also sensitive to rifampicin treatment. We conclude that a transcription-mediated form of DNA replication priming, which may be downregulated during gametogenesis, is indispensable for robust maintenance of PS+ genomes. These results highlight the potential for chloroplast genome copy number regulation through alternative replication strategies.

The essential JIL-1 histone H3S10 kinase is a key regulator of chromatin structure that functions to maintain euchromatic domains while counteracting heterochromatization and gene silencing. In the absence of the JIL-1 kinase, two of the major heterochromatin markers H3K9me2 and HP1a spread in tandem to ectopic locations on the chromosome arms. Here we address the role of the third major heterochromatin component, the zinc-finger protein Su(var)3-7. We show that the lethality but not the chromosome morphology defects associated with the null JIL-1 phenotype to a large degree can be rescued by reducing the dose of the Su(var)3-7 gene and that Su(var)3-7 and JIL-1 loss-of-function mutations have an antagonistic and counterbalancing effect on position-effect variegation (PEV). Furthermore, we show that in the absence of JIL-1 kinase activity, Su(var)3-7 gets redistributed and upregulated on the chromosome arms. Reducing the dose of the Su(var)3-7 gene dramatically decreases this redistribution; however, the spreading of H3K9me2 to the chromosome arms was unaffected, strongly indicating that ectopic Su(var)3-9 activity is not a direct cause of lethality. These observations suggest a model where Su(var)3-7 functions as an effector downstream of Su(var)3-9 and H3K9 dimethylation in heterochromatic spreading and gene silencing that is normally counteracted by JIL-1 kinase activity.

A variety of cellular factors affect the movement of the retrovirus-like transposon Ty1. To identify genes involved in Ty1 virus-like particle (VLP) function, the level of the major capsid protein (Gag-p45) and its proteolytic precursor (Gag-p49p) was monitored in a subset of Ty1 cofactor mutants. Twenty-nine of 87 mutants contained alterations in the level of Gag; however, only bud22 showed a striking defect in Gag processing. BUD22 affected the +1 translational frameshifting event required to express the Pol proteins protease, integrase, and reverse transcriptase. Therefore, it is possible that the bud22 mutant may not produce enough functional Ty1 protease to completely process Gag-p49 to p45. Furthermore, BUD22 is required for 18S rRNA processing and 40S subunit biogenesis and influences polysome density. Together our results suggest that BUD22 is involved in a step in ribosome biogenesis that not only affects general translation, but also may alter the frameshifting efficiency of ribosomes, an event central to Ty1 retrotransposition.

Carbon dioxide (CO₂) sensing and metabolism via carbonic anhydrases (CAs) play pivotal roles in survival and proliferation of pathogenic fungi infecting human hosts from natural environments due to the drastic difference in CO₂ levels. In Cryptococcus neoformans, which causes fatal fungal meningoencephalitis, the Can2 CA plays essential roles during both cellular growth in air and sexual differentiation of the pathogen. However the signaling networks downstream of Can2 are largely unknown. To address this question, the present study employed comparative transcriptome DNA microarray analysis of a C. neoformans strain in which CAN2 expression is artificially controlled by the CTR4 (copper transporter) promoter. The P_CTR4::CAN2 strain showed growth defects in a CO₂-dependent manner when CAN2 was repressed but resumed normal growth when CAN2 was overexpressed. The Can2-dependent genes identified by the transcriptome analysis include FAS1 (fatty acid synthase 1) and GPB1 (G-protein β subunit), supporting the roles of Can2 in fatty acid biosynthesis and sexual differentiation. Cas3, a capsular structure designer protein, was also discovered to be Can2-dependent and yet was not involved in CO₂-mediated capsule induction. Most notably, a majority of Can2-dependent genes were environmental stress-regulated (ESR) genes. Supporting this, the CAN2 overexpression strain was hypersensitive to oxidative and genotoxic stress as well as antifungal drugs, such as polyene and azole drugs, potentially due to defective membrane integrity. Finally, an oxidative stress-responsive Atf1 transcription factor was also found to be Can2-dependent. Atf1 not only plays an important role in diverse stress responses, including thermotolerance and antifungal drug resistance, but also represses melanin and capsule production in C. neoformans. In conclusion, this study provides insights into the comprehensive signaling networks orchestrated by CA/CO₂-sensing pathways in pathogenic fungi.

Numerous studies have shown that the clinical antidepressant sertraline (Zoloft) is biologically active in model systems, including fungi, which do not express its putative protein target, the serotonin/5-HT transporter, thus demonstrating the existence of one or more secondary targets. Here we show that in the absence of its putative protein target, sertraline targets phospholipid membranes that comprise the acidic organelles of the intracellular vesicle transport system by a mechanism consistent with the bilayer couple hypothesis. On the basis of a combination of drug-resistance selection and chemical-genomic screening, we hypothesize that loss of vacuolar ATPase activity reduces uptake of sertraline into cells, whereas dysregulation of clathrin function reduces the affinity of membranes for sertraline. Remarkably, sublethal doses of sertraline stimulate growth of mutants with impaired clathrin function. Ultrastructural studies of sertraline-treated cells revealed a phenotype that resembles phospholipidosis induced by cationic amphiphilic drugs in mammalian cells. Using reconstituted enzyme assays, we also demonstrated that sertraline inhibits phospholipase A₁ and phospholipase D, exhibits mixed effects on phospholipase C, and activates phospholipase A₂. Overall, our study identifies two evolutionarily conserved membrane-active processes—vacuolar acidification and clathrin-coat formation—as modulators of sertraline's action at membranes.

The elongator (ELP) complex consisting of Elp1-6p has been indicated to play roles in multiple cellular processes. In yeast, the ELP complex has been shown to genetically interact with Uba4p/Urm1p and Kti11-13p for a function in tRNA modification. Through a Caenorhabditis elegans genetic suppressor screen and positional cloning, we discovered that loss-of-function mutations of moc-3 and dph-3, orthologs of the yeast UBA4 and KTI11, respectively, effectively suppress the Multivulva (Muv) phenotype of the lin-1(e1275, R175Opal) mutation. These mutations do not suppress the Muv phenotype caused by other lin-1 alleles or by gain-of-function alleles of ras or raf that act upstream of lin-1. The suppression can also be reverted by RNA interference of lin-1. Furthermore, we showed that dph-3(lf) also suppressed the defect of lin-1(e1275) in promoting the expression of a downstream target (egl-17). These results indicate that suppression by the moc-3 and dph-3 mutations is due to the elevated activity of lin-1(e1275) itself rather than the altered activity of a factor downstream of lin-1. We further showed that loss-of-function mutations of urm-1 and elpc-1-4, the worm counterparts of URM1 and ELP complex components in yeast, also suppressed lin-1(e1275). We also confirmed that moc-3(lf) and dph-3(lf) have defects in tRNA modifications as do the mutants of their yeast orthologs. These results, together with the observation of a likely readthrough product from a lin-1(e1275)::gfp fusion transgene indicate that the aberrant tRNA modification led to failed recognition of a premature stop codon in lin-1(e1275). Our genetic data suggest that the functional interaction of moc-3/urm-1 and dph-3 with the ELP complex is an evolutionarily conserved mechanism involved in tRNA functions that are important for accurate translation.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 F_ST. Software is available at http://www.eve.ucdavis.edu/gmcoop/.

S₁ 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 S₁ 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 S₁ 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 S₁ was also found. In light of the available data, a model where incompatibilities in epistatic interactions between S₁ 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 S₁. According to our model, the recombination ratio and allelic combinations between these factors would determine the final allelic frequencies observed for a given cross.

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.

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.

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.

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.

Directional epistasis describes a situation in which epistasis consistently increases or decreases the effect of allele substitutions, thereby affecting the amount of additive genetic variance available for selection in a given direction. This study applies a recent parameterization of directionality of epistasis to empirical data. Data stems from a QTL mapping study on an intercross between inbred mouse (Mus musculus) strains LG/J and SM/J, originally selected for large and small body mass, respectively. Results show a negative average directionality of epistasis for body-composition traits, predicting a reduction in additive allelic effects and in the response to selection for increased size. Focusing on average modification of additive effect of single loci, we find a more complex picture, whereby the effects of some loci are enhanced consistently across backgrounds, while effects of other loci are decreased, potentially contributing to either enhancement or reduction of allelic effects when selection acts at single loci. We demonstrate and discuss how the interpretation of the overall measurement of directionality depends on the complexity of the genotype–phenotype map. The measure of directionality changes with the power of scale in a predictable way; however, its expected effect with respect to the modification of additive genetic effects remains constant.

Transposons of the Mutator (Mu) superfamily have been shown to play a critical role in the evolution of plant genomes. However, the identification of Mutator transposons in other eukaryotes has been quite limited. Here we describe a previously uncharacterized group of DNA transposons designated Phantom identified in the genomes of a wide range of eukaryotic taxa, including many animals, and provide evidence for its inclusion within the Mutator superfamily. Interestingly three Phantom proteins were also identified in two insect viruses and phylogenetic analysis suggests horizontal movement from insect to virus, providing a new line of evidence for the role of viruses in the horizontal transfer of DNA transposons in animals. Many of the Phantom transposases are predicted to harbor a FLYWCH domain in the amino terminus, which displays a WRKY–GCM1 fold characteristic of the DNA binding domain (DBD) of Mutator transposases and of several transcription factors. While some Phantom elements have terminal inverted repeats similar in length and structure to Mutator elements, some display subterminal inverted repeats (sub-TIRs) and others have more complex termini reminiscent of so-called Foldback (FB) transposons. The structural plasticity of Phantom and the distant relationship of its encoded protein to known transposases may have impeded the discovery of this group of transposons and it suggests that structure in itself is not a reliable character for transposon classification.

The distal arm of the fourth ("dot") chromosome of Drosophila melanogaster is unusual in that it exhibits an amalgamation of heterochromatic properties (e.g., dense packaging, late replication) and euchromatic properties (e.g., gene density similar to euchromatic domains, replication during polytenization). To examine the evolution of this unusual domain, we undertook a comparative study by generating high-quality sequence data and manually curating gene models for the dot chromosome of D. virilis (Tucson strain 15010–1051.88). Our analysis shows that the dot chromosomes of D. melanogaster and D. virilis have higher repeat density, larger gene size, lower codon bias, and a higher rate of gene rearrangement compared to a reference euchromatic domain. Analysis of eight "wanderer" genes (present in a euchromatic chromosome arm in one species and on the dot chromosome in the other) shows that their characteristics are similar to other genes in the same domain, which suggests that these characteristics are features of the domain and are not required for these genes to function. Comparison of this strain of D. virilis with the strain sequenced by the Drosophila 12 Genomes Consortium (Tucson strain 15010–1051.87) indicates that most genes on the dot are under weak purifying selection. Collectively, despite the heterochromatin-like properties of this domain, genes on the dot evolve to maintain function while being responsive to changes in their local environment.