Novel rice iron biofortification approaches using expression of ZmYS1 and OsTOM1 controlled by tissue-specific promoters.

Intrinsic improvement of iron (Fe) concentration in rice grains, called rice Fe biofortification, is a promising countermeasure against widespread human Fe deficiency. In this study, two novel rice Fe biofortification approaches are reported. The first approach (Y approach) involved the expression of maize YELLOW STRIPE 1 controlled by the HEAVY METAL ATPASE 2 promoter. The Y approach increased the polished grain Fe concentrations up to 4.8-fold compared with the non-transgenic (NT) line. The second approach (T approach) involved the expression of rice TRANSPORTER OF MUGINEIC ACID 1 controlled by the FERRIC REDUCTASE DEFECTIVE LIKE 1 promoter. The T approach increased the polished grain Fe concentrations by up to 3.2-fold. No synergistic increases in the polished grain Fe concentrations were observed when Y and T approaches were combined (YT approach). However, the polished grain Fe concentrations further increased by 5.1- to 9.3-fold compared with the NT line, when YT approach was combined with the endosperm-specific expression of FERRITIN (YTF approach), or when YTF approach was combined with the constitutive expression of NICOTIANAMINE SYNTHASE (YTFN approach). Total grain weight per plant in most Y, T, YT, and YTFN lines was comparable to that in the NT line, while it was significantly decreased in most YTF lines. The novel approaches reported in this study expand the portfolio of genetic engineering strategies that can be used for Fe biofortification in rice.

The Cassava Source-Sink project: opportunities and challenges for crop improvement by metabolic engineering.

Cassava (Manihot esculenta Crantz) is one of the important staple foods in Sub-Saharan Africa. It produces starchy storage roots that provide food and income for several hundred million people, mainly in tropical agriculture zones. Increasing cassava storage root and starch yield is one of the major breeding targets with respect to securing the future food supply for the growing population of Sub-Saharan Africa. The Cassava Source-Sink (CASS) project aims to increase cassava storage root and starch yield by strategically integrating approaches from different disciplines. We present our perspective and progress on cassava as an applied research organism and provide insight into the CASS strategy, which can serve as a blueprint for the improvement of other root and tuber crops. Extensive profiling of different field-grown cassava genotypes generates information for leaf, phloem, and root metabolic and physiological processes that are relevant for biotechnological improvements. A multi-national pipeline for genetic engineering of cassava plants covers all steps from gene discovery, cloning, transformation, molecular and biochemical characterization, confined field trials, and phenotyping of the seasonal dynamics of shoot traits under field conditions. Together, the CASS project generates comprehensive data to facilitate conventional breeding strategies for high-yielding cassava genotypes. It also builds the foundation for genome-scale metabolic modelling aiming to predict targets and bottlenecks in metabolic pathways. This information is used to engineer cassava genotypes with improved source-sink relations and increased yield potential.

Diurnal dynamics of the Arabidopsis rosette proteome and phosphoproteome.

Plant growth depends on the diurnal regulation of cellular processes, but it is not well understood if and how transcriptional regulation controls diurnal fluctuations at the protein level. Here, we report a high-resolution Arabidopsis thaliana (Arabidopsis) leaf rosette proteome acquired over a 12 hr light:12 hr dark diurnal cycle and the phosphoproteome immediately before and after the light-to-dark and dark-to-light transitions. We quantified nearly 5,000 proteins and 800 phosphoproteins, of which 288 fluctuated in their abundance and 226 fluctuated in their phosphorylation status. Of the phosphoproteins, 60% were quantified for changes in protein abundance. This revealed six proteins involved in nitrogen and hormone metabolism that had concurrent changes in both protein abundance and phosphorylation status. The diurnal proteome and phosphoproteome changes involve proteins in key cellular processes, including protein translation, light perception, photosynthesis, metabolism and transport. The phosphoproteome at the light-dark transitions revealed the dynamics at phosphorylation sites in either anticipation of or response to a change in light regime. Phosphorylation site motif analyses implicate casein kinase II and calcium/calmodulin-dependent kinases among the primary light-dark transition kinases. The comparative analysis of the diurnal proteome and diurnal and circadian transcriptome established how mRNA and protein accumulation intersect in leaves during the diurnal cycle of the plant.

A new full-length circular DNA sequencing method for viral-sized genomes reveals that RNAi transgenic plants provoke a shift in geminivirus populations in the field.

We present a new method, CIDER-Seq (Circular DNA Enrichment sequencing) for the unbiased enrichment and long-read sequencing of viral-sized circular DNA molecules. We used CIDER-Seq to produce single-read full-length virus genomes for the first time. CIDER-Seq combines PCR-free virus enrichment with Single Molecule Real Time sequencing and a new sequence de-concatenation algorithm. We apply our technique to produce >1200 full-length, highly accurate geminivirus genomes from RNAi-transgenic and control plants in a field trial in Kenya. Using CIDER-Seq we can demonstrate for the first time that the expression of antiviral double-stranded RNA (dsRNA) in transgenic plants causes a consistent shift in virus populations towards species sharing low homology to the transgene derived dsRNA. Our method and its application in an economically important crop plant opens new possibilities in periodic virus sequence surveillance and accurate profiling of diverse circular DNA elements.

Diurnal changes in concerted plant protein phosphorylation and acetylation in Arabidopsis organs and seedlings.

Protein phosphorylation and acetylation are the two most abundant post-translational modifications (PTMs) that regulate protein functions in eukaryotes. In plants, these PTMs have been investigated individually; however, their co-occurrence and dynamics on proteins is currently unknown. Using Arabidopsis thaliana, we quantified changes in protein phosphorylation, acetylation and protein abundance in leaf rosettes, roots, flowers, siliques and seedlings at the end of day (ED) and at the end of night (EN). This identified 2549 phosphorylated and 909 acetylated proteins, of which 1724 phosphorylated and 536 acetylated proteins were also quantified for changes in PTM abundance between ED and EN. Using a sequential dual-PTM workflow, we identified significant PTM changes and intersections in these organs and plant developmental stages. In particular, cellular process-, pathway- and protein-level analyses reveal that the phosphoproteome and acetylome predominantly intersect at the pathway- and cellular process-level at ED versus EN. We found 134 proteins involved in core plant cell processes, such as light harvesting and photosynthesis, translation, metabolism and cellular transport, that were both phosphorylated and acetylated. Our results establish connections between PTM motifs, PTM catalyzing enzymes and putative substrate networks. We also identified PTM motifs for further characterization of the regulatory mechanisms that control cellular processes during the diurnal cycle in different Arabidopsis organs and seedlings. The sequential dual-PTM analysis expands our understanding of diurnal plant cell regulation by PTMs and provides a useful resource for future analyses, while emphasizing the importance of analyzing multiple PTMs simultaneously to elucidate when, where and how they are involved in plant cell regulation.

Enhancement of vitamin B6 levels in rice expressing Arabidopsis vitamin B6 biosynthesis de novo genes.

Vitamin B6 (pyridoxine) is vital for key metabolic reactions and reported to have antioxidant properties in planta. Therefore, enhancement of vitamin B6 content has been hypothesized to be a route to improve resistance to biotic and abiotic stresses. Most of the current studies on vitamin B6 in plants are on eudicot species, with monocots remaining largely unexplored. In this study, we investigated vitamin B6 biosynthesis in rice, with a view to examining the feasibility and impact of enhancing vitamin B6 levels. Constitutive expression in rice of two Arabidopsis thaliana genes from the vitamin B6 biosynthesis de novo pathway, AtPDX1.1 and AtPDX2, resulted in a considerable increase in vitamin B6 in leaves (up to 28.3-fold) and roots (up to 12-fold), with minimal impact on general growth. Rice lines accumulating high levels of vitamin B6 did not display enhanced tolerance to abiotic stress (salt) or biotic stress (resistance to Xanthomonas oryzae infection). While a significant increase in vitamin B6 content could also be achieved in rice seeds (up to 3.1-fold), the increase was largely due to its accumulation in seed coat and embryo tissues, with little enhancement observed in the endosperm. However, seed yield was affected in some vitamin B6 -enhanced lines. Notably, expression of the transgenes did not affect the expression of the endogenous rice PDX genes. Intriguingly, despite transgene expression in leaves and seeds, the corresponding proteins were only detectable in leaves and could not be observed in seeds, possibly pointing to a mode of regulation in this organ.

Haplotype-resolved genomes of geminivirus-resistant and geminivirus-susceptible African cassava cultivars.

BACKGROUND: Cassava is an important food crop in tropical and sub-tropical regions worldwide. In Africa, cassava production is widely affected by cassava mosaic disease (CMD), which is caused by the African cassava mosaic geminivirus that is transmitted by whiteflies. Cassava breeders often use a single locus, CMD2, for introducing CMD resistance into susceptible cultivars. The CMD2 locus has been genetically mapped to a 10-Mbp region, but its organization and genes as well as their functions are unknown. RESULTS: We report haplotype-resolved de novo assemblies and annotations of the genomes for the African cassava cultivar TME (tropical Manihot esculenta), which is the origin of CMD2, and the CMD-susceptible cultivar 60444. The assemblies provide phased haplotype information for over 80% of the genomes. Haplotype comparison identified novel features previously hidden in collapsed and fragmented cassava genomes, including thousands of allelic variants, inter-haplotype diversity in coding regions, and patterns of diversification through allele-specific expression. Reconstruction of the CMD2 locus revealed a highly complex region with nearly identical gene sets but limited microsynteny between the two cultivars. CONCLUSIONS: The genome maps of the CMD2 locus in both 60444 and TME3, together with the newly annotated genes, will help the identification of the causal genetic basis of CMD2 resistance to geminiviruses. Our de novo cassava genome assemblies will also facilitate genetic mapping approaches to narrow the large CMD2 region to a few candidate genes for better informed strategies to develop robust geminivirus resistance in susceptible cassava cultivars.

Targeting intracellular transport combined with efficient uptake and storage significantly increases grain iron and zinc levels in rice.

Rice, a staple food for more than half of the world population, is an important target for iron and zinc biofortification. Current strategies mainly focus on the expression of genes for efficient uptake, long-distance transport and storage. Targeting intracellular iron mobilization to increase grain iron levels has not been reported. Vacuole is an important cell compartment for iron storage and the NATURAL RESISTANCE ASSOCIATED MACROPHAGE PROTEIN (NRAMP) family of transporters export iron from vacuoles to cytosol when needed. We developed transgenic Nipponbare rice lines expressing AtNRAMP3 under the control of the UBIQUITIN or rice embryo/aleurone-specific 18-kDa Oleosin (Ole18) promoter together with NICOTIANAMINE SYNTHASE (AtNAS1) and FERRITIN (PvFER), or expressing only AtNRAMP3 and PvFER together. Iron and zinc were increased close to recommended levels in polished grains of the transformed lines, with maximum levels when AtNRAMP3, AtNAS1 and PvFER were expressed together (12.67 µg/g DW iron and 45.60 µg/g DW zinc in polished grains of line NFON16). Similar high iron and zinc levels were obtained in transgenic Indica IR64 lines expressing the AtNRAMP3, AtNAS1 and PvFER cassette (13.65 µg/g DW iron and 48.18 µg/g DW zinc in polished grains of line IR64_1), equalling more than 90% of the recommended iron increase in rice endosperm. Our results demonstrate that targeting intracellular iron stores in combination with iron and zinc transport and endosperm storage is an effective strategy for iron biofortification. The increases achieved in polished IR64 grains are of dietary relevance for human health and a valuable nutrition trait for breeding programmes.

Photoperiodic control of the Arabidopsis proteome reveals a translational coincidence mechanism.

Plants respond to seasonal cues such as the photoperiod, to adapt to current conditions and to prepare for environmental changes in the season to come. To assess photoperiodic responses at the protein level, we quantified the proteome of the model plant Arabidopsis thaliana by mass spectrometry across four photoperiods. This revealed coordinated changes of abundance in proteins of photosynthesis, primary and secondary metabolism, including pigment biosynthesis, consistent with higher metabolic activity in long photoperiods. Higher translation rates in the day than the night likely contribute to these changes, via an interaction with rhythmic changes in RNA abundance. Photoperiodic control of protein levels might be greatest only if high translation rates coincide with high transcript levels in some photoperiods. We term this proposed mechanism "translational coincidence", mathematically model its components, and demonstrate its effect on the Arabidopsis proteome. Datasets from a green alga and a cyanobacterium suggest that translational coincidence contributes to seasonal control of the proteome in many phototrophic organisms. This may explain why many transcripts but not their cognate proteins exhibit diurnal rhythms.

Symplasmic phloem unloading and radial post-phloem transport via vascular rays in tuberous roots of Manihot esculenta.

Cassava (Manihot esculenta) is one of the most important staple food crops worldwide. Its starchy tuberous roots supply over 800 million people with carbohydrates. Yet, surprisingly little is known about the processes involved in filling of those vital storage organs. A better understanding of cassava carbohydrate allocation and starch storage is key to improving storage root yield. Here, we studied cassava morphology and phloem sap flow from source to sink using transgenic pAtSUC2::GFP plants, the phloem tracers esculin and 5(6)-carboxyfluorescein diacetate, as well as several staining techniques. We show that cassava performs apoplasmic phloem loading in source leaves and symplasmic unloading into phloem parenchyma cells of tuberous roots. We demonstrate that vascular rays play an important role in radial transport from the phloem to xylem parenchyma cells in tuberous roots. Furthermore, enzymatic and proteomic measurements of storage root tissues confirmed high abundance and activity of enzymes involved in the sucrose synthase-mediated pathway and indicated that starch is stored most efficiently in the outer xylem layers of tuberous roots. Our findings form the basis for biotechnological approaches aimed at improved phloem loading and enhanced carbohydrate allocation and storage in order to increase tuberous root yield of cassava.

BRR2a Affects Flowering Time via FLC Splicing.

Several pathways control time to flowering in Arabidopsis thaliana through transcriptional and posttranscriptional gene regulation. In recent years, mRNA processing has gained interest as a critical regulator of flowering time control in plants. However, the molecular mechanisms linking RNA splicing to flowering time are not well understood. In a screen for Arabidopsis early flowering mutants we identified an allele of BRR2a. BRR2 proteins are components of the spliceosome and highly conserved in eukaryotes. Arabidopsis BRR2a is ubiquitously expressed in all analyzed tissues and involved in the processing of flowering time gene transcripts, most notably FLC. A missense mutation of threonine 895 in BRR2a caused defects in FLC splicing and greatly reduced FLC transcript levels. Reduced FLC expression increased transcription of FT and SOC1 leading to early flowering in both short and long days. Genome-wide experiments established that only a small set of introns was not correctly spliced in the brr2a mutant. Compared to control introns, retained introns were often shorter and GC-poor, had low H3K4me1 and CG methylation levels, and were often derived from genes with a high-H3K27me3-low-H3K36me3 signature. We propose that BRR2a is specifically needed for efficient splicing of a subset of introns characterized by a combination of factors including intron size, sequence and chromatin, and that FLC is most sensitive to splicing defects.

Genome-scale analysis of regulatory protein acetylation enzymes from photosynthetic eukaryotes.

BACKGROUND: Reversible protein acetylation occurring on Lys-Ne has emerged as a key regulatory post-translational modification in eukaryotes. It is mediated by two groups of enzymes: lysine acetyltransferases (KATs) and lysine deacetylases (KDACs) that catalyze the addition and removal of acetyl groups from target proteins. Estimates indicate that protein acetylation is second to protein phosphorylation in abundance, with thousands of acetylated sites now identified in different subcellular compartments. Considering the important regulatory role of protein phosphorylation, elucidating the diversity of KATs and KDACs across photosynthetic eukaryotes is essential in furthering our understanding of the impact of reversible protein acetylation on plant cell processes. RESULTS: We report a genome-scale analysis of lysine acetyltransferase (KAT)- and lysine deacetylase (KDAC)-families from 53 photosynthetic eukaryotes. KAT and KDAC orthologs were identified in sequenced genomes ranging from glaucophytes and algae to land plants and then analyzed for evolutionary relationships. Based on consensus molecular phylogenetic and subcellular localization data we found new sub-classes of enzymes in established KAT- and KDAC-families. Specifically, we identified a non-photosynthetic origin of the HD-tuin family KDACs, a new monocot-specific Class I HDA-family sub-class, and a phylogenetically distinct Class II algal/heterokont sub-class which maintains an ankyrin domain not conserved in land plant Class II KDACs. Protein structure analysis showed that HDA- and SRT-KDACs exist as bare catalytic subunits with highly conserved median protein length, while all KATs maintained auxiliary domains, with CBP- and TAFII250-KATs displaying protein domain gain and loss over the course of photosynthetic eukaryote evolution in addition to variable protein length. Lastly, promoter element enrichment analyses across species revealed conserved cis-regulatory sequences that support KAT and KDAC involvement in the regulation of plant development, cold/drought stress response, as well as cellular processes such as the circadian clock. CONCLUSIONS: Our results reveal new evolutionary, structural, and biological insights into the KAT- and KDAC-families of photosynthetic eukaryotes, including evolutionary parallels to protein kinases and protein phosphatases. Further, we provide a comprehensive annotation framework through our extensive phylogenetic analysis, from which future research investigating aspects of protein acetylation in plants can use to position new findings in a broader context.

Arabidopsis MSI1 connects LHP1 to PRC2 complexes.

Polycomb group (PcG) proteins form essential epigenetic memory systems for controlling gene expression during development in plants and animals. However, the mechanism of plant PcG protein functions remains poorly understood. Here, we probed the composition and function of plant Polycomb repressive complex 2 (PRC2). This work established the fact that all known plant PRC2 complexes contain MSI1, a homologue of Drosophila p55. While p55 is not essential for the in vitro enzymatic activity of PRC2, plant MSI1 was required for the functions of the EMBRYONIC FLOWER and the VERNALIZATION PRC2 complexes including trimethylation of histone H3 Lys27 (H3K27) at the target chromatin, as well as gene repression and establishment of competence to flower. We found that MSI1 serves to link PRC2 to LIKE HETEROCHROMATIN PROTEIN 1 (LHP1), a protein that binds H3K27me3 in vitro and in vivo and is required for a functional plant PcG system. The LHP1-MSI1 interaction forms a positive feedback loop to recruit PRC2 to chromatin that carries H3K27me3. Consequently, this can provide a mechanism for the faithful inheritance of local epigenetic information through replication.

H3K36ac Is an Evolutionary Conserved Plant Histone Modification That Marks Active Genes.

In eukaryotic cells, histones are subject to a large number of posttranslational modifications whose sequential or combinatorial action affects chromatin structure and genome function. We identified acetylation at Lys-36 in histone H3 (H3K36ac) as a new chromatin modification in plants. The H3K36ac modification is evolutionary conserved in seed plants, including the gymnosperm Norway spruce (Picea abies) and the angiosperms rice (Oryza sativa), tobacco (Nicotiana tabacum), and Arabidopsis (Arabidopsis thaliana). In Arabidopsis, H3K36ac is highly enriched in euchromatin but not in heterochromatin. Genome-wide chromatin immunoprecipitation sequencing experiments revealed that H3K36ac peaks at the 5' end of genes, mainly on the two nucleosomes immediately distal to the transcription start site, independently of gene length. H3K36ac overlaps with H3K4me3 and the H2A.Z histone variant. The histone acetyl transferase GCN5 and the histone deacetylase HDA19 are required for H3K36ac homeostasis. H3K36ac and H3K36me3 show negative crosstalk, which is mediated by GCN5 and the histone methyl transferase SDG8. Although H3K36ac is associated with gene activity, we did not find a linear relationship between H3K36ac and transcript levels, suggesting that H3K36ac is a binary indicator of transcription.

Vitamin B1 diversity and characterization of biosynthesis genes in cassava.

Vitamin B1, which consists of the vitamers thiamin and its phosphorylated derivatives, is an essential micronutrient for all living organisms because it is required as a metabolic cofactor in several enzymatic reactions. Genetic diversity of vitamin B1 biosynthesis and accumulation has not been investigated in major crop species other than rice and potato. We analyzed cassava germplasm for accumulation of B1 vitamers. Vitamin B1 content in leaves and roots of 41 cassava accessions showed significant variation between accessions. HPLC analyses of B1 vitamers revealed distinct profiles in cassava leaves and storage roots, with nearly equal relative levels of thiamin pyrophosphate and thiamin monophosphate in leaves, but mostly thiamin pyrophosphate in storage roots. Unusually, the cassava genome has two genes encoding the 4-amino-2-methyl-5-hydroxymethylpyrimidine phosphate synthase, THIC (MeTHIC1 and MeTHIC2), both of which carry a riboswitch in the 3'-UTR, as well as the adenylated thiazole synthase, THI1 (MeTHI1a and MeTHI1b). The THIC and THI1 genes are expressed at very low levels in storage roots compared with the accumulation of vitamin B1, indicating only limited biosynthesis de novo therein. In leaves, vitamin B1 content is negatively correlated with THIC and THI1 expression levels, suggesting post-transcriptional regulation of THIC by the riboswitch present in the 3'-UTR of the THIC mRNA and regulation of THI1 by promoter activity or alternative post-transcriptional mechanisms.

Large-Scale Proteomics of the Cassava Storage Root and Identification of a Target Gene to Reduce Postharvest Deterioration.

Cassava (Manihot esculenta) is the most important root crop in the tropics, but rapid postharvest physiological deterioration (PPD) of the root is a major constraint to commercial cassava production. We established a reliable method for image-based PPD symptom quantification and used label-free quantitative proteomics to generate an extensive cassava root and PPD proteome. Over 2600 unique proteins were identified in the cassava root, and nearly 300 proteins showed significant abundance regulation during PPD. We identified protein abundance modulation in pathways associated with oxidative stress, phenylpropanoid biosynthesis (including scopoletin), the glutathione cycle, fatty acid α-oxidation, folate transformation, and the sulfate reduction II pathway. Increasing protein abundances and enzymatic activities of glutathione-associated enzymes, including glutathione reductases, glutaredoxins, and glutathione S-transferases, indicated a key role for ascorbate/glutathione cycles. Based on combined proteomics data, enzymatic activities, and lipid peroxidation assays, we identified glutathione peroxidase as a candidate for reducing PPD. Transgenic cassava overexpressing a cytosolic glutathione peroxidase in storage roots showed delayed PPD and reduced lipid peroxidation as well as decreased H2O2 accumulation. Quantitative proteomics data from ethene and phenylpropanoid pathways indicate additional gene candidates to further delay PPD. Cassava root proteomics data are available at www.pep2pro.ethz.ch for easy access and comparison with other proteomics data.

Rice NICOTIANAMINE SYNTHASE 2 expression improves dietary iron and zinc levels in wheat.

KEY MESSAGE: Iron and zinc deficiencies negatively impact human health worldwide. We developed wheat lines that meet or exceed recommended dietary target levels for iron and zinc in the grains. These lines represent useful germplasm for breeding new wheat varieties that can reduce iron and zinc deficiency-associated health burdens in the affected populations. Micronutrient deficiencies, including iron and zinc deficiencies, have negative impacts on human health globally. Iron-deficiency; anemia affects nearly two billion people worldwide and is the cause of reduced cognitive development, fatigue and overall low productivity. Similarly, zinc deficiency causes stunted growth, decreased immunity and increased risk of respiratory infections. Biofortification of staple crops is a sustainable and effective approach to reduce the burden of health problems associated with micronutrient deficiencies. Here, we developed wheat lines expressing rice NICOTIANAMINE SYNTHASE 2 (OsNAS2) and bean FERRITIN (PvFERRITIN) as single genes as well as in combination. NAS catalyzes the biosynthesis of nicotianamine (NA), which is a precursor of the iron chelator deoxymugeneic acid (DMA) required for long distance iron translocation. FERRITIN is important for iron storage in plants because it can store up to 4500 iron ions. We obtained significant increases of iron and zinc content in wheat grains of plants expressing either OsNAS2 or PvFERRTIN, or both genes. In particular, wheat lines expressing OsNAS2 greatly surpass the HarvestPlus recommended target level of 30 % dietary estimated average requirement (EAR) for iron, and 40 % of EAR for zinc, with lines containing 93.1 µg/g of iron and 140.6 µg/g of zinc in the grains. These wheat lines with dietary significant levels of iron and zinc represent useful germplasm for breeding new wheat varieties that can reduce micronutrient deficiencies in affected populations.

Characterization of Brown Streak Virus-Resistant Cassava.

Cassava brown streak disease (CBSD) has become a major constraint to cassava production in East and Central Africa. The identification of new sources of CBSD resistance is essential to deploy CBSD mitigation strategies, as the disease is progressing westwards to new geographical areas. A stringent infection method based on top cleft-grafting combined with precise virus titer quantitation was utilized to screen 14 cassava cultivars and elite breeding lines. When inoculated with mixed infections of Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), the scions of elite breeding lines KBH 2006/18 and KBH 2006/26 remained symptom-free during a 16-week period of virus graft inoculation, while susceptible varieties displayed typical CBSD infection symptoms at 4 weeks after grafting. The identified CBSD resistance was stable under the coinoculation of CBSV and UCBSV with cassava geminiviruses. Double-grafting experiments revealed that transmission of CBSV and UCBSV to CBSD-susceptible top scions was delayed when using intermediate scions of elite breeding lines KBH 2006/18 and KBH 2006/26. Nonetheless, comparison of virus systemic movement using scions from KBH2006/18 and a transgenic CBSD resistant 60444 line (60444-Hp9 line) showed that both CBSV and UCBSV move at undetectable levels through the stems. Further, protoplast-based assays of virus titers showed that the replication of CBSV is inhibited in the resistant line KBH2006/18, suggesting that the identified CBSD resistance is at least partially based on inhibition of virus replication. Our molecular characterization of CBSD resistance in cassava offers a robust virus-host system to further investigate the molecular determinants of CBSD resistance.