An initial study on gene cloning and characterization of unicellular green microalga Ankistrodesmus convolutus was carried out to isolate and characterize the full-length cDNA of ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit (RbcS) as a first step towards elucidating the structure of A. convolutus RbcS gene. The full-length of A. convolutus RbcS cDNA (AcRbcS) contained 28 bp of 5' untranslated region (UTR), 225 bp of 3' non-coding region, and an open reading frame of 165 amino acids consisting of a chloroplast transit peptide with 24 amino acids and a mature protein of 141 amino acids. The amino acid sequence has high identity to those of other green algae RbcS genes. The AcRbcS contained a few conserved domains including protein kinase C phosphorylation site, tyrosine kinase phosphorylation site and N-myristoylation sites. The AcRbcS was successfully expressed in Escherichia coli and a ~21 kDa of anticipated protein band was observed on SDS-PAGE. From the phylogenetic analysis of RbcS protein sequences, it was found that the RbcS of A. convolutus has closer genetic relationship with green microalgae species compared to those of green seaweed and green macroalgae species. Southern hybridization analysis revealed that the AcRbcS is a member of a small multigene family comprising of two to six members in A. convolutus genome. Under different illumination conditions, RT-PCR analysis showed that AcRbcS transcription was reduced in the dark, and drastically recovered in the light condition. Results presented in this paper established a good foundation for further study on the photosynthetic process of A. convolutus and other green algae species where little information is known on Rubisco small subunit.
We have characterized an oil palm (Elaeis guineensis Jacq.) constitutive promoter that is derived from a translationally control tumor protein (TCTP) gene. The TCTP promoter was fused transcriptionally with the gusA reporter gene and transferred to monocot and dicot systems in order to study its regulatory role in a transient expression study. It was found that the 5' region of TCTP was capable of driving the gusA expression in all the oil palm tissues tested, including immature embryo, embryogenic callus, embryoid, young leaflet from mature palm, green leaf, mesocarp and stem. It could also be used in dicot systems as it was also capable of driving gusA expression in tobacco leaves. The results indicate that the TCTP promoter could be used for the production of recombinant proteins that require constitutive expression in the plant system.
Lignification of the plant cell wall could serve as the first line of defense against pathogen attack, but the molecular mechanisms of virulence and disease between oil palm and Ganoderma boninense are poorly understood. This study presents the biochemical, histochemical, enzymology and gene expression evidences of enhanced lignin biosynthesis in young oil palm as a response to G. boninense (GBLS strain). Comparative studies with control (T1), wounded (T2) and infected (T3) oil palm plantlets showed significant accumulation of total lignin content and monolignol derivatives (syringaldehyde and vanillin). These derivatives were deposited on the epidermal cell wall of infected plants. Moreover, substantial differences were detected in the activities of enzyme and relative expressions of genes encoding phenylalanine ammonia lyase (EC 4.3.1.24), cinnamate 4-hydroxylase (EC 1.14.13.11), caffeic acid O-methyltransferase (EC 2.1.1.68) and cinnamyl alcohol dehydrogenase (CAD, EC 1.1.1.195). These enzymes are key intermediates dedicated to the biosynthesis of lignin monomers, the guaicyl (G), syringyl (S) and ρ-hydroxyphenyl (H) subunits. Results confirmed an early, biphasic and transient positive induction of all gene intermediates, except for CAD enzyme activities. These differences were visualized by anatomical and metabolic changes in the profile of lignin in the oil palm plantlets such as low G lignin, indicating a potential mechanism for enhanced susceptibility toward G. boninense infection.
CBF/DREB1 is a group of transcription factors that are mainly involved in abiotic stress tolerance in plants. They belong to the AP2/ERF superfamily of plant-specific transcription factors. A gene encoding a new member of this group was isolated from ripening oil palm fruit and designated as EgCBF3. The oil palm fruit demonstrates the characteristics of a climacteric fruit like tomato, in which ethylene has a major impact on the ripening process. A transgenic approach was used for functional characterization of the EgCBF3, using tomato as the model plant. The effects of ectopic expression of EgCBF3 were analyzed based on expression profiling of the ethylene biosynthesis-related genes, anti-freeze proteins (AFPs), abiotic stress tolerance and plant growth and development. The EgCBF3 tomatoes demonstrated altered phenotypes compared to the wild type tomatoes. Delayed leaf senescence and flowering, increased chlorophyll content and abnormal flowering were the consequences of overexpression of EgCBF3 in the transgenic tomatoes. The EgCBF3 tomatoes demonstrated enhanced abiotic stress tolerance under in vitro conditions. Further, transcript levels of ethylene biosynthesis-related genes, including three SlACSs and two SlACOs, were altered in the transgenic plants' leaves and roots compared to that in the wild type tomato plant. Among the eight AFPs studied in the wounded leaves of the EgCBF3 tomato plants, transcript levels of SlOSM-L, SlNP24, SlPR5L and SlTSRF1 decreased, while expression of the other four, SlCHI3, SlPR1, SlPR-P2 and SlLAP2, were up-regulated. These findings indicate the possible functions of EgCBF3 in plant growth and development as a regulator of ethylene biosynthesis-related and AFP genes, and as a stimulator of abiotic stress tolerance.
Rice is one of the most important cereal crops with great potential for biotechnology progress. In transformation method, antibiotic resistance genes are routinely used as powerful markers for selecting transformed cells from surrounding non-transformed cells. In this study, the toxicity level of hygromycin was optimized for two selected mutant rice lines, MR219 line 4 and line 9. The mature embryos were isolated and cultured on an MS medium with different hygromycin concentrations (0, 20, 40, 60, 80 and 100 mg L(-1)). Evidently, above 60 mg L(-1) was effective for callus formation and observed completely dead. Further there were tested for specific concentration (0-60). Although, 21.28% calli survived on the medium containing 45 mg L(-1) hygromycin, it seemed suitable for the identification of putative transformants. These findings indicated that a system for rice transformation in a relatively high frequency and the transgenes are stably expressed in the transgenic plants. Green shoots were regenerated from the explant under hygromycin stress. RT-PCR using hptII and gus sequence specific primer and Southern blot analysis were used to confirm the presence of the transgene and to determine the transformation efficiency for their stable integration in regenerated plants. This study demonstrated that the hygromycin resistance can be used as an effective marker for rice transformation.
Silicon (Si) plays an important role in reducing plant susceptibility against a variety of different biotic and abiotic stresses; and also has an important regulatory role in soil to avoid heavy metal toxicity and providing suitable growing conditions for plants. A full-length cDNAs of 696bp of serine-rich protein was cloned from mangrove plant (Rhizophora apiculata) by amplification of cDNA ends from an expressed sequence tag homologous to groundnut (Arachis hypogaea), submitted to NCBI (KF211374). This serine-rich protein gene encodes a deduced protein of 223 amino acids. The transcript titre of the serine-rich protein was found to be strongly enriched in roots compared with the leaves of two month old mangrove plants and expression level of this serine-rich protein was found to be strongly induced when the mangrove seedlings were exposed to SiO2. Expression of the serine-rich protein transgenic was detected in transgenic Arabidopsis thaliana, where the amount of serine increased from 1.02 to 37.8mg/g. The same trend was also seen in Si content in the roots (14.3%) and leaves (7.4%) of the transgenic A. thaliana compared to the wild-type plants under Si treatment. The biological results demonstrated that the accumulation of the serine amino acid in the vegetative tissues of the transgenic plants enhanced their ability to absorb and accumulate more Si in the roots and leaves and suggests that the serine-rich protein gene has potential for use in genetic engineering of different stress tolerance characteristics.
Panduratin A extracted from Boesenbergia rotunda is a flavonoid reported to possess a range of medicinal indications which include anti-dengue, anti-HIV, anti-cancer, antioxidant and anti-inflammatory properties. Boesenbergia rotunda is a plant from the Zingiberaceae family commonly used as a food ingredient and traditional medicine in Southeast Asia and China. Reports on the health benefits of secondary metabolites extracted from Boesenbergia rotunda over the last few years has resulted in rising demands for panduratin A. However large scale extraction has been hindered by the naturally low abundance of the compound and limited knowledge of its biosynthetic pathway.
A key event in the domestication and breeding of the oil palm Elaeis guineensis was loss of the thick coconut-like shell surrounding the kernel. Modern E. guineensis has three fruit forms, dura (thick-shelled), pisifera (shell-less) and tenera (thin-shelled), a hybrid between dura and pisifera. The pisifera palm is usually female-sterile. The tenera palm yields far more oil than dura, and is the basis for commercial palm oil production in all of southeast Asia. Here we describe the mapping and identification of the SHELL gene responsible for the different fruit forms. Using homozygosity mapping by sequencing, we found two independent mutations in the DNA-binding domain of a homologue of the MADS-box gene SEEDSTICK (STK, also known as AGAMOUS-LIKE 11), which controls ovule identity and seed development in Arabidopsis. The SHELL gene is responsible for the tenera phenotype in both cultivated and wild palms from sub-Saharan Africa, and our findings provide a genetic explanation for the single gene hybrid vigour (or heterosis) attributed to SHELL, via heterodimerization. This gene mutation explains the single most important economic trait in oil palm, and has implications for the competing interests of global edible oil production, biofuels and rainforest conservation.
The polymorphisms of Waxy (Wx) microsatellite and G-T single-nucleotide polymorphism (SNP) in the Wx gene region were analyzed using simplified techniques in fifteen rice varieties. A rapid and reliable electrophoresis method, MetaPhor agarose gel electrophoresis (MAGE), was effectively employed as an alternative to polyacrylamide gel electrophoresis (PAGE) for separating Wx microsatellite alleles. The amplified products containing the Wx microsatellite ranged from 100 to 130 bp in length. Five Wx microsatellite alleles, namely (CT)(10), (CT)(11), (CT)(16), (CT)(17), and (CT)(18) were identified. Of these, (CT)(11) and (CT)(17) were the predominant classes among the tested varieties. All varieties with an apparent amylose content higher than 24% were associated with the shorter repeat alleles; (CT)(10) and (CT)(11), while varieties with 24% or less amylose were associated with the longer repeat alleles. All varieties with intermediate and high amylose content had the sequence AGGTATA at the 5'-leader intron splice site, while varieties with low amylose content had the sequence AGTTATA. The G-T polymorphism was further verified by the PCR-AccI cleaved amplified polymorphic sequence (CAPS) method, in which only genotypes containing the AGGTATA sequence were cleaved by AccI. Hence, varieties with desirable amylose levels can be developed rapidly using the Wx microsatellite and G-T SNP, along with MAGE.
Improving the quality of the non-climacteric fruit, pineapple, is possible with information on the expression of genes that occur during the process of fruit ripening. This can be made known though the generation of partial mRNA transcript sequences known as expressed sequence tags (ESTs). ESTs are useful not only for gene discovery but also function as a resource for the identification of molecular markers, such as simple sequence repeats (SSRs). This paper reports on firstly, the construction of a normalized library of the mature green pineapple fruit and secondly, the mining of EST-SSRs markers using the newly obtained pineapple ESTs as well as publically available pineapple ESTs deposited in GenBank. Sequencing of the clones from the EST library resulted in 282 good sequences. Assembly of sequences generated 168 unique transcripts (UTs) consisting of 34 contigs and 134 singletons with an average length of ≈500 bp. Annotation of the UTs categorized the known proteins transcripts into the three ontologies as: molecular function (34.88%), biological process (38.43%), and cellular component (26.69%). Approximately 7% (416) of the pineapple ESTs contained SSRs with an abundance of trinucleotide SSRs (48.3%) being identified. This was followed by dinucleotide and tetranucleotide SSRs with frequency of 46 and 57%, respectively. From these EST-containing SSRs, 355 (85.3%) matched to known proteins while 133 contained flanking regions for primer design. Both the ESTs were sequenced and the mined EST-SSRs will be useful in the understanding of non-climacteric ripening and the screening of biomarkers linked to fruit quality traits.
MicroRNAs (miRNAs) are a class of small, endogenous non-coding RNAs that negatively regulate gene expression, resulting in the silencing of target mRNA transcripts through mRNA cleavage or translational inhibition. MiRNAs play significant roles in various biological and physiological processes in plants. However, the miRNA-mediated gene regulatory network in pineapple, the model tropical non-climacteric fruit, remains largely unexplored. Here, we report a complete list of pineapple mature miRNAs obtained from high-throughput small RNA sequencing and precursor miRNAs (pre-miRNAs) obtained from ESTs. Two small RNA libraries were constructed from pineapple fruits and leaves, respectively, using Illumina's Solexa technology. Sequence similarity analysis using miRBase revealed 579,179 reads homologous to 153 miRNAs from 41 miRNA families. In addition, a pineapple fruit transcriptome library consisting of approximately 30,000 EST contigs constructed using Solexa sequencing was used for the discovery of pre-miRNAs. In all, four pre-miRNAs were identified (MIR156, MIR399, MIR444 and MIR2673). Furthermore, the same pineapple transcriptome was used to dissect the function of the miRNAs in pineapple by predicting their putative targets in conjunction with their regulatory networks. In total, 23 metabolic pathways were found to be regulated by miRNAs in pineapple. The use of high-throughput sequencing in pineapples to unveil the presence of miRNAs and their regulatory pathways provides insight into the repertoire of miRNA regulation used exclusively in this non-climacteric model plant.
Suppression subtractive hybridization (SSH) is an effective method to identify different genes with different expression levels involved in a variety of biological processes. This method has often been used to study molecular mechanisms of plants in complex relationships with different pathogens and a variety of biotic stresses. Compared to other techniques used in gene expression profiling, SSH needs relatively smaller amounts of the initial materials, with lower costs, and fewer false positives present within the results. Extraction of total RNA from plant species rich in phenolic compounds, carbohydrates, and polysaccharides that easily bind to nucleic acids through cellular mechanisms is difficult and needs to be considered. Remarkable advancement has been achieved in the next-generation sequencing (NGS) field. As a result of progress within fields related to molecular chemistry and biology as well as specialized engineering, parallelization in the sequencing reaction has exceptionally enhanced the overall read number of generated sequences per run. Currently available sequencing platforms support an earlier unparalleled view directly into complex mixes associated with RNA in addition to DNA samples. NGS technology has demonstrated the ability to sequence DNA with remarkable swiftness, therefore allowing previously unthinkable scientific accomplishments along with novel biological purposes. However, the massive amounts of data generated by NGS impose a substantial challenge with regard to data safe-keeping and analysis. This review examines some simple but vital points involved in preparing the initial material for SSH and introduces this method as well as its associated applications to detect different novel genes from different plant species. This review evaluates general concepts, basic applications, plus the probable results of NGS technology in genomics, with unique mention of feasible potential tools as well as bioinformatics.
Papaya is one of the most nutritional fruits, rich in vitamins, carotenoids, flavonoids and other antioxidants. Previous studies showed phytonutrient improvement without affecting quality in tomato fruit and rapeseed through the suppression of DE-ETIOLATED-1 (DET1), a negative regulator in photomorphogenesis. This study is conducted to study the effects of DET1 gene suppression in papaya embryogenic callus. Immature zygotic embryos were transformed with constitutive expression of a hairpin DET1 construct (hpDET1). PCR screening of transformed calli and reverse transcription quantitative PCR (RT-qPCR) verified that DET1 gene downregulation in two of the positive transformants. High-throughput cDNA 3' ends sequencing on DET1-suppressed and control calli for transcriptomic analysis of global gene expression identified a total of 452 significant (FDR genes (DEGs) upon DET1 suppression. The 123 upregulated DEGs were mainly involved in phenylpropanoid biosynthesis and stress responses, compared to 329 downregulated DEGs involved in developmental processes, lipid metabolism, and response to various stimuli. This is the first study to demonstrate transcriptome-wide relationship between light-regulated pathway and secondary metabolite biosynthetic pathways in papaya. This further supports that the manipulation of regulatory gene involved in light-regulated pathway is possible for phytonutrient improvement of tropical fruit crops.
Abiotic stress reduces plant growth and crop productivity. However, the mechanism underlying posttranscriptional regulations of stress response remains elusive. Herein, we report the posttranscriptional mechanism of nucleocytoplasmic RNA transport of stress-responsive transcripts mediated by EgRBP42, a heterogeneous nuclear ribonucleoprotein-like RNA-binding protein from oil palm, which could be necessary for rapid protein translation to confer abiotic stress tolerance in plants. Transgenic Arabidopsis overexpressing EgRBP42 showed early flowering through alteration of gene expression of flowering regulators and exhibited tolerance towards heat, cold, drought, flood, and salinity stresses with enhanced poststress recovery response by increasing the expression of its target stress-responsive genes. EgRBP42 harbours nucleocytoplasmic shuttling activity mediated by the nuclear localization signal and the M9-like domain of EgRBP42 and interacts directly with regulators in the nucleus, membrane, and the cytoplasm. EgRBP42 regulates the nucleocytoplasmic RNA transport of target stress-responsive transcripts through direct binding to their AG-rich motifs. Additionally, EgRBP42 transcript and protein induction by environmental stimuli are regulated at the transcriptional and posttranscriptional levels. Taken together, the posttranscriptional regulation of RNA transport mediated by EgRBP42 may change the stress-responsive protein profiles under abiotic stress conditions leading to a better adaptation of plants to environmental changes.
Morus alba is an important medicinal plant that is used to treat human diseases. The leaf, branch, and root of Morus can be applied as antidiabetic, antioxidant, and anti-inflammatory medicines, respectively. To explore the molecular mechanisms underlying the various pharmacological functions within different parts of Morus, organ-specific proteomics were performed. Protein profiles of the Morus leaf, branch, and root were determined using a gel-free/label-free proteomic technique. In the Morus leaf, branch, and root, a total of 492, 414, and 355 proteins were identified, respectively, including 84 common proteins. In leaf, the main function was related to protein degradation, photosynthesis, and redox ascorbate/glutathione metabolism. In branch, the main function was related to protein synthesis/degradation, stress, and redox ascorbate/glutathione metabolism. In root, the main function was related to protein synthesis/degradation, stress, and cell wall. Additionally, organ-specific metabolites and antioxidant activities were analyzed. These results revealed that flavonoids were highly accumulated in Morus root compared with the branch and leaf. Accordingly, two root-specific proteins named chalcone flavanone isomerase and flavonoid 3,5-hydroxylase were accumulated in the flavonoid pathway. Consistent with this finding, the content of the total flavonoids was higher in root compared to those detected in branch and leaf. These results suggest that the flavonoids in Morus root might be responsible for its biological activity and the root is the main part for flavonoid biosynthesis in Morus.
KEY MESSAGE: This study highlights the changes in umami-related nucleotide and glutamate levels when the AMP deaminase gene was elevated in transgenic tomato. Taste is perceived as one of a combination of five sensations, sweet, sour, bitter, salty, and umami. The umami taste is best known as a savoury sensation and plays a central role in food flavour, palatability, and eating satisfaction. Umami flavour can be imparted by the presence of glutamate and is greatly enhanced by the addition of ribonucleotides, such as inosine monophosphate (IMP) and guanosine monophosphate (GMP). The production of IMP is regulated by the enzyme adenosine monophosphate (AMP) deaminase which functions to convert AMP into IMP. We have generated transgenic tomato (Solanum lycopersicum) lines over expressing AMP deaminase under the control of a fruit-specific promoter. The transgenic lines showed substantially enhanced levels of AMP deaminase expression in comparison to the wild-type control. Elevated AMP deaminase levels resulted in the reduced accumulation of glutamate and increased levels of the umami nucleotide GMP. AMP concentrations were unchanged. The effects on the levels of glutamate and GMP were unexpected and are discussed in relation to the metabolite flux within this pathway.
DNA mismatch repair proteins play an essential role in maintaining genomic integrity during replication and genetic recombination. We successfully isolated a full length MSH2 and partial MSH7 cDNAs from tomato, based on sequence similarity between MutS and plant MSH homologues. Semi-quantitative RT-PCR reveals higher levels of mRNA expression of both genes in young leaves and floral buds. Genetic mapping placed MSH2 and MSH7 on chromosomes 6 and 7, respectively, and indicates that these genes exist as single copies in the tomato genome. Analysis of protein sequences and phylogeny of the plant MSH gene family show that these proteins are evolutionarily conserved, and follow the classical model of asymmetric protein evolution. Genetic manipulation of the expression of these MSH genes in tomato will provide a potentially useful tool for modifying genetic recombination and hybrid fertility between wide crosses.
Black pepper is an important commodity crop in Malaysia that generates millions of annual revenue for the country. However, black pepper yield is affected by slow decline disease caused by a soil-borne fungus Fusarium solani. RNA sequencing transcriptomics approach has been employed in this study to explore the differential gene expression in susceptible Piper nigrum L. and resistant Piper colubrinum Link. Gene expression comparative analysis of the two pepper species has yielded 2,361 differentially expressed genes (DEGs). Among them, higher expression of 1,426 DEGs was detected in resistant plant. These DEGs practically demonstrated the major branches of plant-pathogen interaction pathway (Path: ko04626). We selected five groups of defence-related DEGs for downstream qRT-PCR analysis. Cf-9, the gene responsible for recognizing fungal avirulence protein activity was found inexpressible in susceptible plant. However, this gene exhibited promising expression in resistant plant. Inactivation of Cf-9 could be the factor that causes susceptible plant fail in recognition of F. solani and subsequently delay activation of adaptive response to fungal invasion. This vital study advance the understanding of pepper plant defence in response to F. solani and aid in identifying potential solution to manage slow decline disease in black pepper cultivation.
Transcription factors (TFs) form the major class of regulatory genes and play key roles in multiple plant stress responses. In most eukaryotic plants, transcription factor (TF) families (WRKY, MADS-box and MYB) activate unique cellular-level abiotic and biotic stress-responsive strategies, which are considered as key determinants for defense and developmental processes. Arabidopsis and rice are two important representative model systems for dicot and monocot plants, respectively. A comprehensive comparative study on 101 OsWRKY, 34 OsMADS box and 122 OsMYB genes (rice genome) and, 71 AtWRKY, 66 AtMADS box and 144 AtMYB genes (Arabidopsis genome) showed various relationships among TFs across species. The phylogenetic analysis clustered WRKY, MADS-box and MYB TF family members into 10, 7 and 14 clades, respectively. All clades in WRKY and MYB TF families and almost half of the total number of clades in the MADS-box TF family are shared between both species. Chromosomal and gene structure analysis showed that the Arabidopsis-rice orthologous TF gene pairs were unevenly localized within their chromosomes whilst the distribution of exon-intron gene structure and motif conservation indicated plausible functional similarity in both species. The abiotic and biotic stress-responsive cis-regulatory element type and distribution patterns in the promoter regions of Arabidopsis and rice WRKY, MADS-box and MYB orthologous gene pairs provide better knowledge on their role as conserved regulators in both species. Co-expression network analysis showed the correlation between WRKY, MADs-box and MYB genes in each independent rice and Arabidopsis network indicating their role in stress responsiveness and developmental processes.
A cross between IR64 (high-yielding but drought-susceptible) and Aday Sel (drought-tolerant) rice cultivars yielded a stable line with enhanced grain yield under drought screening field trials at International Rice Research Institute. The major effect qDTY4.1 drought tolerance and yield QTL was detected in the IR77298-14-1-2-10 Backcrossed Inbred Line (BIL) and its IR87705-7-15-B Near Isogenic Line (NIL) with 93.9% genetic similarity to IR64. Although rice yield is extremely susceptible to water stress at reproductive stage, currently, there is only one report on the detection of drought-responsive microRNAs in inflorescence tissue of a Japonica rice line. In this study, more drought-responsive microRNAs were identified in the inflorescence tissues of IR64, IR77298-14-1-2-10 and IR87705-7-15-B via next-generation sequencing. Among the 32 families of inflorescence-specific non-conserved microRNAs that were identified, 22 families were up-regulated in IR87705-7-15-B. Overall 9 conserved and 34 non-conserved microRNA families were found as drought-responsive in rice inflorescence with 5 conserved and 30 non-conserved families induced in the IR87705-7-15-B. The observation of more drought-responsive non-conserved microRNAs may imply their prominence over conserved microRNAs in drought response mechanisms of rice inflorescence. Gene Ontology annotation analysis on the target genes of drought-responsive microRNAs identified in IR87705-7-15-B revealed over-representation of biological processes including development, signalling and response to stimulus. Particularly, four inflorescence-specific microRNAs viz. osa-miR5485, osa-miR5487, osa-miR5492 and osa-miR5517, and two non-inflorescence specific microRNAs viz. osa-miR169d and osa-miR169f.2 target genes that are involved in flower or embryonic development. Among them, osa-miR169d, osa-miR5492 and osa-miR5517 are related to flowering time control. It is also worth mentioning that osa-miR2118 and osa-miR2275, which are implicated in the biosynthesis of rice inflorescence-specific small interfering RNAs, were induced in IR87705-7-15-B but repressed in IR77298-14-1-2-10. Further, gene search within qDTY4.1 QTL region had identified multiple copies of NBS-LRR resistance genes (potential target of osa-miR2118), subtilisins and genes implicated in stomatal movement, ABA metabolism and cuticular wax biosynthesis.