The chromosome-scale reference genome of mirid bugs (Adelphocoris suturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation.

BACKGROUND: Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. RESULTS: Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. CONCLUSIONS: The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution.

BACKGROUND: Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism. RESULTS: Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa. CONCLUSIONS: We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.

An enhanced photosynthesis and carbohydrate metabolic capability contributes to heterosis of the cotton (Gossypium hirsutum) hybrid 'Huaza Mian H318', as revealed by genome-wide gene expression analysis.

BACKGROUND: Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton. RESULTS: The elite cotton hybrid variety 'Huaza Mian H318' (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4-5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318. CONCLUSION: The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

Orchestration of plant development and defense by indirect crosstalk of salicylic acid and brassinosteorid signaling via transcription factor GhTINY2.

Salicylic acid (SA) and brassinosteroids (BRs) are well known to regulate diverse processes of plant development and stress responses, but the mechanisms by which these phytohormones mediate the growth and defense trade-off are largely unclear. In addition, little is known about the roles of DEHYDRATION RESPONSIVE ELEMENT BINDING transcription factors, especially in biotic stress and plant growth. Here, we identified a cotton (Gossypium hirsutum) APETALA2/ETHYLENE RESPONSIVE FACTOR gene GhTINY2 that is strongly induced by Verticillium dahliae. Overexpression of GhTINY2 in cotton and Arabidopsis enhanced tolerance to V. dahliae, while knockdown of expression increased the susceptibility of cotton to the pathogen. GhTINY2 was found to promote SA accumulation and SA signaling transduction by directly activating expression of WRKY51. Moreover, GhTINY2-overexpressing cotton and Arabidopsis showed retardation of growth, increased sensitivity to inhibitors of BR biosynthesis, down-regulation of several BR-induced genes, and up-regulation of BR-repressed genes, while GhTINY2-RNAi cotton showed the opposite effects. We further determined that GhTINY2 negatively regulates BR signaling by interacting with BRASSINAZOLE-RESISTANT 1 (BZR1) and restraining its transcriptional activation of the expression of INDOLE-3-ACETIC ACID INDUCIBLE 19 (IAA19). These findings indicate that GhTINY2 fine-tunes the trade-off between immunity and growth via indirect crosstalk between WRKY51-mediated SA biosynthesis and BZR1-IAA19-regulated BR signaling.

Multi-omics maps of cotton fibre reveal epigenetic basis for staged single-cell differentiation.

Epigenetic modifications are highlighted for their great importance in regulating plant development, but their function associated with single-cell differentiation remains undetermined. Here, we used the cotton fibre, which is the epidermal hair on the cotton ovule, as a model to investigate the regulatory role of DNA methylation in cell differentiation. The level of CHH (H = A, T, or C) DNA methylation level was found to increase during fibre development, accompanied by a decrease in RNA-directed DNA methylation (RdDM). Examination of nucleosome positioning revealed a gradual transition from euchromatin to heterochromatin for chromatin dynamics in developing fibres, which could shape the DNA methylation landscape. The observed increase in DNA methylation in fibres, compared with other ovule tissue, was demonstrated to be mediated predominantly by an active H3K9me2-dependent pathway rather than the RdDM pathway, which was inactive. Furthermore, integrated multi-omics analyses revealed that dynamic DNA methylation played a role in the regulation of lipid biosynthesis and spatio-temporal modulation of reactive oxygen species during fibre differentiation. Our study illustrates two divergent pathways mediating a continuous increase of DNA methylation and also sheds further light on the epigenetic basis for single-cell differentiation in plants. These data and analyses are made available to the wider research community through a comprehensive web portal.

GbEXPATR, a species-specific expansin, enhances cotton fibre elongation through cell wall restructuring.

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre-specific α-expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the DT genome) is a classical α-expansin, while its homoeolog, GbEXPATR (AT genome), encodes a truncated protein lacking the normal C-terminal polysaccharide-binding domain of other α-expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall-associated genes, especially chitinase-like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α-expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.

Long noncoding RNAs and their proposed functions in fibre development of cotton (Gossypium spp.).

Long noncoding RNAs (lncRNAs) are transcripts of at least 200 bp in length, possess no apparent coding capacity and are involved in various biological regulatory processes. Until now, no systematic identification of lncRNAs has been reported in cotton (Gossypium spp.). Here, we describe the identification of 30 550 long intergenic noncoding RNA (lincRNA) loci (50 566 transcripts) and 4718 long noncoding natural antisense transcript (lncNAT) loci (5826 transcripts). LncRNAs are rich in repetitive sequences and preferentially expressed in a tissue-specific manner. The detection of abundant genome-specific and/or lineage-specific lncRNAs indicated their weak evolutionary conservation. Approximately 76% of homoeologous lncRNAs exhibit biased expression patterns towards the At or Dt subgenomes. Compared with protein-coding genes, lncRNAs showed overall higher methylation levels and their expression was less affected by gene body methylation. Expression validation in different cotton accessions and coexpression network construction helped to identify several functional lncRNA candidates involved in cotton fibre initiation and elongation. Analysis of integrated expression from the subgenomes of lncRNAs generating miR397 and its targets as a result of genome polyploidization indicated their pivotal functions in regulating lignin metabolism in domesticated tetraploid cotton fibres. This study provides the first comprehensive identification of lncRNAs in Gossypium.

Cotton GhCKI disrupts normal male reproduction by delaying tapetum programmed cell death via inactivating starch synthase.

Anther infertility under high temperature (HT) conditions is a critical factor contributing to yield loss in cotton (Gossypium hirsutum). Using large-scale expression profile sequencing, we studied the effect of HT on cotton anther development. Our analysis revealed that altered carbohydrate metabolism or disrupted tapetal programmed cell death (PCD) underlie anther sterility. Expression of the Gossypium hirsutum casein kinase I (GhCKI) gene, which encodes a homolog of casein kinase I (CKI), was induced in an HT-sensitive cotton line after exposure to HT. As mammalian homologs of GhCKI are involved in inactivation of glycogen synthase and the regulation of apoptosis, GhCKI may be considered a target gene for improving anther fertility under HT conditions. Our studies suggest that GhCKI exhibits starch synthase kinase activity, increases glucose content in early-stage buds and activates the accumulation of abscisic acid, thereby disturbing the balance of reactive oxygen species and eventually disrupting tapetal PCD, leading to anther abortion or indehiscence. These results indicate that GhCKI may be a key regulator of tapetal PCD and anther dehiscence, with the potential to facilitate regulation of HT tolerance in crops.

Comparative transcriptome analysis between somatic embryos (SEs) and zygotic embryos in cotton: evidence for stress response functions in SE development.

As a product of asexual reproduction in plants, the somatic embryo (SE) differentiates into a new plantlet via a zygotic embryogenesis-like process. Here, we present the phenotypic and cellular differences between SEs and zygotic embryos (ZEs) revealed by histological section scanning using three parallel development stages of the two types of embryos of cotton (Gossypium hirsutum cv. YZ1), including globular, torpedo and cotyledonary-stages. To identify the molecular characteristics of SE development in cotton, the digital gene expression system was used to profile the genes active during SE and ZE development. A total of 4242 differentially expressed genes (DEGs) were identified in at least one developmental stage. Expression pattern and functional classification analysis based on these DEGs reveals that SE development exhibits a transcriptional activation of stress responses. RT-PCR analysis further confirmed enhanced expression levels of stress-related genes in SEs than in ZEs. Experimental stress treatment, induced by NaCl and ABA, accelerated SE development and increased the transcription of genes related to stress response, in parallel with decelerated proliferation of embryogenic calluses under stress treatment. Our data reveal that SE development involves the activation of stress responses, which we suggest may regulate the balance between cell proliferation and differentiation. These results provide new insight into the molecular mechanisms of SE development and suggest strategies that can be used for regulating the developmental processes of somatic embryogenesis.

Functional characterization of cotton genes responsive to Verticillium dahliae through bioinformatics and reverse genetics strategies.

Verticillium wilt causes dramatic cotton yield loss in China. Although some genes or biological processes involved in the interaction between cotton and Verticillium dahliae have been identified, the molecular mechanism of cotton resistance to this disease is still poorly understood. The basic innate immune response for defence is somewhat conserved among plant species to defend themselves in complex environments, which makes it possible to characterize genes involved in cotton immunity based on information from model plants. With the availability of Arabidopsis databases, a data-mining strategy accompanied by virus-induced gene silencing (VIGS) and heterologous expression were adopted in cotton and tobacco, respectively, for global screening and gene function characterization. A total of 232 Arabidopsis genes putatively involved in basic innate immunity were screened as candidate genes, and bioinformatic analysis suggested a role of these genes in the immune response. In total, 38 homologous genes from cotton were singled out to characterize their response to V. dahliae and methyl jasmonate treatment through quantitative real-time PCR. The results revealed that 24 genes were differentially regulated by pathogen inoculation, and most of these genes responded to both Verticillium infection and jasmonic acid stimuli. Furthermore, the efficiency of the strategy was illustrated by the functional identification of six candidate genes via heterologous expression in tobacco or a knock-down approach using VIGS in cotton. Functional categorization of these 24 differentially expressed genes as well as functional analysis suggest that reactive oxygen species, salicylic acid- and jasmonic acid-signalling pathways are involved in the cotton disease resistance response to V. dahliae. Our data demonstrate how information from model plants can allow the rapid translation of information into non-model species without complete genome sequencing, via high-throughput screening and functional identification of target genes based on data-mining and VIGS.

Small RNA and degradome sequencing reveal complex miRNA regulation during cotton somatic embryogenesis.

MicroRNAs (miRNAs) are endogenous non-coding ~21 nucleotide RNAs that regulate gene expression at the transcriptional and post-transcriptional levels in plants and animals. They play an important role in development, abiotic stress, and pathogen responses. miRNAs with their targets have been widely studied in model plants, but limited knowledge is available on the small RNA population of cotton (Gossypium hirsutum)-an important economic crop, and global identification of related targets through degradome sequencing has not been developed previously. In this study, small RNAs and their targets were identified during cotton somatic embryogenesis (SE) through high-throughput small RNA and degradome sequencing, comparing seedling hypocotyl and embryogenic callus (EC) of G. hirsutum YZ1. A total of 36 known miRNA families were found to be differentially expressed, of which 19 miRNA families were represented by 29 precursors. Twenty-five novel miRNAs were identified. A total of 234 transcripts in EC and 322 transcripts in control (CK) were found to be the targets of 23 and 30 known miRNA families, respectively, and 16 transcripts were targeted by eight novel miRNAs. Interestingly, four trans-acting small interfering RNAs (tas3-siRNAs) were also found in degradome libraries, three of which perfectly matched their precursors. Several targets were further validated via RNA ligase-mediated rapid amplification of 5' cDNA ends (RLM 5'-RACE). The profiling of the miRNAs and their target genes provides new information on the miRNAs network during cotton SE.

Genome and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense.

INTRODUCTION: Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare. OBJECTIVES: We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton. METHODS: Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton. RESULTS: One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian. CONCLUSION: Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Single-cell resolution analysis reveals the preparation for reprogramming the fate of stem cell niche in cotton lateral meristem.

BACKGROUND: Somatic embryogenesis is a major process for plant regeneration. However, cell communication and the gene regulatory network responsible for cell reprogramming during somatic embryogenesis are still largely unclear. Recent advances in single-cell technologies enable us to explore the mechanism of plant regeneration at single-cell resolution. RESULTS: We generate a high-resolution single-cell transcriptomic landscape of hypocotyl tissue from the highly regenerable cotton genotype Jin668 and the recalcitrant TM-1. We identify nine putative cell clusters and 23 cluster-specific marker genes for both cultivars. We find that the primary vascular cell is the major cell type that undergoes cell fate transition in response to external stimulation. Further developmental trajectory and gene regulatory network analysis of these cell clusters reveals that a total of 41 hormone response-related genes, including LAX2, LAX1, and LOX3, exhibit different expression patterns in the primary xylem and cambium region of Jin668 and TM-1. We also identify novel genes, including CSEF, PIS1, AFB2, ATHB2, PLC2, and PLT3, that are involved in regeneration. We demonstrate that LAX2, LAX1 and LOX3 play important roles in callus proliferation and plant regeneration by CRISPR/Cas9 editing and overexpression assay. CONCLUSIONS: This study provides novel insights on the role of the regulatory network in cell fate transition and reprogramming during plant regeneration driven by somatic embryogenesis.

Regulatory controls of duplicated gene expression during fiber development in allotetraploid cotton.

Polyploidy complicates transcriptional regulation and increases phenotypic diversity in organisms. The dynamics of genetic regulation of gene expression between coresident subgenomes in polyploids remains to be understood. Here we document the genetic regulation of fiber development in allotetraploid cotton Gossypium hirsutum by sequencing 376 genomes and 2,215 time-series transcriptomes. We characterize 1,258 genes comprising 36 genetic modules that control staged fiber development and uncover genetic components governing their partitioned expression relative to subgenomic duplicated genes (homoeologs). Only about 30% of fiber quality-related homoeologs show phenotypically favorable allele aggregation in cultivars, highlighting the potential for subgenome additivity in fiber improvement. We envision a genome-enabled breeding strategy, with particular attention to 48 favorable alleles related to fiber phenotypes that have been subjected to purifying selection during domestication. Our work delineates the dynamics of gene regulation during fiber development and highlights the potential of subgenomic coordination underpinning phenotypes in polyploid plants.

Transcript profiling reveals complex auxin signalling pathway and transcription regulation involved in dedifferentiation and redifferentiation during somatic embryogenesis in cotton.

BACKGROUND: Somatic embryogenesis (SE), by which somatic cells of higher plants can dedifferentiate and reorganize into new plants, is a notable illustration of cell totipotency. However, the precise molecular mechanisms regulating SE remain unclear. To characterize the molecular events of this unique process, transcriptome analysis, in combination with biochemical and histological approaches, were conducted in cotton, a typical plant species in SE. Genome-wide profiling of gene expression allowed the identification of novel molecular markers characteristic of this developmental process. RESULTS: RNA-Seq was used to identify 5,076 differentially expressed genes during cotton SE. Expression profile and functional assignments of these genes indicated significant transcriptional complexity during this process, associated with morphological, histological changes and endogenous indole-3-acetic acid (IAA) alteration. Bioinformatics analysis showed that the genes were enriched for basic processes such as metabolic pathways and biosynthesis of secondary metabolites. Unigenes were abundant for the functions of protein binding and hydrolase activity. Transcription factor-encoding genes were found to be differentially regulated during SE. The complex pathways of auxin abundance, transport and response with differentially regulated genes revealed that the auxin-related transcripts belonged to IAA biosynthesis, indole-3-butyric acid (IBA) metabolism, IAA conjugate metabolism, auxin transport, auxin-responsive protein/indoleacetic acid-induced protein (Aux/IAA), auxin response factor (ARF), small auxin-up RNA (SAUR), Aux/IAA degradation, and other auxin-related proteins, which allow an intricate system of auxin utilization to achieve multiple purposes in SE. Quantitative real-time PCR (qRT-PCR) was performed on selected genes with different expression patterns and functional assignments were made to demonstrate the utility of RNA-Seq for gene expression profiles during cotton SE. CONCLUSION: We report here the first comprehensive analysis of transcriptome dynamics that may serve as a gene expression profile blueprint in cotton SE. Our main goal was to adapt the RNA-Seq technology to this notable development process and to analyse the gene expression profile. Complex auxin signalling pathway and transcription regulation were highlighted. Together with biochemical and histological approaches, this study provides comprehensive gene expression data sets for cotton SE that serve as an important platform resource for further functional studies in plant embryogenesis.

Increasing cotton genome coverage with polymorphic SSRs as revealed by SSCP.

Simple sequence repeat (SSR) markers are widely used in plant genetics and breeding. However, there are many SSR markers that do not reveal polymorphism in cotton. Traditional SSR genotyping methods only provide information on product sizes. This leaves many marker polymorphism undetected, thus, lowering the utility of SSRs. In the present study, monomorphic SSRs between two mapping parents, 'Emian22' and 3-79, were subjected to single-strand conformation polymorphism (SSCP) analysis to reveal polymorphism. Of the 4194 monomorphic SSR primer pairs, 158 pairs (3.77%) showed polymorphism and revealed 174 polymorphic loci. Sequence analysis showed that the differences in PCR products between the mapping parents were solely due to base transition or transversion, which was in agreement with SSCP principles. SSCP also revealed SSRs with motifs of AT/TA and GAA/CTT were more polymorphic in dinucleotides and trinucleotides, respectively. Genetic mapping integrated 160 loci into our interspecific BC(1) linkage map, 5 of which associated with QTLs related to cotton fiber quality. The technique discussed in the present study enables us to detect polymorphism of monomorphic SSRs, and increase the utilization efficiency of the existing SSR primers.

Dual Domestication, Diversity, and Differential Introgression in Old World Cotton Diploids.

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 6 × 10-3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

Highly Efficient Genome Editing Using Geminivirus-Based CRISPR/Cas9 System in Cotton Plant.

Upland cotton (Gossypium hirsutum), an allotetraploid, contains At- and Dt- subgenome and most genes have multiple homologous copies, which pose a huge challenge to investigate genes' function due to the functional redundancy. Therefore, it is of great significance to establish effective techniques for the functional genomics in cotton. In this study, we tested two novel genome editing vectors and compared them with the CRISPR/Cas9 system (pRGEB32-GhU6.7) developed in our laboratory previously. In the first new vector, the sgRNA transcription unite was constructed into the replicon (LIR-Donor-SIR-Rep-LIR) of the bean yellow dwarf virus (BeYDV) and named as pBeYDV-Cas9-KO and in the second vector, the ubiquitin promoter that drives Cas9 protein was replaced with a constitutive CaMV 35S promoter and defined as pRGEB32-35S. The results from transgenic cotton calli/plants revealed that pBeYDV-Cas9-KO vector showed the highest editing efficiency of GhCLA1 in At and Dt subgenomes edited simultaneously up to 73.3% compared to the 44.6% of pRGEB32-GhU6.7 and 51.2% of pRGEB32-35S. The editing efficiency of GhCLA1 in At and Dt subgenome by pBeYDV-Cas9-KO was 85.7% and 97.2%, respectively, whereas the efficiency by pRGEB32-GhU6.7 and pRGEB32-35S vectors was 67.7%, 86.5%, 84%, and 87.2%, respectively. The editing profile of pBeYDV-Cas9-KO was mainly composed of fragment deletion, accounting for 84.0% and ranging 1-10 bp in length. The main editing sites are located at positions 11-17 upstream of PAM site. The off-target effects were not detected in all potential off-target sites. Taken together, the pBeYDV-Cas9-KO system has high editing efficiency and specificity with wide editing range than the traditional CRISPR/Cas9 system, which provides a powerful tool for cotton functional genomics research and molecular breeding.