Fast isogenic mapping-by-sequencing of ethyl methanesulfonate-induced mutant bulks

Agile Genetics: Single gene resolution without the fuss

Gene discovery reveals new biology, expands the utility of marker-assisted selection, and enables targeted mutagenesis. Still, such discoveries can take over a decade. We present a general strategy, "Agile Genetics," that uses nested, structured populations to overcome common limits on gene resolution. Extensive simulation work on realistic genetic architectures shows that, at population sizes of >5000 samples, single gene-resolution can be achieved using bulk segregant pools. At this scale, read depth and technical replication become major drivers of resolution. Emerging enrichment methods to address coverage are on the horizon; we describe one possibility - iterative depth sequencing (ID-seq). In addition, graph-based pangenomics in experimental populations will continue to maximize accuracy and improve interpretation. Based on this merger of agronomic scale with molecular and bioinformatic innovation, we predict a new age of rapid gene discovery.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Reducing PHYTOENE SYNTHASE activity fine-tunes the abundance of a cis-carotene-derived signal that regulates the PIF3/HY5 module and plastid biogenesis

PHYTOENE SYNTHASE (PSY) is a rate-limiting enzyme catalysing the first committed step of carotenoid biosynthesis, and changes in PSY gene expression and/or protein activity alter carotenoid composition and plastid differentiation in plants. Four genetic variants of PSY (psy-4, psy-90, psy-130, and psy-145) were identified using a forward genetics approach that rescued leaf virescence phenotypes and plastid abnormalities displayed by the Arabidopsis CAROTENOID ISOMERASE (CRTISO) mutant ccr2 (carotenoid and chloroplast regulation 2) when grown under a shorter photoperiod. The four non-lethal mutations affected alternative splicing, enzyme-substrate interactions, and PSY:ORANGE multi-enzyme complex binding, constituting the dynamic post-transcriptional fine-tuning of PSY levels and activity without changing localization to the stroma and protothylakoid membranes. psy genetic variants did not alter total xanthophyll or β-carotene accumulation in ccr2, yet they reduced specific acyclic linear cis-carotenes linked to the biosynthesis of a currently unidentified apocarotenoid signal regulating plastid biogenesis, chlorophyll biosynthesis, and photomorphogenic regulation. ccr2 psy variants modulated the PHYTOCHROME-INTERACTING FACTOR 3/ELONGATED HYPOCOTYL 5 (PIF3/HY5) ratio, and displayed a normal prolamellar body formation in etioplasts and chlorophyll accumulation during seedling photomorphogenesis. Thus, suppressing PSY activity and impairing PSY:ORANGE protein interactions revealed how cis-carotene abundance can be fine-tuned through holoenzyme-metabolon interactions to control plastid development.

Synthesis and import of GDP-l-fucose into the Golgi affect plant-water relations

Land plants evolved multiple adaptations to restrict transpiration. However, the underlying molecular mechanisms are not sufficiently understood. We used an ozone-sensitivity forward genetics approach to identify Arabidopsis thaliana mutants impaired in gas exchange regulation. High water loss from detached leaves and impaired decrease of leaf conductance in response to multiple stomata-closing stimuli were identified in a mutant of MURUS1 (MUR1), an enzyme required for GDP-l-fucose biosynthesis. High water loss observed in mur1 was independent from stomatal movements and instead could be linked to metabolic defects. Plants defective in import of GDP-l-Fuc into the Golgi apparatus phenocopied the high water loss of mur1 mutants, linking this phenotype to Golgi-localized fucosylation events. However, impaired fucosylation of xyloglucan, N-linked glycans, and arabinogalactan proteins did not explain the aberrant water loss of mur1 mutants. Partial reversion of mur1 water loss phenotype by borate supplementation and high water loss observed in boron uptake mutants link mur1 gas exchange phenotypes to pleiotropic consequences of l-fucose and boron deficiency, which in turn affect mechanical and morphological properties of stomatal complexes and whole-plant physiology. Our work emphasizes the impact of fucose metabolism and boron uptake on plant-water relations.

Mutation of Arabidopsis SME1 and Sm core assembly improves oxidative stress resilience

Alternative splicing is a key posttranscriptional gene regulatory process, acting in diverse adaptive and basal plant processes. Splicing of precursor-messenger RNA (pre-mRNA) is catalyzed by a dynamic ribonucleoprotein complex, designated the spliceosome. In a suppressor screen, we identified a nonsense mutation in the Smith (Sm) antigen protein SME1 to alleviate photorespiratory H₂O₂-dependent cell death in catalase deficient plants. Similar attenuation of cell death was observed upon chemical inhibition of the spliceosome, suggesting pre-mRNA splicing inhibition to be responsible for the observed cell death alleviation. Furthermore, the sme1-2 mutants showed increased tolerance to the reactive oxygen species inducing herbicide methyl viologen. Both an mRNA-seq and shotgun proteomic analysis in sme1-2 mutants displayed a constitutive molecular stress response, together with extensive alterations in pre-mRNA splicing of transcripts encoding metabolic enzymes and RNA binding proteins, even under unstressed conditions. Using SME1 as a bait to identify protein interactors, we provide experimental evidence for almost 50 homologs of the mammalian spliceosome-associated protein to reside in the Arabidopsis thaliana spliceosome complexes and propose roles in pre-mRNA splicing for four uncharacterized plant proteins. Furthermore, as for sme1-2, a mutant in the Sm core assembly protein ICLN resulted in a decreased sensitivity to methyl viologen. Taken together, these data show that both a perturbed Sm core composition and assembly results in the activation of a defense response and in enhanced resilience to oxidative stress.

Versatile mapping-by-sequencing with Easymap v.2

Mapping-by-sequencing combines Next Generation Sequencing (NGS) with classical genetic mapping by linkage analysis to establish gene-to-phenotype relationships. Although numerous tools have been developed to analyze NGS datasets, only a few are available for mapping-by-sequencing. One such tool is Easymap, a versatile, easy-to-use package that performs automated mapping of point mutations and large DNA insertions. Here, we describe Easymap v.2, which also maps small insertion/deletions (InDels), and includes workflows to perform QTL-seq and variant density mapping analyses. Each mapping workflow can accommodate different experimental designs, including outcrossing and backcrossing, F₂, M₂, and M₃ mapping populations, chemically induced mutation and natural variant mapping, input files containing single-end or paired-end reads of genomic or complementary DNA sequences, and alternative control sample files in FASTQ and VCF formats. Easymap v.2 can also be used as a variant analyzer in the absence of a mapping algorithm and includes a multi-threading option.

EXTRA LARGE G-PROTEIN2 mediates cell death and hyperimmunity in the chitin elicitor receptor kinase 1-4 mutant

Heterotrimeric G-proteins are signal transduction complexes that comprised three subunits, Gα, Gβ, and Gγ, and are involved in many aspects of plant life. The noncanonical Gα subunit EXTRA LARGE G-PROTEIN2 (XLG2) mediates pathogen-associated molecular pattern (PAMP)-induced reactive oxygen species (ROS) generation and immunity downstream of pattern recognition receptors. A mutant of the chitin receptor component CHITIN ELICITOR RECEPTOR KINASE1 (CERK1), cerk1-4, maintains normal chitin signaling capacity but shows excessive cell death upon infection with powdery mildew fungi. We identified XLG2 mutants as suppressors of the cerk1-4 phenotype. Mutations in XLG2 complex partners ARABIDOPSIS Gβ1 (AGB1) and Gγ1 (AGG1) have a partial cerk1-4 suppressor effect. Contrary to its role in PAMP-induced immunity, XLG2-mediated control of ROS production by RESPIRATORY BURST OXIDASE HOMOLOGUE D (RBOHD) is not critical for cerk1-4-associated cell death and hyperimmunity. The cerk1-4 phenotype is also independent of the co-receptor/adapter kinases BRI1-ASSOCIATED RECEPTOR KINASE 1 (BAK1) and SUPPRESSOR OF BIR1 1 (SOBIR1), but requires the E3 ubiquitin ligase PLANT U-BOX 2 (PUB2). XLG2 localizes to both the cell periphery and nucleus, and the cerk1-4 cell death phenotype is mediated by the cell periphery pool of XLG2. Integrity of the XLG2 N-terminal domain, but not its phosphorylation, is essential for correct XLG2 localization and formation of the cerk1-4 phenotype. Our results support a model in which XLG2 acts downstream of an unknown cell surface receptor that activates an NADPH oxidase-independent cell death pathway in Arabidopsis (Arabidopsis thaliana).

inGAP-family: Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering Out Artificial Variants due to Genome Complexities

Accurately identifying DNA polymorphisms can bridge the gap between phenotypes and genotypes and is essential for molecular marker assisted genetic studies. Genome complexities, including large-scale structural variations, bring great challenges to bioinformatic analysis for obtaining high-confidence genomic variants, as sequence differences between non-allelic loci of two or more genomes can be misinterpreted as polymorphisms. It is important to correctly filter out artificial variants to avoid false genotyping or estimation of allele frequencies. Here, we present an efficient and effective framework, inGAP-family, to discover, filter, and visualize DNA polymorphisms and structural variants (SVs) from alignment of short reads. Applying this method to polymorphism detection on real datasets shows that elimination of artificial variants greatly facilitates the precise identification of meiotic recombination points as well as causal mutations in mutant genomes or quantitative trait loci. In addition, inGAP-family provides a user-friendly graphical interface for detecting polymorphisms and SVs, further evaluating predicted variants and identifying mutations related to genotypes. It is accessible at https://sourceforge.net/projects/ingap-family/.

Systematic prediction of EMS-induced mutations in a sorghum mutant population

The precise detection of causal DNA mutations (deoxyribonucleic acid) is very crucial for forward genetic studies. Several sources of errors contribute to false-positive detections by current variant-calling algorithms, which impact associating phenotypes with genotypes. To improve the accuracy of mutation detection, we implemented a binning method for the accurate detection of likely ethyl methanesulfonate (EMS)-induced mutations in a sequenced mutant population. We also implemented a clustering algorithm for detecting likely false negatives with high accuracy. Sorghum bicolor is a very valuable crop species with tremendous potential for uncovering novel gene functions associated with highly desirable agronomical traits. We demonstrate the precision of the described approach in the detection of likely EMS-induced mutations from the publicly available low-cost sequencing of the M₃ generation from 600 sorghum BTx623 mutants. The approach detected 3,274,606 single nucleotide polymorphisms (SNPs), of which 96% (3,141,908) were G/C to A/T DNA substitutions, as expected by EMS-mutagenesis mode of action. We demonstrated the general applicability of the described method and showed a high concordance, 94% (3,074,759) SNPs overlap between SAMtools-based and GATK-based variant-calling algorithms. Our clustering algorithm uncovered evidence for an additional 223,048 likely false-negative shared EMS-induced mutations. The final 3,497,654 SNPs represent an 87% increase in SNPs detected from the previous analysis of the mutant population, with an average of one SNP per 125 kb in the sorghum genome. Annotation of the final SNPs revealed 10,263 high-impact and 136,639 moderate-impact SNPs, including 7217 stop-gained mutations, which averages 12 stop-gained mutations per mutant, and four high- or medium-impact SNPs per sorghum gene. We have implemented a public search database for this new genetic resource of 30,285 distinct sorghum genes containing medium- or high-impact EMS-induced mutations. Seedstock for a select 486 of the 600 described mutants are publicly available in the Germplasm Resources Information Network (GRIN) database.

Broadening the horizon of crop research: a decade of advancements in plant molecular genetics to divulge phenotype governing genes

Advancements in sequencing, genotyping, and computational technologies during the last decade (2011-2020) enabled new forward-genetic approaches, which subdue the impediments of precise gene mapping in varied crops. The modern crop improvement programs rely heavily on two major steps-trait-associated QTL/gene/marker's identification and molecular breeding. Thus, it is vital for basic and translational crop research to identify genomic regions that govern the phenotype of interest. Until the advent of next-generation sequencing, the forward-genetic techniques were laborious and time-consuming. Over the last 10 years, advancements in the area of genome assembly, genotyping, large-scale data analysis, and statistical algorithms have led faster identification of genomic variations regulating the complex agronomic traits and pathogen resistance. In this review, we describe the latest developments in genome sequencing and genotyping along with a comprehensive evaluation of the last 10-year headways in forward-genetic techniques that have shifted the focus of plant research from model plants to diverse crops. We have classified the available molecular genetic methods under bulk-segregant analysis-based (QTL-seq, GradedPool-Seq, QTG-Seq, Exome QTL-seq, and RapMap), target sequence enrichment-based (RenSeq, AgRenSeq, and TACCA), and mutation-based groups (MutMap, NIKS algorithm, MutRenSeq, MutChromSeq), alongside improvements in classical mapping and genome-wide association analyses. Newer methods for outcrossing, heterozygous, and polyploid plant genetics have also been discussed. The use of k-mers has enriched the nature of genetic variants which can be utilized to identify the phenotype-causing genes, independent of reference genomes. We envisage that the recent methods discussed herein will expand the repertoire of useful alleles and help in developing high-yielding and climate-resilient crops.

The Sorghum (Sorghum bicolor) Brown Midrib 30 Gene Encodes a Chalcone Isomerase Required for Cell Wall Lignification

In sorghum (Sorghum bicolor) and other C4 grasses, brown midrib (bmr) mutants have long been associated with plants impaired in their ability to synthesize lignin. The brown midrib 30 (Bmr30) gene, identified using a bulk segregant analysis and next-generation sequencing, was determined to encode a chalcone isomerase (CHI). Two independent mutations within this gene confirmed that loss of its function was responsible for the brown leaf midrib phenotype and reduced lignin concentration. Loss of the Bmr30 gene function, as shown by histochemical staining of leaf midrib and stalk sections, resulted in altered cell wall composition. In the bmr30 mutants, CHI activity was drastically reduced, and the accumulation of total flavonoids and total anthocyanins was impaired, which is consistent with its function in flavonoid biosynthesis. The level of the flavone lignin monomer tricin was reduced 20-fold in the stem relative to wild type, and to undetectable levels in the leaf tissue of the mutants. The bmr30 mutant, therefore, harbors a mutation in a phenylpropanoid biosynthetic gene that is key to the interconnection between flavonoids and monolignols, both of which are utilized for lignin synthesis in the grasses.

Separable regulation of POW1 in grain size and leaf angle development in rice

Leaf angle is one of the key factors that determines rice plant architecture. However, the improvement of leaf angle erectness is often accompanied by unfavourable changes in other traits, especially grain size reduction. In this study, we identified the pow1 (put on weight 1) mutant that leads to increased grain size and leaf angle, typical brassinosteroid (BR)-related phenotypes caused by excessive cell proliferation and cell expansion. We show that modulation of the BR biosynthesis genes OsDWARF4 (D4) and D11 and the BR signalling gene D61 could rescue the phenotype of leaf angle but not grain size in the pow1 mutant. We further demonstrated that POW1 functions in grain size regulation by repressing the transactivation activity of the interacting protein TAF2, a highly conserved member of the TFIID transcription initiation complex. Down-regulation of TAF2 rescued the enlarged grain size of pow1 but had little effect on the increased leaf angle phenotype of the mutant. The separable functions of the POW1-TAF2 and POW1-BR modules in grain size and leaf angle control provide a promising strategy for designing varieties with compact plant architecture and increased grain size, thus promoting high-yield breeding in rice.

Sorghum genetic, genomic, and breeding resources

Sorghum research has entered an exciting and fruitful era due to the genetic, genomic, and breeding resources that are now available to researchers and plant breeders. As the world faces the challenges of a rising population and a changing global climate, new agricultural solutions will need to be developed to address the food and fiber needs of the future. To that end, sorghum will be an invaluable crop species as it is a stress-resistant C₄ plant that is well adapted for semi-arid and arid regions. Sorghum has already remained as a staple food crop in many parts of Africa and Asia and is critically important for animal feed and niche culinary applications in other regions, such as the United States. In addition, sorghum has begun to be developed into a promising feedstock for forage and bioenergy production. Due to this increasing demand for sorghum and its potential to address these needs, the continuous development of powerful community resources is required. These resources include vast collections of sorghum germplasm, high-quality reference genome sequences, sorghum association panels for genome-wide association studies of traits involved in food and bioenergy production, mutant populations for rapid discovery of causative genes for phenotypes relevant to sorghum improvement, gene expression atlas, and online databases that integrate all resources and provide the sorghum community with tools that can be used in breeding and genomic studies. Used in tandem, these valuable resources will ensure that the rate, quality, and collaborative potential of ongoing sorghum improvement efforts is able to rival that of other major crops.

One factor, many systems: the floral homeotic protein AGAMOUS and its epigenetic regulatory mechanisms

Tissue-specific transcription factors allow cells to specify new fates by exerting control over gene regulatory networks and the epigenetic landscape of a cell. However, our knowledge of the molecular mechanisms underlying cell fate decisions is limited. In Arabidopsis, the MADS-box transcription factor AGAMOUS (AG) plays a central role in regulating reproductive organ identity and meristem determinacy during flower development. During the vegetative phase, AG transcription is repressed by Polycomb complexes and intronic noncoding RNA. Once AG is transcribed in a spatiotemporally regulated manner during the reproductive phase, AG functions with chromatin regulators to change the chromatin structure at key target gene loci. The concerted actions of AG and the transcription factors functioning downstream of AG recruit general transcription machinery for proper cell fate decision. In this review, we describe progress in AG research that has provided important insights into the regulatory and epigenetic mechanisms underlying cell fate determination in plants.

Mapping of a major QTL controlling plant height using a high-density genetic map and QTL-seq methods based on whole-genome resequencing in Brassica napus

Plant height is a crucial element related to plant architecture that influences the seed yield of oilseed rape (Brassica napus L.). In this study, we isolated a natural B. napus mutant, namely a semi-dwarf mutant (sdw-e), which exhibits a 30% reduction in plant height compared with Zhongshuang 11-HP (ZS11-HP). Quantitative trait locus sequencing (QTL-seq) was conducted using two extreme DNA bulks in F₂ populations in Wuchang-2017 derived from ZS11-HP × sdw-e to identify QTLs associated with plant height. The result suggested that two QTL intervals were located on chromosome A10. The F₂ population consisting of 200 individuals in Yangluo-2018 derived from ZS11-HP × sdw-e was used to construct a high-density linkage map using whole-genome resequencing. The high-density linkage map harbored 4323 bin markers and covered a total distance of 2026.52 cM with an average marker interval of 0.47 cM. The major QTL for plant height named qPHA10 was identified on linkage group A10 by interval mapping and composite interval mapping methods. The major QTL qPHA10 was highly consistent with the QTL-seq results. And then, we integrated the variation sites and expression levels of genes in the major QTL interval to predict the candidate genes. Thus, the identified QTL and candidate genes could be used in marker-assisted selection for B. napus breeding in the future.

Easymap: A User-Friendly Software Package for Rapid Mapping-by-Sequencing of Point Mutations and Large Insertions

Mapping-by-sequencing strategies combine next-generation sequencing (NGS) with classical linkage analysis, allowing rapid identification of the causal mutations of the phenotypes exhibited by mutants isolated in a genetic screen. Computer programs that analyze NGS data obtained from a mapping population of individuals derived from a mutant of interest to identify a causal mutation are available; however, the installation and usage of such programs requires bioinformatic skills, modifying or combining pieces of existing software, or purchasing licenses. To ease this process, we developed Easymap, an open-source program that simplifies the data analysis workflows from raw NGS reads to candidate mutations. Easymap can perform bulked segregant mapping of point mutations induced by ethyl methanesulfonate (EMS) with DNA-seq or RNA-seq datasets, as well as tagged-sequence mapping for large insertions, such as transposons or T-DNAs. The mapping analyses implemented in Easymap have been validated with experimental and simulated datasets from different plant and animal model species. Easymap was designed to be accessible to all users regardless of their bioinformatics skills by implementing a user-friendly graphical interface, a simple universal installation script, and detailed mapping reports, including informative images and complementary data for assessment of the mapping results. Easymap is available at http://genetics.edu.umh.es/resources/easymap; its Quickstart Installation Guide details the recommended procedure for installation.

Beyond flowering time: diverse roles of an APETALA2-like transcription factor in shoot architecture and perennial traits

Polycarpic perennials maintain vegetative growth after flowering. PERPETUAL FLOWERING 1 (PEP1), the orthologue of FLOWERING LOCUS C (FLC) in Arabis alpina regulates flowering and contributes to polycarpy in a vernalisation-dependent pathway. pep1 mutants do not require vernalisation to flower and have reduced return to vegetative growth as all of their axillary branches become reproductive. To identify additional genes that regulate flowering and contribute to perennial traits we performed an enhancer screen of pep1. Using mapping-by-sequencing, we cloned a mutant (enhancer of pep1-055, eop055), performed transcriptome analysis and physiologically characterised the role it plays on perennial traits in an introgression line carrying the eop055 mutation and a functional PEP1 wild-type allele. eop055 flowers earlier than pep1 and carries a lesion in the A. alpina orthologue of the APETALA2 (AP2)-like gene, TARGET OF EAT2 (AaTOE2). AaTOE2 is a floral repressor and acts upstream of SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE 5 (AaSPL5). In the wild-type background, which requires cold treatment to flower, AaTOE2 regulates the age-dependent response to vernalisation. In addition, AaTOE2 ensures the maintenance of vegetative growth by delaying axillary meristem initiation and repressing flowering of axillary buds before and during cold exposure. We conclude that AaTOE2 is instrumental in fine tuning different developmental traits in the perennial life cycle of A. alpina.

An extreme-phenotype genome-wide association study identifies candidate cannabinoid pathway genes in Cannabis

Cannabis produces a class of isoprenylated resorcinyl polyketides known as cannabinoids, a subset of which are medically important and exclusive to this plant. The cannabinoid alkyl group is a critical structural feature that governs therapeutic activity. Genetic enhancement of the alkyl side-chain could lead to the development of novel chemical phenotypes (chemotypes) for pharmaceutical end-use. However, the genetic determinants underlying in planta variation of cannabinoid alkyl side-chain length remain uncharacterised. Using a diversity panel derived from the Ecofibre Cannabis germplasm collection, an extreme-phenotype genome-wide association study (XP-GWAS) was used to enrich for alkyl cannabinoid polymorphic regions. Resequencing of chemotypically extreme pools revealed a known cannabinoid synthesis pathway locus as well as a series of chemotype-associated genomic regions. One of these regions contained a candidate gene encoding a β-keto acyl carrier protein (ACP) reductase (BKR) putatively associated with polyketide fatty acid starter unit synthesis and alkyl side-chain length. Association analysis revealed twenty-two polymorphic variants spanning the length of this gene, including two nonsynonymous substitutions. The success of this first reported application of XP-GWAS for an obligate outcrossing and highly heterozygote plant genus suggests that this approach may have generic application for other plant species.

Discovery of genomic regions and candidate genes for grain weight employing next generation sequencing based QTL-seq approach in rice (Oryza sativa L.)

Rice (Oryza sativa L.) yield enhancement is one of the prime objectives of plant breeders. Elucidation of the inheritance of grain weight, a key yield component trait, is of paramount importance for raising the yield thresholds in rice. In the present investigation, we employed Next-Generation Sequencing based QTL-seq approach to identify major genomic regions associated with grain weight using mapping populations derived from a cross between BPT5204 and MTU3626. QTL-seq analysis identified three grain weight quantitative trait loci (QTL) viz., qGW1 (35-40 Mb), qGW7 (10-18 Mb), and qGW8 (2-5 Mb) on chromosomes 1, 7 and 8, respectively and all are found to be novel. Further, qGW8 was confirmed through conventional QTL mapping in F₂, F₃ and BC₁F₂ populations and found to explain the phenotypic variance of 17.88%, 16.70% and 15.00%, respectively, indicating a major QTL for grain weight. Based on previous reports, two candidate genes in the qGW8 QTL were predicted i.e., LOC_Os08g01490 (Cytochrome P450), and LOC_Os08g01680 (WD domain, G-beta repeat domain containing protein) and through in silico analysis they were found to be highly expressed in reproductive organs during different stages of grain development. Here, we have demonstrated that QTL-seq is one of the rapid approaches to uncover novel QTLs controlling complex traits. The candidate genes identified in the present study undoubtedly enhance our understanding of the mechanism and inheritance of the grain weight. These candidate genes can be exploited for yield enhancement after confirmation through complementary studies.

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.

The domesticated transposase ALP2 mediates formation of a novel Polycomb protein complex by direct interaction with MSI1, a core subunit of Polycomb Repressive Complex 2 (PRC2)

A large fraction of plant genomes is composed of transposable elements (TE), which provide a potential source of novel genes through “domestication”–the process whereby the proteins encoded by TE diverge in sequence, lose their ability to catalyse transposition and instead acquire novel functions for their hosts. In Arabidopsis, ANTAGONIST OF LIKE HETEROCHROMATIN PROTEIN 1 (ALP1) arose by domestication of the nuclease component of Harbinger class TE and acquired a new function as a component of POLYCOMB REPRESSIVE COMPLEX 2 (PRC2), a histone H3K27me3 methyltransferase involved in regulation of host genes and in some cases TE. It was not clear how ALP1 associated with PRC2, nor what the functional consequence was. Here, we identify ALP2 genetically as a suppressor of Polycomb-group (PcG) mutant phenotypes and show that it arose from the second, DNA binding component of Harbinger transposases. Molecular analysis of PcG compromised backgrounds reveals that ALP genes oppose silencing and H3K27me3 deposition at key PcG target genes. Proteomic analysis reveals that ALP1 and ALP2 are components of a variant PRC2 complex that contains the four core components but lacks plant-specific accessory components such as the H3K27me3 reader LIKE HETEROCHROMATION PROTEIN 1 (LHP1). We show that the N-terminus of ALP2 interacts directly with ALP1, whereas the C-terminus of ALP2 interacts with MULTICOPY SUPPRESSOR OF IRA1 (MSI1), a core component of PRC2. Proteomic analysis reveals that in alp2 mutant backgrounds ALP1 protein no longer associates with PRC2, consistent with a role for ALP2 in recruitment of ALP1. We suggest that the propensity of Harbinger TE to insert in gene-rich regions of the genome, together with the modular two component nature of their transposases, has predisposed them for domestication and incorporation into chromatin modifying complexes.

A large part of the genomes of plants and animals consists of transposable elements (TE), which are usually considered as selfish or parasitic as they encode proteins (transposases) which promote TE proliferation but not functions useful for their hosts. As a result, hosts have evolved ways of reducing TE proliferation, usually by modifying the DNA or chromatin of TE so that their transposases are no longer produced. Once the TE are inactivated they can no longer proliferate and over time they accumulate mutations and can evolve new functions, often beneficial for their hosts. This process is known as domestication and is increasingly recognised as a potent source of evolutionary novelty. For example, the CRISPR/Cas system that has provided the basis for a revolution in genetic engineering (“genome editing”) has evolved via domestication of transposons in bacteria. We have identified the ALP proteins, two domesticated transposases which function as components of an enzyme complex (PRC2) involved in modifying chromatin and regulating host gene activity in plants. Here we show how ALPs contact PRC2 and direct formation of a novel complex that lacks several of the usual components. The ALPs and related proteins will provide valuable tools for manipulating plant chromatin.

Next-generation forward genetic screens: using simulated data to improve the design of mapping-by-sequencing experiments in Arabidopsis

Forward genetic screens have successfully identified many genes and continue to be powerful tools for dissecting biological processes in Arabidopsis and other model species. Next-generation sequencing technologies have revolutionized the time-consuming process of identifying the mutations that cause a phenotype of interest. However, due to the cost of such mapping-by-sequencing experiments, special attention should be paid to experimental design and technical decisions so that the read data allows to map the desired mutation. Here, we simulated different mapping-by-sequencing scenarios. We first evaluated which short-read technology was best suited for analyzing gene-rich genomic regions in Arabidopsis and determined the minimum sequencing depth required to confidently call single nucleotide variants. We also designed ways to discriminate mutagenesis-induced mutations from background Single Nucleotide Polymorphisms in mutants isolated in Arabidopsis non-reference lines. In addition, we simulated bulked segregant mapping populations for identifying point mutations and monitored how the size of the mapping population and the sequencing depth affect mapping precision. Finally, we provide the computational basis of a protocol that we already used to map T-DNA insertions with paired-end Illumina-like reads, using very low sequencing depths and pooling several mutants together; this approach can also be used with single-end reads as well as to map any other insertional mutagen. All these simulations proved useful for designing experiments that allowed us to map several mutations in Arabidopsis.

Mutagenesis of a Quintuple Mutant Impaired in Environmental Responses Reveals Roles for CHROMATIN REMODELING4 in the Arabidopsis Floral Transition

Several pathways conferring environmental flowering responses in Arabidopsis (Arabidopsis thaliana) converge on developmental processes that mediate the floral transition in the shoot apical meristem. Many characterized mutations disrupt these environmental responses, but downstream developmental processes have been more refractory to mutagenesis. Here, we constructed a quintuple mutant impaired in several environmental pathways and showed that it possesses severely reduced flowering responses to changes in photoperiod and ambient temperature. RNA-sequencing (RNA-seq) analysis of the quintuple mutant showed that the expression of genes encoding gibberellin biosynthesis enzymes and transcription factors involved in the age pathway correlates with flowering. Mutagenesis of the quintuple mutant generated two late-flowering mutants, quintuple ems1 (qem1) and qem2 The mutated genes were identified by isogenic mapping and transgenic complementation. The qem1 mutant is an allele of the gibberellin 20-oxidase gene ga20ox2, confirming the importance of gibberellin for flowering in the absence of environmental responses. By contrast, qem2 is impaired in CHROMATIN REMODELING4 (CHR4), which has not been genetically implicated in floral induction. Using co-immunoprecipitation, RNA-seq, and chromatin immunoprecipitation sequencing, we show that CHR4 interacts with transcription factors involved in floral meristem identity and affects the expression of key floral regulators. Therefore, CHR4 mediates the response to endogenous flowering pathways in the inflorescence meristem to promote floral identity.

A cis-carotene derived apocarotenoid regulates etioplast and chloroplast development

Carotenoids are a core plastid component and yet their regulatory function during plastid biogenesis remains enigmatic. A unique carotenoid biosynthesis mutant, carotenoid chloroplast regulation 2 (ccr2), that has no prolamellar body (PLB) and normal PROTOCHLOROPHYLLIDE OXIDOREDUCTASE (POR) levels, was used to demonstrate a regulatory function for carotenoids and their derivatives under varied dark-light regimes. A forward genetics approach revealed how an epistatic interaction between a ζ-carotene isomerase mutant (ziso-155) and ccr2 blocked the biosynthesis of specific cis-carotenes and restored PLB formation in etioplasts. We attributed this to a novel apocarotenoid retrograde signal, as chemical inhibition of carotenoid cleavage dioxygenase activity restored PLB formation in ccr2 etioplasts during skotomorphogenesis. The apocarotenoid acted in parallel to the repressor of photomorphogenesis, DEETIOLATED1 (DET1), to transcriptionally regulate PROTOCHLOROPHYLLIDE OXIDOREDUCTASE (POR), PHYTOCHROME INTERACTING FACTOR3 (PIF3) and ELONGATED HYPOCOTYL5 (HY5). The unknown apocarotenoid signal restored POR protein levels and PLB formation in det1, thereby controlling plastid development.

Next-generation sequencing identified genomic region and diagnostic markers for resistance to bacterial wilt on chromosome B02 in peanut (Arachis hypogaea L.)

Bacterial wilt, caused by Ralstonia solanacearum, is a devastating disease affecting over 350 plant species. A few peanut cultivars were found to possess stable and durable bacterial wilt resistance (BWR). Genomics-assisted breeding can accelerate the process of developing resistant cultivars by using diagnostic markers. Here, we deployed sequencing-based trait mapping approach, QTL-seq, to discover genomic regions, candidate genes and diagnostic markers for BWR in a recombination inbred line population (195 progenies) of peanut. The QTL-seq analysis identified one candidate genomic region on chromosome B02 significantly associated with BWR. Mapping of newly developed single nucleotide polymorphism (SNP) markers narrowed down the region to 2.07 Mb and confirmed its major effects and stable expressions across three environments. This candidate genomic region had 49 nonsynonymous SNPs affecting 19 putative candidate genes including seven putative resistance genes (R-genes). Two diagnostic markers were successfully validated in diverse breeding lines and cultivars and could be deployed in genomics-assisted breeding of varieties with enhanced BWR.

Combining next-generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M₂ ) family derived from a heterozygous (M₁ ) parent were identified using whole-genome shotgun (WGS) sequencing of a small number of (M₂ ) individuals with mutant and wild-type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self-pollinated (M₃ ) offspring. Finally, we filtered for segregating EMS-induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M₂ ) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non-synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.

artMAP: A user-friendly tool for mapping ethyl methanesulfonate-induced mutations in Arabidopsis

Mapping-by-sequencing is a rapid method for identifying both natural as well as induced variations in the genome. However, it requires extensive bioinformatics expertise along with the computational infrastructure to analyze the sequencing data and these requirements have limited its widespread adoption. In the current study, we develop an easy to use tool, artMAP, to discover ethyl methanesulfonate (EMS) induced mutations in the Arabidopsis genome. The artMAP pipeline consists of well-established tools including TrimGalore, BWA, BEDTools, SAMtools, and SnpEff which were integrated in a Docker container. artMAP provides a graphical user interface and can be run on a regular laptop and desktop, thereby limiting the bioinformatics expertise required. artMAP can process input sequencing files generated from single or paired-end sequencing. The results of the analysis are presented in interactive graphs which display the annotation details of each mutation. Due to its ease of use, artMAP makes the identification of EMS-induced mutations in Arabidopsis possible with only a few mouse click. The source code of artMAP is available on Github (https://github.com/RihaLab/artMAP).

Comparative Subsequence Sets Analysis (CoSSA) is a robust approach to identify haplotype specific SNPs; mapping and pedigree analysis of a potato wart disease resistance gene Sen3

BACKGROUND

Standard strategies to identify genomic regions involved in a specific trait variation are often limited by time and resource consuming genotyping methods. Other limiting pre-requisites are the phenotyping of large segregating populations or of diversity panels and the availability and quality of a closely related reference genome. To overcome these limitations, we designed efficient Comparative Subsequence Sets Analysis (CoSSA) workflows to identify haplotype specific SNPs linked to a trait of interest from Whole Genome Sequencing data.

RESULTS

As a model, we used the resistance to Synchytrium endobioticum pathotypes 2, 6 and 18 that co-segregated in a tetraploid full sib population. Genomic DNA from both parents, pedigree genotypes, unrelated potato varieties lacking the wart resistance traits and pools of resistant and susceptible siblings were sequenced. Set algebra and depth filtering of subsequences (k-mers) were used to delete unlinked and common SNPs and to enrich for SNPs from the haplotype(s) harboring the resistance gene(s). Using CoSSA, we identified a major and a minor effect locus. Upon comparison to the reference genome, it was inferred that the major resistance locus, referred to as Sen3, was located on the north arm of chromosome 11 between 1,259,552 and 1,519,485 bp. Furthermore, we could anchor the unanchored superscaffold DMB734 from the potato reference genome to a synthenous interval. CoSSA was also successful in identifying Sen3 in a reference genome independent way thanks to the de novo assembly of paired end reads matching haplotype specific k-mers. The de novo assembly provided more R haplotype specific polymorphisms than the reference genome corresponding region. CoSSA also offers possibilities for pedigree analysis. The origin of Sen3 was traced back until Ora. Finally, the diagnostic power of the haplotype specific markers was shown using a panel of 56 tetraploid varieties.

CONCLUSIONS

CoSSA is an efficient, robust and versatile set of workflows for the genetic analysis of a trait of interest using WGS data. Because the WGS data are used without intermediate reads mapping, CoSSA does not require the use of a reference genome. This approach allowed the identification of Sen3 and the design of haplotype specific, diagnostic markers.

Whole Genome Resequencing from Bulked Populations as a Rapid QTL and Gene Identification Method in Rice

Most Quantitative Trait Loci (QTL) and gene isolation approaches, such as positional- or map-based cloning, are time-consuming and low-throughput methods. Understanding and detecting the genetic material that controls a phenotype is a key means to functionally analyzing genes as well as to enhance crop agronomic traits. In this regard, high-throughput technologies have great prospects for changing the paradigms of DNA marker revealing, genotyping, and for discovering crop genetics and genomic study. Bulk segregant analysis, based on whole genome resequencing approaches, permits the rapid isolation of the genes or QTL responsible for the causative mutation of the phenotypes. MutMap, MutMap Gap, MutMap+, modified MutMap, and QTL-seq methods are among those approaches that have been confirmed to be fruitful gene mapping approaches for crop plants, such as rice, irrespective of whether the characters are determined by polygenes. As a result, in the present study we reviewed the progress made by all these methods to identify QTL or genes in rice.

Consensus Coexpression Network Analysis Identifies Key Regulators of Flower and Fruit Development in Wild Strawberry

The diploid strawberry, Fragaria vesca, is a developing model system for the economically important Rosaceae family. Strawberry fleshy fruit develops from the floral receptacle and its ripening is nonclimacteric. The external seed configuration of strawberry fruit facilitates the study of seed-to-fruit cross tissue communication, particularly phytohormone biosynthesis and transport. To investigate strawberry fruit development, we previously generated spatial and temporal transcriptome data profiling F. vesca flower and fruit development pre- and postfertilization. In this study, we combined 46 of our existing RNA-seq libraries to generate coexpression networks using the Weighted Gene Co-Expression Network Analysis package in R. We then applied a post-hoc consensus clustering approach and used bootstrapping to demonstrate consensus clustering's ability to produce robust and reproducible clusters. Further, we experimentally tested hypotheses based on the networks, including increased iron transport from the receptacle to the seed postfertilization and characterized a F. vesca floral mutant and its candidate gene. To increase their utility, the networks are presented in a web interface (www.fv.rosaceaefruits.org) for easy exploration and identification of coexpressed genes. Together, the work reported here illustrates ways to generate robust networks optimized for the mining of large transcriptome data sets, thereby providing a useful resource for hypothesis generation and experimental design in strawberry and related Rosaceae fruit crops.

Fine mapping and identification of a novel locus qGL12.2 control grain length in wild rice (Oryza rufipogon Griff.)

A wild rice QTL qGL12.2 for grain length was fine mapped to an 82-kb interval in chromosome 12 containing six candidate genes and none was reported previously. Grain length is an important trait for yield and commercial value in rice. Wild rice seeds have a very slender shape and have many desirable genes that have been lost in cultivated rice during domestication. In this study, we identified a quantitative trait locus, qGL12.2, which controls grain length in wild rice. First, a wild rice chromosome segment substitution line, CSSL41, was selected that has longer glume and grains than does the Oryza sativa indica cultivar, 9311. Next, an F₂ population was constructed from a cross between CSSL41 and 9311. Using the next-generation sequencing combined with bulked-segregant analysis and F₃ recombinants analysis, qGL12.2 was finally fine mapped to an 82-kb interval in chromosome 12. Six candidate genes were found, and no reported grain length genes were found in this interval. Using scanning electron microscopy, we found that CSSL41 cells are significantly longer than those of 9311, but there is no difference in cell widths. These data suggest that qGL12.2 is a novel gene that controls grain cell length in wild rice. Our study provides a new genetic resource for rice breeding and a starting point for functional characterization of the wild rice GL gene.

INCURVATA11 and CUPULIFORMIS2 Are Redundant Genes That Encode Epigenetic Machinery Components in Arabidopsis

All critical developmental and physiological events in a plant's life cycle depend on the proper activation and repression of specific gene sets, and this often involves epigenetic mechanisms. Some Arabidopsis thaliana mutants with disorders of the epigenetic machinery exhibit pleiotropic defects, including incurved leaves and early flowering, due to the ectopic and heterochronic derepression of developmental regulators. Here, we studied one such mutant class, the incurvata11 (icu11) loss-of-function mutants. We have identified ICU11 as the founding member of a small gene family that we have named CUPULIFORMIS (CP). This family is part of the 2-oxoglutarate/Fe(II)-dependent dioxygenase superfamily. ICU11 and its closest paralog, CP2, have unequally redundant functions: although cp2 mutants are phenotypically wild type, icu11 cp2 double mutants skip vegetative development and flower upon germination. This phenotype is reminiscent of loss-of-function mutants of the Polycomb-group genes EMBRYONIC FLOWER1 (EMF1) and EMF2 Double mutants harboring icu11 alleles and loss-of-function alleles of genes encoding components of the epigenetic machinery exhibit synergistic, severe phenotypes, and some are similar to those of emf mutants. Hundreds of genes are misexpressed in icu11 plants, including SEPALLATA3 (SEP3), and derepression of SEP3 causes the leaf phenotype of icu11 ICU11 and CP2 are nucleoplasmic proteins that act as epigenetic repressors through an unknown mechanism involving histone modification, but not DNA methylation.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis

KEY MESSAGE

GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus

To unlock the power of next generation sequencing-based bulked segregant analysis in allele discovery in out-crossing woody species, and to understand the genetic control of the weeping trait, an F1 population from the cross 'Cheal's Weeping' × 'Evereste' was used to create two genomic DNA pools 'weeping' (17 progeny) and 'standard' (16 progeny). Illumina pair-end (2 × 151 bp) sequencing of the pools to a 27.1× (weeping) and a 30.4× (standard) genome (742.3 Mb) coverage allowed detection of 84562 DNA variants specific to 'weeping', 92148 specific to 'standard', and 173169 common to both pools. A detailed analysis of the DNA variant genotypes in the pools predicted three informative segregation types of variants: (type I) in weeping pool-specific variants, and (type II) and (type III) in variants common to both pools, where the first allele is assumed to be weeping linked and the allele shown in bold is a variant in relation to the reference genome. Conducting variant allele frequency and density-based mappings revealed four genomic regions with a significant association with weeping: a major locus, Weeping (W), on chromosome 13 and others on chromosomes 10 (W2), 16 (W3), and 5 (W4). The results from type I variants were noisier and less certain than those from type II and type III variants, demonstrating that although type I variants are often the first choice, type II and type III variants represent an important source of DNA variants that can be exploited for genetic mapping in out-crossing woody species. Confirmation of the mapping of W and W2, investigation into their genetic interactions, and identification of expressed genes in the W and W2 regions provided insight into the genetic control of weeping and its expressivity in Malus.

Sequencing of bulks of segregants allows dissection of genetic control of amylose content in rice

Amylose content (AC) is a key quality trait in rice. A cross between Oryza glaberrima (African rice) and Oryza sativa (Asian rice) segregating for AC was analysed by sequencing bulks of individuals with high and low AC. SNP associated with the granule bound starch synthase (GBSS1) locus on chromosome 6 were polymorphic between the bulks. In particular, a G/A SNP that would result in an Asp to Asn mutation was identified. This amino acid substitution may be responsible for differences in GBSS activity as it is adjacent to a disulphide linkage conserved in all grass GBSS proteins. Other polymorphisms in genomic regions closely surrounding this variation may be the result of linkage drag. In addition to the variant in the starch biosynthesis gene, SNP on chromosomes 1 and 11 linked to AC was also identified. SNP was found in the genes encoding the NAC and CCAAT-HAP5 transcription factors that have previously been linked to starch biosynthesis. This study has demonstrated that the approach of sequencing bulks was able to identify genes on different chromosomes associated with this complex trait.

Mutation of Arabidopsis SMC4 identifies condensin as a corepressor of pericentromeric transposons and conditionally expressed genes

In this study, Wang et al. perform genome-wide analyses that implicate condensin in the suppression of hundreds of loci, acting in both DNA methylation-dependent and methylation-independent pathways. They show that silencing of transposons in the pericentromeric heterochromatin of Arabidopsis thaliana requires SMC4, a core subunit of condensins I and II, acting in conjunction with CG methylation by MET1, CHG methylation by CMT3, the chromatin remodeler DDM1, and histone modifications, including H3K27me1, imparted by ATXR5 and ATXR6.

In eukaryotes, transcriptionally inactive loci are enriched within highly condensed heterochromatin. In plants, as in mammals, the DNA of heterochromatin is densely methylated and wrapped by histones displaying a characteristic subset of post-translational modifications. Growing evidence indicates that these chromatin modifications are not sufficient for silencing. Instead, they are prerequisites for further assembly of higher-order chromatin structures that are refractory to transcription but not fully understood. We show that silencing of transposons in the pericentromeric heterochromatin of Arabidopsis thaliana requires SMC4, a core subunit of condensins I and II, acting in conjunction with CG methylation by MET1 (DNA METHYLTRANSFERASE 1), CHG methylation by CMT3 (CHROMOMETHYLASE 3), the chromatin remodeler DDM1 (DECREASE IN DNA METHYLATION 1), and histone modifications, including histone H3 Lys 27 monomethylation (H3K27me1), imparted by ATXR5 and ATXR6. SMC4/condensin also acts within the mostly euchromatic chromosome arms to suppress conditionally expressed genes involved in flowering or DNA repair, including the DNA glycosylase ROS1, which facilitates DNA demethylation. Collectively, our genome-wide analyses implicate condensin in the suppression of hundreds of loci, acting in both DNA methylation-dependent and methylation-independent pathways.

Indel-seq: a fast-forward genetics approach for identification of trait-associated putative candidate genomic regions and its application in pigeonpea (Cajanus cajan)

Identification of candidate genomic regions associated with target traits using conventional mapping methods is challenging and time-consuming. In recent years, a number of single nucleotide polymorphism (SNP)-based mapping approaches have been developed and used for identification of candidate/putative genomic regions. However, in the majority of these studies, insertion-deletion (Indel) were largely ignored. For efficient use of Indels in mapping target traits, we propose Indel-seq approach, which is a combination of whole-genome resequencing (WGRS) and bulked segregant analysis (BSA) and relies on the Indel frequencies in extreme bulks. Deployment of Indel-seq approach for identification of candidate genomic regions associated with fusarium wilt (FW) and sterility mosaic disease (SMD) resistance in pigeonpea has identified 16 Indels affecting 26 putative candidate genes. Of these 26 affected putative candidate genes, 24 genes showed effect in the upstream/downstream of the genic region and two genes showed effect in the genes. Validation of these 16 candidate Indels in other FW- and SMD-resistant and FW- and SMD-susceptible genotypes revealed a significant association of five Indels (three for FW and two for SMD resistance). Comparative analysis of Indel-seq with other genetic mapping approaches highlighted the importance of the approach in identification of significant genomic regions associated with target traits. Therefore, the Indel-seq approach can be used for quick and precise identification of candidate genomic regions for any target traits in any crop species.

Large-scale SNP screenings identify markers linked with GCRV resistant traits through transcriptomes of individuals and cell lines in Ctenopharyngodon idella

Grass carp (Ctenopharyngodon idella) is an important economic species in freshwater aquaculture and its industry has been confined due to variety degeneration and frequent diseases. Marker-assisted selection is a feasible method for selective breeding of new varieties. Transcriptome data have greatly facilitated high-throughput single nucleotide polymorphism (SNP) marker discovery and phenotype association study. In this study, we gained a total of 25,981 and 5,775 high quality SNPs in two transcriptomes from individuals and cell lines, respectively. Comparative transcriptome analysis identified 413 and 832 grass carp reovirus (GCRV)-resistant-association SNPs as well as 1,381 and 1,606 GCRV-susceptible-association SNPs in individuals and cell lines, respectively. Integrated analysis indicated 22 genes with single SNP share common resistant/susceptible traits in two transcriptomes. Furthermore, we infected grass carp with GCRV, genotyping and association analyses were performed, and 9 in 22 SNPs were confirmed by PCR-RFLP. Meanwhile, mRNA expression profiles of 6 genes containing confirmed SNPs were examined by qRT-PCR. The results demonstrated that mRNA expressions were significant differences in resistant/susceptible individuals and cell lines. The present study develops an important strategy for high throughput screening of phenotype association genetic markers and the results will serve in grass carp breeding for GCRV resistance.

Ctf4-related protein recruits LHP1-PRC2 to maintain H3K27me3 levels in dividing cells in Arabidopsis thaliana

Polycomb Repressive Complex (PRC) 2 catalyzes the H3K27me3 modification that warrants inheritance of a repressive chromatin structure during cell division, thereby assuring stable target gene repression in differentiated cells. It is still under investigation how H3K27me3 is passed on from maternal to filial strands during DNA replication; however, cell division can reinforce H3K27me3 coverage at target regions. To identify novel factors involved in the Polycomb pathway in plants, we performed a forward genetic screen for enhancers of the like heterochromatin protein 1 (lhp1) mutant, which shows relatively mild phenotypic alterations compared with other plant PRC mutants. We mapped enhancer of lhp1 (eol) 1 to a gene related to yeast Chromosome transmission fidelity 4 (Ctf4) based on phylogenetic analysis, structural similarities, physical interaction with the CMG helicase component SLD5, and an expression pattern confined to actively dividing cells. A combination of eol1 with the curly leaf (clf) allele, carrying a mutation in the catalytic core of PRC2, strongly enhanced the clf phenotype; furthermore, H3K27me3 coverage at target genes was strongly reduced in eol1 clf double mutants compared with clf single mutants. EOL1 physically interacted with CLF, its partially redundant paralog SWINGER (SWN), and LHP1. We propose that EOL1 interacts with LHP1-PRC2 complexes during replication and thereby participates in maintaining the H3K27me3 mark at target genes.

Governing the Silencing State of Chromatin: The Roles of Polycomb Repressive Complex 1 in Arabidopsis

Polycomb group proteins form multiple protein complexes such as Polycomb Repressive Complex (PRC) 1 and PRC2, which repress the expression of thousands of genes. PRC1 and PRC2 are essential for normal development in Arabidopsis. Recently, significant progress has been made in understanding the functions and regulatory mechanisms of PRC1. In this review, we focus on the discovery of the composition of PRC1, functions of its components, the recruitment of PRC1 to target genes and the control of PRC1 function in Arabidopsis. Perspectives on dissecting the roles of PRC1 in plant gene expression and development are also given.

Efficient Identification of Causal Mutations through Sequencing of Bulked F 2 from Two Allelic Bloomless Mutants of Sorghum bicolor

Sorghum (Sorghum bicolor Moench, L.) plant accumulates copious layers of epi-cuticular wax (EW) on its aerial surfaces, to a greater extent than most other crops. EW provides a vapor barrier that reduces water loss, and is therefore considered to be a major determinant of sorghum's drought tolerance. However, little is known about the genes responsible for wax accumulation in sorghum. We isolated two allelic mutants, bloomless40-1 (bm40-1) and bm40-2, from a mutant library constructed from ethyl methane sulfonate (EMS) treated seeds of an inbred, BTx623. Both mutants were nearly devoid of the EW layer. Each bm mutant was crossed to the un-mutated BTx623 to generated F₂ populations that segregated for the bm phenotype. Genomic DNA from 20 bm F₂ plants from each population was bulked for whole genome sequencing. A single gene, Sobic.001G228100, encoding a GDSL-like lipase/acylhydrolase, had unique homozygous mutations in each bulked F₂ population. Mutant bm40-1 harbored a missense mutation in the gene, whereas bm40-2 had a splice donor site mutation. Our findings thus provide strong evidence that mutation in this GDSL-like lipase gene causes the bm phenotype, and further demonstrate that this approach of sequencing two independent allelic mutant populations is an efficient method for identifying causal mutations. Combined with allelic mutants, MutMap provides powerful method to identify all causal genes for the large collection of bm mutants in sorghum, which will provide insight into how sorghum plants accumulate such abundant EW on their aerial surface. This knowledge may facilitate the development of tools for engineering drought-tolerant crops with reduced water loss.

Development and utilization of a new chemically-induced soybean library with a high mutation density

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.

Next-Generation Sequencing from Bulked-Segregant Analysis Accelerates the Simultaneous Identification of Two Qualitative Genes in Soybean

Next-generation sequencing (NGS)-based bulked-segregant analysis (BSA) approaches have been proven successful for rapidly mapping genes in plant species. However, most such methods are based on mutants and usually only one gene controlling the mutant phenotype is identified. In this study, NGS-based BSA was employed to map simultaneously two qualitative genes controlling cotyledon color of seed in soybean. Yellow-cotyledon (YC) and green-cotyledon (GC) bulks from progenies of a biparental population (Zhonghuang 30 × Jiyu 102) were sequenced. The SNP-index of each SNP locus in YC and GC bulks was calculated and two genomic regions on chromosomes 1 and 11 harboring, respectively, loci qCC1 and qCC2 were identified by Δ(SNP-index) analysis. These two BSA-seq-derived loci were further validated with SSR markers and fine-mapped. qCC1 was mapped to a 30.7-kb region containing four annotated genes and qCC2 was mapped to a 67.7-kb region with nine genes. These two regions contained, respectively, genes D1 and D2, which had previously been identified by homology-based cloning as being associated with cotyledon color. Sequence analysis of the NGS data also identified a frameshift deletion in the coding region of D1. These results suggested that BSA-seq could accelerate the mapping of loci controlling qualitative traits, even if a trait is controlled by more than one locus.

From Mendel's discovery on pea to today's plant genetics and breeding : Commemorating the 150th anniversary of the reading of Mendel's discovery

KEY MESSAGE

This work discusses several selected topics of plant genetics and breeding in relation to the 150th anniversary of the seminal work of Gregor Johann Mendel. In 2015, we celebrated the 150th anniversary of the presentation of the seminal work of Gregor Johann Mendel. While Darwin's theory of evolution was based on differential survival and differential reproductive success, Mendel's theory of heredity relies on equality and stability throughout all stages of the life cycle. Darwin's concepts were continuous variation and "soft" heredity; Mendel espoused discontinuous variation and "hard" heredity. Thus, the combination of Mendelian genetics with Darwin's theory of natural selection was the process that resulted in the modern synthesis of evolutionary biology. Although biology, genetics, and genomics have been revolutionized in recent years, modern genetics will forever rely on simple principles founded on pea breeding using seven single gene characters. Purposeful use of mutants to study gene function is one of the essential tools of modern genetics. Today, over 100 plant species genomes have been sequenced. Mapping populations and their use in segregation of molecular markers and marker-trait association to map and isolate genes, were developed on the basis of Mendel's work. Genome-wide or genomic selection is a recent approach for the development of improved breeding lines. The analysis of complex traits has been enhanced by high-throughput phenotyping and developments in statistical and modeling methods for the analysis of phenotypic data. Introgression of novel alleles from landraces and wild relatives widens genetic diversity and improves traits; transgenic methodologies allow for the introduction of novel genes from diverse sources, and gene editing approaches offer possibilities to manipulate gene in a precise manner.

Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single-nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild-type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.

Isolation of New Gravitropic Mutants under Hypergravity Conditions

Forward genetics is a powerful approach used to link genotypes and phenotypes, and mutant screening/analysis has provided deep insights into many aspects of plant physiology. Gravitropism is a tropistic response in plants, in which hypocotyls and stems sense the direction of gravity and grow upward. Previous studies of gravitropic mutants have suggested that shoot endodermal cells in Arabidopsis stems and hypocotyls are capable of sensing gravity (i.e., statocytes). In the present study, we report a new screening system using hypergravity conditions to isolate enhancers of gravitropism mutants, and we also describe a rapid and efficient genome mapping method, using next-generation sequencing (NGS) and single nucleotide polymorphism (SNP)-based markers. Using the endodermal-amyloplast less 1 (eal1) mutant, which exhibits defective development of endodermal cells and gravitropism, we found that hypergravity (10 g) restored the reduced gravity responsiveness in eal1 hypocotyls and could, therefore, be used to obtain mutants with further reduction in gravitropism in the eal1 background. Using the new screening system, we successfully isolated six ene (enhancer of eal1) mutants that exhibited little or no gravitropism under hypergravity conditions, and using NGS and map-based cloning with SNP markers, we narrowed down the potential causative genes, which revealed a new genetic network for shoot gravitropism in Arabidopsis.

Exploring Genetic Diversity in Plants Using High-Throughput Sequencing Techniques

Food security has emerged as an urgent concern because of the rising world population. To meet the food demands of the near future, it is required to improve the productivity of various crops, not just of staple food crops. The genetic diversity among plant populations in a given species allows the plants to adapt to various environmental conditions. Such diversity could therefore yield valuable traits that could overcome the food-security challenges. To explore genetic diversity comprehensively and to rapidly identify useful genes and/or allele, advanced high-throughput sequencing techniques, also called next-generation sequencing (NGS) technologies, have been developed. These provide practical solutions to the challenges in crop genomics. Here, we review various sources of genetic diversity in plants, newly developed genetic diversity-mining tools synergized with NGS techniques, and related genetic approaches such as quantitative trait locus analysis and genome-wide association study.

Mapping-by-sequencing in complex polyploid genomes using genic sequence capture: a case study to map yellow rust resistance in hexaploid wheat

Previously we extended the utility of mapping-by-sequencing by combining it with sequence capture and mapping sequence data to pseudo-chromosomes that were organized using wheat-Brachypodium synteny. This, with a bespoke haplotyping algorithm, enabled us to map the flowering time locus in the diploid wheat Triticum monococcum L. identifying a set of deleted genes (Gardiner et al., 2014). Here, we develop this combination of gene enrichment and sliding window mapping-by-synteny analysis to map the Yr6 locus for yellow stripe rust resistance in hexaploid wheat. A 110 MB NimbleGen capture probe set was used to enrich and sequence a doubled haploid mapping population of hexaploid wheat derived from an Avalon and Cadenza cross. The Yr6 locus was identified by mapping to the POPSEQ chromosomal pseudomolecules using a bespoke pipeline and algorithm (Chapman et al., 2015). Furthermore the same locus was identified using newly developed pseudo-chromosome sequences as a mapping reference that are based on the genic sequence used for sequence enrichment. The pseudo-chromosomes allow us to demonstrate the application of mapping-by-sequencing to even poorly defined polyploidy genomes where chromosomes are incomplete and sub-genome assemblies are collapsed. This analysis uniquely enabled us to: compare wheat genome annotations; identify the Yr6 locus - defining a smaller genic region than was previously possible; associate the interval with one wheat sub-genome and increase the density of SNP markers associated. Finally, we built the pipeline in iPlant, making it a user-friendly community resource for phenotype mapping.