Allelic variation in PtGA20Ox associates with growth and wood properties in Populus spp

Treatment of Parkinson's disease with biologics that penetrate the blood-brain barrier via receptor-mediated transport

Parkinson's disease (PD) is characterized by neurodegeneration of nigral-striatal neurons in parallel with the formation of intra-neuronal α-synuclein aggregates, and these processes are exacerbated by neuro-inflammation. All 3 components of PD pathology are potentially treatable with biologics. Neurotrophins, such as glial derived neurotrophic factor or erythropoietin, can promote neural repair. Therapeutic antibodies can lead to disaggregation of α-synuclein neuronal inclusions. Decoy receptors can block the activity of pro-inflammatory cytokines in brain. However, these biologic drugs do not cross the blood-brain barrier (BBB). Biologics can be made transportable through the BBB following the re-engineering of the biologic as an IgG fusion protein, where the IgG domain targets an endogenous receptor-mediated transcytosis (RMT) system within the BBB, such as the insulin receptor or transferrin receptor. The receptor-specific antibody domain of the fusion protein acts as a molecular Trojan horse to ferry the biologic into brain via the BBB RMT pathway. This review describes the re-engineering of all 3 classes of biologics (neurotrophins, decoy receptor, therapeutic antibodies) for BBB delivery and treatment of PD. Targeting the RMT pathway at the BBB also enables non-viral gene therapy of PD using lipid nanoparticles (LNP) encapsulated with plasmid DNA encoding therapeutic genes. The surface of the lipid nanoparticle is conjugated with a receptor-specific IgG that triggers RMT of the LNP across the BBB in vivo.

Genome-Wide Association Study With Growth-Related Traits and Secondary Metabolite Contents in Red- and White-Heart Chinese Fir

Chinese fir [Cunninghamia lanceolata (Lamb.) Hook] is an important evergreen coniferous tree species that is widely distributed in many southern provinces of China and has important economic value. The Chinese fir accounts for 1/4 and 1/3 of the total artificial forest area and stock volume, respectively. Red-heart Chinese fir is popular in the market because of its high density and red heartwood. The long-growth cycle hindered the breeding process of Chinese fir, while molecular marker-assisted breeding could accelerate it. However, Chinese fir, a perennial conifer species, has a large genome, which has not yet been published. In this study, the growth-related traits and secondary metabolite contents of red- and white-heart Chinese fir were measured and found to be different between them. There are extremely significant differences among growth-related traits (p < 0.001), but secondary metabolite contents have different correlations due to differences in chemical structure. Moreover, genotype effect analysis of the substantially correlated single nucleotide polymorphisms (SNPs) revealed that most of the loci related to each growth-related traits were different from each other, indicating a type specificity of the genes regulated different growth-related traits. Furthermore, among the loci related to secondary metabolite contents, nine loci associated with multiple metabolite phenotypes such as Marker21022_4, Marker21022_172, Marker24559_31, Marker27425_37, Marker20748_85, Marker18841_115, Marker18841_198, Marker65846_146, and Marker21486_163, suggesting the presence of pleiotropic genes. This study identified the potential SNP markers associated with secondary metabolites in Chinese fir, thus setting the basis for molecular marker-assisted selection.

Variation in genomic methylation in natural populations of Populus simonii is associated with leaf shape and photosynthetic traits

DNA methylation, one of the best-studied types of chromatin modification, suppresses the expression of transposable elements, pseudogenes, repetitive sequences, and individual genes. However, the extent and variation of genome-wide DNA methylation in natural populations of plants remain relatively unknown. To investigate variation in DNA methylation and whether this variation associates with important plant traits, including leaf shape and photosynthesis, 20 413 DNA methylation sites were examined in a poplar association population (505 individuals) using methylation-sensitive amplification polymorphism (MSAP) technology. Calculation of epi-population structure and kinships assigned individuals into subsets (K=3), revealing that the natural population of P. simonii consists of three subpopulations. Population epigenetic distance and geographic distance showed a significant correlation (r=0.4688, P<0.001), suggesting that environmental factors may affect epigenetics. Single-marker approaches were also used to identify significant marker-trait associations, which found 1087 high-confidence DNA methylation markers associated with different phenotypic traits explaining ~5-15% of the phenotypic variance. Among these loci, 147 differentially methylated fragments were obtained by sequencing, representing 130 candidate genes. Expression analysis of six candidate genes indicated that these genes might play important roles in leaf development and regulation of photosynthesis. This study provides association analysis to study the effects of DNA methylation on plant development and these data indicate that epigenetics bridges environmental and genetic factors in affecting plant growth and development.

Association genetics in Populus reveals the interactions between Pto-miR160a and its target Pto-ARF16

MicroRNAs (miRNAs) play important roles in the regulation of gene expression in various biological processes. However, the interactions between miRNAs and their targets are largely unknown in plants. As a powerful tool for identification of variation associated with traits, association genetics provides another strategy for exploration of interactions between miRNAs and their targets. Here, we conducted expression analysis and association mapping to evaluate the interaction between Pto-miR160a and its target Pto-ARF16 in Populus tomentosa. By examining the expression patterns of Pto-MIR160a and Pto-ARF16, we identified a significant, negative correlation between their expression levels, indicating that Pto-miR160a may affect the expression of Pto-ARF16. Among the single nucleotide polymorphisms (SNPs) identified in this study, one common SNP in the pre-miRNA region of Pto-miR160a altered its predicted secondary structure while another common SNP in the predicted miRNA target site changed the binding affinity of Pto-miR160a. Linkage disequilibrium (LD) analysis revealed low LD levels of Pto-MIR160a and Pto-ARF16, indicating that they are suitable for candidate gene-based association analysis. Single SNP-based association analysis identified 19 SNPs (false discovery rate Q < 0.05) in Pto-MIR160a and Pto-ARF16 associated with three phenotypic traits. Epistasis analysis further identified 36 SNP-SNP interactions between SNPs in Pto-MIR160a and SNPs in Pto-ARF16, reflecting the possible genetic interaction of Pto-miR160a and Pto-ARF16. Taking these results together, our study identified SNPs in Pto-MIR160a and Pto-ARF16 associated with tree growth and wood properties, providing SNPs with potential applications in marker-assisted breeding and evidence for the genetic interaction of Pto-miR160a and Pto-ARF16.

Association Studies in Populus tomentosa Reveal the Genetic Interactions of Pto-MIR156c and Its Targets in Wood Formation

MicroRNAs (miRNAs) regulate gene expression in many biological processes, but the significance of the interaction between a miRNA and its targets in perennial trees remains largely unknown. Here, we employed transcript profiling and association studies in Populus tomentosa (Pto) to decipher the effect of genetic variation and interactions between Pto-miR156c and its potential targets (Pto-SPL15, Pto-SPL20, and Pto-SPL25) in 435 unrelated individuals from a natural population of P. tomentosa. Single-SNP (single-nucleotide polymorphism) based association studies with analysis of the underlying additive and dominant effects identified 69 significant associations (P < 0.01), representing 51 common SNPs (minor allele frequency > 0.05) from Pto-MIR156c and its three potential targets, with six wood and growth traits, revealing their common roles in wood formation. Epistasis analysis uncovered 129 significant SNP-SNP associations with ten traits, indicating the potential genetic interactions of Pto-MIR156c and its three putative targets. Interestingly, expression analysis in stem (phloem, cambium, and xylem) revealed that Pto-miR156c expression showed strong negative correlations with Pto-SPL20 (r = -0.90, P < 0.01) and Pto-SPL25 (r = -0.65, P < 0.01), and a positive correlation with Pto-SPL15 (r = 0.40, P < 0.01), which also indicated the putative interactions of Pto-miR156c and its potential targets and their common roles in wood formation. Thus, our study provided an alternative approach to decipher the interaction between miRNAs and their targets and to dissect the genetic architecture of complex traits in trees.

The genetic regulatory network centered on Pto-Wuschela and its targets involved in wood formation revealed by association studies

Transcription factors (TFs) regulate gene expression and can strongly affect phenotypes. However, few studies have examined TF variants and TF interactions with their targets in plants. Here, we used genetic association in 435 unrelated individuals of Populus tomentosa to explore the variants in Pto-Wuschela and its targets to decipher the genetic regulatory network of Pto-Wuschela. Our bioinformatics and co-expression analysis identified 53 genes with the motif TCACGTGA as putative targets of Pto-Wuschela. Single-marker association analysis showed that Pto-Wuschela was associated with wood properties, which is in agreement with the observation that it has higher expression in stem vascular tissues in Populus. Also, SNPs in the 53 targets were associated with growth or wood properties under additive or dominance effects, suggesting these genes and Pto-Wuschela may act in the same genetic pathways that affect variation in these quantitative traits. Epistasis analysis indicated that 75.5% of these genes directly or indirectly interacted Pto-Wuschela, revealing the coordinated genetic regulatory network formed by Pto-Wuschela and its targets. Thus, our study provides an alternative method for dissection of the interactions between a TF and its targets, which will strength our understanding of the regulatory roles of TFs in complex traits in plants.

Exogenous GA₃ Application Enhances Xylem Development and Induces the Expression of Secondary Wall Biosynthesis Related Genes in Betula platyphylla

Gibberellin (GA) is a key signal molecule inducing differentiation of tracheary elements, fibers, and xylogenesis. However the molecular mechanisms underlying the effect of GA on xylem elongation and secondary wall development in tree species remain to be determined. In this study, Betula platyphylla (birch) seeds were treated with 300 ppm GA₃ and/or 300 ppm paclobutrazol (PAC), seed germination was recorded, and transverse sections of hypocotyls were stained with toluidine blue; the two-month-old seedlings were treated with 50 μM GA₃ and/or 50 μM PAC, transverse sections of seedling stems were stained using phloroglucinol–HCl, and secondary wall biosynthesis related genes expression was analyzed by real-time quantitative PCR. Results indicated that germination percentage, energy and time of seeds, hypocotyl height and seedling fresh weight were enhanced by GA₃, and reduced by PAC; the xylem development was wider in GA₃-treated plants than in the control; the expression of NAC and MYB transcription factors, CESA, PAL, and GA oxidase was up-regulated during GA₃ treatment, suggesting their role in GA₃-induced xylem development in the birch. Our results suggest that GA₃ induces the expression of secondary wall biosynthesis related genes to trigger xylogenesis in the birch plants.

Association mapping in Populus reveals the interaction between Pto-miR530a and its target Pto-KNAT1

We used transcript profiling and multi-SNP association to investigate the genetic regulatory relationship between miRNA Pto-miR530a and its target Pto-KNAT1, identifying additive, dominant, and epistatic effects. MicroRNAs (miRNAs) play crucial roles in the post-transcriptional regulation of plant growth and development; indeed, many studies have described the importance of miRNA-target interactions in herbaceous species. However, elucidation of the miRNA-target interactions in trees may require novel strategies. In the present study, we describe a strategy combining expression profiling by reverse transcription quantitative PCR (RT-qPCR) and association mapping with multiple single nucleotide polymorphisms (SNPs) to evaluate the interaction between Pto-miR530a and its target Pto-KNAT1 in Populus tomentosa. RT-qPCR analysis showed a negative correlation (r = -0.62, P < 0.05) between expression levels of Pto-miR530a and Pto-KNAT1 in eight tissues. We used a Bayesian hierarchical model to identify allelic variants of Pto-miR530a and Pto-KNAT1 that associated with eight traits related to growth and wood properties, in a population of 460 unrelated individuals of P. tomentosa. This analysis identified 27 associations, with the proportions of phenotypic variance (R (2)) contributed by each SNP ranging of 0.82-15.81 %, the additive effects of each SNP ranging of 0.16-18.09, and the dominant effects ranging from -14.09 to 19.00. Epistatic interaction models showed a strong interaction among SNPs in the miRNA target with R (2) of 0.1-3.56 %, and information gain of significant SNP pairs of -3.09 to 0.93 %, representing the regulatory interactions between the miRNA and the mRNA. Thus, we used a new strategy that combines association genetics and expression profiling based on SNPs to study the regulatory relationship between this miRNA and its target mRNA, thereby providing novel advances in our understanding of the genetic architecture of important traits.

Identification of additive, dominant, and epistatic variation conferred by key genes in cellulose biosynthesis pathway in Populus tomentosa†

Economically important traits in many species generally show polygenic, quantitative inheritance. The components of genetic variation (additive, dominant and epistatic effects) of these traits conferred by multiple genes in shared biological pathways remain to be defined. Here, we investigated 11 full-length genes in cellulose biosynthesis, on 10 growth and wood-property traits, within a population of 460 unrelated Populus tomentosa individuals, via multi-gene association. To validate positive associations, we conducted single-marker analysis in a linkage population of 1,200 individuals. We identified 118, 121, and 43 associations (P< 0.01) corresponding to additive, dominant, and epistatic effects, respectively, with low to moderate proportions of phenotypic variance (R(2)). Epistatic interaction models uncovered a combination of three non-synonymous sites from three unique genes, representing a significant epistasis for diameter at breast height and stem volume. Single-marker analysis validated 61 associations (false discovery rate, Q ≤ 0.10), representing 38 SNPs from nine genes, and its average effect (R(2) = 3.8%) nearly 2-fold higher than that identified with multi-gene association, suggesting that multi-gene association can capture smaller individual variants. Moreover, a structural gene-gene network based on tissue-specific transcript abundances provides a better understanding of the multi-gene pathway affecting tree growth and lignocellulose biosynthesis. Our study highlights the importance of pathway-based multiple gene associations to uncover the nature of genetic variance for quantitative traits and may drive novel progress in molecular breeding.

Single nucleotide polymorphisms in a cellulose synthase gene (PtoCesA3) are associated with growth and wood properties in Populus tomentosa

In plants, the composition and organization of the cell wall determine cell shape, enable cell expansion, and affect the properties of woody tissues. Cellulose synthase (CesA) genes encode the enzymes involved in the synthesis of cellulose which is the major component of plant primary and secondary cell walls. Here, we isolated a full-length PtoCesA3 cDNA from the stem cambium tissue of Populus tomentosa. Tissue-specific expression profiling showed that PtoCesA3 is highly expressed during primary cell wall formation. Estimation of single nucleotide polymorphism (SNP) diversity and linkage disequilibrium (LD) revealed that PtoCesA3 harbors high SNP diversity (π(T) = 0.00995 and θ(w) = 0.0102) and low LD (r(2) ≥ 0.1, within 1,280 bp). Association analysis in a P. tomentosa association population (460 individuals) showed that seven SNPs (false discovery rate Q < 0.10) and five haplotypes (Q < 0.10) were significantly associated with growth and wood properties, explaining 4.09-7.02% of the phenotypic variance. All significant marker-trait associations were validated in at least one of the three smaller subsets (climatic regions) while five associations were repeated in the linkage population. Variation in RNA transcript abundance among genotypic classes of significant loci was also confirmed in the association or linkage populations. Identification of PtoCesA3 and examining its allelic polymorphisms using association studies open an avenue to understand the mechanism of cellulose synthesis in the primary cell wall and its effects on the properties of woody tissues.

Wood biosynthesis and typologies: a molecular rhapsody

Wood represents one of the most important renewable commodities for humanity and plays a crucial role in terrestrial ecosystem carbon-cycling. Wood formation is the result of a multitude of events that require the concerted action of endogenous and exogenous factors under the influence of photoperiod, for instance genes and plant growth regulators. Beyond providing mechanical support and being responsible for the increase in stem radial diameter, woody tissues constitute the vascular system of trees and are capable of reacting to environmental stimuli, and as such are therefore quite plastic and responsive. Despite the ecological and economic importance of wood, not all aspects of its formation have been unveiled. Many gaps in our knowledge are still present, which hinder the maximal exploitation of this precious bioresource. This review aims at surveying the current knowledge of wood formation and the available molecular data addressing the relationship between wood production and environmental factors, which have crucial influences on the rhythmic regulation of cambial activity and exert profound effects on tree stem growth, wood yield and properties. We will here go beyond wood sensu stricto, i.e., secondary xylem, and extend our survey to other tissues, namely vascular cambium, phloem and fibres. The purpose is to provide the reader with an overview of the complexity of the topic and to highlight the importance of progressing in the future towards an integrated knowledge on the subject.

Allelic variation within the S-adenosyl-L-homocysteine hydrolase gene family is associated with wood properties in Chinese white poplar (Populus tomentosa)

Background

S-adenosyl-l-homocysteine hydrolase (SAHH) is the only eukaryotic enzyme capable of S-adenosyl-l-homocysteine (SAH) catabolism for the maintenance of cellular transmethylation potential. Recently, biochemical and genetic studies in herbaceous species have obtained important discoveries in the function of SAHH, and an extensive characterization of SAHH family in even one tree species is essential, but currently lacking.

Results

Here, we first identified the SAHH family from Populus tomentosa using molecular cloning method. Phylogenetic analyses of 28 SAHH proteins from dicotyledons, monocotyledons, and lower plants revealed that the sequences formed two monophyletic groups: the PtrSAHHA with PtoSAHHA and PtrSAHHB with PtoSAHHB. Examination of tissue-specific expression profiles of the PtoSAHH family revealed similar expression patterns; high levels of expression in xylem were found. Nucleotide diversity and linkage disequilibrium (LD) in the PtoSAHH family, sampled from P. tomentosa natural distribution, revealed that PtoSAHH harbors high single-nucleotide polymorphism (SNP) diversity (π=0.01059±0.00122 and 0.00930±0.00079,respectively) and low LD (r² > 0.1, within 800 bp and 2,200 bp, respectively). Using an LD-linkage analysis approach, two noncoding SNPs (PtoSAHHB_1065 and PtoSAHHA_2203) and the corresponding haplotypes were found to significantly associate with α-cellulose content, and a nonsynonymous SNP (PtoSAHHB_410) within the SAHH signature motifs showed significant association with fiber length, with an average of 3.14% of the phenotypic variance explained.

Conclusions

The present study demonstrates that PtoSAHHs were split off prior to the divergence of interspecies in Populus, and SAHHs may play a key role promoting transmethylation reactions in the secondary cell walls biosynthesis in trees. Hence, our findings provide insights into SAHH function and evolution in woody species and also offer a theoretical basis for marker-aided selection breeding to improve the wood quality of Populus.

Expression and molecular evolution of two DREB1 genes in black poplar (Populus nigra)

Environmental stresses such as low temperature, drought, and high salinity significantly affect plant growth and yield. As selective forces, these adverse factors play essential roles in shaping phenotypic variation in plant populations. Black poplar (Populus nigra) is an economically and ecologically important forest tree species with widely distributed populations and is thus suitable for experiments detecting evolutionary footprints left by stress. Here, we performed expression and evolutionary analysis of two duplicated DREB A1-subgroup (DREB1) genes, PnDREB68 and PnDREB69, encoding transcription factors that are involved in stress responses. The two genes showed partially overlapping but distinct expression patterns in response to stresses. These genes were strongly and rapidly induced by cold stress in leaves, stems, and roots. In leaf tissue, dehydration stress induced the expression of PnDREB68 but not PnDREB69. PnDREB69 displayed more rapid responses and longer expression durations than PnDREB68 under salt and ABA stress, respectively. Based on single nucleotide polymorphism (SNP) analysis, we found significant population genetic differentiation, with a greater F_ST value (0.09189) for PnDREB69 than for PnDREB68 (0.07743). Nucleotide diversity analysis revealed a two-fold higher π_T for PnDREB68 than for PnDREB69 (0.00563 vs. 0.00243), reflecting strong purifying selection acting on the former. The results suggest that positive selection acted on PnDREB69, as evidenced by neutral testing using Tajima’s D statistic. The distinct selective forces to which each of the genes was subjected may be associated with expression divergence. Linkage disequilibrium (LD) was low for the sequenced region, with a higher level for PnDREB68 than for PnDREB69. Additionally, analysis of the relationship among carbon isotope ratios, SNP classes and gene expression, together with motif and domain analysis, suggested that 14 polymorphisms within the two genes may be candidates for an association study of important traits such as water use efficiency/drought tolerance in black poplar.

Single-nucleotide polymorphisms in PtoCesA7 and their association with growth and wood properties in Populus tomentosa

Cellulose synthase (CesA) genes encode the enzymes that synthesize cellulose; therefore, CesAs play central roles in plant development and affect the yield and quality of wood, essential properties for industrial applications of plant biomass. To effectively manipulate wood biosynthesis in trees and improve wood quality, we thus require a better understanding of the natural variation in CesAs. Association studies have emerged as a powerful tool for identification of variation associated with quantitative traits. Here, we used a candidate gene-based association mapping approach to identify PtoCesA7 allelic variants that associate with growth and wood quality traits in Populus tomentosa. We isolated a full-length PtoCesA7 cDNA and observed high PtoCesA7 expression in xylem, consistent with the xylem-specific expression of CesA7. Nucleotide diversity and linkage disequilibrium (LD) in PtoCesA7, sampled from the P. tomentosa natural distribution, revealed that PtoCesA7 harbors high nucleotide diversity (π(T) = 0.0091) and low LD (r(2) ≥ 0.1, within 800 bp). By association analysis, we identified seven single-nucleotide polymorphisms (SNPs) (false discovery rate Q < 0.10) and 12 haplotypes (Q < 0.10) that associated with growth and wood properties, explaining 3.62-10.59 % of the phenotypic variance. We also validated 9 of the 10 significant marker-trait associations in at least one of three smaller subsets (climatic regions) or in a linkage-mapping population. Thus, our study identified functional PtoCesA7 allelic variants associated with growth and wood quality traits, giving new insights into genes affecting wood quality and quantity. From an applied perspective, the SNPs revealed in this study have potential applications in marker-assisted breeding.

Allelic variation in a cellulose synthase gene (PtoCesA4) associated with growth and wood properties in Populus tomentosa

Lignocellulosic biomass from trees provides a renewable feedstock for biofuels, lumber, pulp, paper, and other uses. Dissecting the mechanism underlying natural variation of the complex traits controlling growth and lignocellulose biosynthesis in trees can enable marker-assisted breeding to improve wood quality and yield. Here, we combined linkage disequilibrium (LD)-based association analysis with traditional linkage analysis to detect the genetic effect of a Populus tomentosa cellulose synthase gene, PtoCesA4. PtoCesA4 is strongly expressed in developing xylem and leaves. Nucleotide diversity and LD in PtoCesA4, sampled from the P. tomentosa natural distribution, revealed that PtoCesA4 harbors high single nucleotide polymorphism (SNP) diversity (πT = 0.0080 and θw = 0.0098) and low LD (r(2) ≥ 0.1, within 1400 bp), demonstrating that the potential of a candidate-gene-based LD approach in understanding the molecular basis underlying quantitative variation in this species. By combining single SNP, multi-SNP, and haplotype-based associations in an association population of 460 individuals with single SNP linkage analysis in a family-based linkage populations (1200 individuals), we identified three strong associations (false discovery rate Q < 0.05) in both populations. These include two nonsynonymous markers (SNP49 associated with α-cellulose content and SNP59 associated with fiber width) and a noncoding marker (SNP18 associated with α-cellulose content). Variation in RNA transcript abundance among genotypic classes of SNP49 was confirmed in these two populations. Therefore, combining different methods allowed us to examine functional PtoCesA4 allelic variation underlying natural variation in complex quantitative traits related to growth and lignocellulosic biosynthesis.

The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics

In woody crop plants, the oligosaccharide components of the cell wall are essential for important traits such as bioenergy content, growth, and structural wood properties. UDP-glucuronate decarboxylase (UXS) is a key enzyme in the synthesis of UDP-xylose for the formation of xylans during cell wall biosynthesis. Here, we isolated a multigene family of seven members (PtUXS1-7) encoding UXS from Populus tomentosa, the first investigation of UXSs in a tree species. Analysis of gene structure and phylogeny showed that the PtUXS family could be divided into three groups (PtUXS1/4, PtUXS2/5, and PtUXS3/6/7), consistent with the tissue-specific expression patterns of each PtUXS. We further evaluated the functional consequences of nucleotide polymorphisms in PtUXS1. In total, 243 single-nucleotide polymorphisms (SNPs) were identified, with a high frequency of SNPs (1/18 bp) and nucleotide diversity (π_T = 0.01033, θ_w =0.01280). Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r²<0.1, P<0.001, within 700 bp). SNP- and haplotype-based association analysis showed that nine SNPs (Q <0.10) and 12 haplotypes (P<0.05) were significantly associated with growth and wood property traits in the association population (426 individuals), with 2.70% to 12.37% of the phenotypic variation explained. Four significant single-marker associations (Q <0.10) were validated in a linkage mapping population of 1200 individuals. Also, RNA transcript accumulation varies among genotypic classes of SNP10 was further confirmed in the association population. This is the first comprehensive study of the UXS gene family in woody plants, and lays the foundation for genetic improvements of wood properties and growth in trees using genetic engineering or marker-assisted breeding.