De novo Transcriptome Assembly of Common Wild Rice (Oryza rufipogon Griff.) and Discovery of Drought-Response Genes in Root Tissue Based on Transcriptomic Data

Comparative Cytological and Gene Expression Analysis Reveals That a Common Wild Rice Inbred Line Showed Stronger Drought Tolerance Compared with the Cultivar Rice

Common wild rice (Oryza rufipogon Griff.) is an important germplasm resource containing valuable genes. Our previous analysis reported a stable wild rice inbred line, Huaye3, which derives from the common wild rice of Guangdong Province. However, there was no information about its drought tolerance ability. Here, we assessed the germination characteristics and seedling growth between the Dawennuo and Huaye3 under five concentrations of PEG6000 treatment (0, 5%, 10%, 15%, and 20%). Huaye3 showed a stronger drought tolerance ability, and its seed germination rate still reached more than 52.50% compared with Dawennuo, which was only 25.83% under the 20% PEG6000 treatment. Cytological observations between the Dawennuo and Huaye3 indicated the root tip elongation zone and buds of Huaye3 were less affected by the PEG6000 treatment, resulting in a lower percentage of abnormalities of cortical cells, stele, and shrinkage of epidermal cells. Using the re-sequencing analysis, we detected 13,909 genes that existed in the genetic variation compared with Dawennuo. Of these genes, 39 were annotated as drought stress-related genes and their variance existed in the CDS region. Our study proved the strong drought stress tolerance ability of Huaye3, which provides the theoretical basis for the drought resistance germplasm selection in rice.

Physiological and Multi-Omics Approaches for Explaining Drought Stress Tolerance and Supporting Sustainable Production of Rice

Drought differs from other natural disasters in several respects, largely because of the complexity of a crop's response to it and also because we have the least understanding of a crop's inductive mechanism for addressing drought tolerance among all abiotic stressors. Overall, the growth and productivity of crops at a global level is now thought to be an issue that is more severe and arises more frequently due to climatic change-induced drought stress. Among the major crops, rice is a frontline staple cereal crop of the developing world and is critical to sustaining populations on a daily basis. Worldwide, studies have reported a reduction in rice productivity over the years as a consequence of drought. Plants are evolutionarily primed to withstand a substantial number of environmental cues by undergoing a wide range of changes at the molecular level, involving gene, protein and metabolite interactions to protect the growing plant. Currently, an in-depth, precise and systemic understanding of fundamental biological and cellular mechanisms activated by crop plants during stress is accomplished by an umbrella of -omics technologies, such as transcriptomics, metabolomics and proteomics. This combination of multi-omics approaches provides a comprehensive understanding of cellular dynamics during drought or other stress conditions in comparison to a single -omics approach. Thus a greater need to utilize information (big-omics data) from various molecular pathways to develop drought-resilient crop varieties for cultivation in ever-changing climatic conditions. This review article is focused on assembling current peer-reviewed published knowledge on the use of multi-omics approaches toward expediting the development of drought-tolerant rice plants for sustainable rice production and realizing global food security.

Systems biology for crop improvement

In recent years, generation of large-scale data from genome, transcriptome, proteome, metabolome, epigenome, and others, has become routine in several plant species. Most of these datasets in different crop species, however, were studied independently and as a result, full insight could not be gained on the molecular basis of complex traits and biological networks. A systems biology approach involving integration of multiple omics data, modeling, and prediction of the cellular functions is required to understand the flow of biological information that underlies complex traits. In this context, systems biology with multiomics data integration is crucial and allows a holistic understanding of the dynamic system with the different levels of biological organization interacting with external environment for a phenotypic expression. Here, we present recent progress made in the area of various omics studies-integrative and systems biology approaches with a special focus on application to crop improvement. We have also discussed the challenges and opportunities in multiomics data integration, modeling, and understanding of the biology of complex traits underpinning yield and stress tolerance in major cereals and legumes.

Transcriptomic analysis of Eruca vesicaria subs. sativa lines with contrasting tolerance to polyethylene glycol-simulated drought stress

BACKGROUND

Eruca vesicaria subsp. sativa is one of the Cruciferae species most tolerant to drought stress. In our previous study some extremely drought-tolerant/sensitive Eruca lines were obtained. However little is known about the mechanism for drought tolerance in Eruca.

METHODS

In this study two E. vesicaria subs. sativa lines with contrasting drought tolerance were treated with liquid MS/PEG solution. Total RNA was isolated from 7-day old whole seedlings and then applied to Illumina sequencing platform for high-throughput transcriptional sequencing.

RESULTS

KEGG pathway analysis indicated that differentially expressed genes (DEGs) involved in alpha-Linolenic acid metabolism, Tyrosine metabolism, Phenylalanine, Tyrosine and tryptophan biosynthesis, Galactose metabolism, Isoquinoline alkaloid biosynthesis, Tropane, Piperidine and pyridine alkaloid biosynthesis, Mineral absorption, were all up-regulated specifically in drought-tolerant (DT) Eruca line under drought stress, while DEGs involved in ribosome, ribosome biogenesis, Pyrimidine metabolism, RNA degradation, Glyoxylate and dicarboxylate metabolism, Aminoacyl-tRNA biosynthesis, Citrate cycle, Methane metabolism, Carbon fixation in photosynthetic organisms, were all down-regulated. 51 DEGs were found to be most significantly up-regulated (log2 ratio ≥ 8) specifically in the DT line under PEG treatment, including those for ethylene-responsive transcription factors, WRKY and bHLH transcription factors, calmodulin-binding transcription activator, cysteine-rich receptor-like protein kinase, mitogen-activated protein kinase kinase, WD repeat-containing protein, OPDA reductase, allene oxide cyclase, aquaporin, O-acyltransferase WSD1, C-5 sterol desaturase, sugar transporter ERD6-like 12, trehalose-phosphate phosphatase and galactinol synthase 4. Eight of these 51 DEGs wre enriched in 8 COG and 17 KEGG pathways.

CONCLUSIONS

DEGs that were found to be most significantly up-regulated specifically in the DT line under PEG treatment, up-regulation of DEGs involved in Arginine and proline metabolism, alpha-linolenic acid metabolism and down-regulation of carbon fixation and protein synthesis might be critical for the drought tolerance in Eruca. These results will be valuable for revealing mechanism of drought tolerance in Eruca and also for genetic engineering to improve drought tolerance in crops.

RNA-Seq analysis of soft rush (Juncus effusus): transcriptome sequencing, de novo assembly, annotation, and polymorphism identification

BACKGROUND

Juncus effusus L. (family: Juncaceae; order: Poales) is a helophytic rush growing in temperate damp or wet terrestrial habitats and is of almost cosmopolitan distribution. The species has been studied intensively with respect to its interaction with co-occurring plants as well as microbes being involved in major biogeochemical cycles. J. effusus has biotechnological value as component of Constructed Wetlands where the plant has been employed in phytoremediation of contaminated water. Its genome has not been sequenced.

RESULTS

In this study we carried out functional annotation and polymorphism analysis of de novo assembled RNA-Seq data from 18 genotypes using 249 million paired-end Illumina HiSeq reads and 2.8 million 454 Titanium reads. The assembly comprised 158,591 contigs with a mean contig length of 780 bp. The assembly was annotated using the dammit! annotation pipeline, which queries the databases OrthoDB, Pfam-A, Rfam, and runs BUSCO (Benchmarking Single-Copy Ortholog genes). In total, 111,567 contigs (70.3%) were annotated with functional descriptions, assigned gene ontology terms, and conserved protein domains, which resulted in 30,932 non-redundant gene sequences. Results of BUSCO and KEGG pathway analyses were similar for J. effusus as for the well-studied members of the Poales, Oryza sativa and Sorghum bicolor. A total of 566,433 polymorphisms were identified in transcribed regions with an average frequency of 1 polymorphism in every 171 bases.

CONCLUSIONS

The transcriptome assembly was of high quality and genome coverage was sufficient for global analyses. This annotated knowledge resource can be utilized for future gene expression analysis, genomic feature comparisons, genotyping, primer design, and functional genomics in J. effusus.

Development of introgression lines of AA genome Oryza species, O. glaberrima, O. rufipogon, and O. nivara, in the genetic background of O. sativa L. cv. Taichung 65

To evaluate and utilize potentially valuable quantitative trait loci or genes of wild relatives in the genetic background of domesticated crop species, chromosome segment substitution lines (CSSLs) are a valuable tool. CSSLs can be constructed through the exchange of chromosome segments of AA genome species of the genus Oryza with cultivated rice, Oryza sativa L. Here we report the development of three sets of CSSLs carrying segments of AA genome species closely related to Oryza sativa-O. glaberrima (IRGC 103777 from Mali), O. rufipogon (W1962 from China), and O. nivara (IRGC 105715 from Cambodia)-in the genetic background of ssp. japonica cultivar Taichung 65 through the use of 101 to 121 simple-sequence-repeat markers in whole-genome genotyping and marker-assisted selection. The materials are available via the National Bioresource Project (Rice) Oryzabase Web page.

Leaf Transcriptome Assembly of Protium copal (Burseraceae) and Annotation of Terpene Biosynthetic Genes

Plants in the Burseraceae are globally recognized for producing resins and essential oils with medicinal properties and have economic value. In addition, most of the aromatic and non-aromatic components of Burseraceae resins are derived from a variety of terpene and terpenoid chemicals. Although terpene genes have been identified in model plant crops (e.g., Citrus, Arabidopsis), very few genomic resources are available for non-model groups, including the highly diverse Burseraceae family. Here we report the assembly of a leaf transcriptome of Protium copal, an aromatic tree that has a large distribution in Central America, describe the functional annotation of putative terpene biosynthetic genes and compare terpene biosynthetic genes found in P. copal with those identified in other Burseraceae taxa. The genomic resources of Protium copal can be used to generate novel sequencing markers for population genetics and comparative phylogenetic studies, and to investigate the diversity and evolution of terpene genes in the Burseraceae.

Enriched networks 'nucleoside/nucleotide and ribonucleoside/ribonucleotide metabolic processes' and 'response to stimulus' potentially conferred to drought adaptation of the epiphytic orchid Dendrobium wangliangii

Dendrobium wangliangii is an endangered and epiphytic orchid with tolerance to seasonally extreme arid conditions and occurs exclusively in the hot-dry valley area of southwestern China. To reveal its molecular basis responsible for ecological adaptation, large-scale transcriptome sequencing was performed using Illumina sequencing with pooled mRNA extracted from whole plants and pseudobulbs during drought and rainy seasons. Based on the target transcript selection, the differentially expressed genes were related to 8 well-known drought-tolerant categories, and to morphological traits in resistance to water stress including pseudobulbs and roots. Further gene ontology enrichment analysis revealed that 'nucleoside/nucleotide and ribonucleoside/ribonucleotide metabolic processes' and 'response to stimulus' were the two most important aspects in resistance to drought stress with respect to the whole plant. In addition, the difference in the number and category of differentially expressed genes in whole plant and stem suggested the involvement of genes specifically localized in the stem, such as GTP-binding protein, lipases, signaling related transcripts and those involved in the ATP metabolic process. The comprehensive analysis of the epiphytic orchid in response to water deprivation indicates that integral tactics lead to active adaptation as a basal defense response to drought stress by the endangered epiphyte, including the collaboration of metabolic processes, responses to a various stimulus and other candidate genes contribute to its extreme drought tolerance. Insights from this study can be further utilized to understand stress-responsive genes in other medicinally important species and to improve the drought tolerance of food crops.

Full-length RNA sequencing reveals unique transcriptome composition in bermudagrass

Bermudagrass [Cynodon dactylon (L.) Pers.] is an important perennial warm-season turfgrass species with great economic value. However, the reference genome and transcriptome information are still deficient in bermudagrass, which severely impedes functional and molecular breeding studies. In this study, through analyzing a mixture sample of leaves, stolons, shoots, roots and flowers with single-molecule long-read sequencing technology from Pacific Biosciences (PacBio), we reported the first full-length transcriptome dataset of bermudagrass (C. dactylon cultivar Yangjiang) comprising 78,192 unigenes. Among the unigenes, 66,409 were functionally annotated, whereas 27,946 were found to have two or more isoforms. The annotated full-length unigenes provided many new insights into gene sequence characteristics and systematic phylogeny of bermudagrass. By comparison with transcriptome dataset in nine grass species, KEGG pathway analyses further revealed that C₄ photosynthesis-related genes, notably the phosphoenolpyruvate carboxylase and pyruvate, phosphate dikinase genes, are specifically enriched in bermudagrass. These results not only explained the possible reason why bermudagrass flourishes in warm areas but also provided a solid basis for future studies in this important turfgrass species.

Isolation and Characterization of a Green-Tissue Promoter from Common Wild Rice (Oryza rufipogon Griff.)

Promoters play a very important role in the initiation and regulation of gene transcription. Green-tissue promoter is of great significance to the development of genetically modified crops. Based on RNA-seq data and RT-PCR expression analysis, this study screened a gene, OrGSE (GREEN SPECIAL EXPRESS), which is expressed specifically in green tissues. The study also isolated the promoter of the OrGSE gene (OrGSEp), and predicted many cis-acting elements, such as the CAAT-Box and TATA-Box, and light-responding elements, including circadian, G-BOX and GT1 CONSENSUS. Histochemical analysis and quantification of GUS activity in transgenic Arabidopsis thaliana plants expressing GUS under the control of OrGSEp revealed that this promoter is not only green tissue-specific, but also light-inducible. The ability of a series of 5’-deletion fragments of OrGSEp to drive GUS expression in Arabidopsis was also evaluated. We found that the promoter region from −54 to −114 is critical for the promoter function, and the region from −374 to −114 may contain core cis-elements involved in light response. In transgenic rice expressing GUS under the control of OrGSEp, visualization and quantification of GUS activity showed that GUS was preferentially expressed in green tissues and not in endosperm. OrGSEp is a useful regulatory element for breeding pest-resistant crops.

Early responses to dehydration in contrasting wild Arachis species

Wild peanut relatives (Arachis spp.) are genetically diverse and were selected throughout evolution to a range of environments constituting, therefore, an important source of allelic diversity for abiotic stress tolerance. In particular, A. duranensis and A. stenosperma, the parents of the reference Arachis A-genome genetic map, show contrasting transpiration behavior under limited water conditions. This study aimed to build a comprehensive gene expression profile of these two wild species under dehydration stress caused by the withdrawal of hydroponic nutrient solution. For this purpose, roots of both genotypes were collected at seven time-points during the early stages of dehydration and used to construct cDNA paired-end libraries. Physiological analyses indicated initial differences in gas exchange parameters between the drought-tolerant genotype of A. duranensis and the drought-sensitive genotype of A. stenosperma. High-quality Illumina reads were mapped against the A. duranensis reference genome and resulted in the identification of 1,235 and 799 Differentially Expressed Genes (DEGs) that responded to the stress treatment in roots of A. duranensis and A. stenosperma, respectively. Further analysis, including functional annotation and identification of biological pathways represented by these DEGs confirmed the distinct gene expression behavior of the two contrasting Arachis species genotypes under dehydration stress. Some species-exclusive and common DEGs were then selected for qRT-PCR analysis, which corroborated the in silico expression profiling. These included genes coding for regulators and effectors involved in drought tolerance responses, such as activation of osmosensing molecular cascades, control of hormone and osmolyte content, and protection of macromolecules. This dataset of transcripts induced during the dehydration process in two wild Arachis genotypes constitute new tools for the understanding of the distinct gene regulation processes in these closely related species but with contrasting drought responsiveness. In addition, our findings provide insights into the nature of drought tolerance in wild germoplasm, which might be explored as novel sources of diversity and useful wild alleles to develop climate-resilient crop varieties.

De novo leaf and root transcriptome analysis to identify putative genes involved in triterpenoid saponins biosynthesis in Hedera helix L

Hedera helix L. is an important traditional medicinal plant in Europe. The main active components are triterpenoid saponins, but none of the potential enzymes involved in triterpenoid saponins biosynthesis have been discovered and annotated. Here is reported the first study of global transcriptome analyses using the Illumina HiSeq™ 2500 platform for H. helix. In total, over 24 million clean reads were produced and 96,333 unigenes were assembled, with an average length of 1385 nt; more than 79,085 unigenes had at least one significant match to an existing gene model. Differentially Expressed Gene analysis identified 6,222 and 7,012 unigenes which were expressed either higher or lower in leaf samples when compared with roots. After functional annotation and classification, two pathways and 410 unigenes related to triterpenoid saponins biosynthesis were discovered. The accuracy of these de novo sequences was validated by RT-qPCR analysis and a RACE clone. These data will enrich our knowledge of triterpenoid saponin biosynthesis and provide a theoretical foundation for molecular research on H. helix.

Selection and Validation of Appropriate Reference Genes for qRT-PCR Analysis in Isatis indigotica Fort

Due to its sensitivity and specificity, real-time quantitative PCR (qRT-PCR) is a popular technique for investigating gene expression levels in plants. Based on the Minimum Information for Publication of Real-Time Quantitative PCR Experiments (MIQE) guidelines, it is necessary to select and validate putative appropriate reference genes for qRT-PCR normalization. In the current study, three algorithms, geNorm, NormFinder, and BestKeeper, were applied to assess the expression stability of 10 candidate reference genes across five different tissues and three different abiotic stresses in Isatis indigotica Fort. Additionally, the IiYUC6 gene associated with IAA biosynthesis was applied to validate the candidate reference genes. The analysis results of the geNorm, NormFinder, and BestKeeper algorithms indicated certain differences for the different sample sets and different experiment conditions. Considering all of the algorithms, PP2A-4 and TUB4 were recommended as the most stable reference genes for total and different tissue samples, respectively. Moreover, RPL15 and PP2A-4 were considered to be the most suitable reference genes for abiotic stress treatments. The obtained experimental results might contribute to improved accuracy and credibility for the expression levels of target genes by qRT-PCR normalization in I. indigotica.

Effects of drought stress on global gene expression profile in leaf and root samples of Dongxiang wild rice (Oryza rufipogon)

Drought is a serious constraint to rice production throughout the world, and although Dongxiang wild rice (Oryza rufipogon, DXWR) possesses a high degree of drought resistance, the underlying mechanisms of this trait remains unclear. In the present study, cDNA libraries were constructed from the leaf and root tissues of drought-stressed and untreated DXWR seedlings, and transcriptome sequencing was performed with the goal of elucidating the molecular mechanisms involved in drought-stress response. The results indicated that 11231 transcripts were differentially expressed in the leaves (4040 up-regulated and 7191 down-regulated) and 7025 transcripts were differentially expressed in the roots (3097 up-regulated and 3928 down-regulated). Among these differentially expressed genes (DEGs), the detection of many transcriptional factors and functional genes demonstrated that multiple regulatory pathways were involved in drought resistance. Meanwhile, the DEGs were also annotated with gene ontology (GO) terms and key pathways via functional classification and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway mapping, respectively. A set of the most interesting candidate genes was then identified by combining the DEGs with previously identified drought-resistant quantitative trait loci (QTL). The present work provides abundant genomic information for functional dissection of the drought resistance of DXWR, and findings will further help the current understanding of the biological regulatory mechanisms of drought resistance in plants and facilitate the breeding of new drought-resistant rice cultivars.

A walk on the wild side: Oryza species as source for rice abiotic stress tolerance

Oryza sativa, the common cultivated rice, is one of the most important crops for human consumption, but production is increasingly threatened by abiotic stresses. Although many efforts have resulted in breeding rice cultivars that are relatively tolerant to their local environments, climate changes and population increase are expected to soon call for new, fast generation of stress tolerant rice germplasm, and current within-species rice diversity might not be enough to overcome such needs. The Oryza genus contains other 23 wild species, with only Oryza glaberrima being also domesticated. Rice domestication was performed with a narrow genetic diversity, and the other Oryza species are a virtually untapped genetic resource for rice stress tolerance improvement. Here we review the origin of domesticated Oryza sativa from wild progenitors, the ecological and genomic diversity of the Oryza genus, and the stress tolerance variation observed for wild Oryza species, including the genetic basis underlying the tolerance mechanisms found. The summary provided here is important to indicate how we should move forward to unlock the full potential of these germplasms for rice improvement.

Identification of microRNAs in Response to Drought in Common Wild Rice (Oryza rufipogon Griff.) Shoots and Roots

Background

Drought is the most important factor that limits rice production in drought-prone environments. Plant microRNAs (miRNAs) are involved in biotic and abiotic stress responses. Common wild rice (Oryza rufipogon Griff.) contains abundant drought-resistant genes, which provide an opportunity to explore these excellent resources as contributors to improve rice resistance, productivity, and quality.

Results

In this study, we constructed four small RNA libraries, called CL and CR from PEG6000-free samples and DL and DR from PEG6000-treated samples, where ‘R’ indicates the root tissue and ‘L’ indicates the shoot tissue. A total of 200 miRNAs were identified to be differentially expressed under the drought-treated conditions (16% PEG6000 for 24 h), and the changes in the miRNA expression profile of the shoot were distinct from those of the root. At the miRNA level, 77 known miRNAs, which belong to 23 families, including 40 up-regulated and 37 down-regulated in the shoot, and 85 known miRNAs in 46 families, including 65 up-regulated and 20 down-regulated in the root, were identified as differentially expressed. In addition, we predicted 26 new miRNA candidates from the shoot and 43 from the root that were differentially expressed during the drought stress. The quantitative real-time PCR analysis results were consistent with high-throughput sequencing data. Moreover, 88 miRNAs that were differentially-expressed were predicted to match with 197 targets for drought-stress.

Conclusion

Our results suggest that the miRNAs of O. rufipogon are responsive to drought stress. The differentially expressed miRNAs that are tissue-specific under drought conditions could play different roles in the regulation of the auxin pathway, the flowering pathway, the drought pathway, and lateral root formation. Thus, the present study provides an account of tissue-specific miRNAs that are involved in the drought adaption of O. rufipogon.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock

BACKGROUND

Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited.

RESULTS

We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes.

CONCLUSIONS

This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.

Genetic determination of the enhanced drought resistance of rice maintainer HuHan2B by pedigree breeding

The ongoing deficit of fresh water resource in rice growing regions has made the selection of water-saving and drought-resistance rice (WDR) a crucial factor in developing sustainable cultivation. HuHan2B, a new japonica maintainer for WDR breeding, had the same yield potential as recurrent parent HanFengB but showed improved drought resistance in fields. We investigated the genomic content accumulation and candidate genes passed from parent to offspring using the genomic and transcriptomic approaches. The genomic constitution indicated that the genetic similarity was 84% between HuHan2B and HanFengB; additionally, 7,256 genes with specific alleles were inherited by HuHan2B from parents other than HanFengB. The differentially expressed genes (DEGs) under drought stress showed that biological function was significantly enriched for transcript regulation in HuHan2B, while the oxidation-reduction process was primarily enriched in HanFengB. Furthermore, 36 DEGs with specific inherited alleles in HuHan2B were almost involved in the regulatory network of TFs and target genes. These findings suggested that major-effect genes were congregated and transformed into offspring in manner of interacting network by breeding. Thus, exploiting the potential biological function of allelic-influencing DEGs would be of great importance for improving selection efficiency and the overall genetic gain of multiple complex traits.