Isabelle M. Henry

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

The authors combined high-throughput global sequencing targeted to genes and custom-designed bioinformatics tools to catalog and characterize chemically induced mutations in rice and wheat. They demonstrate that this approach can be used to develop large-scale induced mutation resources with relatively small investments and is applicable to resource-limited and polyploid organisms. 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.

A Y-chromosome–encoded small RNA acts as a sex determinant in persimmons

Y male plants affect female development Although most plants have both male and female organs within a single flower, some produce separate male and female plants. In some cases, such as persimmons, males are determined by a Y chromosome. Akagi et al. examined the gene transcript differences between male and female persimmons. A gene on the Y chromosome regulated a non–sex chromosome–linked small RNA that suppresses female organ development. This small RNA was localized to male flowers and could affect female development in other plant species. The evolutionary history of these genes suggests that they are tied to the origin of the separation of sexes in the persimmon family. Science, this issue p. 646 The Y chromosome in date plum encodes a small RNA involved in male sex determination. In plants, multiple lineages have evolved sex chromosomes independently, providing a powerful comparative framework, but few specific determinants controlling the expression of a specific sex have been identified. We investigated sex determinants in the Caucasian persimmon, Diospyros lotus, a dioecious plant with heterogametic males (XY). Male-specific short nucleotide sequences were used to define a male-determining region. A combination of transcriptomics and evolutionary approaches detected a Y-specific sex-determinant candidate, OGI, that displays male-specific conservation among Diospyros species. OGI encodes a small RNA targeting the autosomal MeGI gene, a homeodomain transcription factor regulating anther fertility in a dosage-dependent fashion. This identification of a feminizing gene suppressed by a Y-chromosome–encoded small RNA contributes to our understanding of the evolution of sex chromosome systems in higher plants.

A genome-wide analysis of Cas9 binding specificity using ChIP-seq and targeted sequence capture

Clustered regularly interspaced short palindromic repeat (CRISPR) RNA-guided nucleases have gathered considerable excitement as a tool for genome engineering. However, questions remain about the specificity of target site recognition. Cleavage specificity is typically evaluated by low throughput assays (T7 endonuclease I assay, target amplification followed by high-throughput sequencing), which are limited to a subset of potential off-target sites. Here, we used ChIP-seq to examine genome-wide CRISPR binding specificity at gRNA-specific and gRNA-independent sites for two guide RNAs. RNA-guided Cas9 binding was highly specific to the target site while off-target binding occurred at much lower intensities. Cas9-bound regions were highly enriched in NGG sites, a sequence required for target site recognition by Streptococcus pyogenes Cas9. To determine the relationship between Cas9 binding and endonuclease activity, we applied targeted sequence capture, which allowed us to survey 1200 genomic loci simultaneously including potential off-target sites identified by ChIP-seq and by computational prediction. A high frequency of indels was observed at both target sites and one off-target site, while no cleavage activity could be detected at other ChIP-bound regions. Our results confirm the high-specificity of CRISPR endonucleases and demonstrate that sequence capture can be used as a high-throughput genome-wide approach to identify off-target activity.

Highly active zinc finger nucleases by extended modular assembly

Zinc-finger nucleases (ZFNs) are important tools for genome engineering. Despite intense interest by many academic groups, the lack of robust noncommercial methods has hindered their widespread use. The modular assembly (MA) of ZFNs from publicly available one-finger archives provides a rapid method to create proteins that can recognize a very broad spectrum of DNA sequences. However, three- and four-finger arrays often fail to produce active nucleases. Efforts to improve the specificity of the one-finger archives have not increased the success rate above 25%, suggesting that the MA method might be inherently inefficient due to its insensitivity to context-dependent effects. Here we present the first systematic study on the effect of array length on ZFN activity. ZFNs composed of six-finger MA arrays produced mutations at 15 of 21 (71%) targeted loci in human and mouse cells. A novel drop-out linker scheme was used to rapidly assess three- to six-finger combinations, demonstrating that shorter arrays could improve activity in some cases. Analysis of 268 array variants revealed that half of MA ZFNs of any array composition that exceed an ab initio B-score cutoff of 15 were active. These results suggest that, when used appropriately, MA ZFNs are able to target more DNA sequences with higher success rates than other current methods.

Phenotypic Consequences of Aneuploidy in Arabidopsis thaliana

Aneuploid cells are characterized by incomplete chromosome sets. The resulting imbalance in gene dosage has phenotypic consequences that are specific to each karyotype. Even in the case of Down syndrome, the most viable and studied form of human aneuploidy, the mechanisms underlying the connected phenotypes remain mostly unclear. Because of their tolerance to aneuploidy, plants provide a powerful system for a genome-wide investigation of aneuploid syndromes, an approach that is not feasible in animal systems. Indeed, in many plant species, populations of aneuploid individuals can be easily obtained from triploid individuals. We phenotyped a population of Arabidopsis thaliana aneuploid individuals containing 25 different karyotypes. Even in this highly heterogeneous population, we demonstrate that certain traits are strongly associated with the dosage of specific chromosome types and that chromosomal effects can be additive. Further, we identified subtle developmental phenotypes expressed in the diploid progeny of aneuploid parent(s) but not in euploid controls from diploid lineages. These results indicate long-term phenotypic consequences of aneuploidy that can persist after chromosomal balance has been restored. We verified the diploid nature of these individuals by whole-genome sequencing and discuss the possibility that trans-generational phenotypic effects stem from epigenetic modifications passed from aneuploid parents to their diploid progeny.

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 F1. Haploid plants produced from an F1 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.

Selection and validation of reference genes for quantitative RT-PCR expression studies of the non-model crop Musa

Gene expression analysis by reverse transcriptase real-time or quantitative polymerase chain reaction (RT-qPCR) is becoming widely used for non-model plant species. Given the high sensitivity of this method, normalization using multiple housekeeping or reference genes is critical, and careful selection of these reference genes is one of the most important steps to obtain reliable results. In this study, reference genes commonly used for other plant species were investigated to identify genes displaying highly uniform expression patterns in different varieties, tissues, developmental stages, fungal infection, and osmotic stress conditions for the non-model crop Musa (banana and plantains). The expression stability of six candidate reference genes was tested on six different sample sets, and the results were analyzed using the publicly available algorithms geNorm and NormFinder. Our results show that variety, plant material, primer set, and gene identity can all influence the robustness and outcome of RT-qPCR analysis. In the case of Musa, a combination of three reference genes (EF1, TUB and ACT) can be used for normalization of gene expression data from greenhouse leaf samples. In the case of shoot meristem cultures, numerous combinations can be used because the investigated reference genes exhibited limited variability. In contrast, variability in expression of the reference genes was much larger among leaf samples from plants grown in vitro, for which the best combination of reference genes (L2 and ACT genes) is still suboptimal. Overall, our data confirm that the stability of candidate reference genes should be thoroughly investigated for each experimental condition under investigation.

Genetic Basis for Dosage Sensitivity in Arabidopsis thaliana

Aneuploidy, the relative excess or deficiency of specific chromosome types, results in gene dosage imbalance. Plants can produce viable and fertile aneuploid individuals, while most animal aneuploids are inviable or developmentally abnormal. The swarms of aneuploid progeny produced by Arabidopsis triploids constitute an excellent model to investigate the mechanisms governing dosage sensitivity and aneuploid syndromes. Indeed, genotype alters the frequency of aneuploid types within these swarms. Recombinant inbred lines that were derived from a triploid hybrid segregated into diploid and tetraploid individuals. In these recombinant inbred lines, a single locus, which we call SENSITIVE TO DOSAGE IMBALANCE (SDI), exhibited segregation distortion in the tetraploid subpopulation only. Recent progress in quantitative genotyping now allows molecular karyotyping and genetic analysis of aneuploid populations. In this study, we investigated the causes of the ploidy-specific distortion at SDI. Allele frequency was distorted in the aneuploid swarms produced by the triploid hybrid. We developed a simple quantitative measure for aneuploidy lethality and using this measure demonstrated that distortion was greatest in the aneuploids facing the strongest viability selection. When triploids were crossed to euploids, the progeny, which lack severe aneuploids, exhibited no distortion at SDI. Genetic characterization of SDI in the aneuploid swarm identified a mechanism governing aneuploid survival, perhaps by buffering the effects of dosage imbalance. As such, SDI could increase the likelihood of retaining genomic rearrangements such as segmental duplications. Additionally, in species where triploids are fertile, aneuploid survival would facilitate gene flow between diploid and tetraploid populations via a triploid bridge and prevent polyploid speciation. Our results demonstrate that positional cloning of loci affecting traits in populations containing ploidy and chromosome number variants is now feasible using quantitative genotyping approaches.

Structure and regulation of the Asr gene family in banana

Abscisic acid, stress, ripening proteins (ASR) are a family of plant-specific small hydrophilic proteins. Studies in various plant species have highlighted their role in increased resistance to abiotic stress, including drought, but their specific function remains unknown. As a first step toward their potential use in crop improvement, we investigated the structure and regulation of the Asr gene family in Musa species (bananas and plantains). We determined that the MusaAsr gene family contained at least four members, all of which exhibited the typical two exons, one intron structure of Asr genes and the “ABA/WDS” (abscisic acid/water deficit stress) domain characteristic of Asr genes. Phylogenetic analyses determined that the MusaAsr genes were closely related to each other, probably as the product of recent duplication events. For two of the four members, two versions corresponding to the two sub-genomes of Musa, acuminata and balbisiana were identified. Gene expression and protein analyses were performed and Asr expression could be detected in meristem cultures, root, pseudostem, leaf and cormus. In meristem cultures, mAsr1 and mAsr3 were induced by osmotic stress and wounding, while mAsr3 and mAsr4 were induced by exposure to ABA. mASR3 exhibited the most variation both in terms of amino acid sequence and expression pattern, making it the most promising candidate for further functional study and use in crop improvement.

Significant enhancement of fatty acid composition in seeds of the allohexaploid, Camelina sativa, using CRISPR/Cas9 gene editing

Summary The CRISPR/Cas9 nuclease system is a powerful and flexible tool for genome editing, and novel applications of this system are being developed rapidly. Here, we used CRISPR/Cas9 to target the FAD2 gene in Arabidopsis thaliana and in the closely related emerging oil seed plant, Camelina sativa, with the goal of improving seed oil composition. We successfully obtained Camelina seeds in which oleic acid content was increased from 16% to over 50% of the fatty acid composition. These increases were associated with significant decreases in the less desirable polyunsaturated fatty acids, linoleic acid (i.e. a decrease from ~16% to <4%) and linolenic acid (a decrease from ~35% to <10%). These changes result in oils that are superior on multiple levels: they are healthier, more oxidatively stable and better suited for production of certain commercial chemicals, including biofuels. As expected, A. thaliana T2 and T3 generation seeds exhibiting these types of altered fatty acid profiles were homozygous for disrupted FAD2 alleles. In the allohexaploid, Camelina, guide RNAs were designed that simultaneously targeted all three homoeologous FAD2 genes. This strategy that significantly enhanced oil composition in T3 and T4 generation Camelina seeds was associated with a combination of germ‐line mutations and somatic cell mutations in FAD2 genes in each of the three Camelina subgenomes.