Discovery of chemically induced mutations in rice by TILLING

Background  

Rice is both a food source for a majority of the world's population and an important model system. Available functional genomics resources include targeted insertion mutagenesis and transgenic tools. While these can be powerful, a non-transgenic, unbiased targeted mutagenesis method that can generate a range of allele types would add considerably to the analysis of the rice genome. TILLING (Targeting Induced Local Lesions in Genomes), a general reverse genetic technique that combines traditional mutagenesis with high throughput methods for mutation discovery, is such a method.     

Discovery of rare mutations in populations: TILLING by sequencing

Discovery of rare mutations in populations requires methods, such as TILLING (for Targeting Induced Local Lesions in Genomes), for processing and analyzing many individuals in parallel. Previous TILLING protocols employed enzymatic or physical discrimination of heteroduplexed from homoduplexed target DNA. Using mutant populations of rice (Oryza sativa) and wheat (Triticum durum), we developed a method based on Illumina sequencing of target genes amplified from multidimensionally pooled templates representing 768 individuals per experiment. Parallel processing of sequencing libraries was aided by unique tracer sequences and barcodes allowing flexibility in the number and pooling arrangement of targeted genes, species, and pooling scheme. Sequencing reads were processed and aligned to the reference to identify possible single-nucleotide changes, which were then evaluated for frequency, sequencing quality, intersection pattern in pools, and statistical relevance to produce a Bayesian score with an associated confidence threshold. Discovery was robust both in rice and wheat using either bidimensional or tridimensional pooling schemes. The method compared favorably with other molecular and computational approaches, providing high sensitivity and specificity.     

TILLING to detect induced mutations in soybean

Background

Potato is a staple food in the diet of the world's population and also being used as animal feed. Compared to other crops, however, potato tubers are relatively poor in the essential amino acid, methionine. Our aim was to increase the methionine content of tubers by co-expressing a gene involved in methionine synthesis with a gene encoding a methionine-rich storage protein in potato plants.

Results

In higher plants, cystathionine γ-synthase (CgS) is the first enzyme specific to methionine biosynthesis. We attempted to increase the methionine content of tubers by expressing the deleted form of theArabidopsis CgS (CgS _Δ90), which is not regulated by methionine, in potato plants. To increase the incorporation of free methionine into a storage protein theCgS _Δ90was co-transformed with the methionine-rich15-kD β-zein. Results demonstrated a 2- to 6-fold increase in the free methionine content and in the methionine content of the zein-containing protein fraction of the transgenic tubers. In addition, in line with higher methionine content, the amounts of soluble isoleucine and serine were also increased. However, all of the lines with high level of CgS_Δ90 expression were phenotypically abnormal showing severe growth retardation, changes in leaf architecture and 40- to 60% reduction in tuber yield. Furthermore, the colour of the transgenic tubers was altered due to the reduced amounts of anthocyanin pigments. The mRNA levels of phenylalanine ammonia-lyase (PAL), the enzyme catalysing the first step of anthocyanin synthesis, were decreased.

Conclusion

Ectopic expression of CgS_Δ90 increases the methionine content of tubers, however, results in phenotypic aberrations in potato. Co-expression of the 15-kD β-zein with CgS_Δ90 results in elevation of protein-bound methionine content of tubers, but can not overcome the phenotypical changes caused by CgS_Δ90 and can not significantly improve the nutritional value of tubers. The level ofPAL mRNA and consequently the amount of anthocyanin pigments are reduced in the CgS_Δ90 transgenic tubers suggesting that methionine synthesis and production of anthocyanins is linked.     

TILLING by sequencing to identify induced mutations in stress resistance genes of peanut (Arachis hypogaea)

Background

Targeting Induced Local Lesions in Genomes (TILLING) is a powerful reverse genetics approach for functional genomics studies. We used high-throughput sequencing, combined with a two-dimensional pooling strategy, with either minimum read percentage with non-reference nucleotide or minimum variance multiplier as mutation prediction parameters, to detect genes related to abiotic and biotic stress resistances. In peanut, lipoxygenase genes were reported to be highly induced in mature seeds infected with Aspergillus spp., indicating their importance in plant-fungus interactions. Recent studies showed that phospholipase D (PLD) expression was elevated more quickly in drought sensitive lines than in drought tolerant lines of peanut. A newly discovered lipoxygenase (LOX) gene in peanut, along with two peanut PLD genes from previous publications were selected for TILLING. Additionally, two major allergen genes Ara h 1 and Ara h 2, and fatty acid desaturase AhFAD2, a gene which controls the ratio of oleic to linoleic acid in the seed, were also used in our study. The objectives of this research were to develop a suitable TILLING by sequencing method for this allotetraploid, and use this method to identify mutations induced in stress related genes.

Results

We screened a peanut root cDNA library and identified three candidate LOX genes. The gene AhLOX7 was selected for TILLING due to its high expression in seeds and roots. By screening 768 M2 lines from the TILLING population, four missense mutations were identified for AhLOX7, three missense mutations were identified for AhPLD, one missense and two silent mutations were identified for Ara h 1.01, three silent and five missense mutations were identified for Ara h 1.02, one missense mutation was identified for AhFAD2B, and one silent mutation was identified for Ara h 2.02. The overall mutation frequency was 1 SNP/1,066 kb. The SNP detection frequency for single copy genes was 1 SNP/344 kb and 1 SNP/3,028 kb for multiple copy genes.

Conclusions

Our TILLING by sequencing approach is efficient to identify mutations in single and multi-copy genes. The mutations identified in our study can be used to further study gene function and have potential usefulness in breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1348-0) contains supplementary material, which is available to authorized users.     

A modified TILLING approach to detect induced mutations in tetraploid and hexaploid wheat

Background  

Wheat (Triticum ssp.) is an important food source for humans in many regions around the world. However, the ability to understand and modify gene function for crop improvement is hindered by the lack of available genomic resources. TILLING is a powerful reverse genetics approach that combines chemical mutagenesis with a high-throughput screen for mutations. Wheat is specially well-suited for TILLING due to the high mutation densities tolerated by polyploids, which allow for very efficient screens. Despite this, few TILLING populations are currently available. In addition, current TILLING screening protocols require high-throughput genotyping platforms, limiting their use.     

Multiple pyruvate decarboxylase genes in maize are induced by hypoxia

Two cDNA clones corresponding to anaerobically induced maize mRNAs were found to have homology to a previously identified maize pyruvate decarboxylase gene. DNA sequencing and RFLP mapping indicate that these cDNAs represent two additional maize pdc genes. Each of the clones is approximately 85% identical in predicted amino acid sequence to the other two. All three clones are induced by hypoxic stress, but with different levels and kinetics of induction.     

Analysis of point mutations induced by ultraviolet light in human cells

Mutations induced in cultured human cells by 254-nm UV light were analyzed within exon 3 of the hypoxanthine guanine phosphoribosyl transferase (HPRT) gene. Five large independent cultures of human lymphoblastoid cells, line TK6, were exposed to 4 J/m2 of 254-nm UV light and mutants at the HPRT locus were selected en masse by 6-thioguanine (6TG) resistance. Exon 3 of the HPRT gene was amplified from the mutant cells by polymerase chain reaction (PCR) using modified T7 DNA polymerase. Denaturing gradient gel electrophoresis (DGGE) was used to separate the mutant sequences from the wild type as mutant/wild-type heteroduplexes. Individual mutant bands were isolated from the gel and the nature of the mutations was determined by direct sequencing. Eight predominant mutations were detected in the 184-bp exon 3 sequence. Of these, 3 transition, including 2 G-C to A-T and 1 A-T to G-C and 2 A-T to C-G transversions, appeared in all 5 UV-treated cultures but not in untreated cultures and were thus considered to be mutational hotspots. These observations are similar in nature to those previously reported in bacterial and rodent cells. A single G deletion, a tandem substitution of CpT for TpA, and a tandem triple substitution of GpGpA for ApApG were also observed but in only 2, 2 and 3 of the 5 UV-treated cultures, respectively. Numerical analysis of the mutant fractions of these 8 mutations indicated that each of them was distributed as a set of non-random and independent events, i.e., a mutational hotspot.     

Identification of novel alleles induced by EMS-mutagenesis in key genes of kernel hardness and starch biosynthesis in wheat by TILLING

To identify novel allelic variations in key genes of wheat quality, the present study used the targeting induced local lesions in genomes platform to detect point mutations in target genes. The wheat variety Longfumai 17 was treated by the mutagen ethyl methanesulfonate to produce a bulk M₂ generation, and the population included 1122 plants. A total length of 3906.80 kb nucleotides was analyzed, and the average mutation density was 1/244.17 kb. The identified mutations included G>A substitutions (43.75%), C>T substitutions (31.25%), A insertions (12.50%), T insertions (6.25%), and deletions (6.25%). These point mutations led to changes in amino acids and thus the encoded protein sequences, ultimately producing 18.75% of missense mutations, 12.50% of frame shift mutations, 6.25% of nonsense mutations, 25.00% of silent mutations and 37.50% of non-coding region mutations. In the kernel hardness gene Pinb and 3 starch synthesis genes waxy, Agp2 and SSIIa-A, we detected 16 different point mutations in 25 mutant lines. The Pinb gene harbored two missense mutations and a nonsense mutation; the C>T missense mutation resulted in a novel allele, this novel allele and the nonsense mutation alerted protein 3D structure; the waxy gene presented missense and frame shift mutations; the Agp2 gene carried a missense mutation; the SSIIa-A incurred a missense mutation and a frame shift mutation that resulted in premature protein termination. All the frame shift mutations, nonsense mutations and the Pinb novel allele resulted in allelic variation of their corresponding genes, which in turn affected their gene functions. The identified mutant lines can be used as intermediate materials in wheat quality improvement schemes.     

Characterisation of alleles of tomato light signalling genes generated by TILLING

Targeting Induced Local Lesions IN Genomes (TILLING) combines chemical mutagenesis with high throughput screening to allow the generation of alleles of selected genes. In this study, TILLING has been applied to produce a series of mutations in genes encoding essential components of the tomato light signal transduction pathway in an attempt to enhance fruit nutritional quality. Point mutations to DEETIOLATED1 (DET1), which is responsible for the high pigment2 (hp2) tomato mutant, resulted in elevated levels of both carotenoid and phenylpropanoid phytonutrients in ripe fruit, whilst immature fruit showed increased chlorophyll content, photosynthetic capacity and altered fruit morphology. Furthermore, genotypes with mutations to the UV-DAMAGED DNA BINDING PROTEIN 1 (DDB1), COP1 and COP1like were also characterised. These genotypes largely did not display phenotypes characteristic of mutation to light signalling components but their characterisation has enabled interrogation of structure function relationships of the mutated genes.     

Precise excision of transposons and point mutations induced by chemicals.

The ability of 23 chemicals (carcinogens and non-carcinogens) to induce precise excision of Tn10 and point mutations was studied in experiments with a single strain. The mutation assay was shown to detect a wider spectrum of genotoxic agents than the assay of Tn10 precise excision. The latter was induced only by potent SOS mutagens, which is in accordance with data on the SOS dependence of the induction of precise excision of Tn10. The precise excision assay as an additional test contributing to the knowledge of particular features of the action of a tested mutagen is discussed. The induction of precise excision of Tn10 by pyrene (and its failure to induce point mutations in this strain) demonstrates the value of using the transposon excision assay in cases of 'problem' mutagens.     

Helicobacter pylori genes involved in avoidance of mutations induced by 8-oxoguanine

Chromosomal rearrangements and base substitutions contribute to the large intraspecies genetic diversity of Helicobacter pylori. Here we explored the base excision repair pathway for the highly mutagenic 8-oxo-7,8-dihydroguanine (8-oxoG), a ubiquitous form of oxidized guanine. In most organisms, 8-oxoG is removed by a specific DNA glycosylase (Fpg in bacteria or OGG1 in eukaryotes). In the case where replication of the lesion yields an A/8-oxoG base pair, a second DNA glycosylase (MutY) can excise the adenine and thus avoid the fixation of the mutation in the next round of replication. In a genetic screen for H. pylori genes complementing the hypermutator phenotype of an Escherichia coli fpg mutY strain, open reading frame HP0142, a putative MutY coding gene, was isolated. Besides its capacity to complement E. coli mutY strains, HP0142 expression resulted in a strong adenine DNA glycosylase activity in E. coli mutY extracts. Consistently, the purified protein also exhibited such an activity. Inactivation of HP0142 in H. pylori resulted in an increase in spontaneous mutation frequencies. An Mg-dependent AP (abasic site) endonuclease activity, potentially allowing the processing of the abasic site resulting from H. pylori MutY activity, was detected in H. pylori cell extracts. Disruption of HP1526, a putative xth homolog, confirmed that this gene is responsible for the AP endonuclease activity. The lack of evidence for an Fpg/OGG1 functional homolog is also discussed.     

Long-culm mutations with dominant genes are induced by mPing transposon in rice

Development of semidwarf rice cultivars contributed to the epoch of high yielding crops called the 'Green Revolution'. However, over-reliance on semidwarf rice also has intrinsic limitations to supply food for an ever expanding world population. As a solution to the food supply problem, we propose the development of 'tall dwarf' rice cultivars that are characterized by increased biomass with long culms or large grains. However, genetic studies on the elongation of rice culms have remained scarce. This study seeks to analyze mutant genes involved in culm elongation in long-culm mutants induced by the MITE transposon mPing, which has been shown to be active in the japonica cultivar Gimbozu. Through analysis of the experimental results, we have confirmed that the three mutant long-culm genes exhibit genetic dominance. These represent rare cases of artificially induced dominant mutations. It is very likely that the mPing transposons played an important role in inducing the dominant mutations and also play an evolutionary interesting role.     

A functional genomics resource for Brassica napus: development of an EMS mutagenized population and discovery of FAE1 point mutations by TILLING

Two ethylmethanesulfonate (EMS) mutant populations of the semi-winter rapeseed cv. Ningyou7 were constructed with high mutant load, to provide a TILLING platform for functional genomics in Brassica napus, and for introduction of novel allelic variation in rapeseed breeding. Forward genetic screening of mutants from the M2 populations resulted in identification of a large number of novel phenotypes. Reverse genetic screening focused on the potentially multi-paralogous gene FAE1 (fatty acid elongase1), which controls seed erucic acid synthesis in rapeseed. A B. napus BAC library was screened, and loci in a reference mapping population (TNDH) were mapped to conclude that there are two paralogous copies of FAE1, one on each of the B. napus A and C genomes. A new procedure is demonstrated to identify novel mutations in situations where two or more very similar paralogous gene copies exist in a genome. The procedure involves TILLING of single plants, using existing SNPs as a positive control, and is able to distinguish novel mutations based on primer pairs designed to amplify both FAE1 paralogues simultaneously. The procedure was applied to 1344 M2 plants, with 19 mutations identified, of which three were functionally compromised with reduced seed erucic acid content.     

Stochastic traits of molecular evolution--acceptance of point mutations in native actin genes

A stochastic matrix of nucleotide mutation probabilities is derived by counting differences and identities in alignments of native actin genes, with the aim of obtaining a more reliable data base for regular modes of molecular evolution. The evolution of DNA sequences is thereby considered as a Markov process consisting of events (point mutations) characterized by a stochastic matrix for codon-codon interchanges. The genetic distance is set to 1 PAM (percentage of accepted point mutations). The results can be reproduced by Monte Carlo simulations which are subjected to selective constraints. The latter are observed as nonrandom codon usage and ratios of silent to recognizable point mutations. Specific patterns within the matrix of mutation probabilities attest to preferences of natural selection in the evolution of a specific protein.     

Stability of deletion,insertion and point mutations at the bronze locus in maize

Summary Phenotypic revertants from several kinds of mutations, including deletions, have been detected by pollen analysis at the wx and Adh loci in maize. Mutations in these genes give phenotypic revertants with median frequencies of 0.7 and 0.5×10^–5, respectively. However, the nature of such revertants can only be analyzed following their recovery from conventional matings. In the current study large seed populations derived from crosses involving several bz (bronze) locus mutations in maize were examined for reversion to a Bz (purple) expression. Deletion, insertion and point mutations were included in the study. Principally, over 2 million gametes of the bz-R mutation, which is shown here to be associated with a 340 base pair deletion within the transcribed region of the gene, have been screened for reversion. No revertants from it or any of the other bz mutations have been recovered, even though a total of almost 5 million gametes from homoallelic crosses have been examined to date. Results from seed analysis are discussed in reference to those from pollen analysis in maize.     

Computational analysis of the quaternary structural changes induced by point mutations in human UDP-glucose dehydrogenase

UDP-glucose dehydrogenase (UGDH) is an enzyme catalyzing the conversion of UDP-glucose to UDP-glucuronic acid. Site-directed mutagenesis studies have revealed that human UGDH (hUGDH) has distinct oligomeric states that vary with different point mutations. In this study we have investigated how the changes in the oligomer-forming propensity may be involved in the thermal motion of wild-type hUGDH and its mutants, using normal mode analysis (NMA). Our results show that the perturbation caused by the mutation of a residue at a considerably distant location from the oligomeric interfaces is preferentially distributed throughout specific sites, especially the large flexible regions in the hUGDH structure, thereby changing the motional fluctuation pattern at the oligomeric interfaces. A large-magnitude cooperative motion at the oligomeric interfaces is a critical factor in interfering with the hexamer formation of the enzyme. In particular, structural stability at the dimeric interface is necessary to retain the hexameric structure of hUGDH.