Next-generation sequencing as a tool to quickly identify causative EMS-generated mutations

The advent of next generation sequencing has influenced every aspect of biological research. Many labs are now using whole genome sequencing in Arabidopsis thaliana as a means to quickly identify EMS-generated mutations present in isolated mutants. Following identification of these mutations, examination of T-DNA insertional alleles defective in candidate genes or complementation of the mutant phenotype with a wild type copy of candidate genes can be used to verify which mutation is causative for the phenotype of interest. Here, we discuss the benefits and pitfalls of using this method to identify mutations underlying phenotypes.     

Assignment of 44 Ds Insertions to the Linkage Map of Arabidopsis

Transposon tagging is a useful tool for biological studies. Transposon insertions have been used to obtain new mutants which are extremely helpful in understanding gene function. These insertions immediately provide a means to isolate the corresponding genes. Transposon tagging has also been used to clone genes previously defined by point mutations. In addition, transposon insertions into cloned genes that lack mutations can be generated to facilitate functional analysis. The maize Ac/Ds transposon elements are known to transpose to local sites with high frequencies and have been shown to function in several dicots. To generate a collection of Ds elements for the purpose of targeted insertional mutagenesis of mapped genes in Arabidopsis, we have mapped 44 Ds insertions by simple sequence length polymorphism (SSLP). Because the Arabidopsis genome project is advancing rapidly, many genes will be discovered whose functions are unknown. The mapped 44 Ds insertions will be a useful resource for post-genome analysis of gene functions in Arabidopsis.     

Testing hypotheses regarding the genetics of adaptation

Many of the hypotheses regarding the genetics of adaptation require that one know specific details about the genetic basis of complex traits, such as the number and effects of the loci involved. Developments in molecular biology have made it possible to create relatively dense maps of markers that can potentially be used to map genes underlying specific traits. However, there are a number of reasons to doubt that such mapping will provide the level of resolution necessary to specifically address many evolutionary questions. Moreover, evolutionary change is built upon the substitution of individual mutations, many of which may now be cosegregating in the same allele. In order for this developing area not to become a mirage that traps the efforts of an entire field, the genetic dissection of adaptive traits should be conducted within a strict hypothesis-testing framework and within systems that promise a reasonable chance of identifying the specific genetic changes of interest. Continuing advances in molecular technology may lead the way here, but some form of genetic testing is likely to be forever required.     

Dwarfism in Dexter cattle is not caused by the mutations in FGFR3 responsible for achondroplasia in humans

Dexter cattle carry a genetic defect causing a dwarf phenotype in the heterozygotes (Dx +/–), while homozygotes (Dx +/+) are stillborn with extreme shortening of limbs and gross craniofacial defects and are described as 'bulldog' calves. The heterozygous phenotype has been likened to achondroplastic dwarfism in humans (ACH), which has recently been shown to be the result of mutations in the transmembrane region of the fibroblast growth factor receptor 3 (FGFR3) gene. We have sequenced the transmembrane region of bovine FGFR3 from normal Dexter cattle (Dx -/-) and bulldog calves (Dx +/+). The sequence from both is identical and therefore excludes mutations in the transmembrane region of FGFR3 as the cause of Dexter dwarfism.     

Mapping endpoints of partial proteolysis fragments from regulatory subunit of type I cyclic AMP-dependent protein kinase

Methods for mapping endpoints of partial proteolysis fragments from regulatory subunit of type I cyclic AMP-dependent protein kinase are described with a view to using such data for fine-structure analysis of mutations and/or modifications affecting the protein's electrostatic charge. Peptides generated from [35S]methionine-labeled regulatory subunit were separated by high-resolution two-dimensional gel electrophoresis. Sites of papain cleavage in denatured regulatory subunit were deduced from the kinetics of the appearance, molecular weights, and relative isoelectric points of the fragments produced. These sites and sites of chymotrypsin digestion in the native protein were confirmed by studying peptide overlaps. Carboxy-terminal peptides were identified both by overlaps with cyclic AMP-binding chymotryptic fragments and by their preferential labeling during polysome runoff mediated by pactamycin, an inhibitor of protein initiation. Since peptides containing modifications or mutations that alter protein charge can be identified by shifts in first-dimension isoelectric focusing gel positions, knowledge of fragment endpoints will permit rapid mapping of sites of such alterations by two-dimensional gel analysis of partial proteolytic digests. Such a mapping procedure is inexpensive, can be applied to partially purified proteins or to proteins eluted from polyacrylamide gels, requires only nanogram amounts of the protein of interest, and does not require sequence data to determine relative positions of peptides. Therefore, it provides an attractive alternative to more classical peptide analysis for studying point mutations in cellular proteins of low abundance.     

Fast, accurate and simulation-free stochastic mapping

Mapping evolutionary trajectories of discrete traits onto phylogenies receives considerable attention in evolutionary biology. Given the trait observations at the tips of a phylogenetic tree, researchers are often interested where on the tree the trait changes its state and whether some changes are preferential in certain parts of the tree. In a model-based phylogenetic framework, such questions translate into characterizing probabilistic properties of evolutionary trajectories. Current methods of assessing these properties rely on computationally expensive simulations. In this paper, we present an efficient, simulation-free algorithm for computing two important and ubiquitous evolutionary trajectory properties. The first is the mean number of trait changes, where changes can be divided into classes of interest (e.g. synonymous/non-synonymous mutations). The mean evolutionary reward, accrued proportionally to the time a trait occupies each of its states, is the second property. To illustrate the usefulness of our results, we first employ our simulation-free stochastic mapping to execute a posterior predictive test of correlation between two evolutionary traits. We conclude by mapping synonymous and non-synonymous mutations onto branches of an HIV intrahost phylogenetic tree and comparing selection pressure on terminal and internal tree branches.     

高通量测序技术结合正向遗传学手段在基因定位研究中的应用

传统的利用正向遗传学方法的基因定位一般是通过构建遗传连锁图谱进行的,该过程步骤繁琐、耗时耗力,很多情形下定位精确度低、区间大。随着高通量测序技术的快速发展以及测序成本的不断降低,多种简单快捷的利用测序手段定位基因的方法被开发出来,包括对突变体基因组直接测序定位、突变体材料构建混池测序定位和遗传分离群体测序构建图谱定位等,还可以对转录组和部分基因组进行测序定位。这些方法可以在核苷酸水平鉴定突变位点,并已推广到复杂的遗传背景中。近期报道的一些测序定位甚至是在不依赖于参考基因组序列、遗传杂交和连锁信息的情况下完成的,这使得很多非模式物种也能开展正向遗传学研究。本文就这些新技术及其在基因定位中的应用进行了综述。     

Indel arrays: an affordable alternative for genotyping

Natural variation and induced mutations are important resources for gene discovery and the elucidation of genetic circuits. Mapping such polymorphisms requires rapid and cost-efficient methods for genome-wide genotyping. Here we report the development of a microarray-based method that assesses 240 unique markers in a single hybridization experiment at a cost of less than US$50 in materials per line. Our genotyping array is built with 70-mer oligonucleotide elements representing insertion/deletion (indel) polymorphisms between the Arabidopsis thaliana accessions Columbia-0 (Col) and Landsberg erecta (Ler). These indel polymorphisms are recognized with great precision by comparative genomic hybridization, eliminating the need for array replicates and complex statistical analysis. Markers are present genome-wide, with an average spacing of approximately 500 kb. PCR primer information is provided for all array indels, allowing rapid single-locus inquiries. Multi-well chips allow groups of 16 lines to be genotyped in a single experiment. We demonstrate the utility of the array for accurately mapping recessive mutations, RIL populations and mixed genetic backgrounds from accessions other than Col and Ler. Given the ease of use of shotgun sequencing to generate partial genomic sequences of unsequenced species, this approach is readily transferable to non-model organisms.     

Traffic Lines: New Tools for Genetic Analysis in Arabidopsis thaliana

Genetic analysis requires the ability to identify the genotypes of individuals in a segregating population. This task is straightforward if each genotype has a distinctive phenotype, but is difficult if these genotypes are phenotypically similar or identical. We show that Arabidopsis seeds homozygous or heterozygous for a mutation of interest can be identified in a segregating family by placing the mutation in trans to a chromosome carrying a pair of seed-expressed green and red fluorescent transgenes (a “traffic line”) that flank the mutation. Nonfluorescent seeds in the self-pollinated progeny of such a heterozygous plant are usually homozygous for the mutation, whereas seeds with intermediate green and red fluorescence are typically heterozygous for the mutation. This makes it possible to identify seedlings homozygous for mutations that lack an obvious seedling phenotype, and also facilitates the analysis of lethal or sterile mutations, which must be propagated in heterozygous condition. Traffic lines can also be used to identify progeny that have undergone recombination within a defined region of the genome, facilitating genetic mapping and the production of near-isogenic lines. We produced 488 transgenic lines containing single genome-mapped insertions of NAP:dsRED and NAP:eGFP in Columbia (330 lines) and Landsberg erecta (158 lines) and generated sets of traffic lines that span most regions of the Arabidopsis genome. We demonstrated the utility of these lines for identifying seeds of a specific genotype and for generating near-isogenic lines using mutations of WUSCHEL and SHOOTMERISTEMLESS. This new resource significantly decreases the effort and cost of genotyping segregating families and increases the efficiency of experiments that rely on the ability to detect recombination in a defined chromosomal segment.     

Exploring variation in the d(N)/d(S) ratio among sites and lineages using mutational mappings: applications to the influenza virus

We use a likelihood-based method for mapping mutations on a phylogeny in a way that allows for both site-specific and lineage-specific variation in selection intensity. The method accounts for many of the potential sources of bias encountered in mapping of mutations on trees while still being computationally efficient. We apply the method to a previously published influenza data set to investigate hypotheses about changes in selection intensity in influenza strains. Influenza virus is sometimes propagated in chicken cells for several generations before sequencing, a process that has been hypothesized to induce mutations adapting the virus to the lab medium. Our analysis suggests that there are approximately twice as many replacement substitutions in lineages propagated in chicken eggs as in lineages that are not. Previous studies have attempted to predict which viral strains future epidemics may arise from using inferences regarding positive selection. The assumption is that future epidemics are more likely to arise from the strains in which positive selection on the so-called “trunk lineages” of the evolutionary tree is most pervasive. However, we find no difference in the strength of selection in the trunk lineages versus other evolutionary lineages. Our results suggest that it may be more difficult to use inferences regarding the strength of selection on mutations to make predictions regarding viral epidemics than previously thought. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Willie Swanson     

Exhaustive mutational analysis of severe acute respiratory syndrome coronavirus 2 ORF3a: An essential component in the pathogen's infectivity cycle

Detailed knowledge of a protein's key residues may assist in understanding its function and designing inhibitors against it. Consequently, such knowledge of one of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)'s proteins is advantageous since the virus is the etiological agent behind one of the biggest health crises of recent times. To that end, we constructed an exhaustive library of bacteria differing from each other by the mutated version of the virus's ORF3a viroporin they harbor. Since the protein is harmful to bacterial growth due to its channel activity, genetic selection followed by deep sequencing could readily identify mutations that abolish the protein's function. Our results have yielded numerous mutations dispersed throughout the sequence that counteract ORF3a's ability to slow bacterial growth. Comparing these data with the conservation pattern of ORF3a within the coronavirinae provided interesting insights: Deleterious mutations obtained in our study corresponded to conserved residues in the protein. However, despite the comprehensive nature of our mutagenesis coverage (108 average mutations per site), we could not reveal all of the protein's conserved residues. Therefore, it is tempting to speculate that our study unearthed positions in the protein pertinent to channel activity, while other conserved residues may correspond to different functionalities of ORF3a. In conclusion, our study provides important information on a key component of SARS-CoV-2 and establishes a procedure to analyze other viroporins comprehensively.     

Fine‐mapping natural alleles: quantitative complementation to the rescue

Mapping the genes responsible for natural variation and divergence is a challenging task. Many studies have mapped genes to genomic regions or generated lists of candidates, but few studies have implicated specific genes with a high standard of evidence. I propose that combining recent advances in genomic engineering with a modified version of the quantitative complementation test will help turn candidate genes into causal genes. By creating loss‐of‐function mutations in natural strains, and using these mutations to quantitatively fail‐to‐complement natural alleles, fine mapping should be greatly facilitated. As an example, I propose that the CRISPR/Cas9 system could be combined with the FLP/FRT system to fine‐map genes in the numerous systems where inversions have frustrated these efforts.     

World-wide survey of an Accord insertion and its association with DDT resistance in Drosophila melanogaster

Previous work showed that insecticide resistance in Drosophila melanogaster is correlated with the insertion of an Accord-like element into the 5' region of the cytochrome P450 gene, Cyp6g1. Here, we study the distribution of the Accord-like element in 673 recently collected D. melanogaster lines from 34 world-wide populations. We also examine the extent of microsatellite variability along a 180-kilobase (kb) genomic region of chromosome II encompassing the resistance gene. We confirm a 100% correlation of the Accord insertion with insecticide resistance and a significant reduction in variability extending at least 20 kb downstream of the Cyp6g1 gene. The frequency of the Accord insertion differs significantly between East African (32-55%) and nonAfrican (85-100%) populations. This pattern is consistent with a selective sweep driving the Accord insertion close to fixation in nonAfrican populations as a result of the insecticide resistance phenotype it confers. This study confirms that hitchhiking mapping can be used to identify beneficial mutations in natural populations.     

The Human Obesity Gene Map: The 1998 Update

PÉRUSSE, LOUIS, YVON C. CHAGNON, JOHN WEISNAGEL, AND CLAUDE BOUCHARD. The human obesity gene map: the 1998 update. Obes Res. 1999;7:111–129. An update of the human obesity gene map incorporating published results up to the end of October 1998 is presented. Evidence from the human obesity cases caused by single gene mutations; other Mendelian disorders exhibiting obesity as a clinical feature; quantitative trait loci uncovered in human genome-wide scans and in crossbreeding experiments with mouse, rat, and pig models; association and case-control studies with candidate genes; and linkage studies with genes and other markers is reviewed. The most noticeable changes from the 1997 update is the number of obesity cases due to single gene mutations that increased from three cases due to mutations in two genes to 25 cases due to 12 mutations in seven genes. A look at the obesity gene map depicted in Figure 1 reveals that putative loci affecting obesity-related phenotypes are found on all but chromosome Y of the human chromosomes. Some chromosomes show at least three putative loci related to obesity on both arms (1, 2, 3, 6, 7, 8, 9, 11, 17, 19, 20, and X) and several on one chromosome arm only (4q, 5q, 10q, 12q, 13q, 15q, 16p, and 22q). The number of genes and other markers that have been associated or linked with human obesity phenotypes is increasing very rapidly and now approaches 27.     

Isolation and Chromosomal Localization of a Cornea-Specific Human Keratin 12 Gene and Detection of Four Mutations in Meesmann Corneal Epithelial Dystrophy

Keratin 12 (K12) is an intermediate-filament protein expressed specifically in corneal epithelium. Recently, we isolated K12 cDNA from a human corneal epithelial cDNA library and determined its full sequence. Herein, we present the exon-intron boundary structure and chromosomal localization of human K12. In addition, we report four K12 mutations in Meesmann corneal epithelial dystrophy (MCD), an autosomal dominant disorder characterized by intraepithelial microcysts and corneal epithelial fragility in which mutations in keratin 3 (K3) and K12 have recently been implicated. In the human K12 gene, we identified seven introns, defining eight individual exons that cover the coding sequence. Together the exons and introns span approximately 6 kb of genomic DNA. Using FISH, we found that the K12 gene mapped to 17q12, where a type I keratin cluster exists. In this study, four new K12 mutations (Arg135Gly, Arg135Ile, Tyr429Asp, and Leu140Arg) were identified in three unrelated MCD pedigrees and in one individual with MCD. All mutations were either in the highly conserved alpha-helix-initiation motif of rod domain 1A or in the alpha-helix-termination motif of rod domain 2B. These sites are essential for keratin filament assembly, suggesting that the mutations described above may be causative for MCD. Of particular interest, one of these mutations (Tyr429Asp), detected in both affected individuals in one of our pedigrees, is the first mutation to be identified within the alpha-helix-termination motif in type I keratin.     

Genomic organization and localization of the human CRMP-1 gene.

The Collapsin Response Mediator Protein-1 (CRMP-1) is a brain specific protein considered to be involved in the collapsin-induced growth cone collapse during neural development. CRMP-1 belongs to the Unc-33 gene family. Here we report the genomic structure and the localization of the human CRMP-1 gene to chromosome 4p16.1. Sequence analysis revealed that the human CRMP-1 gene consists of 14 exons. We have also established sequencing assays for all its coding exons. This should permit the rapid screening for mutations to assess CRMP-1 role in genetic disorders mapped in the 4p16.1 region.