Heterosis in early maize ear inflorescence development: a genome-wide transcription analysis for two maize inbred lines and their hybrid

Mapping and Functional Analysis of QTL for Kernel Number per Row in Tropical and Temperate-Tropical Introgression Lines of Maize (Zea mays L.)

Kernel number per row (KNR) is an essential component of maize (Zea mays L.) grain yield (GY), and understanding its genetic mechanism is crucial to improve GY. In this study, two F₇ recombinant inbred line (RIL) populations were created using a temperate-tropical introgression line TML418 and a tropical inbred line CML312 as female parents and a backbone maize inbred line Ye107 as the common male parent. Bi-parental quantitative trait locus (QTL) mapping and genome-wide association analysis (GWAS) were then performed on 399 lines of the two maize RIL populations for KNR in two different environments using 4118 validated single nucleotide polymorphism (SNP) markers. This study aimed to: (1) detect molecular markers and/or the genomic regions associated with KNR; (2) identify the candidate genes controlling KNR; and (3) analyze whether the candidate genes are useful in improving GY. The authors reported a total of 7 QTLs tightly linked to KNR through bi-parental QTL mapping and identified 21 SNPs significantly associated with KNR through GWAS. Among these, a highly confident locus qKNR7-1 was detected at two locations, Dehong and Baoshan, with both mapping approaches. At this locus, three novel candidate genes (Zm00001d022202, Zm00001d022168, Zm00001d022169) were identified to be associated with KNR. These candidate genes were primarily involved in the processes related to compound metabolism, biosynthesis, protein modification, degradation, and denaturation, all of which were related to the inflorescence development affecting KNR. These three candidate genes were not reported previously and are considered new candidate genes for KNR. The progeny of the hybrid Ye107 × TML418 exhibited strong heterosis for KNR, which the authors believe might be related to qKNR7-1. This study provides a theoretical foundation for future research on the genetic mechanism underlying KNR in maize and the use of heterotic patterns to develop high-yielding hybrids.

Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines

Heterosis is the phenomenon in which some hybrid traits are superior to those of their parents. Most studies have analyzed the heterosis of agronomic traits of crops; however, heterosis of the panicles can improve yield and is important for crop breeding. Therefore, a systematic study of panicle heterosis is needed, especially during the reproductive stage. RNA sequencing (RNA Seq) and transcriptome analysis are suitable for further study of heterosis. Using the Illumina Nova Seq platform, the transcriptome of ZhongZheYou 10 (ZZY10), an elite rice hybrid, the maintainer line ZhongZhe B (ZZB), and the restorer line Z7-10 were analyzed at the heading date in Hangzhou, 2022. 581 million high-quality short reads were obtained by sequencing and were aligned against the Nipponbare reference genome. A total of 9000 differential expression genes were found between the hybrids and their parents (DG_HP). Of the DG_HP, 60.71% were up-regulated and 39.29% were down-regulated in the hybrid. Comparative transcriptome analysis revealed that 5235 and 3765 DG_HP were between ZZY10 and ZhongZhe B and between ZZY10 and Z7-10, respectively. This result is consistent with the transcriptome profile of ZZY10 and was similar to Z7-10. The expression patterns of DG_HP mainly exhibited over-dominance, under-dominance, and additivity. Among the DG_HP-involved GO terms, pathways such as photosynthesis, DNA integration, cell wall modification, thylakoid, and photosystem were significant. 21 DG_HP, which were involved in photosynthesis, and 17 random DG_HP were selected for qRT-PCR validation. The up-regulated PsbQ and down-regulated subunits of PSI and PSII and photosynthetic electron transport in the photosynthesis pathway were observed in our study. Extensive transcriptome data were obtained by RNA-Seq, providing a comprehensive overview of panicle transcriptomes at the heading stage in a heterotic hybrid.

Transcriptome Analysis Reveals a Major Gene Expression Pattern and Important Metabolic Pathways in the Control of Heterosis in Chinese Cabbage

Although heterosis is commonly used in Chinese cabbage, its molecular basis is poorly understood. In this study, 16Chinese cabbage hybrids were utilized as test subjects to explore the potential molecular mechanism of heterosis. RNA sequencing revealed 5815-10,252 differentially expressed genes (DEGs) (female parent vs. male parent), 1796-5990 DEGs (female parent-vs-hybrid), and 2244-7063 DEGs (male parent vs. hybrid) in 16 cross combinations at the middle stage of heading. Among of them, 72.83-84.20% DEGs conformed to the dominant expression pattern, which is the predominant expression pattern in hybrids. There were 13 pathways in which DEGs were significantly enriched in most cross combinations. Among them, the plant-pathogen interaction (ko04626) and circadian rhythm-plant (ko04712)were significantly enriched by DEGs in strong heterosis hybrids. WGCNA also proved that the two pathways were significantly related to heterosis in Chinese cabbage.

Comparative transcriptomics analysis at the key stage of maize ear development dissect heterosis

Important traits related to maize (Zea mays L.) grain yield, such as kernel row number, ear length, kernel number per row, are determined during the development of female inflorescence. There is a significant positive correlation between yield component and the activity of inflorescence meristem (IM). To find the key stage of heterosis in the development of the ear, immature ears (from the IM stage until the end of the floral meristem [FM] stage) of Yudan888 and its parent lines were sampled to assay phenotype and for comparative transcriptomics analysis. The immature ear length of Yudan888 at the IM stage fitted an additive (mid-parental) model, but it showed high parental dominance at the spikelet-pair meristem (SPM) stage. Comparative analysis of transcriptomes suggested significant differences between additive and nonadditive expression patterns for different developmental stages. The number of distinct maternal or paternal genes (DMP) (genes expressed only in one parental line and their hybrid but silenced in another line) was greater than ABF1 (genes expressed in both parental lines but silenced in hybrid) at each stage. Gene Ontology (GO) enrichment suggested that the cell redox homeostasis genes with overdominance expression patterns in hybrids have an important contribution to heterosis. According to our research, an ear length heterosis network was established. The discovery of the inflection point for ear length heterosis allows us for inferring that the transition state of IM to SPM may be the starting point of ear length heterosis. These findings improved the understanding of maize ear length heterosis.

Dynamic patterns of gene expression and regulatory variation in the maize seed coat

BACKGROUND

Seed size is an important factor contributing to maize yield, but its molecular mechanism remains unclear. The seed coat, which serves as one of the three components of the maize grain, determines seed size to a certain extent. The seed coat also shares the maternal genotype and is an ideal material for studying heterosis.

RESULTS

In this study, the self-pollinated seeds of the maize hybrid Yudan888 and its parental lines were continuously collected from 0 day after pollination (DAP) to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP and 8 DAP are the best time points to study maize seed coat heterosis. Cytological observations indicated that maize seed coat heterosis might be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results showed that the nonadditive genes changed significantly between 3 and 8 DAP. However, the number of genes expressed additively was not significantly different. Our findings suggest that seed coat heterosis in hybrid is the result of nonadditive expression caused by dynamic changes in genes at different time points during seed expansion and seed coat development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment indicated that genes related to DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation contributed significantly to hybrid seed coat heterosis.

CONCLUSION

Maize seed coat phenotyping allowed us to infer that 3 DAP and 8 DAP are important time points in the study of seed coat heterosis. Our findings provide evidence for genes involved in DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation in hybrid with high or low parental expression as major contributors to hybrid seed coat heterosis.

Integrated Analysis of mRNA and Non-coding RNA Transcriptome in Pepper (Capsicum chinense) Hybrid at Seedling and Flowering Stages

Pepper is an important vegetable in the world. In this work, mRNA and ncRNA transcriptome profiles were applied to understand the heterosis effect on the alteration in the gene expression at the seedling and flowering stages between the hybrid and its parents in Capsicum chinense. Our phenotypic data indicated that the hybrid has dominance in leaf area, plant scope, plant height, and fruit-related traits. Kyoto Encyclopedia of Genes and Genomes analysis showed that nine members of the plant hormone signal transduction pathway were upregulated in the seedling and flowering stages of the hybrid, which was supported by weighted gene coexpression network analysis and that BC332_23046 (auxin response factor 8), BC332_18317 (auxin-responsive protein IAA20), BC332_13398 (ethylene-responsive transcription factor), and BC332_27606 (ethylene-responsive transcription factor WIN1) were candidate hub genes, suggesting the important potential role of the plant hormone signal transduction in pepper heterosis. Furthermore, some transcription factor families, including bHLH, MYB, and HSF were greatly over-dominant. We also identified 2,525 long ncRNAs (lncRNAs), 47 micro RNAs (miRNAs), and 71 circle RNAs (circRNAs) in the hybrid. In particular, downregulation of miR156, miR169, and miR369 in the hybrid suggested their relationship with pepper growth vigor. Moreover, we constructed some lncRNA–miRNA–mRNA regulatory networks that showed a multi-dimension to understand the ncRNA relationship with heterosis. These results will provide guidance for a better understanding of the molecular mechanism involved in pepper heterosis.

An enhanced photosynthesis and carbohydrate metabolic capability contributes to heterosis of the cotton (Gossypium hirsutum) hybrid 'Huaza Mian H318', as revealed by genome-wide gene expression analysis

Background

Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton.

Results

The elite cotton hybrid variety ‘Huaza Mian H318’ (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4–5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318.

Conclusion

The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07580-8.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

Robust non-syntenic gene expression patterns in diverse maize hybrids during root development

Distantly related maize (Zea mays L.) inbred lines exhibit an exceptional degree of structural genomic diversity, which is probably unique among plants. This study systematically investigated the developmental and genotype-dependent regulation of the primary root transcriptomes of a genetically diverse panel of maize F1-hybrids and their parental inbred lines. While we observed substantial transcriptomic changes during primary root development, we demonstrated that hybrid-associated gene expression patterns, including differential, non-additive, and allele-specific transcriptome profiles, are particularly robust to these developmental fluctuations. For instance, differentially expressed genes with preferential expression in hybrids were highly conserved during development in comparison to their parental counterparts. Similarly, in hybrids a major proportion of non-additively expressed genes with expression levels between the parental values were particularly conserved during development. Importantly, in these expression patterns non-syntenic genes that evolved after the separation of the maize and sorghum lineages were systemically enriched. Furthermore, non-syntenic genes were substantially linked to the conservation of all surveyed gene expression patterns during primary root development. Among all F1-hybrids, between ~40% of the non-syntenic genes with unexpected allelic expression ratios and ~60% of the non-syntenic differentially and non-additively expressed genes were conserved and therefore robust to developmental changes. Hence, the enrichment of non-syntenic genes during primary root development might be involved in the developmental adaptation of maize roots and thus the superior performance of hybrids.

Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea

Heterosis refers to the superior performance of F1 hybrids over their respective parental inbred lines. Although the genetic and expression basis of heterosis have been previously investigated, the metabolic basis for this phenomenon is poorly understood. In a preliminary morphological study in Brassica juncea, we observed significant heterosis at the 50% flowering stage, wherein both the growth and reproduction of F1 reciprocal hybrids were greater than that of their parents. To identify the possible metabolic causes or consequences of this heterosis, we carried out targeted LC-MS analysis of 48 primary (amino acids and sugars) and secondary metabolites (phytohormones, glucosinolates, flavonoids, and phenolic esters) in five developmental tissues at 50% flowering in hybrids and inbred parents. Principal component analysis (PCA) of metabolites clearly separated inbred lines from their hybrids, particularly in the bud tissues. In general, secondary metabolites displayed more negative heterosis values in comparison to primary metabolites. The tested primary and secondary metabolites displayed both additive and non-additive modes of inheritance in F1 hybrids, wherein the number of metabolites showing an additive mode of inheritance were higher in buds and siliques (52.77-97.14%) compared to leaf tissues (47.37-80%). Partial least regression (PLS) analysis further showed that primary metabolites, in general, displayed higher association with morphological parameters in F1 hybrids. Overall, our results are consistent with a resource-cost model for heterosis in B. juncea, where metabolite allocation in hybrids appears to favor growth, at the expense of secondary metabolism.

Genome-wide proteomic profiling reveals the role of dominance protein expression in heterosis in immature maize ears

Heterosis refers to the phenomenon in which hybrid progeny show superior performance relative to their parents. Early maize ear development shows strong heterosis in ear architecture traits and greatly affects grain yield. To explore the underlying molecular mechanisms, genome-wide proteomics of immature ears of maize hybrid ZD909 and its parents were analyzed using tandem mass tag (TMT) technology. A total of 9,713 proteins were identified in all three genotypes. Among them, 3,752 (38.6%) proteins were differentially expressed between ZD909 and its parents. Multiple modes of protein action were discovered in the hybrid, while dominance expression patterns accounted for 63.6% of the total differentially expressed proteins (DEPs). Protein pathway enrichment analysis revealed that high parent dominance proteins mainly participated in carbon metabolism and nitrogen assimilation processes. Our results suggested that the dominant expression of favorable alleles related to C/N metabolism in the hybrid may be essential for ZD909 ear growth and heterosis formation. Integrated analysis of proteomic and quantitative trait locus (QTL) data further support our DEP identification and provide useful information for the discovery of genes associated with ear development. Our study provides comprehensive insight into the molecular mechanisms underlying heterosis in immature maize ears from a proteomic perspective.

QTL for Maize Midparent Heterosis in the Heterotic Pattern American Dent × European Flint under Corn Borer Pressure

Despite the importance of heterosis and the efforts to comprehend this phenomenon, its molecular bases are still unknown. In this study, we intended to detect Quantitative trait loci (QTL) for mid-parent heterosis under infestation with the Mediterranean corn borer (MCB, Sesamia nonagrioides Lef.) using a North Carolina design III approach with a RIL population derived from a European flint inbred (EP42) × American dent inbred (A637) cross. QTL for heterosis of kernel yield have been positioned in regions corresponding to previously identified QTL for the same trait in different backgrounds. These results reinforce the high congruency of genes controlling heterosis across populations, even when populations have been developed from different heterotic patterns. A high percentage of genetic variation for mid-parent heterosis (Z2) for kernel yield could not be explained. Furthermore, genomic regions involved in heterosis for yield and plant height were not found despite the high genetic correlation between Z2 transformations for kernel yield and plant height. The moderate power in detecting QTL for mid-parent heterosis suggests that many genes with low augmented dominance effects contribute to the genetic architecture of mid-parent heterosis; dominance and additive-additive epistatic effects could also contribute to heterosis. However, results from this and previous studies suggest that the region 8.03-8.05 deserves special attention in future works in order to fine map loci involved in mid-parent heterosis for yield.

Transcriptome analysis of neo-tetraploid rice reveals specific differential gene expressions associated with fertility and heterosis

Polyploid rice hybrids have a powerful biological and yield potential that may become a new way for rice breeding; however, low fertility is major hindrance in commercial utilization. Here, we developed a neo-tetraploid rice that could overcome the sterility of autotetraploid rice and produce high heterosis. Transcriptome analysis of F₁ hybrid developed by crossing neo-tetraploid with autotetraploid rice displayed 807, 663 and 866 differentially expressed genes that uniquely associated with F₁ and specific to (DEGFu-sp) anther, ovary and leaf, respectively. Of the DEGFu-sp, 1224 genes displayed nonadditive expression; 44 and 10 genes were annotated as TFs and methyltransferase or hydroxymethyltransferase, respectively. Gene ontology enrichment and co-expression analysis revealed specific differential gene expressions in the DEGFu-sp to leaf, anther and ovary, such as genes related to photosynthesis, metabolic process and transport, and co-expression network including fertility, resistance and epigenetic elements. Of the DEGFu-sp to anther, 42 meiosis stage-specific genes, eight meiosis-related genes, such as RAD51 and SMC2, were identified. We identified 38 miRNAs from DEGFu-sp to anther, and their targets were associated with pollen fertility and retrotransposon protein. Our study provides new germplasm for polyploid rice breeding, and revealed complex regulatory mechanisms that might be associated with heterosis and fertility.

Nonsyntenic Genes Drive Tissue-Specific Dynamics of Differential, Nonadditive, and Allelic Expression Patterns in Maize Hybrids

Distantly related maize (Zea mays) inbred lines display an exceptional degree of genomic diversity. F1 progeny of such inbred lines are often more vigorous than their parents, a phenomenon known as heterosis. In this study, we investigated how the genetic divergence of the maize inbred lines B73 and Mo17 and their F1 hybrid progeny is reflected in differential, nonadditive, and allelic expression patterns in primary root tissues. In pairwise comparisons of the four genotypes, the number of differentially expressed genes between the two parental inbred lines significantly exceeded those of parent versus hybrid comparisons in all four tissues under analysis. No differentially expressed genes were detected between reciprocal hybrids, which share the same nuclear genome. Moreover, hundreds of nonadditive and allelic expression ratios that were different from the expression ratios of the parents were observed in the reciprocal hybrids. The overlap of both nonadditive and allelic expression patterns in the reciprocal hybrids significantly exceeded the expected values. For all studied types of expression - differential, nonadditive, and allelic - substantial tissue-specific plasticity was observed. Significantly, nonsyntenic genes that evolved after the last whole genome duplication of a maize progenitor from genes with synteny to sorghum (Sorghum bicolor) were highly overrepresented among differential, nonadditive, and allelic expression patterns compared with the fraction of these genes among all expressed genes. This observation underscores the role of nonsyntenic genes in shaping the transcriptomic landscape of maize hybrids during the early developmental manifestation of heterosis in root tissues of maize hybrids.

Comparative transcriptome analysis among parental inbred and crosses reveals the role of dominance gene expression in heterosis in Drosophila melanogaster

We observed heteroses for body weight in Drosophila melanogaster after generating hybrids from three inbred lines. To better understand the mechanism for this phenomenon at the mRNA level, we compared the mRNA profiles of the parental and hybrid lines using high-throughput RNA-seq. A total of 5877 differentially expressed genes (DEGs) were found and about 92% of these exhibited parental expression level dominance. Genes in the dominance category were functionally characterized using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the gene classifications offered by the Gene Ontology (GO) Consortium. The analysis identified genes associated with crucial processes such as development and growth in all three crosses. Functional assignments involving aminoglycan metabolism, starch and sucrose metabolism, and galactose metabolism are significantly overrepresented amongst the 215 common dominance DEGs. We conclude that dominance DEGs are important in heteroses in Drosophila melanogaster and contribute specifically to body weight heterosis.

Transcriptome Analysis of Interspecific Hybrid between Brassica napus and B. rapa Reveals Heterosis for Oil Rape Improvement

The hybrid between Brassica napus and B. rapa displays obvious heterosis in both growth performance and stress tolerances. A comparative transcriptome analysis for B. napus (AⁿAⁿCC genome), B. rapa (A^rA^r genome), and its hybrid F1 (AⁿA^rC genome) was carried out to reveal the possible molecular mechanisms of heterosis at the gene expression level. A total of 40,320 nonredundant unigenes were identified using B. rapa (AA genome) and B. oleracea (CC genome) as reference genomes. A total of 6,816 differentially expressed genes (DEGs) were mapped in the A and C genomes with 4,946 DEGs displayed nonadditively by comparing the gene expression patterns among the three samples. The coexistence of nonadditive DEGs including high-parent dominance, low-parent dominance, overdominance, and underdominance was observed in the gene action modes of F1 hybrid, which were potentially related to the heterosis. The coexistence of multiple gene actions in the hybrid was observed and provided a list of candidate genes and pathways for heterosis. The expression bias of transposable element-associated genes was also observed in the hybrid compared to their parents. The present study could be helpful for the better understanding of the determination and regulation of mechanisms of heterosis to aid Brassica improvement.

Hybrid mimics and hybrid vigor in Arabidopsis

F1 hybrids can outperform their parents in yield and vegetative biomass, features of hybrid vigor that form the basis of the hybrid seed industry. The yield advantage of the F1 is lost in the F2 and subsequent generations. In Arabidopsis, from F2 plants that have a F1-like phenotype, we have by recurrent selection produced pure breeding F5/F6 lines, hybrid mimics, in which the characteristics of the F1 hybrid are stabilized. These hybrid mimic lines, like the F1 hybrid, have larger leaves than the parent plant, and the leaves have increased photosynthetic cell numbers, and in some lines, increased size of cells, suggesting an increased supply of photosynthate. A comparison of the differentially expressed genes in the F1 hybrid with those of eight hybrid mimic lines identified metabolic pathways altered in both; these pathways include down-regulation of defense response pathways and altered abiotic response pathways. F6 hybrid mimic lines are mostly homozygous at each locus in the genome and yet retain the large F1-like phenotype. Many alleles in the F6 plants, when they are homozygous, have expression levels different to the level in the parent. We consider this altered expression to be a consequence of transregulation of genes from one parent by genes from the other parent. Transregulation could also arise from epigenetic modifications in the F1. The pure breeding hybrid mimics have been valuable in probing the mechanisms of hybrid vigor and may also prove to be useful hybrid vigor equivalents in agriculture.

Genome-wide identification and analysis of drought-responsive genes and microRNAs in tobacco

Drought stress response is a complex trait regulated at transcriptional and post-transcriptional levels in tobacco. Since the 1990s, many studies have shown that miRNAs act in many ways to regulate target expression in plant growth, development and stress response. The recent draft genome sequence of Nicotiana benthamiana has provided a framework for Digital Gene Expression (DGE) and small RNA sequencing to understand patterns of transcription in the context of plant response to environmental stress. We sequenced and analyzed three Digital Gene Expression (DGE) libraries from roots of normal and drought-stressed tobacco plants, and four small RNA populations from roots, stems and leaves of control or drought-treated tobacco plants, respectively. We identified 276 candidate drought responsive genes (DRGs) with sequence similarities to 64 known DRGs from other model plant crops, 82 were transcription factors (TFs) including WRKY, NAC, ERF and bZIP families. Of these tobacco DRGs, 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. An integrated gene regulatory network is proposed for the molecular mechanisms of tobacco root response to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and their target transcripts. This network analysis serves as a reference for future studies on tobacco response stresses such as drought, cold and heavy metals.