The Role of Transmembrane Proteins in Plant Growth, Development, and Stress Responses

Genome-wide exploration: Evolution, structural characterization, molecular docking, molecular dynamics simulation and expression analysis of sugar transporter (ST) gene family in potato (Solanum tuberosum)

Sugars are the basic structural components of carbohydrates. Sugar transport is crucial for plants to ensure their optimal growth and development. Long-distance sugar transport occurs through either diffusion-based passive or active transport mediated by transporter proteins. In potatoes, STs play a vital role in sugar transport and total sugar accumulation. To better understand the roles of these transporters, in-depth structural, protein characterization, and tissue-specific expression analysis were performed. A total of 61 StSTs were identified and classified into eight sub-families (STP, PLT, ERD6L, INT, TMT, pGlcT, SUC, and VGT). The majority of StSTs were found in the plasma membrane, and all of them were dispersed throughout the 12 chromosomes. Exon and motif counts ranged from 1-18 and 1-10, respectively. In synteny analysis with four plant genomes, the highest (38) orthologous gene pair was found with S. lycopersicum (tomato). In 3D protein modeling, the alpha helix and transmembrane helices range varied from 32 % to 78 % and 53 %-57 %, respectively. During molecular docking analysis, the lowest binding energy was observed for Glu-StINT1 (ΔG: - 6.6 kcal/mol), Fru-StVGT1 (ΔG: - 6.1 kcal/mol), Gal-StSTP10 (ΔG: - 6.5 kcal/mol), and Suc-StINT2 (ΔG: - 7.5 kcal/mol), among 244 docking results. These complexes showed significant hydrogen and hydrophobic interactions, due to having significant amino acid residues. The molecular dynamics (MD) simulation of four complexes (Glu-StINT1, Fru-StVGT1, Gal-StSTP10, and Suc-StINT2) validated the ligand's stable attachment to the intended target proteins and it can be predicted that these complexes are the best sugar transporters of potato. In RNA-seq mediated expression analysis, StSTP12, StERD6L-6, 12, StpGlcT3, StVGT1, and StVGT2, were significantly upregulated in vegetative tissues/organs, revealing their significant role in vegetative organ development. In addition, stu-miRNA395 was the largest family interacting with StERD6L-1 and StTMT2 genes, demonstrating their significant role in sulfate metabolism. The detection and visualization of potential transcription factors (TFs) like ERF, Dof, MYB, BBR-BPC, LBD, and NAC in conjunction with the StSTs gene indicate their significant contribution to stress tolerance and DNA conversion and transcription into RNA. A significant interaction of StSTs in the PPI network might be due to their cumulative role in the same signaling pathways. The integration of these findings will guide the development of programming-based sugar transporter-mediated genetic circuits to improve the sugar accumulation in potatoes using synthetic biology approaches.

Identification and functional analyses of drought stress resistance genes by transcriptomics of the Mongolian grassland plant Chloris virgata

BACKGROUND

Mongolian grasslands, including the Gobi Desert, have been exposed to drought conditions with few rains. In such harsh environments, plants with highly resistant abilities against drought stress survive over long periods. We hypothesized that these plants could harbor novel and valuable genes for enhancing drought stress resistance.

RESULTS

In this study, we identified Chloris virgata, a Mongolian grassland plant with strong drought resistance. RNA-seq-based transcriptome analysis was performed to uncover genes associated with drought stress resistance in C. virgata. De novo transcriptome assembly revealed 25,469 protein-coding transcripts and 1,219 upregulated genes after 3- and 6-hr drought stress treatments. Analysis by homology search and Gene Ontology (GO) enrichment indicated that abscisic acid (ABA)- and drought stress-related GO terms were enriched. Among the highly induced genes, ten candidate cDNAs were selected and overexpressed in Arabidopsis. When subjected to drought stress, three of these genes conferred strong drought resistance in the transgenic plants. We named these genes Mongolian Grassland plant Drought-stress resistance genes 1, 2, and 3 (MGD1, MGD2, and MGD3). Gene expression analyses in the transformants suggested that MGD1, MGD2, and MGD3 may activate drought stress-related signalling pathways.

CONCLUSION

This study highlighted the drought resistance of C. virgata and identified three novel genes that enhance drought stress resistance.

Transmembrane proteins in grape immunity: current knowledge and methodological advances

Transmembrane proteins (TMPs) are pivotal components of plant defence mechanisms, serving as essential mediators in the response to biotic stresses. These proteins are among the most complex and diverse within plant cells, making their study challenging. In spite of this, relatively few studies have focused on the investigation and characterization of TMPs in plants. This is particularly true for grapevine. This review aims to provide a comprehensive overview of TMP-encoding genes involved in grapevine immunity. These genes include Lysin Motif Receptor-Like Kinases (LysM-RLKs), which are involved in the recognition of pathogens at the apoplastic level, Plant Respiratory Burst Oxidase Homologs (Rbohs), which generate reactive oxygen species (ROS) for host defense, and Sugars Will Eventually be Exported Transporters (SWEETs), which play a role in nutrient allocation and stress responses. Furthermore, the review discusses the methodologies employed to study TMPs, including in vivo, in vitro and in silico approaches, highlighting their strengths and limitations. In vivo studies include the assessment of TMP function in whole plants or plant tissues, while in vitro experiments focus on isolating and characterizing either specific TMPs or their components. In silico analyses utilize computational tools to predict protein structure, function, and interactions. By identifying and characterizing genes encoding TMPs involved in grapevine immunity, researchers can develop strategies to enhance grapevine resilience and lead to more sustainable viticulture.

Dissecting genomic regions and underlying candidate genes in groundnut MAGIC population for drought tolerance

BACKGROUND

Groundnut is mainly grown in the semi-arid tropic (SAT) regions worldwide, where abiotic stress like drought is persistent. However, a major research gap exists regarding exploring the genetic and genomic underpinnings of tolerance to drought. In this study, a multi-parent advanced generation inter-cross (MAGIC) population was developed and evaluated for five seasons at two locations for three consecutive years (2018-19, 2019-20 and 2020-21) under drought stress and normal environments.

RESULTS

Phenotyping data of drought tolerance related traits, combined with the high-quality 10,556 polymorphic SNPs, were used to perform multi-locus model genome-wide association study (GWAS) analysis. We identified 37 significant marker-trait associations (MTAs) (Bonferroni-corrected) accounting, 0.91- 9.82% of the phenotypic variance. Intriguingly, 26 significant MTAs overlap on four chromosomes (Ah03, Ah07, Ah10 and Ah18) (harboring 70% of MTAs), indicating genomic hotspot regions governing drought tolerance traits. Furthermore, important candidate genes associated with leaf senescence (NAC transcription factor), flowering (B3 domain-containing transcription factor, Ulp1 protease family, and Ankyrin repeat-containing protein), involved in chlorophyll biosynthesis (FAR1 DNA-binding domain protein), stomatal regulation (Rop guanine nucleotide exchange factor; Galacturonosyltransferases), and associated with yield traits (Fasciclin-like arabinogalactan protein 11 and Fasciclin-like arabinogalactan protein 21) were found in the vicinity of significant MTAs genomic regions.

CONCLUSION

The findings of our investigation have the potential to provide a basis for significant MTAs validation, gene discovery and development of functional markers, which could be employed in genomics-assisted breeding to develop climate-resilient groundnut varieties.

GWAS for Drought Resilience Traits in Red Clover (Trifolium pratense L.)

Red clover (Trifolium pratense L.) is a well-appreciated grassland crop in temperate climates but suffers from increasingly frequent and severe drought periods. Molecular markers for drought resilience (DR) would benefit breeding initiatives for red clover, as would a better understanding of the genes involved in DR. Two previous studies, as follows, have: (1) identified phenotypic DR traits in a diverse set of red clover accessions; and (2) produced genotypic data using a pooled genotyping-by-sequencing (GBS) approach in the same collection. In the present study, we performed genome-wide association studies (GWAS) for DR using the available phenotypic and genotypic data. Single nucleotide polymorphism (SNP) calling was performed using GBS data and the following two red clover genome assemblies: the recent HEN-17 assembly and the Milvus assembly. SNP positions with significant associations were used to delineate flanking regions in both genome assemblies, while functional annotations were retrieved from Medicago truncatula orthologs. GWAS revealed 19 significant SNPs in the HEN-17-derived SNP set, explaining between 5.3 and 23.2% of the phenotypic variation per SNP-trait combination for DR traits. Among the genes in the SNP-flanking regions, we identified candidate genes related to cell wall structuring, genes encoding sugar-modifying proteins, an ureide permease gene, and other genes linked to stress metabolism pathways. GWAS revealed 29 SNPs in the Milvus-derived SNP set that explained substantially more phenotypic variation for DR traits, between 5.3 and 42.3% per SNP-trait combination. Candidate genes included a DEAD-box ATP-dependent RNA helicase gene, a P-loop nucleoside triphosphate hydrolase gene, a Myb/SANT-like DNA-binding domain protein, and an ubiquitin-protein ligase gene. Most accessions in this study are genetically more closely related to the Milvus genotype than to HEN-17, possibly explaining how the Milvus-derived SNP set yielded more robust associations. The Milvus-derived SNP set pinpointed 10 genomic regions that explained more than 25% of the phenotypic variation for DR traits. A possible next step could be the implementation of these SNP markers in practical breeding programs, which would help to improve DR in red clover. Candidate genes could be further characterized in future research to unravel drought stress resilience in red clover in more detail.

Spontaneous Dimerization and Distinct Packing Modes of Transmembrane Domains in Receptor Tyrosine Kinases

The insulin receptor (IR) and the insulin-like growth factor-1 receptor (IGF1R) are homodimeric transmembrane glycoproteins that transduce signals across the membrane on binding of extracellular peptide ligands. The structures of IR/IGF1R fragments in apo and liganded states have revealed that the extracellular subunits of these receptors adopt Λ-shaped configurations to which are connected the intracellular tyrosine kinase (TK) domains. The binding of peptide ligands induces structural transitions in the extracellular subunits leading to potential dimerization of transmembrane domains (TMDs) and autophosphorylation in TKs. However, the activation mechanisms of IR/IGF1R, especially the role of TMDs in coordinating signal-inducing structural transitions, remain poorly understood, in part due to the lack of structures of full-length receptors in apo or liganded states. While atomistic simulations of IR/IGF1R TMDs showed that these domains can dimerize in single component membranes, spontaneous unbiased dimerization in a plasma membrane having a physiologically representative lipid composition has not been observed. We address this limitation by employing coarse-grained (CG) molecular dynamics simulations to probe the dimerization propensity of IR/IGF1R TMDs. We observed that TMDs in both receptors spontaneously dimerized independent of their initial orientations in their dissociated states, signifying their natural propensity for dimerization. In the dimeric state, IR TMDs predominantly adopted X-shaped configurations with asymmetric helical packing and significant tilt relative to the membrane normal, while IGF1R TMDs adopted symmetric V-shaped or parallel configurations with either no tilt or a small tilt relative to the membrane normal. Our results suggest that IR/IGF1R TMDs spontaneously dimerize and adopt distinct dimerized configurations.

Cross-Species Insights into PR Proteins: A Comprehensive Study of Arabidopsis thaliana, Solanum lycopersicum, and Solanum tuberosum

This study provides a comprehensive analysis of pathogenesis-related (PR) proteins, focusing on PR1, PR5, and PR10, in three plant species: Arabidopsis thaliana (At), Solanum lycopersicum (Sl), and Solanum tuberosum (St). We investigated various physico-chemical properties, including protein length, molecular weight, isoelectric point (pI), hydrophobicity, and structural characteristics, such as RMSD, using state-of-the-art tools like AlphaFold and PyMOL. Our analysis found that the SlPR10-StPR10 protein pair had the highest sequence identity (80.00%), lowest RMSD value (0.307 Å), and a high number of overlapping residues (160) among all other protein pairs, indicating their remarkable similarity. Additionally, we used bioinformatics tools such as Cello, Euk-mPLoc 2.0, and Wolfpsort to predict subcellular localization, with AtPR1, AtPR5, and SlPR5 proteins predicted to be located in the extracellular space in both Arabidopsis and S. lycopersicum, while AtPR10 was predicted to be located in the cytoplasm. This comprehensive analysis, including the use of cutting-edge structural prediction and subcellular localization tools, enhances our understanding of the structural, functional, and localization aspects of PR proteins, shedding light on their roles in plant defense mechanisms across different plant species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-024-01343-1.

Spontaneous Dimerization and Distinct Packing Modes of Transmembrane Domains in Receptor Tyrosine Kinases

The insulin receptor (IR) and the insulin-like growth factor-1 receptor (IGF1R) are homodimeric transmembrane glycoproteins that transduce signals across the membrane on binding of extracellular peptide ligands. The structures of IR/IGF1R fragments in apo and liganded states have revealed that the extracellular subunits of these receptors adopt Λ-shaped configurations to which are connected the intracellular tyrosine kinase (TK) domains. The binding of peptide ligands induces structural transitions in the extracellular subunits leading to potential dimerization of transmembrane domains (TMDs) and autophosphorylation in TKs. However, the activation mechanisms of IR/IGF1R, especially the role of TMDs in coordinating signal-inducing structural transitions, remain poorly understood, in part due to the lack of structures of full-length receptors in apo or liganded states. While atomistic simulations of IR/IGF1R TMDs showed that these domains can dimerize in single component membranes, spontaneous unbiased dimerization in a plasma membrane having physiologically representative lipid composition has not been observed. We address this limitation by employing coarse-grained (CG) molecular dynamics simulations to probe the dimerization propensity of IR/IGF1R TMDs. We observed that TMDs in both receptors spontaneously dimerized independent of their initial orientations in their dissociated states, signifying their natural propensity for dimerization. In the dimeric state, IR TMDs predominantly adopted X-shaped configurations with asymmetric helical packing and significant tilt relative to the membrane normal, while IGF1R TMDs adopted symmetric V-shaped or parallel configurations with either no tilt or a small tilt relative to the membrane normal. Our results suggest that IR/IGF1R TMDs spontaneously dimerize and adopt distinct dimerized configurations.

Up-regulation of growth-related gene expression in tobacco by volatile compounds released by Bacillus velezensis WSW007

In order to investigate the mechanism of plant growth promoting (PGP) effects of strain Bacillus velezensis WSW007, its PGP traits and production of volatile organic compounds (VOCs) were tested. The effects of VOCs produced by strain WSW007 on plant growth were observed by co-culturing this strain with tobacco seedlings in I-plates. Meanwhile, the effects of VOCs on tobacco gene expression were analysed by a transcriptome analysis and VOCs were identified by solid phase micro extraction coupled with gas chromatography-mass spectrometry (SPME-GC-MS) analysis. As results, strains WSW007 produced acetic acid and siderophore, and could solubilize phosphate; while it also significantly increased the fresh weight of tobacco seedlings via production of VOCs. In transcriptome analysis, plants co-cultured with strain WSW007 presented the highest up-regulated expression for the genes involved in plant growth and development processes, implying that the bacterial VOCs played a role as regulator of plant gene expression. Conclusively, the up-regulation in expression of growth- and development-related genes via VOCs production is an important PGP mechanism in strain B. velezensis WSW007.

Expression and functional analysis of a recombinant aquaporin Z from Antarctic Pseudomonas sp. AMS3

Aquaporin (AQP) is a water channel protein from the family of transmembrane proteins which facilitates the movement of water across the cell membrane. It is ubiquitous in nature, however the understanding of the water transport mechanism, especially for AQPs in microbes adapted to low temperatures, remains limited. AQP also has been recognized for its ability to be used for water filtration, but knowledge of the biochemical features necessary for its potential applications in industrial processes has been lacking. Therefore, this research was conducted to express, extract, solubilize, purify, and study the functional adaptations of the aquaporin Z family from Pseudomonas sp. AMS3 via molecular approaches. In this study, AqpZ1 AMS3 was successfully subcloned and expressed in E. coli BL21 (DE3) as a recombinant protein. The AqpZ1 AMS3 gene was expressed under optimized conditions and the best optimized condition for the AQP was in 0.5 mM IPTG incubated at 25°C for 20 h induction time. A zwitterionic mild detergent [(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate was the suitable surfactant for the protein solubilization. The protein was then purified via affinity chromatography. Liposome and proteoliposome was reconstituted to determine the particle size using dynamic light scattering. This information obtained from this psychrophilic AQP identified provides new insights into the structural adaptation of this protein at low temperatures and could be useful for low temperature application and molecular engineering purposes in the future.

Genome assembly and microsatellite marker development using Illumina and PacBio sequencing in Persicaria maackiana (Polygonaceae) from Korea

BACKGROUND

Persicaria maackiana (Regel) is a potential medicinal plant that exerts anti-diabetic effects. However, the lack of genomic information on P. maackiana hinders research at the molecular level.

OBJECTIVE

Herein, we aimed to construct a draft genome assembly and obtain comprehensive genomic information on P. maackiana using high-throughput sequencing tools PacBio Sequel II and Illumina.

METHODS

Persicaria maackiana samples from three natural populations in Gaecheon, Gichi, and Uiryeong reservoirs in South Korea were used to generate genomic DNA libraries, perform genome de novo assembly, gene ontology analysis, phylogenetic tree analysis, genotyping, and identify microsatellite markers.

RESULTS

The assembled P. maackiana genome yielded 32,179 contigs. Assessment of assembly integrity revealed 1503 (93.12%) complete Benchmarking Universal Single-Copy Orthologs. A total of 64,712 protein-coding genes were predicted and annotated successfully in the protein database. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs, 13,778 genes were annotated into 18 categories. Genes that activated AMPK were identified in the KEGG pathway. A total of 316,992 microsatellite loci were identified, and primers targeting the flanking regions were developed for 292,059 microsatellite loci. Of these, 150 primer sets were randomly selected for amplification, and 30 of these primer sets were identified as polymorphic. These primers amplified 3-9 alleles. The mean observed and expected heterozygosity were 0.189 and 0.593, respectively. Polymorphism information content values of the markers were 0.361-0.754.

CONCLUSION

Collectively, our study provides a valuable resource for future comparative genomics, phylogeny, and population studies of P. maackiana.

Genome-wide identification of potassium channels in maize showed evolutionary patterns and variable functional responses to abiotic stresses

Potassium (K) channels are essential components of plant biology, mediating not only K ion (K+) homeostasis but also regulating several physiological processes and stress tolerance. In the current investigation, we identified 27 K+ channels in maize and deciphered the evolution and divergence pattern with four monocots and five dicot species. Chromosomal localization and expansion of K+ channel genes showed uneven distribution and were independent of genome size. The dispersed duplication is the major force in expanding K+ channels in the target genomes. The mean Ka/Ks ratio of <0.5 in paralogs and orthologs indicates horizontal and vertical expansions of K+ channel genes under strong purifying selection. The one-to-one K+ channel orthologs were prominent among the closely related species, with higher synteny between maize and the rest of the monocots. Comprehensive K+ channels promoter analysis revealed various cis-regulatory elements mediating stress tolerance with the predominance of MYB and STRE binding sites. The regulatory network showed AP2-EREBP TFs, miR164 and miR399 are prominent regulatory elements of K+ channels. The qRT-PCR analysis of K+ channels and regulatory miRNAs showed significant expressions in response to drought and waterlogging stresses. The present study expanded the knowledge on K+ channels in maize and will serve as a basis for an in-depth functional analysis.

Decoding VaCOLD1 Function in Grapevines: A Membrane Protein Enhancing Cold Stress Tolerance

In globally cultivated grapevines, low-temperature stress poses a persistent challenge. Although COLD1 is recognized as a cold receptor in rice, its function in grapevine cold signaling is unclear. Here, we identified VaCOLD1, a transmembrane protein from the cold-tolerant Vitis amurensis Rupr, which is primarily located on plasma and endoplasmic reticulum membranes. Broadly expressed across multiple tissues, VaCOLD1 responds to various environmental stresses, particularly to cold. Its promoter contains distinct hormone- and stress-responsive elements, with GUS assays confirming widespread expression in Arabidopsis thaliana. Validation of interaction between VaCOLD1 and VaGPA1, together with their combined expression in yeast and grape calli, notably improved cold endurance. Overexpression of VaCOLD1 enhances cold tolerance in Arabidopsis by strengthening the CBF-COR signaling pathway. This is achieved through shielding against osmotic disturbances and modifying the expression of ABA-mediated genes. These findings emphasize the critical role of the VaCOLD1-VaGPA1 complex in mediating the response to cold stress via the CBF-COR pathway.

Sequence Characterization of Extra-Chromosomal Circular DNA Content in Multiple Blackgrass (Alopecurus myosuroides) Populations

Alopecurus myosuroides (blackgrass) is a problematic weed of Western European winter wheat, and its success is largely due to widespread multiple-herbicide resistance. Previous analysis of F2 seed families derived from two distinct blackgrass populations exhibiting equivalent non-target site resistance (NTSR) phenotypes shows resistance is polygenic and evolves from standing genetic variation. Using a CIDER-seq pipeline, we show that herbicide-resistant (HR) and herbicide-sensitive (HS) F3 plants from these F2 seed families as well as the parent populations they were derived from carry extra-chromosomal circular DNA (eccDNA). We identify the similarities and differences in the coding structures within and between resistant and sensitive populations. Although the numbers and size of detected eccDNAs varied between the populations, comparisons between the HR and HS blackgrass populations identified shared and unique coding content, predicted genes, and functional protein domains. These include genes related to herbicide detoxification such as Cytochrome P450s, ATP-binding cassette transporters, and glutathione transferases including AmGSTF1. eccDNA content was mapped to the A. myosuroides reference genome, revealing genomic regions at the distal end of chromosome 5 and the near center of chromosomes 1 and 7 as regions with a high number of mapped eccDNA gene density. Mapping to 15 known herbicide-resistant QTL regions showed that the eccDNA coding sequences matched twelve, with four QTL matching HS coding sequences; only one region contained HR coding sequences. These findings establish that, like other pernicious weeds, blackgrass has eccDNAs that contain homologs of chromosomal genes, and these may contribute genetic heterogeneity and evolutionary innovation to rapidly adapt to abiotic stresses, including herbicide treatment.

Genome-Wide Association Study of Seed Morphology Traits in Senegalese Sorghum Cultivars

Sorghum is considered the fifth most important crop in the world. Despite the potential value of Senegalese germplasm for various traits, such as resistance to fungal diseases, there is limited information on the study of sorghum seed morphology. In this study, 162 Senegalese germplasms were evaluated for seed area size, length, width, length-to-width ratio, perimeter, circularity, the distance between the intersection of length & width (IS) and center of gravity (CG), and seed darkness and brightness by scanning and analyzing morphology-related traits with SmartGrain software at the USDA-ARS Plant Science Research Unit. Correlations between seed morphology-related traits and traits associated with anthracnose and head smut resistance were analyzed. Lastly, genome-wide association studies were performed on phenotypic data collected from over 16,000 seeds and 193,727 publicly available single nucleotide polymorphisms (SNPs). Several significant SNPs were found and mapped to the reference sorghum genome to uncover multiple candidate genes potentially associated with seed morphology. The results indicate clear correlations among seed morphology-related traits and potential associations between seed morphology and the defense response of sorghum. GWAS analysis listed candidate genes associated with seed morphologies that can be used for sorghum breeding in the future.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.