Global insights into homoeolog gene expression in upland cotton under intraspecific hybridization

doi:10.21203/rs.3.rs-234267/v1

Download PDF

Research

Global insights into homoeolog gene expression in upland cotton under intraspecific hybridization

https://doi.org/10.21203/rs.3.rs-234267/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Background: Hybridization is useful to enhance yield potential of agronomic crops in the world. Cotton has genome doubling due to alloteraploidy process and hybridization process in coordinate with duplicated genome can produce more yield and adaptability. Therefore, expression of homoeolog gene pairs between hybrids and inbred parents are vital to characterize genetic source of heterosis in cotton.

Results: Investigation results of homoeolog gene pairs between two contrasting hybrids and their respective inbred parents identified 36853 homoeolog genes in hybrids. It was observed both high and low hybrids had similar trends in homoeolog gene expression patterns in each tissue under study. An average of 96% of homoeolog genes had no biased expression and their expressions were derived from the equal contribution of both parents. Besides, very few homoeolog genes (An average of 1%) showed no biased or novel expression in both hybrids. The functional analysis described secondary metabolic pathways had a majority of novel biased homoeolog genes in hybrids.

Conclusions: These results contribute preliminary knowledge about how hybridization affects expression patterns of homoeolog gene pairs in upland cotton hybrids. Our study also highlights the functional genomics of metabolic genes to explore the genetic mechanism of heterosis in cotton.

Molecular Genetics

Structural Biology

Agricultural Engineering

Cotton

Heterosis

Homoeolog expression

RNA-Seq

Secondary metabolism

Hybridization or polyploidy is a vital process, produces exceptional phenotypes through the interaction of two distinct genomes, and has a substantial influence on the plant genome (Dufresne et al. 1994; Yoo et al. 2013). The novel phenotypes in plants can be consequent of many modifications not only at genomic and epigenetic levels but also in gene expression patterns (Doyle et al. 2008; Kashkush et al. 2003). Several examples in cotton (Flagel et al. 2010) Brassica napus (Gaeta et al. 2009), and Spartina (Chelaifa et al. 2010) stated that hybridization (genome merger) modified more transcriptome than polyploidy. About 15% of flowering plants are polyploidy (Wood et al. 2009) and divided into categories of autopolyploid and allopolyploid. Autopolyploid is evolved from a genome merger of the same or similar species but allopolyploids contain a genomic combination of two different species (Doyle et al. 2008). Allopolyploids have enhanced growth and wider adaptations relative to autopolyploid. Many important crops are allopolyploids together with upland cotton (Gossypium hirsutum). In allopolyploids, hybridization or polyploidy is allied with various kinds of exchanges and interactions between homoeolog genes (duplicate genes) (Kovarik et al. 2008; Salmon et al. 2010; Szadkowski et al. 2010). Further modifications are consisting of epigenetic (Doyle et al. 2008; Madlung et al. 2002), gene activation (Kashkush et al. 2003; Kraitshtein et al. 2010), and sequence removal (Anssour et al. 2009; Han et al. 2005). The result of different genomic interactions by hybridization or polyploidy produced several expression patterns of duplicate or genes with similar functions. For instance, these may be gene silencing, gene activation, novel expression, and altered parental contribution (Doyle et al. 2008; Li et al. 2014; Yoo et al. 2013). However, the direction and magnitude of duplicate gene expression are varied with crop, tissue, and organ (Adams et al. 2003; Bottley et al. 2006; Buggs et al. 2010; Chaudhary et al. 2009).

Previous studies have identified four different expression forms of duplicated genes in plants e.g., one duplicate gene copy may be lost its expression pattern as compared to other (De Smet et al. 2013), duplicate gene copies may maintain original expression (McGrath et al. 2014), one duplicate gene copy may adopt new expression (Hughes et al. 2014; Liu et al. 2010), and sometimes original expression is the equal contribution of duplicates (Kleinjan et al. 2008). Importantly, the duplicate gene pairs may show homoeolog expression bias in allopolyploid (Grover et al. 2012). Where, one homoeolog display preferential expression to the transcriptome as compared to the other. The phenomena of homoeolog expression have already been stated for allopolyploids like cotton (Flagel et al. 2008; Hovav et al. 2008; Samuel Yang et al. 2006), wheat (Bottley et al. 2006; Mochida et al. 2004), Arabidopsis (Chang et al. 2010; Wang et al. 2004), and brassica ( Auger et al. 2009; Gaeta et al. 2007). The analysis results in cotton showed that hybridization caused more homoeolog expression changes and fewer changes in gene expressions (Yoo et al. 2013). In contrast, polyploidy formation is linked with more gene expression modifications and fewer changes in homoeolog expression. A study in cyprinid fish also reported hybridization brought more homoeolog expression changes than polyploidy (Ren et al. 2016). Homoeolog expression can be either balanced or unbalanced (Grover et al. 2012). When duplicate gene pairs display biased expression towards one parental genome is almost equal to the biased expression towards other parental genomes. It shows balanced homoeolog expression. On the other hand, an unbalanced homoeolog expression change states to a situation where more duplicate gene pairs demonstrate expression twisted toward one parental genome. An unbalanced form of homoeolog expression is mostly reported in allopolyploids (Akhunova et al. 2010; Schnable et al. 2011; Wang et al. 2006).

Upland cotton is a prime fiber, an important source of edible oil, and cash crops in the world (Shahzad et al. 2019a). The demand for cotton fiber in the textile industry is increasing. While supply is reducing due to low area of cultivation and staggered yield potential of cultivars. Heterosis is and will continue to be a basic breeding tool to produce superior genotypes (Birchler et al. 2003). Many researchers reported meaningful heterosis in seed cotton yield and fiber quality traits in upland cotton (Ali et al. 2016; Li et al. 2019). Further, cotton genotypes produced through the utilization of heterosis exhibit more stability and adaptability to changing the environment (Shahzad et al. 2019b). The genetic mechanism of heterosis is a perturbed and so-called genetic model of dominance (Xiao et al. 1995), overdominance (Shahzad et al. 2020a), and epistasis (Yu et al. 1997) does not reveal it. However, progress has already been achieved in the recent few years with advanced molecular research. For instance, high-throughput transcriptome sequencing is a useful technology to get accurate knowledge about global gene expressions. This also creates a possibility to discriminate the expression of homologous genes (Shahzad et al. 2020b; Wang et al. 2009; Wu et al. 2018). To date, many previous studies in agronomic crops analyzed patterns of homoeolog expression in allotetraploids and their diploid progenitors e.g., cotton (Yoo et al. 2013) and brassica (Li et al. 2020). However, information regarding how hybridization effect homoeolog expression changes in newly formed upland cotton hybrids is limited. Thus, we conducted a study to investigate the effect of hybridization on homoeolog expression modifications on a genome-wide scale. The direction and extent of homoeologue expressions between contrasting F₁ hybrids and their reciprocal parents were analyzed by RNA-seq of leaf, root, and flower buds. The result will provide new insights into how homoeologue expressions contribute to growth heterosis in upland cotton and enhance our genetic knowledge for hybridization process.

Plant material and RNA sequencing

The plant material used in this study included two seed cotton yield contrasting F₁ hybrids of upland cotton and their three reciprocal inbred parents named as A, B, and D. The high yielding hybrid (H) was cross of inbred A and B. whereas, a low hybrid was a combination of cross A and D. The seedling was raised in the lab for all these materials, root and leaf samples were picked at true leaf stage, and flower buds were picked at squaring stage from the field in three biological repeats. Total RNA was isolated using the Sigma Spectrum Plant Total RNA kit (Sigma-Aldrich, USA) as described by the manufacturing company. After quality checks and preparations, paired-end sequencing (300 ± 50 bp) on a second-generation Illumina Hiseq 4000 was performed for hybrid H and L. These sequencing details can be seen at NCBI with accession number GSE144676 (Shahzad et al. 2020a) and GSE150052 (Shahzad et al. 2020b). The full-length Pacbio transcriptome sequencing (Unpublished) for parental A, B, and D were used for a better comparison of homoeolog expression patterns between parents and progeny.

Analysis of homoeolog expression changes and functional enrichment

To find homologous sequences between reciprocal parents of hybrid H and L, the full-length Pacbio sequences of the two parents were homologously aligned using the reciprocal BLAST (BLASTN) hit method with screening index of e-20 (Ren et al. 2016). Two sequences were homologous, if each of them was the best hit of the other, and if the sequences were aligned over 300 bp. The total number of homologous gene pairs identified in parent A and B were 45216. While the parents A and D had 44920 homologous gene pairs. These homologous sequences were used as a reference for the already published second-generation Illumina sequencing H and L hybrid samples (Shahzad et al. 2020a; Shahzad et al. 2020b). Each sample was subjected to mutation detection and specific SNPs. The hybrids offspring samples may be homozygous or heterozygous at the specific SNP sites of their respective parents. Later, we compared the SNPs of two parents to select the homoeologue-specific SNPs bins using custom python scripts (Li et al. 2020). Further, genes with homoeologue-specific SNPs were compared on a case by case basis between the parents and then progeny to perform synchronous classification of homoeolog expression bias. The index to determine the significant difference was P-value <0.05. Compared with the diploid progenitors, the homologous gene pairs between both parents of hybrids were divided into 12 main categories with three subcategories, including the parental condition, no bias, and novel bias in progeny as already described (Yoo et al. 2013). Functional analyses including gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) was performed for all homoeolog expression genes that showed novel bias in different tissues of hybrids. GO terms and enriched pathways were retrieved from the cotton functional genomics database (https://cottonfgd.org/analyz) with P-value < 0.05.

Homoeolog expression changes in different tissues of high hybrid

To observe the pattern of homoeolog expression in different tissues of high hybrid, a case by case gene comparison was performed for high hybrid (H) and its parents (maternal=A, paternal=B). The results are detailed in Table 1. Notably, a pattern of homoeolog expression was conserved in this hybrid for the greater number of genes in all tissues under investigation. Fewer homoeolog gene copies reflected novel biased pattern. For example, 36853 total number homoeolog gene pairs were identified in each tissue. We observed that 96% in the root, 97% in leaf, and 97% in flower buds had maintained parental expression patterns in hybrid. The total expression level of a homologous gene pair had no change or equal to that of both parents. An average of 1% homoeolog gene pairs showed no bias in each hybrid tissue. Moreover, an average of 2.1% of genes reflected novel bias expression in each tissue of hybrid. Further analysis indicated that 1.2% homoeolog gene pairs in the root, 1.06% in leaf, and 1.1% in flower buds demonstrated overall maternal parent biased expression in hybrid. Whereas, overall homoeolog gene pairs with paternal parent biased expression in hybrid was 1.1%, 1.04%, and 0.95% in the root, leaf, and flower buds, respectively. The overall statistics of homoeolog expression changes in different tissues of high hybrid seem to show a highly balanced biased expression of homoeolog gene pairs in each tissue of high hybrid. These results suggested homoeolog gene pairs that displayed fewer novel patterns of homoeolog expression may be responsible for the high level of heterosis in upland cotton.

Homoeolog expression changes in different tissues of low hybrid

The results of homoeolog expression analysis in different tissues of low hybrid (L) and its reciprocal parents (A= maternal, B= paternal) are represented in Table 2. Interestingly, the homoeolog expression patterns in hybrid were in the same direction for the majority of homoeolog gene pairs as observed in parents. For instance, 97.3% in the root, 97.2% in leaf, and 97.6% in flower buds out of all homoeolog gene pairs (36853) reflected equal parental expression patterns in this hybrid. These result predicted that expression biases pre-existing in parents was simply conserved in these tissues of hybrid. The group of homoeolog gene pairs with bias expression in parents and reverted to no bias expression in hybrid had an average of 1% in root and leaf. On the other hand, only 0.60% of homologous gene pairs in flower buds had this reversion. According to the results of novel bias expression in hybrid, 1.75, 1.78, and 1.76% homologous gene pairs adopted novel bias in the root, leaf, and flower buds, respectively. Overall, a group of maternal or paternal biased homologous gene pairs had very similar trends in each tissue. There was 0.83% maternal and 0.84% paternal biased homologous gene expression in the root of hybrid, whereas leaf had 0.98 and 0.81% maternal and paternal biased expression, respectively. Flower buds also showed a similar percentage of maternal and parental biased expression. Altogether, these results suggest that the majority of the homoeolog gene exhibited balanced or non-differential homoeolog expression patterns in low hybrid. Further, homologous genes with few novel patterns of homoeolog expression probably contribute to the low level of heterosis in upland cotton.

Functional analysis of homoeolog genes exhibiting novel expression changes in hybrids

As to understand possible functions of homoeolog genes that had novel biased expression under hybridization in cotton, we performed GO and KEGG enrichment analysis in each hybrid. GO enrichment analysis (p-value < 0.05) for high hybrid showed a majority of novel biased expressed homoeolog genes had functional annotation related to catalytic activity, nucleus, metabolic process, transmembrane transport, binding of nucleotide, nucleic acid, protein, ATP, and RNA (Fig. 1a). Genes associated with binding of protein and ATP had the highest portion of novel biased expressed homoeolog genes in this hybrid. Most enriched GO terms for novel biased expressed homoeolog genes in low hybrid were involved in functional processes related to the nucleus, catalytic activity, binding of ATP, protein, and nucleotide (Fig. 1b). Interestingly, genes that performed the binding function of protein and ATP also had the majority of novel biased expressed homoeolog genes in this hybrid.

KEGG analysis revealed that the majority of novel biased expressed homoeolog genes in high hybrid were enriched in pathways such as biosynthesis of amino acids, mRNA surveillance pathway, Glycolysis / Gluconeogenesis, biosynthesis of antibiotics, metabolism of carbon, glutathione, cysteine/ methionine, biosynthesis of secondary metabolites, and microbial metabolism in diverse environments (Fig. 2a). Biosynthesis of secondary metabolites had more significant and gene enrichment as compared to other pathways. It was observed that a major portion of novel biased expressed homoeolog genes in low hybrid had enrichment in mRNA surveillance pathway, RNA degradation, Oocyte meiosis, metabolism of pyruvate, carbon, sugars, and biosynthesis of secondary metabolites (Fig. 2b). In a low hybrid, biosynthesis of secondary metabolites also had a maximum number of novel biased expressed homoeolog genes as compared to other pathways.

Cotton is the backbone of the textile industry because of its unique fiber characteristics. However, the area of cotton cultivation, as well as production, is decreasing in a recent decade all over the world. Hybridization has effectively resolved the problem of decreased crop productivity. For example, maize and rice cultivars developed with exploitation of heterosis have more harvest potential than conservative genotypes (Schnable et al. 2013). Cotton belongs to the genus Gossypium and upland cotton (G. hirsutum) is stated as new world cotton. The hybridization in upland cotton whether it is between species (interspecific) or within species (intraspecific) produced meaningful heterosis in yield-related traits (Li et al. 2019). It was observed that intraspecific crosses and improvement in fiber quality traits require more breeding efforts in cotton. How hybridization increased yield potential in upland cotton needs plenty of investigation at the genomic, transcriptomic, and epigenetic levels. Interestingly, the previous study has shown hybridization itself had sound effects on duplicate genes (homoeolog genes) in allopolyploids (Ren et al. 2016; Yoo et al. 2013). This study analyzed the effects of intraspecific hybridization on naturally occurring homoeolog genes in upland cotton through RNA-seq data of contrasting hybrids and full-length Pacbio sequences of their respective parents. In particular, gene expression data from three tissues (root, leaf, and flower buds) were used to explore changes in homoeolog expression and functional annotation of novel biased genes of hybrids.

Hybridization produced parent legacy of homoeolog expressions in upland cotton

Hybridization is a vital phenomenon to produce genetic improvement in crops. It is a common method in which two genetically different genotypes are crossed to obtain vigor progeny in desired traits. More than two sets of chromosomes in allotetraploid upland cotton caused genetic, gene expression, and epigenetic changes. These effects may result in more vigorous growth or heterosis (Chen 2007). Previous studies in maize (Yao et al. 2013) and alfalfa (Bingham et al. 1994) suggested that allopolyploids exhibit more heterosis as compared to autopolyploid and tetraploid display increased heterosis as compared with diploids. It was found that intraspecific crossing in upland cotton produced progeny that often displays heterosis in early biomass growth and seed cotton yield (Shahzad et al. 2020a, 2019b). However, there has been no study to analyze changes in homoeolog gene expression of inbred parents due to hybridization. The result of our study revealed that intraspecific hybridization in upland cotton produced highly balanced changes in the pattern of homoeolog expression. The pattern was conserved in hybrids for the greater number of genes. For example, both higher and lower hybrids had around 95% homoeolog gene pairs with similar expression pre-existing in parents in all tissues under investigation. Fewer homoeolog gene copies imitated novel homoeolog expression patterns in hybrids. Further, fewest homoeolog gene pairs exhibited maternal or paternal overall biases expression. A previous study in allopolyploid cotton showed that most homoeolog gene pairs reflected parental expression conditions in hybrid and synthetic allopolyploids (Yao et al. 2013). Similar results were reported in hybrid fish and its progenitors (Ren et al. 2016; Ren et al. 2017). Assessment of homoeolog expressions in Brassica napus described that majority of homoeolog gene pairs maintained parental expression patterns in synthetic hybrids and natural allotetraploids (Zhang et al. 2016). The effect of polyploidization on homoeolog expression in hexaploid wheat revealed equal homoeolog expression of parental lines is often conserved in progeny (Li et al. 2014). Altogether, these results indicated that the effect of hybridization or genome doubling on the expression of homoeolog gene pairs in derived hybrids is often inherited from diploid parents. Though the relationship between heterosis and homoeolog gene pairs is still not fully understood, these findings suggest that fewer changes at the transcriptomic level could be associated with heterosis in upland cotton.

Novel homoeolog gene expressions may be regulated heterosis in cotton through secondary metabolism

As already described in the majority of the result of homoeolog gene pairs had a parental legacy in hybrids. Theoretically, these gene pairs do not contribute to the superior or inferior performance of cotton hybrids. Homoeolog gene pairs with a novel or differential expression in hybrids could regulate heterosis. Functional analysis of these novel homoeolog gene pairs revealed that the binding function of protein and ATP through the secondary metabolism pathway plays the role to increase growth potential in offspring relative to their parents. Metabolites are an essential component for adaptation, stress management, and regulation of different biological processes in plants (Kleessen et al. 2014; Miller et al. 2015). Modification of secondary metabolism is useful to growth heterosis in oilseed crop Brassica juncea (Bajpai et al. 2019). Further, these researchers reported an additive mode of inheritance of metabolites in hybrids. Previous metabolic investigation in hybrids of maize (Wang et al. 2014), rice (Ma et al. 2011), and Arabidopsis (Korn et al. 2010) provided preliminary perceptions and enhanced knowledge about the molecular mechanism of plant heterosis. There is plenty of emphasis on genetic, gene expression, and epigenetic related studies to explore the biological basis of heterosis in cotton. Yet, the relationship between heterosis and metabolic changes in hybrids has been poorly known in cotton. Our results underpin the importance of genomic inheritance of secondary metabolites between hybrids and inbred parents in cotton.

In this study, global homoeolog gene expression patterns between intraspecific contrasting cotton hybrids and their reciprocal inbred parents were analyzed using RNA sequencing. The results indicated that hybridization caused no expression biases of homoeolog gene pairs in hybrids. Hybrids showed no bias expression preference towards any parents. Further, functional annotation of novel expressed homoeolog gene pairs explained changed secondary metabolic gene expression between parents and hybrids may regulate heterosis in upland cotton. Together, our results provide a general understanding of homoeolog gene expression patterns and may help to explore the genetic basis of heterosis in cotton.

Acknowledgments

The authors thank Yasir Muhammad for his valuable comments on this paper.

Authors’ contributions

Xing CZ and Wu JY designed this research. Gou LP, Qi TX, Tang HN, Qiao XQ, and Wang HL prepared plant population. Zhang XX, Zhang M, and Zhang BB carried out RNA extraction and performed data analysis. Shahzad K, Zhang XX, Xing CZ, and Wu JY drafted the manuscript. All authors read and approved the final manuscript.

Funding

This work was funded by National Key Research and Development Program of China (2016YFD0101400). The funding organization was not involved in the design of the study, analysis, and writing the manuscript.

Availability of data and materials

All data generated or analyzed in this study included in published article

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Adams KL, Cronn R, Percifield R, et al. Genes duplicated by polyploidy show unequal contributions to the transcriptome and organ-specific reciprocal silencing. Proc Natl Acad Sci. 2003;100(8):4649-54. https://doi.org/10.1073/pnas.0630618100.

Akhunova AR, Matniyazov RT, Liang H, et al. Homoeolog-specific transcriptional bias in allopolyploid wheat. BMC genomics. 2010;11(1):505. https://doi.org/10.1186/1471-2164-11-505.

Ali I, Shakeel A, Saeed A, et al. Combining ability analysis and heterotic studies for within-boll yield components and fibre quality in cotton. J Anim Plant Sci. 2016;26(1):156-62.

Anssour S, Krügel T, Sharbel T, et al. Phenotypic, genetic and genomic consequences of natural and synthetic polyploidization of Nicotiana attenuata and Nicotiana obtusifolia. Ann Bot. 2009;103(8):1207-17. https://doi.org/10.1093/aob/mcp058.

Auger B, Baron C, Lucas M-O, et al. Brassica orthologs from BANYULS belong to a small multigene family, which is involved in procyanidin accumulation in the seed. Planta. 2009;230(6):1167. https://doi.org/10.1007/s00425-009-1017-0.

Bajpai PK, Reichelt M, Augustine R, et al. Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea. Heredity. 2019;123(3):318-36. https://doi.org/10.1038/s41437-019-0213-3.

Bingham E, Groose R, Woodfield D, et al. Complementary gene interactions in alfalfa are greater in autotetraploids than diploids. Crop Sci. 1994;34(4):823-9. https://doi.org/10.2135/cropsci1994.0011183X003400040001x.

Birchler JA, Auger DL, Riddle NC. In search of the molecular basis of heterosis. Plant Cell. 2003;15(10):2236-9. https://doi.org/10.1105/tpc.151030.

Bottley A, Xia G, Koebner R. Homoeologous gene silencing in hexaploid wheat. Plant J. 2006;47(6):897-906. https://doi.org/10.1111/j.1365-313x.2006.02841.x.

Buggs RJ, Elliott NM, Zhang L, et al. Tissue‐specific silencing of homoeologs in natural populations of the recent allopolyploid Tragopogon mirus. New Phytol. 2010;186(1):175-83. https://doi.org/10.1111/j.1469-8137.2010.03205.x.

Chang PL, Dilkes BP, McMahon M, et al. Homoeolog-specific retention and use in allotetraploid Arabidopsis suecica depends on parent of origin and network partners. Genome biol. 2010;11(12):R125. https://doi.org/10.1186/gb-2010-11-12-r125.

Chaudhary B, Flagel L, Stupar RM, et al. Reciprocal silencing, transcriptional bias and functional divergence of homeologs in polyploid cotton (Gossypium). Genetics. 2009;182(2):503-17. https://doi.org/10.1534/genetics.109.102608.

Chelaifa H, Monnier A, Ainouche M. Transcriptomic changes following recent natural hybridization and allopolyploidy in the salt marsh species Spartina× townsendii and Spartina anglica (Poaceae). New Phytol. 2010;186(1):161-74. https://doi.org/10.1111/j.1469-8137.2010.03179.x.

Chen ZJ. Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu Rev Plant Biol. 2007;58:377-406. https://doi.org/10.1146/annurev.arplant.58.032806.103835.

De Smet R, Adams KL, Vandepoele K, et al. Convergent gene loss following gene and genome duplications creates single-copy families in flowering plants. Proc Natl Acad Sci. 2013;110(8):2898-903. https://doi.org/10.1073/pnas.1300127110.

Doyle JJ, Flagel LE, Paterson AH, et al. Evolutionary genetics of genome merger and doubling in plants. Annu Rev Genet. 2008;42:443-61. https://doi.org/10.1146/annurev.genet.42.110807.091524.

Dufresne F, Hebert PD. Hybridization and origins of polyploidy. Proc R Soc Lond. Series B: Biological Sciences. 1994;258(1352):141-6. http://dx.doi.org/10.1098/rspb.1994.0154.

Flagel L, Udall J, Nettleton D, et al. Duplicate gene expression in allopolyploid Gossypiumreveals two temporally distinct phases of expression evolution. BMC Biol. 2008;6(1):16. https://doi.org/10.1186/1741-7007-6-16.

Flagel LE, Wendel JF. Evolutionary rate variation, genomic dominance and duplicate gene expression evolution during allotetraploid cotton speciation. New Phytol. 2010;186(1):184-93. https://doi.org/10.1111/j.1469-8137.2009.03107.x.

Gaeta RT, Pires JC, Iniguez-Luy F, et al. Genomic changes in resynthesized Brassica napus and their effect on gene expression and phenotype. Plant Cell. 2007;19(11):3403-17. https://doi.org/10.1105/tpc.107.054346.

Gaeta RT, Yoo S-Y, Pires J, et al. Analysis of gene expression in resynthesized Brassica napus allopolyploids using Arabidopsis 70mer oligo microarrays. PloS one. 2009;4(3):e4760. https://doi.org/10.1371/journal.pone.0004760.

Grover C, Gallagher J, Szadkowski E, et al. Homoeolog expression bias and expression level dominance in allopolyploids. New Phytol. 2012;196(4):966-71. https://doi.org/10.1111/j.1469-8137.2012.04365.x.

Han F, Fedak G, Guo W, et al. Rapid and repeatable elimination of a parental genome-specific DNA repeat (pGc1R-1a) in newly synthesized wheat allopolyploids. Genetics. 2005;170(3):1239-45. https://dx.doi.org/ 10.1534/genetics.104.039263.

Hovav R, Udall JA, Chaudhary B, et al. Partitioned expression of duplicated genes during development and evolution of a single cell in a polyploid plant. Proc Natl Acad Sci. 2008;105(16):6191-5. https://doi.org/10.1073/pnas.0711569105.

Hughes TE, Langdale JA, Kelly S. The impact of widespread regulatory neofunctionalization on homeolog gene evolution following whole-genome duplication in maize. Genome Res. 2014;24(8):1348-55. https://dx.doi.org/ 10.1101/gr.172684.

Kashkush K, Feldman M, Levy AA. Transcriptional activation of retrotransposons alters the expression of adjacent genes in wheat. Nat Genet. 2003;33(1):102-6. https://doi.org/10.1038/ng1063.

Kleessen S, Laitinen R, Fusari CM, et al. Metabolic efficiency underpins performance trade-offs in growth of Arabidopsis thaliana. Nat Commun. 2014;5(1):1-10. https://doi.org/10.1038/ncomms4537.

Kleinjan DA, Bancewicz RM, Gautier P, et al. Subfunctionalization of duplicated zebrafish pax6 genes by cis-regulatory divergence. PLoS Genet. 2008;4(2):e29. https://doi.org/10.1371/journal.pgen.0040029.

Korn M, Gärtner T, Erban A, et al. Predicting Arabidopsis freezing tolerance and heterosis in freezing tolerance from metabolite composition. Mol Plant. 2010;3(1):224-35. https://doi.org/10.1093/mp/ssp105.

Kovarik A, Dadejova M, Lim YK, et al. Evolution of rDNA in Nicotiana allopolyploids: a potential link between rDNA homogenization and epigenetics. Ann Bot. 2008;101(6):815-23. https://doi.org/10.1093/aob/mcn019.

Kraitshtein Z, Yaakov B, Khasdan V, et al. Genetic and epigenetic dynamics of a retrotransposon after allopolyploidization of wheat. Genetics. 2010;186(3):801-12. https://doi.org/10.1534/genetics.110.120790.

Li A, Liu D, Wu J, et al. mRNA and small RNA transcriptomes reveal insights into dynamic homoeolog regulation of allopolyploid heterosis in nascent hexaploid wheat. Plant Cell. 2014;26(5):1878-900. https://doi.org/10.1105/tpc.114.124388.

Li M, Wang R, Wu X, et al. Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus. BMC genomics. 2020;21:1-15. https://doi.org/10.1186/s12864-020-6747-1.

Li X, Shahzad K, Guo L, et al. Using yield quantitative trait locus targeted SSR markers to study the relationship between genetic distance and yield heterosis in upland cotton (Gossypium hirsutum). Plant Breed. 2019;138(1):105-13. https://doi.org/10.1111/pbr.12668.

Liu S-L, Adams KL. Dramatic change in function and expression pattern of a gene duplicated by polyploidy created a paternal effect gene in the Brassicaceae. Mol Biol Evol. 2010;27(12):2817-28. https://doi.org/10.1093/molbev/msq169.

Ma Q, Hedden P, Zhang Q. Heterosis in rice seedlings: its relationship to gibberellin content and expression of gibberellin metabolism and signaling genes. Plant Physiol. 2011;156(4):1905-20. https://doi.org/10.1104/pp.111.178046.

Madlung A, Masuelli RW, Watson B, et al. Remodeling of DNA methylation and phenotypic and transcriptional changes in synthetic Arabidopsis allotetraploids. Plant Physiol. 2002;129(2):733-46. https://doi.org/10.1104/pp.003095.

McGrath CL, Gout J-F, Johri P, et al. Differential retention and divergent resolution of duplicate genes following whole-genome duplication. Genome Res. 2014;24(10):1665-75. 10.1101/gr.173740.114.

Miller M, Song Q, Shi X, et al. Natural variation in timing of stress-responsive gene expression predicts heterosis in intraspecific hybrids of Arabidopsis. Nat Commun. 2015;6(1):1-13. https://doi.org/10.1038/ncomms8453.

Mochida K, Yamazaki Y, Ogihara Y. Discrimination of homoeologous gene expression in hexaploid wheat by SNP analysis of contigs grouped from a large number of expressed sequence tags. Mol Genet and Genomics. 2004;270(5):371-7. https://doi.org/10.1007/s00438-003-0939-7.

Ren L, Cui J, Wang J, et al. Analyzing homoeolog expression provides insights into the rediploidization event in gynogenetic hybrids of Carassius auratus red var.× Cyprinus carpio. Sci Rep. 2017;7(1):1-11. https://doi.org/10.1038/s41598-017-14084-7.

Ren L, Li W, Tao M, et al. Homoeologue expression insights into the basis of growth heterosis at the intersection of ploidy and hybridity in Cyprinidae. Sci Rep. 2016;6(1):1-12. https://doi.org/10.1038/srep27040.

Salmon A, Flagel L, Ying B, et al. Homoeologous nonreciprocal recombination in polyploid cotton. New Phytol. 2010;186(1):123-34. 10.1111/j.1469-8137.2009.03093.x.

Samuel Yang S, Cheung F, Lee JJ, et al. Accumulation of genome‐specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton. Plant J. 2006;47(5):761-75. 10.1111/j.1365-313X.2006.02829.x.

Schnable JC, Freeling M. Genes identified by visible mutant phenotypes show increased bias toward one of two subgenomes of maize. PloS one. 2011;6(3):e17855. https://doi.org/10.1371/journal.pone.0017855.

Schnable PS, Springer NM. Progress toward understanding heterosis in crop plants. Annu Rev Plant Biol. 2013;64. 10.1146/annurev-arplant-042110-103827.

Shahzad K, Qi T, Guo L, et al. Adaptability and stability comparisons of inbred and hybrid cotton in yield and fiber quality traits. Agronomy. 2019b;9(9):516. https://doi.org/10.3390/agronomy9090516.

Shahzad K, Xue L, Tingxiang Q, et al. Genetic analysis of yield and fiber quality traits in upland cotton (Gossypium hirsutum L.) cultivated in different ecological regions of China. J Cotton Res. 2019a;2(1):14. https://doi.org/10.1186/s42397-019-0031-4.

Shahzad K, Zhang X, Guo L, et al. Comparative transcriptome analysis of inbred lines and contrasting hybrids reveals overdominance mediate early biomass vigor in hybrid cotton. BMC genomics. 2020a;21(1):1-16. https://doi.org/10.1186/s12864-020-6561-9.

Shahzad K, Zhang X, Guo L, et al. Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton. BMC Plant Biol. 2020b;20:1-18. https://doi.org/10.1186/s12870-020-02442-z.

Szadkowski E, Eber F, Huteau V, et al. The first meiosis of resynthesized Brassica napus, a genome blender. New Phytol. 2010;186(1):102-12. 10.1111/j.1469-8137.2010.03182.x.

Wang J, Tian L, Lee H-S, et al. Genomewide nonadditive gene regulation in Arabidopsis allotetraploids. Genetics. 2006;172(1):507-17. 10.1534/genetics.105.047894.

Wang J, Tian L, Madlung A, et al. Stochastic and epigenetic changes of gene expression in Arabidopsis polyploids. Genetics. 2004;167(4):1961-73. 10.1534/genetics.104.027896

Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57-63. https://doi.org/10.1038/nrg2484.

Wang Z, Xue Z, Wang T. Differential analysis of proteomes and metabolomes reveals additively balanced networking for metabolism in maize heterosis. J Proteome Res. 2014;13(9):3987-4001. https://doi.org/10.1021/pr500337j.

Wood TE, Takebayashi N, Barker MS, et al. The frequency of polyploid speciation in vascular plants. Proc Natl Acad Sci. 2009;106(33):13875-9. https://doi.org/10.1073/pnas.0811575106.

Wu J, Lin L, Xu M, et al. Homoeolog expression bias and expression level dominance in resynthesized allopolyploid Brassica napus. BMC genomics. 2018;19(1):1-13. https://doi.org/10.1186/s12864-018-4966-5.

Xiao J, Li J, Yuan L, et al. Dominance is the major genetic basis of heterosis in rice as revealed by QTL analysis using molecular markers. Genetics. 1995;140(2):745-54.

Yao H, Gray AD, Auger DL, et al. Genomic dosage effects on heterosis in triploid maize. Proc Natl Acad Sci. 2013;110(7):2665-9. https://doi.org/10.1073/pnas.1221966110.

Yoo M, Szadkowski E, Wendel J. Homoeolog expression bias and expression level dominance in allopolyploid cotton. Heredity. 2013;110(2):171. https://doi.org/10.1038/hdy.2012.94.

Yu S, Li J, Xu C, et al. Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci. 1997;94(17):9226-31. https://doi.org/10.1073/pnas.94.17.9226.

Zhang D, Pan Q, Tan C, et al. Genome-wide gene expressions respond differently to A-subgenome origins in Brassica napus synthetic hybrids and natural allotetraploid. Front Plant Sci. 2016;7:1508. https://doi.org/10.3389/fpls.2016.01508.

Table 1 Homoeolog expression bias in different tissues of the high F₁upland cotton hybrid

Expression in parents	Expression in progeny	Root (%)	Leaf (%)	Flower buds (%)
A=B	A=B Parental condition	35471 (96.2)	35734 (97.1)	35779 (97.1)
A<B	A<B Parental condition	4 (0.01)	0 (0.0)	1 (0.0)
A>B	A>B Parental condition	5 (0.01)	7 (0.02)	4 (0.01)
A>B	A=B No bias in progeny	190 (0.52)	279 (0.76)	210 (0.57)
A<B	A=B No bias in progeny	345 (0.94)	66 (0.18)	107 (0.29)
A=B	A>B Novel bias in progeny	416 (1.13)	375 (1.02)	391 (1.06)
A<B	A>B Novel bias in progeny	20 (0.05)	10 (0.03)	12 (0.03)
A>B	A<B Novel bias in progeny	10 (0.03)	18 (0.05)	10 (0.03)
A=B	A<B Novel bias in progeny	392 (1.06)	364 (1.0)	339 (0.92)
Overall A-biased in progeny		441 (1.20)	392 (1.06)	407 (1.10)
Overall B-biased in progeny		406 (1.10)	382 (1.04)	350 (0.95)
Total number of genes		36853	36853	36853

Here A= Maternal parent and B= Paternal parent

Table 2 Homoeolog expression bias in different tissues of the low F₁ upland cotton hybrid

Expression in parents	Expression in progeny	Root (%)	Leaf (%)	Flower buds (%)
A=D	A=D Parental condition	35844 (97.3)	35827 (97.2)	35981 (97.6)
A<D	A<D Parental condition	0 (0.0)	0 (0.0)	1 (0.0)
A>D	A>D Parental condition	2 (0.01)	1 (0.0)	1 (0.0)
A>D	A=D No bias in progeny	202 (0.55)	233 (0.63)	107 (0.29)
A<D	A=D No bias in progeny	197 (0.53)	135 (0.37)	116 (0.31)
A=D	A>D Novel bias in progeny	290 (0.89)	352 (0.96)	294 (0.80)
A<D	A>D Novel bias in progeny	12 (0.03)	8 (0.02)	14 (0.04)
A>D	A<D Novel bias in progeny	6 (0.02)	8 (0.02)	8 (0.02)
A=D	A<D Novel bias in progeny	300 (0.81)	288 (0.78)	331 (0.90)
Overall A-biased in progeny		304 (0.83)	361 (0.98)	309 (0.84)
Overall D-biased in progeny		306 (0.84)	297 (0.81)	340 (0.92)
Total number of genes		36853	36853	36853

Here A= Maternal parent and D= Paternal parent

Download PDF

Version 1

posted

You are reading this latest preprint version

Global insights into homoeolog gene expression in upland cotton under intraspecific hybridization

Status:

Version 1

Abstract

Figures

Background

Material And Methods

Results

Discussion

Conclusion

Declarations

References

Tables

Status:

Version 1