Multiple Alleles Research Articles

Simplification of genotyping techniques of the ABO blood type experiment and exploration of population genetics.

The ABO blood type is one of the most common and widely used genetic traits in humans. Three glycosyltransferase-encoding gene alleles, IA, IB and i, produce three red blood cell surface antigens, by which the ABO blood type is classified. By using the ABO blood type experiment as an ideal case for genetics teaching, we can easily introduce to the students several genetic concepts, including multiple alleles, gene interaction, single nucleotide polymorphism (SNP) and gene evolution. Herein we have innovated and integrated our ABO blood type genetics experiments. First, in the section of Molecular Genetics, a new method of ABO blood genotyping was established: specific primers based on SNP sites were designed to distinguish three alleles through quantitative real-time PCR. Next, the experimental teaching method of Gene Evolution was innovated in the Population Genetics section: a gene-evolution software was developed to simulate the evolutionary tendency of the ABO genotype encoding alleles under diverse conditions. Our reform aims to extend the contents of genetics experiments, to provide additional teaching approaches, and to improve the learning efficiency of our students eventually.

Versatile and precise gene-targeting strategies for functional studies in mammalian cell lines.

The advent of programmable nucleases such as ZFNs, TALENs and CRISPR/Cas9 has brought the power of genetic manipulation to widely used model systems. In mammalian cells, nuclease-mediated DNA double strand break is mainly repaired through the error-prone non-homologous end-joining (NHEJ) repair pathway, eventually leading to accumulation of small deletions or insertions (indels) that can inactivate gene function. However, due to the variable size of the indels and the polyploid status of many cell lines (e.g., cancer-derived cells), obtaining a knockout usually requires lengthy screening and characterization procedures. Given the more precise type of modifications that can be introduced upon homology-directed repair (HDR), we have developed HDR-based gene-targeting strategies that greatly facilitate the process of knockout generation in cell lines. To generate reversible knockouts (R-KO), a selectable promoter-less STOP cassette is inserted in an intron, interrupting transcription. Loss-of-function can be validated by RT-qPCR and is removable, enabling subsequent restoration of gene function. A variant of the R-KO procedure can be used to introduce point mutations. To generate constitutive knockouts (C-KO), an exon is targeted, which makes use of HDR-based gene disruption together with NHEJ-induced indels on non-HDR targeted allele(s). Hence the C-KO procedure greatly facilitates simultaneous inactivation of multiple alleles. Overall these genome-editing tools offer superior precision and efficiency for functional genetic approaches. We provide detailed protocols guiding in the design of targeting vectors and in the analysis and validation of gene targeting experiments.

Genetic characterization and diversity among avocado (Persea americana Mill.) genotypes from INIA-CENIAP, Venezuela

Studies on Persea americana have been addressed in different ways with biochemical and molecular techniques. Microsatellites are able to detect multiple alleles for particular loci and are therefore a useful tool to study genealogical relationships, population structures and genetic mapping. Ninety-six samples from 49 cultivars including three horticultural groups and hybrids were collected from the avocado germplasm bank at INIA-CENIAP (Venezuela). A modified DNA extraction protocol was performed. Forty microsatellites were selected from previous references, PCR amplifications were performed, and presence/absence, size, and number of alleles were evaluated on polyacrylamide gels. Attributes for polymorphic alleles were analyzed with POPGENE, and genetic diversity was calculated by effective sample size, number of alleles per locus (Na), effective number of alleles (Ne), Shannon information index (In), observed heterozygosis (H), expected heterozygosity (He), Wright’s fixation index (Fis), and allele frequencies. Only 14 primers were amplified, and AVT106 primer resulted monomorphic. Unique genotypes for each sample were obtained. Nine loci showed allele patterns that can be useful for taxonomic identification of cultivars or varieties. Comparing values of Fis with Ho and He, we found a direct relationship where low heterozygosis alleles identified in the population may affect the expected level. Allele frequencies ranged from 0.5632 to 0.0105. For all loci, at least one rare allele was observed. With the available information from genetic analysis, an identifying system was implemented for selected avocado cultivars maintained at the INIA-CENIAP Venezuelan germplasm bank on the basis of molecular data.

HLA-Bw4 80(T) and multiple HLA-Bw4 copies combined with KIR3DL1 associate with spontaneous clearance of HCV infection in people who inject drugs

Natural killer (NK) cell function is regulated by inhibitory and activating receptors including killer cell immunoglobulin-like receptors (KIRs). Here, we analyzed the impact of different KIR/KIR-ligand genotypes on the outcome of hepatitis C virus (HCV) infection in people who inject drugs (PWID). KIR/KIR-ligand genotypes associated with spontaneous clearance of HCV infection were identified in a cohort of PWID from Germany (n=266) and further validated in a second anti-HCV positive cohort of PWID recruited in North America (n=342). NK cells of PWID and healthy donors were functionally characterized according to their KIR/KIR-ligand genotype by flow cytometry. Multivariate logistic regression analysis revealed that KIR3DL1/HLA-Bw4 80(T) was associated with spontaneous clearance of HCV infection in PWID, which was confirmed in the PWID cohort from North America. Compared with PWID with detectable HCV RNA, the frequency of individuals with multiple HLA-Bw4 alleles was significantly higher in anti-HCV positive PWID with resolved HCV infection (29.7% vs. 15.2%; p=0.0229) and in anti-HCV seronegative PWID (39.2%; p=0.0006). KIR3DL1+ NK cells from HLA-Bw4 80(T)-positive PWID showed superior functionality compared to HLA-Bw4 80(I)-positive PWID. This differential impact was not observed in healthy donors; however, the HLA-Bw4 copy number strongly correlated with the functionality of KIR3DL1+ NK cells. HLA-Bw4-80(T) and multiple HLA-Bw4 copies in combination with KIR3DL1 are associated with protection against chronic hepatitis C in PWID by distinct mechanisms. Better licensing of KIR3DL1+ NK cells in the presence of multiple HLA-Bw4 copies is beneficial prior to seroconversion whereas HLA-Bw4 80(T) may be beneficial during acute hepatitis C. Lay summary: Natural killer (NK) cells are part of the innate immune system and are regulated by a complex network of activating and inhibiting receptors. The regulating receptor-ligand pairs of an individual are genetically determined. Here, we identified a particular set of ligand and receptor genes that are associated with better functionality of NK cells and better outcome upon exposure to HCV in a high-risk group.

The Predicted Cross Value for Genetic Introgression of Multiple Alleles.

We consider the plant genetic improvement challenge of introgressing multiple alleles from a homozygous donor to a recipient. First, we frame the project as an algorithmic process that can be mathematically formulated. We then introduce a novel metric for selecting breeding parents that we refer to as the predicted cross value (PCV). Unlike estimated breeding values, which represent predictions of general combining ability, the PCV predicts specific combining ability. The PCV takes estimates of recombination frequencies as an input vector and calculates the probability that a pair of parents will produce a gamete with desirable alleles at all specified loci. We compared the PCV approach with existing estimated-breeding-value approaches in two simulation experiments, in which 7 and 20 desirable alleles were to be introgressed from a donor line into a recipient line. Results suggest that the PCV is more efficient and effective for multi-allelic trait introgression. We also discuss how operations research can be used for other crop genetic improvement projects and suggest several future research directions.

Maintenance of polygenic sex determination in a fluctuating environment: an individual-based model.

R. A. Fisher predicted that individuals should invest equally in offspring of both sexes, and that the proportion of males and females produced (the primary sex ratio) should evolve towards 1:1 when unconstrained. For many species, sex determination is dependent on sex chromosomes, creating a strong tendency for balanced sex ratios, but in other cases, multiple autosomal genes interact to determine sex. In such cases, the maintenance of multiple sex-determining alleles at multiple loci and the consequent among-family variability in sex ratios presents a puzzle, as theory predicts that such systems should be unstable. Theory also predicts that environmental influences on sex can complicate outcomes of genetic sex determination, and that population structure may play a role. Tigriopus californicus, a copepod that lives in splash-pool metapopulations and exhibits polygenic and environment-dependent sex determination, presents a test case for relevant theory. We use this species as a model for parameterizing an individual-based simulation to investigate conditions that could maintain polygenic sex determination. We find that metapopulation structure can delay the degradation of polygenic sex determination and that periods of alternating frequency-dependent selection, imposed by seasonal fluctuations in environmental conditions, can maintain polygenic sex determination indefinitely.

Introgression of Clearfield™ rice crop traits into weedy red rice outcrosses

Studies were conducted to determine the impact of introgression of crop alleles into weedy rice populations. Seeds of 89 weedy rice accessions collected from Arkansas fields in 11 counties, with a history of Clearfield™ rice production, were planted in 2011 and sprayed twice with imazethapyr (0.071kg ai ha−1). Survivors were genotyped using an allele-specific PCR assay to detect the presence of resistance-endowing als mutations, S653N and G654E, in herbicide-resistant (HR) weedy rice that could have originated from Clearfield™ cultivars planted in the state. Hybridization between Clearfield™ rice and weedy rice was verified using four DNA markers RM 215, RM 220, and RM 234 microsatellites, and the and RID12 indel. Two to five HR plants, representing different plant types per accession per replication (698 plants total) were characterized for 16 morphological and biological traits. Seventy-nine percent of accessions were resistant to imazethapyr. The HR weedy rice plants carried the S653N mutation and were confirmed as outcrosses with Clearfield™ rice using DNA markers. Ten weedy rice accessions (nearly 25% of confirmed outcrosses) harbored multiple alleles from rice cultivars, indicating previous hybridization events with multiple cultivars or outcrossing with weedy rice. The majority of HR weedy rice accessions were similar to cultivated rice in height, culm angle, flag leaf length and width, days to flower, seed shattering and dormancy. Kernels of 20% of the parent weedy rice accessions were segregating for white and red bran color. Two of the parent weedy rice accessions with gold hull were homozygous for white bran and produced progenies with 100% white-bran kernels; thus, weedy rice was no longer exclusively red. Therefore, HR weedy rice outcrosses generally had acquired crop-like traits, which make it easier to deplete from the soil seedbank, but also difficult to scout.

Genome-wide identification of major protein families of cyanobacteria and genomic insight into the circadian rhythm

ABSTRACT Cyanobacteria are among the most abundant organisms present on earth and are considered to be one of the oldest known clades. Cyanobacteria are oxygenic photosynthetic bacteria and are well known as promising renewable sources of energy; therefore, it is important to understand aspects of their genomes in detail across species. Advances in sequencing technology and the availability of several cyanobacterial genomes have provided an excellent opportunity to understand the diversity and evolution of the cyanobacterial genome. Here, we compared the genomes of 62 different phototrophic cyanobacteria. Evaluation of genetic diversity of all the cyanobacteria species studied revealed that evolution from their common ancestors was polyphyletic. In addition, the genomes were very diverse and varied among species, and significant genomic diversity was observed at the species and strain level. Overall, we identified 56 different protein families of cyanobacteria species/strains and found that they varied significantly among strains of a species. The circadian clock proteins KaiA, KaiB and KaiC (KaiABC complex proteins) of cyanobacteria were found to be present and consistent in the majority of cyanobacterial species while absent in a few others. Evolutionary analysis of the KaiABC protein complex showed that the KaiA protein has a high frequency of polymorphism, and multiple alleles were found to be present at high frequency. These results demonstrated that evolution of phosphorylation events occurred via KaiA in the KaiABC complex. Furthermore, multiple sequence alignment showed that KaiA, KaiB and KaiC proteins are highly conserved in nature. Our results provide direct information regarding the presence of different protein or protein families in cyanobacteria. The information presented here will serve as a starting point to explore the genetic diversity of cyanobacteria with the potential to play important roles in biotechnological applications.

Reversion of BRCA1/2 Germline Mutations Detected in Circulating Tumor DNA From Patients With High-Grade Serous Ovarian Cancer.

Purpose Germline BRCA1 or BRCA2 mutations in patients with high-grade serous ovarian cancer (HGSC) are associated with favorable responses to chemotherapy. However, secondary intragenic (reversion) mutations that restore protein function lead to clinically significant rates of acquired resistance. The goal of this study was to determine whether reversion mutations could be found in an unbiased manner in circulating cell-free DNA (cfDNA) to predict treatment response in HGSC. Patients and Methods Plasma and tumor samples were obtained from 30 patients with HGSC with either BRCA1 or BRCA2 germline mutation. Two cohorts were ascertained: patients with a malignancy before undergoing primary HGSC debulking surgery (n = 14) or patients at disease recurrence (n = 16). Paired tumor and plasma samples were available for most patients (24 of 30). Targeted amplicon, next-generation sequencing was performed using primers that flanked germline mutations, whose design did not rely on prior knowledge of reversion sequences. Results Five patients were identified with intragenic mutations predicted to restore BRCA1/2 open reading frames, including two patients with multiple independent reversion alleles. Reversion mutations were only detected in tumor samples from patients with recurrent disease (five of 16) and only in cfDNA from patients with a tumor-detected reversion (three of five). Findings from a rapid autopsy of a patient with multiple independent reversions indicated that reversion-allele frequency in metastatic sites is an important determinant of assay sensitivity. Abundance of tumor-derived DNA in total cell-free DNA, as measured by TP53 mutant allele frequency, also affected assay sensitivity. All patients with reversions detected in tumor-derived DNA were resistant to platin- or poly ADP ribose polymerase inhibitor-based chemotherapy. Conclusion Reversion mutations can be detected in an unbiased analysis of cfDNA, suggesting clinical utility for predicting chemotherapy response in recurrent HGSC.

STELLS2: fast and accurate coalescent-based maximum likelihood inference of species trees from gene tree topologies

It is well known that gene trees and species trees may have different topologies. One explanation is incomplete lineage sorting, which is commonly modeled by the coalescent process. In multispecies coalescent, a gene tree topology is observed with some probability (called the gene tree probability) for a given species tree. Gene tree probability is the main tool for the program STELLS, which finds the maximum likelihood estimate of the species tree from the given gene tree topologies. However, STELLS becomes slow when data size increases. Recently, several fast species tree inference methods have been developed, which can handle large data. However, these methods often do not fully utilize the information in the gene trees. In this paper, we present an algorithm (called STELLS2) for computing the gene tree probability more efficiently than the original STELLS. The key idea of STELLS2 is taking some 'shortcuts' during the computation and computing the gene tree probability approximately. We apply the STELLS2 algorithm in the species tree inference approach in the original STELLS, which leads to a new maximum likelihood species tree inference method (also called STELLS2). Through simulation we demonstrate that the gene tree probabilities computed by STELLS2 and STELLS have strong correlation. We show that STELLS2 is almost as accurate in species tree inference as STELLS. Also STELLS2 is usually more accurate than several existing methods when there is one allele per species, although STELLS2 is slower than these methods. STELLS2 outperforms these methods significantly when there are multiple alleles per species. The program STELLS2 is available for download at: https://github.com/yufengwudcs/STELLS2. yufeng.wu@uconn.edu. Supplementary data are available at Bioinformatics online.

Mass Spectrometry Profiling of HLA-Associated Peptidomes in Mono-allelic Cells Enables More Accurate Epitope Prediction

Identification of human leukocyte antigen (HLA)-bound peptides by liquid chromatography-tandem mass spectrometry (LC-MS/MS) is poised to provide a deep understanding of rules underlying antigen presentation. However, a key obstacle is the ambiguity that arises from the co-expression of multiple HLA alleles. Here, we have implemented a scalable mono-allelic strategy for profiling the HLA peptidome. By using cell lines expressing a single HLA allele, optimizing immunopurifications, and developing an application-specific spectral search algorithm, we identified thousands of peptides bound to 16 different HLA class I alleles. These data enabled the discovery of subdominant binding motifs and an integrative analysis quantifying the contribution of factors critical to epitope presentation, such as protein cleavage and gene expression. We trained neural-network prediction algorithms with our large dataset (>24,000 peptides) and outperformed algorithms trained on datasets of peptides with measured affinities. We thus demonstrate a strategy for systematically learning the rules of endogenous antigen presentation.

Prions.

Infectious proteins (prions) are usually self-templating filamentous protein polymers (amyloids). Yeast prions are genes composed of protein and, like the multiple alleles of DNA-based genes, can have an array of "variants," each a distinct self-propagating amyloid conformation. Like the lethal mammalian prions and amyloid diseases, yeast prions may be lethal, or only mildly detrimental, and show an array of phenotypes depending on the protein involved and the prion variant. Yeast prions are models for both rare mammalian prion diseases and for several very common amyloidoses such as Alzheimer's disease, type 2 diabetes, and Parkinson's disease. Here, we describe their detection and characterization using genetic, cell biological, biochemical, and physical methods.

The pentatricopeptide repeat protein EMP9 is required for mitochondrial ccmB and rps4 transcript editing, mitochondrial complex biogenesis and seed development in maize.

Pentatricopeptide repeat (PPR) proteins comprise a large family of sequence-specific RNA binding proteins in land plants. Because of its large family size and frequent embryo lethality in the mutants, molecular functions and physiological roles of many PPR proteins are unknown. Through characterization of an empty pericarp9 (emp9) mutant in maize (Zea mays), we defined the functions of EMP9 in mitochondrial RNA editing, respiratory complex formation and seed development. Mu insertions in different regions of Emp9 facilitated dissection of the domain functions of the EMP9. Through genetic and functional analyses of multiple alleles, we showed that deletions of two N-terminal PPR motifs and partial E+ domain do not eliminate the editing function of EMP9. Emp9 encodes an E+ subclass PPR protein that is localized in mitochondria. Loss of EMP9 function abolishes the C-to-U editing of ccmB-43 and rps4-335 sites in mitochondria. The loss of editing at ccmB-43 and rps4-335 affects the maturation of cytochrome c and impairs the biogenesis of mitochondrial respiratory complexes, particularly complex III. This work extends our understanding of PPR-E+ protein in editing function and seed development, and provides insights into the molecular function of mitochondrial CcmB protein in higher plants.

Robust and rapid algorithms facilitate large-scale whole genome sequencing downstream analysis in an integrative framework.

Whole genome sequencing (WGS) is a promising strategy to unravel variants or genes responsible for human diseases and traits. However, there is a lack of robust platforms for a comprehensive downstream analysis. In the present study, we first proposed three novel algorithms, sequence gap-filled gene feature annotation, bit-block encoded genotypes and sectional fast access to text lines to address three fundamental problems. The three algorithms then formed the infrastructure of a robust parallel computing framework, KGGSeq, for integrating downstream analysis functions for whole genome sequencing data. KGGSeq has been equipped with a comprehensive set of analysis functions for quality control, filtration, annotation, pathogenic prediction and statistical tests. In the tests with whole genome sequencing data from 1000 Genomes Project, KGGSeq annotated several thousand more reliable non-synonymous variants than other widely used tools (e.g. ANNOVAR and SNPEff). It took only around half an hour on a small server with 10 CPUs to access genotypes of ∼60 million variants of 2504 subjects, while a popular alternative tool required around one day. KGGSeq's bit-block genotype format used 1.5% or less space to flexibly represent phased or unphased genotypes with multiple alleles and achieved a speed of over 1000 times faster to calculate genotypic correlation.

Identification and Characterization of Multiple Intermediate Alleles of the Key Genes Regulating Brassinosteroid Biosynthesis Pathways

Most of the early identified brassinosteroid signaling and biosynthetic mutants are null mutants, exhibiting extremely dwarfed phenotypes and male sterility. These null mutants are usually unable to be directly transformed via a routinely used Agrobacterium-mediated gene transformation system and therefore are less useful for genetic characterization of the brassinosteroid (BR)-related pathways. Identification of intermediate signaling mutants such as bri1–5 and bri1–9 has contributed drastically to the elucidation of BR signaling pathway using both genetic and biochemical approaches. However, intermediate mutants of key genes regulating BR biosynthesis have seldom been reported. Here we report identification of several intermediate BR biosynthesis mutants mainly resulted from leaky transcriptions due to the insertions of T-DNAs in the introns. These mutants are semi-dwarfed and fertile and capable to be transformed. These intermediate mutants could be useful tools for future discovery and analyses of novel components regulating BR biosynthesis and catabolism via genetic modifier screen.

Characterization of the sdw1 semi-dwarf gene in barley

BackgroundThe dwarfing gene sdw1 has been widely used throughout the world to develop commercial barley varieties. There are at least four different alleles at the sdw1 locus.ResultsMutations in the gibberellin 20-oxidase gene (HvGA20ox2) resulted in multiple alleles at the sdw1 locus. The sdw1.d allele from Diamant is due to a 7-bp deletion in exon 1, while the sdw1.c allele from Abed Denso has 1-bp deletion and a 4-bp insertion in the 5’ untranslated region. The sdw1.a allele from Jotun resulted from a total deletion of the HvGA20ox2 gene. The structural changes result in lower gene expression in sdw1.d and lack of expression in sdw1.a. There are three HvGA20ox genes in the barley genome. The partial or total loss of function of the HvGA20ox2 gene could be compensated by enhanced expression of its homolog HvGA20ox1and HvGA20ox3. A diagnostic molecular marker was developed to differentiate between the wild-type, sdw1.d and sdw1.a alleles and another molecular marker for differentiation of sdw1.c and sdw1.a. The markers were further tested in 197 barley varieties, out of which 28 had the sdw1.d allele and two varieties the sdw1.a allele. To date, the sdw1.d and sdw1.a alleles have only been detected in the modern barley varieties and lines.ConclusionsThe results provided further proof that the gibberellin 20-oxidase gene (HvGA20ox2) is the functional gene of the barley sdw1 mutants. Different deletions resulted in different functional alleles for different breeding purposes. Truncated protein could maintain partial function. Partial or total loss of function of the HvGA20ox2 gene could be compensated by enhanced expression of its homolog HvGA20ox1 and HvGA20ox3.

The role of the CNR1 gene in schizophrenia: a systematic review including unpublished data

Objective:Schizophrenia is a multifactorial disorder. It is known that a combination of extensive multiple common alleles may be involved in its etiology, each contributing with a small to moderate effect, and, possibly, some rare alleles with a much larger effect size. We aimed to perform a systematic review of association studies between schizophrenia (and its subphenotypes) and polymorphisms in the CNR1 gene, which encodes cannabinoid receptors classically implicated in schizophrenia pathophysiology, as well as to present unpublished results of an association study in a Brazilian population.Methods:Two reviewers independently searched for eligible studies and extracted outcome data using a structured form. Papers were retrieved from PubMed and ISI Web of Knowledge using the search term schizophrenia in combination with CNR1 or CB1 or cannabinoid receptor. Twenty-four articles met our inclusion criteria. We additionally present data from a study of our own comparing 182 patients with schizophrenia and 244 healthy controls.Results:No consistent evidence is demonstrated.Conclusion:Some seemingly positive association studies stress the need for further investigations of the possible role of endocannabinoid genetics in schizophrenia.

A Drosophila toolkit for defining gene function in spermatogenesis

Expression profiling and genomic sequencing methods enable the accumulation of vast quantities of data that relate to the expression of genes during the maturation of male germ cells from primordial germ cells to spermatozoa and potential mutations that underlie male infertility. However, the determination of gene function in specific aspects of spermatogenesis or linking abnormal gene function with infertility remain rate limiting, as even in an era of CRISPR analysis of gene function in mammalian models, this still requires considerable resources and time. Comparative developmental biology studies have shown the remarkable conservation of spermatogenic developmental processes from insects to vertebrates and provide an avenue of rapid assessment of gene function to inform the potential roles of specific genes in rodent and human spermatogenesis. The vinegar fly, Drosophila melanogaster, has been used as a model organism for developmental genetic studies for over one hundred years, and research with this organism produced seminal findings such as the association of genes with chromosomes, the chromosomal basis for sexual identity, the mutagenic properties of X-irradiation and the isolation of the first tumour suppressor mutations. Drosophila researchers have developed an impressive array of sophisticated genetic techniques for analysis of gene function and genetic interactions. This review focuses on how these techniques can be utilised to study spermatogenesis in an organism with a generation time of 9 days and the capacity to introduce multiple mutant alleles into an individual organism in a relatively short time frame.

Blocking Zebrafish MicroRNAs with Morpholinos.

Antisense morpholino oligonucleotides have been commonly used in zebrafish to inhibit mRNA function, either by inhibiting pre-mRNA splicing or by blocking translation initiation. Even with the advent of genome editing by CRISP/Cas9 technology, morpholinos provide a useful and rapid tool to knockdown gene expression. This is especially true when dealing with multiple alleles and large gene families where genetic redundancy can complicate knockout of all family members. miRNAs are small noncoding RNAs that are often encoded in gene families and can display extensive genetic redundancy. This redundancy, plus their small size which can limit targeting by CRISPR/Cas9, makes morpholino-based strategies particularly attractive for inhibition of miRNA function. We provide the rationale, background, and methods to inhibit miRNA function with antisense morpholinos during early development and in the adult retina in zebrafish.

Rapid Construction of Multiplexed CRISPR-Cas9 Systems for Plant Genome Editing.

Multiplex CRISPR-Cas9 nuclease mediated genome editing is an efficient method for disrupting gene function in plants. Use of CRISPR-Cas9 has escalated rapidly in recent years and is expected to become routine practice in molecular biology and related fields of research. Due to the relatively novel and widespread adoption of this technology, first-time users may not have regular access to experienced guidance or technical support from peers or mentors. Here, we offer guidance and technical support in the form of a detailed and tested protocol for simultaneous targeting of three separate loci on the TRANSPARENT TESTA 4 (TT4) gene in Arabidopsis thaliana using multiplex CRISPR-Cas9. Although we target multiple loci on a single gene in Arabidopsis, the same approach can be used to target multiple genes or alleles in other plant species as well. We recommend the use of a molecular toolkit to streamline the process and make recommendations for this type of approach. The protocol starts with an overview of the reagents and covers designing of gRNAs and assembly of components into a final T-DNA delivery molecule through Golden Gate cloning and Multisite Gateway LR recombination.