Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^

Genome -wide and locus -specific instability

The discovery of expanded simple repeated sequences causing or associated with human disease has lead to a new area of research involved in the elucidation of how the expanded repeat causes disease and how the repeat becomes unstable. ^ To study the genetic basis of the (CTG)n repeat instability in the DMPK gene in myotonic dystrophy (DM1) patients, somatic cell hybrids were constructed between the lymphocytes of DM1 patients and a variety of Chinese hamster ovary (CHO) cell DNA repair gene deficient mutants. By using small pool PCR (SP-PCR), the instability of the (CTG)n can be quantitated for both the frequency and sizes of length change mutations. ^ Additional SP-PCR analysis on 2/11 subclones generated from this original hybrid showed a marked increase in large repeat deletions, ∼50%. A bimodal distribution of repeats was seen around the progenitor allele and at a large deleted product (within the normal range) with no intermediate products present. ^ To determine if the repair capacity of the CHO cell led to a mutator phenotype in the hamster and hybrid clones, SP-PCR was also done on 3 hamster microsatellites in a variety of hamster cell backgrounds. No variant alleles were seen in over 2500 genome equivalents screened. ^ Human-hamster hybrids have long been shown to be chromosomally unstable, yet information about the stability of repeated sequences was not known. To test if repeat instability was associated with either intact or non-intact human chromosomes, more than 300 microsatellite repeats on 13 human chromosomes (intact and non-intact) were analyzed in eight hybrid cells. No variants were seen between the hybrid and patient alleles in the hybrids. ^ To identify whether DM1 patients have a previously undetected level of genome wide instability or if the instability is truly locus specific, SP-PCR was done on 6 human microsatellites within the patient used to make the hybrid cells. No variants were seen in over 1000 genomes screened. ^ These studies show that the somatic cell hybrid approach is a genetically stable system that allows for the determination of factors that could lead to changes in microsatellite instability. It also shows that there is something inherent about the DM1 expanded (CTG)n repeat that it is solely targeted by, as of yet, and unknown mechanism that causes the repeat to be unstable. (Abstract shortened by UMI.)^

Molecular mechanisms involved in the regulation of cell polarity and genome stability by the p21 -activated kinase, Shk1, in the fission yeast, Schizosaccharomyces pombe

The essential p21-activated kinase (PAK), Shk1, is a critical component of a Ras/Cdc42/PAK complex required for cell viability, normal cell polarity, proper regulation of cytoskeletal dynamics, and sexual differentiation in the fission yeast, Schizosaccharomyces pombe. While cellular functions of PAKs have been described in eukaryotes from yeasts to mammals, the molecular mechanisms of PAK regulation and function are poorly understood. This study has characterized a novel Shk1 inhibitor, Skb15, and, in addition, identified the cell polarity regulator, Tea1, as a potential biological substrate of Shk1 in S. pombe. Skb15 is a highly conserved WD repeat protein that was discovered from a two-hybrid screen for proteins that interact with the catalytic domain of Shk1. Molecular data indicate that Skb15 negatively regulates Shk1 kinase activity in S. pombe cells. A null mutation in the skb15 gene is lethal and results in deregulation of actin polymerization and localization, microtubule biogenesis, and the cytokinetic machinery, as well as a substantial uncoupling of these processes from the cell cycle. Loss of Skb15 function is suppressed by partial loss of Shk1, demonstrating that negative regulation of Shk1 by Skb15 is required for proper execution of cytoskeletal remodeling and cytokinetic functions. A mouse homolog of Skb15 can substitute for its counterpart in fission yeast, demonstrating that Skb15 protein function has been substantially conserved through evolution. ^ Our laboratory has recently demonstrated that Shk1, in addition to regulating actin cytoskeletal organization, is required for proper regulation of microtubule dynamics in S. pombe cells. The Shk1 protein localizes to interphase and mitotic microtubules, the septum-forming region, and cell ends. This pattern of localization overlaps with that of the cell polarity regulator, Tea1, in S. pombe cells. The tea1 gene was identified by Paul Nurse's laboratory from a screen for genes involved in the control of cell morphogenesis in S. pombe. In contrast to wild type S. pombe cells, which are rod shaped, tea1 null cells are often bent and/or branched in shape. The Tea1 protein localizes to the cell ends, like Shk1, and the growing tips of interphase microtubules. Thus, experiments were performed to investigate whether Tea1 interacts with Shk1. The tea1 null mutation strongly suppresses the loss of function of Skb15, an essential inhibitor of Shk1 function. All defects associated with the skb15 mutation, including defects in F-actin organization, septation, spindle elongation, and chromosome segregation, are suppressed by tea1Δ, suggesting that Tea1 may function in these diverse processes. Consistent with a role for Tea1 in cytokinesis, tea1Δ cells have a modest cell separation defect that is greatly exacerbated by a shk1 mutation and, like Shk1, Tea1 localizes to the septation site. Molecular analyses showed that Tea1 phosphorylation is significantly dependent on Shk1 function in vivo and that bacterially expressed Tea1 protein is directly phosphorylated by recombinant Shk1 kinase in vitro. Taken together, these results identify Tea1 as a potential biological substrate of Shk1 in S. pombe. ^ In summary, this study provides new insights into a conserved regulatory mechanism for PAKs, and also begins to uncover the molecular mechanisms by which the Ras/Cdc42/PAK complex regulates the microtubule and actin cytoskeletons and cell growth polarization in fission yeast. ^

Genome instability quantified in constitutive tissues as well as putatively stable tumor tissue in hereditary non -polyposis colon cancer patients by small pool PCR

Microsatellite instability (MSI) is a hallmark of the mutator phenotype associated with Hereditary Non-Polyposis Colon Cancer (HNPCC). The MSI-High (MSI-H) HNPCC population has been well characterized, but the microsatellite low and stable (MSI-L/MSS) HNPCC population is much less understood. We hypothesize there are significant levels of MSI in HNPCC DNA classified as MSI-L/MSS, but no single variant allele makes up a sufficient population in the tumor DNA to be detected by standard analysis. Finding variants would suggest there is a mutator phenotype for the MSI-L/MSS HNPCC population that is distinct from the MSI-H HNPCC populations. This study quantified and compared MSI in HNPCC patients previously shown to be MSI-H, MSI-L/MSS and an MSI-H older, sporadic colorectal cancer patient. Small-pool Polymerase Chain Reactions (SP-PCRs) were conducted where the DNAs from each sample and controls are diluted into multiple pools, each containing approximately single genome equivalents. At least 100 alleles/sample were studied at six microsatellite loci. Mutant fragments were identified, quantified, and compared using Poisson statistics. Most of the variants were small deletions or insertions, with more mutants being deletions, as has been previously described in yeast and transgenic mice. SP-PCR, where most of the pools contained only 3 or less fragments, enabled identification of variants too infrequent to be detected by large pool PCR. Mutant fragments in positive control MSI-H tumor samples ranged from 0.26 to 0.68 in at least 4 of the 6 loci tested and were consistent with their MSI-H status. In the so called MSS tumors and constitutive tissues (normal colon tissue, and PBLs) of all the HNPCC patients, low, but significant levels of MSI were seen in at least two of the loci studied. This phenomenon was not seen in the sporadic MSI constitutive tissues nor the normal controls and suggests haploinsufficiency, gain-of-function, or a dominant/negative basis of the instability in HNPCC patients carrying germline mutations for tumor suppressor genes. A different frequency and spectrum of mutant fragments suggests a different genetic basis (other than a major mutation in MLH1 or MSH2) for disease in MSI-L and MSS HNPCC patients. ^

Whole Genome Analysis of Pinkeye Scores in Angus Cattle

Infectious Bovine Keratoconjunctivitis (IBK), known as pinkeye, is a common infectious disease affecting the eyes of cattle. It is characterized by excessive tearing, inflammation of the conjunctiva, and ulceration of the cornea. Although pinkeye is non-fatal, it has a marked economic impact on the cattle industry, due to the decreased performance of infected individuals. Genetic effects on the susceptibility of IBK have been studied and Hereford, Jersey, and Holstein breeds were found to be more susceptible to IBK than Bos Indicus breeds. The objectives of our study were: 1) to estimate genetic parameters of IBK scored in different categories by using genomic threshold model, and 2) to detect markers in linkage disequilibrium with quantitative tract loci (QTL) associated with IBK.

Genome-Wide Analysis of the Response of Dickeya dadantii 3937 to Plant Antimicrobial Peptides

Antimicrobial peptides constitute an important factor in the defense of plants against pathogens, and bacterial resistance to these peptides have previously been shown to be an important virulence factor in Dickeya dadantii, the causal agent of soft-rot disease of vegetables. In order to understand the bacterial response to antimicrobial pep- tides, a transcriptional microarray analysis was performed upon treatment with sub-lethal concentration of thionins, a widespread plant peptide. In all, 36 genes were found to be overexpressed, and were classified according to their deduced function as i) transcriptional regulators, ii) transport, and iii) modification of the bacterial membrane. One gene encoding a uricase was found to be repressed. The majority of these genes are known to be under the control of the PhoP/PhoQ system. Five genes representing the different functions induced were selected for further analysis. The results obtained indicate that the presence of antimicrobial peptides induces a complex response which includes peptide-specific elements and general stress-response elements contributing differentially to the virulence in different hosts.

Whole-genome genotyping of grape using a panel of microsatellite

The use of microsatellite markers in large-scale genetic studies is limited by its low throughput and high cost and labor requirements. Here, we provide a panel of 45 multiplex PCRs for fast and cost-efficient genome-wide fluorescence-based microsatellite analysis in grapevine. The developed multiplex PCRs panel (with up to 15-plex) enables the scoring of 270 loci covering all the grapevine genome (9 to 20 loci/chromosome) using only 45 PCRs and sequencer runs. The 45 multiplex PCRs were validated using a diverse grapevine collection of 207 accessions, selected to represent most of the cultivated Vitis vinifera genetic diversity. Particular attention was paid to quality control throughout the whole process (assay replication, null allele detection, ease of scoring). Genetic diversity summary statistics and features of electrophoretic profiles for each studied marker are provided, as are the genotypes of 25 common cultivars that could be used as references in other studies.

Genomic evolution during a 10,000-generation experiment with bacteria

Molecular methods are used widely to measure genetic diversity within populations and determine relationships among species. However, it is difficult to observe genomic evolution in action because these dynamics are too slow in most organisms. To overcome this limitation, we sampled genomes from populations of Escherichia coli evolving in the laboratory for 10,000 generations. We analyzed the genomes for restriction fragment length polymorphisms (RFLP) using seven insertion sequences (IS) as probes; most polymorphisms detected by this approach reflect rearrangements (including transpositions) rather than point mutations. The evolving genomes became increasingly different from their ancestor over time. Moreover, tremendous diversity accumulated within each population, such that almost every individual had a different genetic fingerprint after 10,000 generations. As has been often suggested, but not previously shown by experiment, the rates of phenotypic and genomic change were discordant, both across replicate populations and over time within a population. Certain pivotal mutations were shared by all descendants in a population, and these are candidates for beneficial mutations, which are rare and difficult to find. More generally, these data show that the genome is highly dynamic even over a time scale that is, from an evolutionary perspective, very brief.

Extent of genomic rearrangement after genome duplication in yeast

Whole-genome duplication approximately 108 years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70–100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150–200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Transgenic DNA integrated into the oat genome is frequently interspersed by host DNA

Integration of transgenic DNA into the plant genome was investigated in 13 transgenic oat (Avena sativa L.) lines produced using microprojectile bombardment with one or two cotransformed plasmids. In all transformation events, the transgenic DNA integrated into the plant genome consisted of intact transgene copies that were accompanied by multiple, rearranged, and/or truncated transgene fragments. All fragments of transgenic DNA cosegregated, indicating that they were integrated at single gene loci. Analysis of the structure of the transgenic loci indicated that the transgenic DNA was interspersed by the host genomic DNA. The number of insertions of transgenic DNA within the transgene loci varied from 2 to 12 among the 13 lines. Restriction endonucleases that do not cleave the introduced plasmids produced restriction fragments ranging from 3.6 to about 60 kb in length hybridizing to a probe comprising the introduced plasmids. Although the size of the interspersing host DNA within the transgene locus is unknown, the sizes of the transgene-hybridizing restriction fragments indicated that the entire transgene locus must be at least from 35–280 kb. The observation that all transgenic lines analyzed exhibited genomic interspersion of multiple clustered transgenes suggests a predominating integration mechanism. We propose that transgene integration at multiple clustered DNA replication forks could account for the observed interspersion of transgenic DNA with host genomic DNA within transgenic loci.

Kinetics of T cell receptor β, γ, and δ rearrangements during adult thymic development: T cell receptor rearrangements are present in CD44+CD25+ Pro-T thymocytes

We performed a comprehensive analysis of T cell receptor (TCR) γ rearrangements in T cell precursors of the mouse adult thymus. Using a sensitive quantitative PCR method, we show that TCRγ rearrangements are present in CD44+CD25+ Pro-T thymocytes much earlier than expected. TCRγ rearrangements increase significantly from the Pro-T to the CD44−CD25+ Pre-T cell transition, and follow different patterns depending on each Vγ gene segment, suggesting that ordered waves of TCRγ rearrangement exist in the adult mouse thymus as has been described in the fetal mouse thymus. Recombinations of TCRγ genes occur concurrently with TCRδ and D-Jβ rearrangements, but before Vβ gene assembly. Productive TCRγ rearrangements do not increase significantly before the Pre-T cell stage and are depleted in CD4+CD8+ double-positive cells from normal mice. In contrast, double-positive thymocytes from TCRδ−/− mice display random proportions of TCRγ rearranged alleles, supporting a role for functional TCRγ/δ rearrangements in the γδ divergence process.

Whole genome analysis: Experimental access to all genome sequenced segments through larger-scale efficient oligonucleotide synthesis and PCR

The recent ability to sequence whole genomes allows ready access to all genetic material. The approaches outlined here allow automated analysis of sequence for the synthesis of optimal primers in an automated multiplex oligonucleotide synthesizer (AMOS). The efficiency is such that all ORFs for an organism can be amplified by PCR. The resulting amplicons can be used directly in the construction of DNA arrays or can be cloned for a large variety of functional analyses. These tools allow a replacement of single-gene analysis with a highly efficient whole-genome analysis.

BAC-VAC, a novel generation of (DNA) vaccines: A bacterial artificial chromosome (BAC) containing a replication-competent, packaging-defective virus genome induces protective immunity against herpes simplex virus 1

This study aimed to exploit bacterial artificial chromosomes (BAC) as large antigen-capacity DNA vaccines (BAC-VAC) against complex pathogens, such as herpes simplex virus 1 (HSV-1). The 152-kbp HSV-1 genome recently has been cloned as an F-plasmid-based BAC in Escherichia coli (fHSV), which can efficiently produce infectious virus progeny upon transfection into mammalian cells. A safe modification of fHSV, fHSVΔpac, does not give rise to progeny virus because the signals necessary to package DNA into virions have been excluded. However, in mammalian cells fHSVΔpac DNA can still replicate, express the HSV-1 genes, cause cytotoxic effects, and produce virus-like particles. Because these functions mimic the lytic cycle of the HSV-1 infection, fHSVΔpac was expected to stimulate the immune system as efficiently as a modified live virus vaccine. To test this hypothesis, mice were immunized with fHSVΔpac DNA applied intradermally by gold-particle bombardment, and the immune responses were compared with those induced by infection with disabled infectious single cycle HSV-1. Immunization with either fHSVΔpac or disabled infectious single cycle HSV-1 induced the priming of HSV-1-specific cytotoxic T cells and the production of virus-specific antibodies and conferred protection against intracerebral injection of wild-type HSV-1 at a dose of 200 LD50. Protection probably was cell-mediated, as transfer of serum from immunized mice did not protect naive animals. We conclude that BAC-VACs per se, or in combination with genetic elements that support replicative amplification of the DNA in the cell nucleus, represent a useful new generation of DNA-based vaccination strategies for many viral and nonviral antigens.

Identification and classification of 16 new kinesin superfamily (KIF) proteins in mouse genome

KIF (kinesin superfamily) proteins are microtubule-dependent molecular motors that play important roles in intracellular transport and cell division. The extent to which KIFs are involved in various transporting phenomena, as well as their regulation mechanism, are unknown. The identification of 16 new KIFs in this report doubles the existing number of KIFs known in the mouse. Conserved nucleotide sequences in the motor domain were amplified by PCR using cDNAs of mouse nervous tissue, kidney, and small intestine as templates. The new KIFs were studied with respect to their expression patterns in different tissues, chromosomal location, and molecular evolution. Our results suggest that (i) there is no apparent tendency among related subclasses of KIFs of cosegregation in chromosomal mapping, and (ii) according to their tissue distribution patterns, KIFs can be divided into two classes–i.e., ubiquitous and specific tissue-dominant. Further characterization of KIFs may elucidate unknown fundamental phenomena underlying intracellular transport. Finally, we propose a straightforward nomenclature system for the members of the mouse kinesin superfamily.