Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing

Genetics of destemming in pepper: A step towards mechanical harvesting

Introduction: The majority of peppers in the US for fresh market and processing are handpicked, and harvesting can account for 20-50% of production costs. Innovation in mechanical harvesting would increase availability; lower the costs of local, healthy vegetable products; and perhaps improve food safety and expand markets. Most processed peppers require removal of pedicels (stem and calyx) from the fruit, but lack of an efficient mechanical process for this operation has hindered adoption of mechanical harvest. In this paper, we present characterization and advancements in breeding green chile peppers for mechanical harvesting. Specifically, we describe inheritance and expression of an easy-destemming trait derived from the landrace UCD-14 that facilitates machine harvest of green chiles. Methods: A torque gauge was used for measuring bending forces similar to those of a harvester and applied to two biparental populations segregating for destemming force and rate. Genotyping by sequencing was used to generate genetic maps for quantitative trait locus (QTL) analyses. Results: A major destemming QTL was found on chromosome 10 across populations and environments. Eight additional population and/or environment-specific QTL were also identified. Chromosome 10 QTL markers were used to help introgress the destemming trait into jalapeño-type peppers. Low destemming force lines combined with improvements in transplant production enabled mechanical harvest of destemmed fruit at a rate of 41% versus 2% with a commercial jalapeńo hybrid. Staining for the presence of lignin at the pedicel/fruit boundary indicated the presence of an abscission zone and homologs of genes known to affect organ abscission were found under several QTL, suggesting that the easy-destemming trait may be due to the presence and activation of a pedicel/fruit abscission zone. Conclusion: Presented here are tools to measure the easy-destemming trait, its physiological basis, possible molecular pathways, and expression of the trait in various genetic backgrounds. Mechanical harvest of destemmed mature green chile fruits was achieved by combining easy-destemming with transplant management.

Novel Plant Breeding Techniques Shake Hands with Cereals to Increase Production

Cereals are the main source of human food on our planet. The ever-increasing food demand, continuously changing environment, and diseases of cereal crops have made adequate production a challenging task for feeding the ever-increasing population. Plant breeders are striving their hardest to increase production by manipulating conventional breeding methods based on the biology of plants, either self-pollinating or cross-pollinating. However, traditional approaches take a decade, space, and inputs in order to make crosses and release improved varieties. Recent advancements in genome editing tools (GETs) have increased the possibility of precise and rapid genome editing. New GETs such as CRISPR/Cas9, CRISPR/Cpf1, prime editing, base editing, dCas9 epigenetic modification, and several other transgene-free genome editing approaches are available to fill the lacuna of selection cycles and limited genetic diversity. Over the last few years, these technologies have led to revolutionary developments and researchers have quickly attained remarkable achievements. However, GETs are associated with various bottlenecks that prevent the scaling development of new varieties that can be dealt with by integrating the GETs with the improved conventional breeding methods such as speed breeding, which would take plant breeding to the next level. In this review, we have summarized all these traditional, molecular, and integrated approaches to speed up the breeding procedure of cereals.

Next Generation Sequencing Based Forward Genetic Approaches for Identification and Mapping of Causal Mutations in Crop Plants: A Comprehensive Review

The recent advancements in forward genetics have expanded the applications of mutation techniques in advanced genetics and genomics, ahead of direct use in breeding programs. The advent of next-generation sequencing (NGS) has enabled easy identification and mapping of causal mutations within a short period and at relatively low cost. Identifying the genetic mutations and genes that underlie phenotypic changes is essential for understanding a wide variety of biological functions. To accelerate the mutation mapping for crop improvement, several high-throughput and novel NGS based forward genetic approaches have been developed and applied in various crops. These techniques are highly efficient in crop plants, as it is relatively easy to grow and screen thousands of individuals. These approaches have improved the resolution in quantitative trait loci (QTL) position/point mutations and assisted in determining the functional causative variations in genes. To be successful in the interpretation of NGS data, bioinformatics computational methods are critical elements in delivering accurate assembly, alignment, and variant detection. Numerous bioinformatics tools/pipelines have been developed for such analysis. This article intends to review the recent advances in NGS based forward genetic approaches to identify and map the causal mutations in the crop genomes. The article also highlights the available bioinformatics tools/pipelines for reducing the complexity of NGS data and delivering the concluding outcomes.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

Genotyping-by-Sequencing and Its Application to Oat Genomic Research

Genotyping-by-sequencing (GBS) has emerged as a useful genomic approach for sampling genome-wide genetic variation, performing genome-wide association mapping, and conducting genomic selection. It is a combined one-step process of SNP marker discovery and genotyping through genome reduction with restriction enzymes and SNP calling with or without a sequenced genome. This approach has the advantage of being rapid, high throughput, cost effective, and applicable to organisms without sequenced genomes. It has been increasingly applied to generate SNP genotype data for plant genetic and genomic studies. To facilitate a wider GBS application, particularly in oat genetic and genomic research, we describe the GBS approach, review the current applications of GBS in plant species, and highlight some applications of GBS to oat research. We also discuss issues in various applications of GBS and provide some perspectives in GBS research. Recent developments of bioinformatics pipelines in high-quality SNP discovery for polyploid crops will enhance the application of GBS to oat genetic and genomic research.

Methods for Genotyping-by-Sequencing

A major goal for biologists is to understand the connection between genes and phenotypic traits, and genetic mapping in experimental populations remains a powerful approach for discovering the causal genes underlying phenotypes. For genetic mapping, the process of genotyping was previously a major rate-limiting step. Modern sequencing technology has greatly improved the resolution and speed of genetic mapping by reducing the time, labor, and cost per genotyping marker. In addition, the ability to perform genotyping-by-sequencing (GBS) has facilitated large-scale population genetic analyses by providing a simpler way to survey segregating genetic variation in natural populations. Here we present two protocols for GBS, using the Illumina platform, that can be applied to a wide range of genotyping projects in different species. The first protocol is for genotyping a subset of marker positions genome-wide using restriction digestion, and the second is for preparing inexpensive paired-end whole-genome libraries. We discuss the suitability of each approach for different genotyping applications and provide notes for adapting these protocols for use with a liquid-handling robot.

A New Advanced Backcross Tomato Population Enables High Resolution Leaf QTL Mapping and Gene Identification

Quantitative Trait Loci (QTL) mapping is a powerful technique for dissecting the genetic basis of traits and species differences. Established tomato mapping populations between domesticated tomato (Solanum lycopersicum) and its more distant interfertile relatives typically follow a near isogenic line (NIL) design, such as the S. pennellii Introgression Line (IL) population, with a single wild introgression per line in an otherwise domesticated genetic background. Here, we report on a new advanced backcross QTL mapping resource for tomato, derived from a cross between the M82 tomato cultivar and S. pennellii This so-called Backcrossed Inbred Line (BIL) population is comprised of a mix of BC2 and BC3 lines, with domesticated tomato as the recurrent parent. The BIL population is complementary to the existing S. pennellii IL population, with which it shares parents. Using the BILs, we mapped traits for leaf complexity, leaflet shape, and flowering time. We demonstrate the utility of the BILs for fine-mapping QTL, particularly QTL initially mapped in the ILs, by fine-mapping several QTL to single or few candidate genes. Moreover, we confirm the value of a backcrossed population with multiple introgressions per line, such as the BILs, for epistatic QTL mapping. Our work was further enabled by the development of our own statistical inference and visualization tools, namely a heterogeneous hidden Markov model for genotyping the lines, and by using state-of-the-art sparse regression techniques for QTL mapping.

Application of genotyping by sequencing technology to a variety of crop breeding programs

Since the Arabidopsis genome was completed, draft sequences or pseudomolecules have been published for more than 100 plant genomes including green algae, in large part due to advances in sequencing technologies. Advanced DNA sequencing technologies have also conferred new opportunities for high-throughput low-cost crop genotyping, based on single-nucleotide polymorphisms (SNPs). However, a recurring complication in crop genotyping that differs from other taxa is a higher level of DNA sequence duplication, noting that all angiosperms are thought to have polyploidy in their evolutionary history. In the current article, we briefly review current genotyping methods using next-generation sequencing (NGS) technologies. We also explore case studies of genotyping-by-sequencing (GBS) applications to several crops differing in genome size, organization and breeding system (paleopolyploids, neo-allopolyploids, neo-autopolyploids). GBS typically shows good results when it is applied to an inbred diploid species with a well-established reference genome. However, we have also made some progress toward GBS of outcrossing species lacking reference genomes and of polyploid populations, which still need much improvement. Regardless of some limitations, low-cost and multiplexed genotyping offered by GBS will be beneficial to breed superior cultivars in many crop species.

High-Throughput Sequencing Reveals Single Nucleotide Variants in Longer-Kernel Bread Wheat

The transcriptomes of bread wheat Yunong 201 and its ethyl methanesulfonate derivative Yunong 3114 were obtained by next-sequencing technology. Single nucleotide variants (SNVs) in the wheat strains were explored and compared. A total of 5907 and 6287 non-synonymous SNVs were acquired for Yunong 201 and 3114, respectively. A total of 4021 genes with SNVs were obtained. The genes that underwent non-synonymous SNVs were significantly involved in ATP binding, protein phosphorylation, and cellular protein metabolic process. The heat map analysis also indicated that most of these mutant genes were significantly differentially expressed at different developmental stages. The SNVs in these genes possibly contribute to the longer kernel length of Yunong 3114. Our data provide useful information on wheat transcriptome for future studies on wheat functional genomics. This study could also help in illustrating the gene functions of the non-synonymous SNVs of Yunong 201 and 3114.

Sequencing-based high throughput mutation detection in bread wheat

Background

Forward genetic approaches have limited use for agronomic traits that can’t be reliably scored on a single plant basis. Thus, mutants in wheat and other crops are more useful for gene function studies by reverse genetic approach. With a long-term goal to develop a sequence-based mutation detection resource in hexaploid wheat, we conducted a feasibility study to accurately differentiate induced mutations from the homoeologs’ sequence variations present among the three wheat genomes.

Results

A reduced representation ApeKI library consisting of 21 Ethylmethane Sulfonate (EMS) induced mutants and two wild type cv. Indian plants was developed using individual barcode adapters and sequenced. A novel bioinformatics pipeline was developed to identify sequence variants using 178,464 wheat unigenes as a reference wheat transcriptome. In total, 14,130 mutational changes [Single Nucleotide Polymorphisms (SNPs) and Insertions/Deletions (INDELs)] and 150,511 homoeologous sequence changes were detected. On an average, 662 SNPs (ranging from 46 to 1,330) and 10 small INDELs (ranging from 0 to 23) were identified for each of the mutants. A mutation frequency of one per 5 Kb was observed with 70 % being transitions and 30 % transversions. The pipeline was tested using the known sequence changes in the three wheat genes. Genes present in the distal regions of the chromosomes were found to be more prone to EMS compared to genes present in the proximal regions. Redefined parameters identified a total of 28,348 mutational changes (1,349/plant).

Conclusions

We conclude that sequencing based mutation detection is a valuable method to identify induced mutations at large.

Electronic supplementary material

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

Domestication selected for deceleration of the circadian clock in cultivated tomato

The circadian clock is a critical regulator of plant physiology and development, controlling key agricultural traits in crop plants. In addition, natural variation in circadian rhythms is important for local adaptation. However, quantitative modulation of circadian rhythms due to artificial selection has not yet been reported. Here we show that the circadian clock of cultivated tomato (Solanum lycopersicum) has slowed during domestication. Allelic variation of the tomato homolog of the Arabidopsis gene EID1 is responsible for a phase delay. Notably, the genomic region harboring EID1 shows signatures of a selective sweep. We find that the EID1 allele in cultivated tomatoes enhances plant performance specifically under long day photoperiods, suggesting that humans selected slower circadian rhythms to adapt the cultivated species to the long summer days it encountered as it was moved away from the equator.

A System for Dosage-Based Functional Genomics in Poplar

Altering gene dosage through variation in gene copy number is a powerful approach to addressing questions regarding gene regulation, quantitative trait loci, and heterosis, but one that is not easily applied to sexually transmitted species. Elite poplar (Populus spp) varieties are created through interspecific hybridization, followed by clonal propagation. Altered gene dosage relationships are believed to contribute to hybrid performance. Clonal propagation allows for replication and maintenance of meiotically unstable ploidy or structural variants and provides an alternative approach to investigating gene dosage effects not possible in sexually propagated species. Here, we built a genome-wide structural variation system for dosage-based functional genomics and breeding of poplar. We pollinated Populus deltoides with gamma-irradiated Populus nigra pollen to produce >500 F1 seedlings containing dosage lesions in the form of deletions and insertions of chromosomal segments (indel mutations). Using high-precision dosage analysis, we detected indel mutations in ∼55% of the progeny. These indels varied in length, position, and number per individual, cumulatively tiling >99% of the genome, with an average of 10 indels per gene. Combined with future phenotype and transcriptome data, this population will provide an excellent resource for creating and characterizing dosage-based variation in poplar, including the contribution of dosage to quantitative traits and heterosis.

Efficient Genome-Wide Detection and Cataloging of EMS-Induced Mutations Using Exome Capture and Next-Generation Sequencing

Chemical mutagenesis efficiently generates phenotypic variation in otherwise homogeneous genetic backgrounds, enabling functional analysis of genes. Advances in mutation detection have brought the utility of induced mutant populations on par with those produced by insertional mutagenesis, but systematic cataloguing of mutations would further increase their utility. We examined the suitability of multiplexed global exome capture and sequencing coupled with custom-developed bioinformatics tools to identify mutations in well-characterized mutant populations of rice (Oryza sativa) and wheat (Triticum aestivum). In rice, we identified ∼18,000 induced mutations from 72 independent M2 individuals. Functional evaluation indicated the recovery of potentially deleterious mutations for >2600 genes. We further observed that specific sequence and cytosine methylation patterns surrounding the targeted guanine residues strongly affect their probability to be alkylated by ethyl methanesulfonate. Application of these methods to six independent M2 lines of tetraploid wheat demonstrated that our bioinformatics pipeline is applicable to polyploids. In conclusion, we provide a method for developing large-scale induced mutation resources with relatively small investments that is applicable to resource-poor organisms. Furthermore, our results demonstrate that large libraries of sequenced mutations can be readily generated, providing enhanced opportunities to study gene function and assess the effect of sequence and chromatin context on mutations.

TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline

Genotyping by sequencing (GBS) is a next generation sequencing based method that takes advantage of reduced representation to enable high throughput genotyping of large numbers of individuals at a large number of SNP markers. The relatively straightforward, robust, and cost-effective GBS protocol is currently being applied in numerous species by a large number of researchers. Herein we describe a bioinformatics pipeline, tassel-gbs, designed for the efficient processing of raw GBS sequence data into SNP genotypes. The tassel-gbs pipeline successfully fulfills the following key design criteria: (1) Ability to run on the modest computing resources that are typically available to small breeding or ecological research programs, including desktop or laptop machines with only 8–16 GB of RAM, (2) Scalability from small to extremely large studies, where hundreds of thousands or even millions of SNPs can be scored in up to 100,000 individuals (e.g., for large breeding programs or genetic surveys), and (3) Applicability in an accelerated breeding context, requiring rapid turnover from tissue collection to genotypes. Although a reference genome is required, the pipeline can also be run with an unfinished “pseudo-reference” consisting of numerous contigs. We describe the tassel-gbs pipeline in detail and benchmark it based upon a large scale, species wide analysis in maize (Zea mays), where the average error rate was reduced to 0.0042 through application of population genetic-based SNP filters. Overall, the GBS assay and the tassel-gbs pipeline provide robust tools for studying genomic diversity.

The BOY NAMED SUE quantitative trait locus confers increased meiotic stability to an adapted natural allopolyploid of Arabidopsis

Whole-genome duplication resulting from polyploidy is ubiquitous in the evolutionary history of plant species. Yet, polyploids must overcome the meiotic challenge of pairing, recombining, and segregating more than two sets of chromosomes. Using genomic sequencing of synthetic and natural allopolyploids of Arabidopsis thaliana and Arabidopsis arenosa, we determined that dosage variation and chromosomal translocations consistent with homoeologous pairing were more frequent in the synthetic allopolyploids. To test the role of structural chromosomal differentiation versus genetic regulation of meiotic pairing, we performed sequenced-based, high-density genetic mapping in F2 hybrids between synthetic and natural lines. This F2 population displayed frequent dosage variation and deleterious homoeologous recombination. The genetic map derived from this population provided no indication of structural evolution of the genome of the natural allopolyploid Arabidopsis suecica, compared with its predicted parents. The F2 population displayed variation in meiotic regularity and pollen viability that correlated with a single quantitative trait locus, which we named BOY NAMED SUE, and whose beneficial allele was contributed by A. suecica. This demonstrates that an additive, gain-of-function allele contributes to meiotic stability and fertility in a recently established allopolyploid and provides an Arabidopsis system to decipher evolutionary and molecular mechanisms of meiotic regularity in polyploids.

Application of genotyping-by-sequencing on semiconductor sequencing platforms: a comparison of genetic and reference-based marker ordering in barley

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

Application of genotyping-by-sequencing on semiconductor sequencing platforms: a comparison of genetic and reference-based marker ordering in barley

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

A combinatorial approach to the restriction of a mouse genome

BACKGROUND

A fragmentation of genomic DNA by restriction digestion is a popular step in many applications. Usually attention is paid to the expected average size of the DNA fragments. Another important parameter, randomness of restriction, is regularly implied but rarely verified. This parameter is crucial to the expectation, that either all fragments made by restriction will be suitable for the method of choice, or only a fraction of those will be effectively used by the method. If only a fraction of the fragments are used, we often should know whether the used fragments are representative of the whole genome. With a modern knowledge of mouse, human and many other genomes, frequencies and distributions of restriction sites and sizes of corresponding DNA fragments can be analyzed in silico. In this manuscript, the mouse genome was systematically scanned for frequencies of complementary 4-base long palindromes.

FINDINGS AND CONCLUSIONS

The study revealed substantial heterogeneity in distribution of those sites genome-wide. Only few palindromes showed close to random pattern of distribution. Overall, the distribution of frequencies for most palindromes is much wider than expected by random occurrence. In practical terms, accessibility of genome upon restriction can be improved by a selective combination of restrictases using a few combinatorial rules. It is recommended to mix at least 3 restrictases, their recognition sequences (palindrome) should be the least similar to each other. Principles of the optimization and optimal combinations of restrictases are provided.

A quantitative genetic basis for leaf morphology in a set of precisely defined tomato introgression lines

Introgression lines (ILs), in which genetic material from wild tomato species is introgressed into a domesticated background, have been used extensively in tomato (Solanum lycopersicum) improvement. Here, we genotype an IL population derived from the wild desert tomato Solanum pennellii at ultrahigh density, providing the exact gene content harbored by each line. To take advantage of this information, we determine IL phenotypes for a suite of vegetative traits, ranging from leaf complexity, shape, and size to cellular traits, such as stomatal density and epidermal cell phenotypes. Elliptical Fourier descriptors on leaflet outlines provide a global analysis of highly heritable, intricate aspects of leaf morphology. We also demonstrate constraints between leaflet size and leaf complexity, pavement cell size, and stomatal density and show independent segregation of traits previously assumed to be genetically coregulated. Meta-analysis of previously measured traits in the ILs shows an unexpected relationship between leaf morphology and fruit sugar levels, which RNA-Seq data suggest may be attributable to genetically coregulated changes in fruit morphology or the impact of leaf shape on photosynthesis. Together, our results both improve upon the utility of an important genetic resource and attest to a complex, genetic basis for differences in leaf morphology between natural populations.

Hierarchical nuclear and cytoplasmic genetic architectures for plant growth and defense within Arabidopsis

To understand how genetic architecture translates between phenotypic levels, we mapped the genetic architecture of growth and defense within the Arabidopsis thaliana Kas × Tsu recombinant inbred line population. We measured plant growth using traditional size measurements and size-corrected growth rates. This population contains genetic variation in both the nuclear and cytoplasmic genomes, allowing us to separate their contributions. The cytoplasmic genome regulated a significant variance in growth but not defense, which was due to cytonuclear epistasis. Furthermore, growth adhered to an infinitesimal model of genetic architecture, while defense metabolism was more of a moderate-effect model. We found a lack of concordance between quantitative trait loci (QTL) regulating defense and those regulating growth. Given the published evidence proving the link between glucosinolates and growth, this is likely a false negative result caused by the limited population size. This size limitation creates an inability to test the entire potential genetic landscape possible between these two parents. We uncovered a significant effect of glucosinolates on growth once we accounted for allelic differences in growth QTLs. Therefore, other growth QTLs can mask the effects of defense upon growth. Investigating direct links across phenotypic hierarchies is fraught with difficulty; we identify issues complicating this analysis.

Identification of SNPs in closely related Temperate Japonica rice cultivars using restriction enzyme-phased sequencing

Very low polymorphism in the germplasm typically used by breeding programs poses a significant bottleneck with regards to molecular breeding and the exploitation of breeding materials for quantitative trait analyses. California rice cultivars, derived from a very small base of temperate japonica germplasm and having a relatively brief breeding history, are a good example. In this study, we employed a reduced representation sequencing approach called Restriction Enzyme Site Comparative Analysis (RESCAN) to simultaneously identify and genotype single nucleotide polymorphisms (SNPs) in forty-five rice cultivars representing the majority of the 100 year-old breeding history in California. Over 20,000 putative SNPs were detected relative to the Nipponbare reference genome which enabled the identification and analysis of inheritance of pedigree haplotypes. Haplotype blocks distinguishing modern California cultivars from each other and from the ancestral short grain temperate japonica cultivars were easily identified. Reduced representation sequencing methods such as RESCAN are a valuable alternative to SNP chip genotyping and low coverage whole genome sequencing.

Molecular analysis of point mutations in a barley genome exposed to MNU and gamma rays

We present studies aimed at determining the types and frequencies of mutations induced in the barley genome after treatment with chemical (N-methyl-N-nitrosourea, MNU) and physical (gamma rays) mutagens. We created M(2) populations of a doubled haploid line and used them for the analysis of mutations in targeted DNA sequences and over an entire barley genome using TILLING (Targeting Induced Local Lesions in Genomes) and AFLP (Amplified Fragment Length Polymorphism) technique, respectively. Based on the TILLING analysis of the total DNA sequence of 4,537,117bp in the MNU population, the average mutation density was estimated as 1/504kb. Only one nucleotide change was found after an analysis of 3,207,444bp derived from the highest dose of gamma rays applied. MNU was clearly a more efficient mutagen than gamma rays in inducing point mutations in barley. The majority (63.6%) of the MNU-induced nucleotide changes were transitions, with a similar number of G>A and C>T substitutions. The similar share of G>A and C>T transitions indicates a lack of bias in the repair of O(6)-methylguanine lesions between DNA strands. There was, however, a strong specificity of the nucleotide surrounding the O(6)-meG at the -1 position. Purines formed 81% of nucleotides observed at the -1 site. Scanning the barley genome with AFLP markers revealed ca. a three times higher level of AFLP polymorphism in MNU-treated as compared to the gamma-irradiated population. In order to check whether AFLP markers can really scan the whole barley genome for mutagen-induced polymorphism, 114 different AFLP products, were cloned and sequenced. 94% of bands were heterogenic, with some bands containing up to 8 different amplicons. The polymorphic AFLP products were characterised in terms of their similarity to the records deposited in a GenBank database. The types of sequences present in the polymorphic bands reflected the organisation of the barley genome.

Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping

Quantitative trait loci (QTL) mapping is a powerful tool for investigating the genetic basis of natural variation. QTL can be mapped using a number of different population designs, but recombinant inbred lines (RILs) are among the most effective. Unfortunately, homozygous RIL populations are time consuming to construct, typically requiring at least six generations of selfing starting from a heterozygous F(1). Haploid plants produced from an F(1) combine the two parental genomes and have only one allele at every locus. Converting these sterile haploids into fertile diploids (termed "doubled haploids," DHs) produces immortal homozygous lines in only two steps. Here we describe a unique technique for rapidly creating recombinant doubled haploid populations in Arabidopsis thaliana: centromere-mediated genome elimination. We generated a population of 238 doubled haploid lines that combine two parental genomes and genotyped them by reduced representation Illumina sequencing. The recombination rate and parental allele frequencies in our population are similar to those found in existing RIL sets. We phenotyped this population for traits related to flowering time and for petiole length and successfully mapped QTL controlling each trait. Our work demonstrates that doubled haploid populations offer a rapid, easy alternative to RILs for Arabidopsis genetic analysis.

Exploring the shallow end; estimating information content in transcriptomics studies

Transcriptomics is a major platform to study organismal biology. The advent of new parallel sequencing technologies has opened up a new avenue of transcriptomics with ever deeper and deeper sequencing to identify and quantify each and every transcript in a sample. However, this may not be the best usage of the parallel sequencing technology for all transcriptomics experiments. I utilized the Shannon Entropy approach to estimate the information contained within a transcriptomics experiment and tested the ability of shallow RNAseq to capture the majority of this information. This analysis showed that it was possible to capture nearly all of the network or genomic information present in a variety of transcriptomics experiments using a subset of the most abundant 5000 transcripts or less within any given sample. Thus, it appears that it should be possible and affordable to conduct large scale factorial analysis with a high degree of replication using parallel sequencing technologies.