Differential gene expression and alternative splicing between diploid and tetraploid watermelon

Friends without benefits: Extensive cytotype sympatry and polyploid persistence in an African geophyte

PREMISE

Polyploidy is a major factor in plant adaptation and speciation. Multiple mechanisms contribute to autopolyploid frequency within populations, but uncertainties remain regarding mechanisms that facilitate polyploid establishment and persistence. Here we aimed to document and predict cytotype distributions of Oxalis obliquifolia Steud. ex A. Rich. across Gauteng, South Africa, and test for evidence of possible mechanisms, including morphological, phenological, and reproductive traits, that may potentially facilitate polyploid persistence.

METHODS

Over 320 O. obliquifolia plants from 25 sites were cytotyped using flow cytometry, and DNA ploidy was confirmed using meiotic chromosome squashes. Cytotypes were mapped and correlations with abiotic variables assessed using ordinations. To assess morphological and phenological associations with cytotype, we grew multiple cytotypes in a common garden, measured phenotypic traits and compared them using linear models and discriminant analyses. Intercytotype reproductive isolation was assessed using crossing experiments, and AMOVAs based on ITS DNA sequences tested for cytogeographic structure.

RESULTS

Six cytotypes were identified, and most sites had multiple cytotypes. Abiotic variables were not predictive of cytotype distribution. A clear gigas effect was present. Differences in flower size and phenology suggested pollinator interactions could play a role in polyploid persistence. Intercytotype crosses produced seed at low frequency. DNA data suggested diploids and polyploids were largely reproductively isolated in situ, and polyploidization events were not frequent enough to explain high cytotype sympatry.

CONCLUSIONS

Diploids and polyploids are behaving as separate species, despite little observable niche differentiation and non-zero potential intercytotype seed set. Tests on biotic interactions and intercytotype F1 fitness may provide insights into diploid and polyploid coexistence.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Whole genome doubling-induced the enrichment of H3K27me3 in genes carrying specific TEs in Aegilops tauschii

Polyploidization plays important roles in the evolution and breeding of the common wheat. Aegilops tauschii, the D-genome progenitor of the common wheat, provides a valuable pool of resistance genes to multiple diseases. Extensive studies focus on the exploration of these genes for wheat improvement. However, few studies have unveiled alternations on genome-wide expression pattern and histone modifications induced by whole-genome doubling (WGD) process. In this study, we conducted transcriptome analysis for the diploid and tetraploid Ae. taushcii lines using the leaf and root tissues. Both lines tend to display similar tissue-specific pattern. Interestingly, we found that TEs located in genic regions were depleted of the repressive histone mark H3K27me3, whereas their adjacent chromatin was enriched with H3K27me3. The tetraploid line exhibited higher levels of H3K27me3 in those regions than the diploid line, particularly for genic regions associated with TEs of the long interspersed nuclear elements (LINEs), CACTA, PIF/Harbinger, Tc1/Mariner and unclassed DNA transposon. Surprisingly, the expression levels of these TEs cognate genes were negatively associated with the levels of H3K27me3 between the tetraploid and diploid lines, suggesting the five types of TEs located within genic regions might be involved in the regulation of the ploidy-related gene expression, possibly through differential enrichment of H3K27me3 in the genic regions. These findings will help to understand the potential role of specific types of TEs on transcription in response to WGD.

Current insights into posttranscriptional regulation of fleshy fruit ripening

Fruit ripening is a complicated process that is accompanied by the formation of fruit quality. It is not only regulated at the transcriptional level via transcription factors or DNA methylation but also fine-tuned after transcription occurs. Here, we review recent advances in our understanding of key regulatory mechanisms of fleshy fruit ripening after transcription. We mainly highlight the typical mechanisms by which fruit ripening is controlled, namely, alternative splicing, mRNA N6-methyladenosine RNA modification methylation, and noncoding RNAs at the posttranscriptional level; regulation of translation efficiency and upstream open reading frame-mediated translational repression at the translational level; and histone modifications, protein phosphorylation, and protein ubiquitination at the posttranslational level. Taken together, these posttranscriptional regulatory mechanisms, along with transcriptional regulation, constitute the molecular framework of fruit ripening. We also critically discuss the potential usage of some mechanisms to improve fruit traits.

Development of homozygous tetraploid potato and whole genome doubling-induced the enrichment of H3K27ac and potentially enhanced resistance to cold-induced sweetening in tubers

Polyploid plants typically display advantages on some agronomically important traits over their diploid counterparts. Extensive studies have shown genetic, transcriptomic, and epigenetic dynamics upon polyploidization in multiple plant species. However, few studies have unveiled those alternations imposed only by ploidy level, without any interference from heterozygosity. Cultivated potato is highly heterozygous. Thus, in this study, we developed two homozygous autotetraploid lines and one homozygous diploid line in parallel from a homozygous diploid potato. We confirmed their ploidy levels using chloroplast counting and karyotyping. Oligo-FISH and genome re-sequencing validated that these potato lines are nearly homozygous. We investigated variations in phenotypes, transcription, and histone modifications between two ploidies. Both autotetraploid lines produced larger but fewer tubers than the diploid line. Interestingly, each autotetraploid line displayed ploidy-related differential expression for various genes. We also discovered a genome-wide enrichment of H3K27ac in genic regions upon whole-genome doubling (WGD). However, such enrichment was not associated with the differential gene expression between two ploidies. The tetraploid lines may exhibit better resistance to cold-induced sweetening (CIS) than the diploid line in tubers, potentially regulated through the expression of CIS-related key genes, which seems to be associated with the levels of H3K4me3 in cold-stored tubers. These findings will help to understand the impacts of autotetraploidization on dynamics of phenotypes, transcription, and histone modifications, as well as on CIS-related genes in response to cold storage.

A Comparative Study of Morphology, Photosynthetic Physiology, and Proteome between Diploid and Tetraploid Watermelon (Citrullus lanatus L.)

Watermelon is an important fruit that is widely distributed around the world. In particular, the production and consumption of watermelon in China ranks first in the world. Watermelon production is severely affected by a variety of biotic and abiotic stresses during cultivation, and polyploidization can promote stress resistance and yield. However, the morphological and physiological characteristics of tetraploid watermelon and the underlying molecular mechanisms are still poorly understood. In this study, we revealed that the leaves, fruits, and seeds of tetraploid watermelon were significantly larger than those of the diploid genotype. Some physiological characteristics, including photosynthetic rate (Pn) and stomatal conductance (Gs), were greater, whereas the intercellular CO₂ concentration (Ci) and transpiration rate (Tr) were lower in tetraploid than in diploid watermelon. Two-dimensional gel electrophoresis combined with tandem mass spectrometry (MALDI-TOF/TOF) was performed to compare proteomic changes between tetraploid and diploid watermelon. A total of 21 differentially expressed proteins were identified; excluding the identical proteins, 8 proteins remained. Among them, four proteins were upregulated and four were downregulated in tetraploid versus diploid genotypes. qRT-PCR results showed inconsistencies in gene expression and protein accumulation, indicating a low correlation between gene expression and protein abundance. Generally, this study extends our understanding of the traits and molecular mechanisms of tetraploid watermelon and provides a theoretical basis for watermelon polyploid breeding.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Genetic and Physiological Responses to Heat Stress in Brassica napus

Given the current rise in global temperatures, heat stress has become a major abiotic challenge affecting the growth and development of various crops and reducing their productivity. Brassica napus, the second largest source of vegetable oil worldwide, experiences a drastic reduction in seed yield and quality in response to heat. This review outlines the latest research that explores the genetic and physiological impact of heat stress on different developmental stages of B. napus with a special attention to the reproductive stages of floral progression, organogenesis, and post flowering. Several studies have shown that extreme temperature fluctuations during these crucial periods have detrimental effects on the plant and often leading to impaired growth and reduced seed production. The underlying mechanisms of heat stress adaptations and associated key regulatory genes are discussed. Furthermore, an overview and the implications of the polyploidy nature of B. napus and the regulatory role of alternative splicing in forming a priming-induced heat-stress memory are presented. New insights into the dynamics of epigenetic modifications during heat stress are discussed. Interestingly, while such studies are scarce in B. napus, opposite trends in expression of key genetic and epigenetic components have been identified in different species and in cultivars within the same species under various abiotic stresses, suggesting a complex role of these genes and their regulation in heat stress tolerance mechanisms. Additionally, omics-based studies are discussed with emphasis on the transcriptome, proteome and metabolome of B. napus, to gain a systems level understanding of how heat stress alters its yield and quality traits. The combination of omics approaches has revealed crucial interactions and regulatory networks taking part in the complex machinery of heat stress tolerance. We identify key knowledge gaps regarding the impact of heat stress on B. napus during its yield determining reproductive stages, where in-depth analysis of this subject is still needed. A deeper knowledge of heat stress response components and mechanisms in tissue specific models would serve as a stepping-stone to gaining insights into the regulation of thermotolerance that takes place in this important crop species and support future breeding of heat tolerant crops.

Mapping and Analysis of a Novel Genic Male Sterility Gene in Watermelon (Citrullus lanatus)

Seed production is critical for watermelon production, which mostly involves first-generation hybrid varieties. However, watermelon hybrid seed production currently requires complex procedures, including artificial isolation and pollination. Therefore, the development and use of a male-sterile system to generate watermelon hybrids can simplify the process. The scarcity of male-sterile watermelon germplasm resources necessitates the use of molecular breeding methods. Unfortunately, the genes responsible for male sterility in watermelon have not been cloned. Thus, the genetic basis of the male sterility remains unknown. In this study, two DNA pools derived from male-sterile and normal plants in the F2 population were used for whole-genome resequencing. The Illumina high-throughput sequencing resulted in 62.99 Gbp clean reads, with a Q30 of 80% after filtering. On the basis of the SNP index association algorithm, eight candidate regions (0.32 Mb) related to specific traits were detected on chromosome 6. Expression pattern analyses and watermelon transformation studies generated preliminary evidence that Cla006625 encodes a pollen-specific leucine-rich repeat protein (ClaPEX1) influencing the male sterility of watermelon. The identification and use of genic male sterility genes will promote watermelon male sterility research and lay the foundation for the efficient application of seed production technology.

Comparative transcriptome analysis reveals the differential response to cadmium stress of two Pleurotus fungi: Pleurotus cornucopiae and Pleurotus ostreatus

Pleurotus has great potential for heavy metal mycoremediation. Using comparative transcriptome analysis, the response of Pleurotus ostreatus and Pleurotus cornucopiae under Cd contamination was evaluated. P. ostreatus and P. cornucopia accumulated 0.34 and 0.46 mg/g Cd in mycelium, respectively. Cd removal elevated with its concentration elevation, which reached 56.47% and 54.60% for P. ostreatus and P. cornucopia with Cd at 20 mg/L. Low-level Cd (≤ 1 mg/L) had no significant influence on either fungus, while varied response was observed under high-level Cd. 705 differentially expressed genes (DEGs) were identified in P. cornucopia at Cd1 and Cd20, whereas 12,551 DEGs in P. ostreatus. Differentially regulated functional categories and pathways were also identified. ATP-binding cassette transporters were involved in Cd transport in P. cornucopia, whereas the endocytosis and phagosome pathways were more enhanced in P. ostreatus. 26 enzymes including peroxisomal enzymes catalase and superoxide dismutase were upregulated in P. ostreatus, whereas only cytosolic catalase was overexpressed in P. cornucopia, suggesting their different Cd detoxification pathways. Also, the mitogen-activated protein kinase signaling pathway involved in Cd resistance in both species instead of glutathione metabolism, although more active in P. ostreatus. These findings provided new insight into the molecular mechanism of mycoremediation and accumulator screening.

The Transcriptional and Splicing Changes Caused by Hybridization Can Be Globally Recovered by Genome Doubling during Allopolyploidization

Polyploidization is a major driving force in plant evolution. Allopolyploidization, involving hybridization and genome doubling, can cause extensive transcriptome reprogramming which confers allopolyploids higher evolutionary potential than their diploid progenitors. To date, little is known about the interplay between hybridization and genome doubling in transcriptome reprogramming. Here, we performed genome-wide analyses of transcriptome reprogramming during allopolyploidization in wheat and brassica lineages. Our results indicated that hybridization-induced transcriptional and splicing changes of genes can be largely recovered to parental levels by genome doubling in allopolyploids. As transcriptome reprogramming is an important contributor to heterosis, our finding updates a longstanding theory that heterosis in interspecific hybrids can be permanently fixed through genome doubling. Our results also indicated that much of the transcriptome reprogramming in interspecific hybrids was not caused by the merging of two parental genomes, providing novel insights into the mechanisms underlying both heterosis and hybrid speciation.

Altered chromatin conformation and transcriptional regulation in watermelon following genome doubling

Polyploidy has played a crucial role in plant evolution, development and function. Synthetic autopolyploid represents an ideal system to investigate the effects of polyploidization on transcriptional regulation. In this study, we deciphered the impact of genome duplication at phenotypic and molecular levels in watermelon. Overall, 88% of the genes in tetraploid watermelon followed a >1:1 dosage effect, and accordingly, differentially expressed genes were largely upregulated. In addition, a great number of hypomethylated regions (1688) were identified in an isogenic tetraploid watermelon. These differentially methylated regions were localized in promoters and intergenic regions and near transcriptional start sites of the identified upregulated genes, which enhances the importance of methylation in gene regulation. These changes were reflected in sophisticated higher-order chromatin structures. The genome doubling caused switching of 108 A and 626 B compartments that harbored genes associated with growth, development and stress responses.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps

BACKGROUND

Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown.

RESULTS

Two high-density genetic maps were constructed by whole-genome resequencing with two F₂ segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci.

CONCLUSION

The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Changes in Alternative Splicing in Response to Domestication and Polyploidization in Wheat

Alternative splicing (AS) occurs extensively in eukaryotes as an important mechanism for regulating transcriptome complexity and proteome diversity, but variation in the AS landscape in response to domestication and polyploidization in crops is unclear. Hexaploid wheat (AABBDD, Triticum aestivum) has undergone two separate allopolyploidization events, providing an ideal model for studying AS changes during domestication and polyploidization events. In this study, we performed high-throughput transcriptome sequencing of roots and leaves from wheat species with varied ploidies, including wild diploids (A^bA^b, Triticum boeoticum) and tetraploids (AABB, Triticum dicoccoides), domesticated diploids (A^mA^m, Triticum monococcum) and tetraploids (AABB, Triticum dicoccum), hexaploid wheat (AABBDD, T aestivum), as well as newly synthesized hexaploids together with their parents. Approximately 22.1% of genes exhibited AS, with the major AS type being intron retention. The number of AS events decreased after domestication in both diploids and tetraploids. Moreover, the frequency of AS occurrence tended to decrease after polyploidization, consistent with the functional sharing model that proposes AS and duplicated genes are complementary in regulating transcriptome plasticity in polyploid crops. In addition, the subgenomes exhibited biased AS responses to polyploidization, and ∼87.1% of homeologs showed AS partitioning in hexaploid wheat. Interestingly, substitution of the D-subgenome modified 42.8% of AS patterns of the A- and B-subgenomes, indicating subgenome interplay reprograms AS profiles at a genome-wide level, although the causal-consequence relationship requires further study. Conclusively, our study shows that AS variation occurs extensively after polyploidization and domestication in wheat species.

Molecular regulatory mechanisms underlying the adaptability of polyploid plants

Polyploidization influences the genetic composition and gene expression of an organism. This multi-level genetic change allows the formation of new regulatory pathways leading to increased adaptability. Although both forms of polyploidization provide advantages, autopolyploids were long thought to have little impact on plant divergence compared to allopolyploids due to their formation through genome duplication only, rather than in combination with hybridization. Recent advances have begun to clarify the molecular regulatory mechanisms such as microRNAs, alternative splicing, RNA-binding proteins, histone modifications, chromatin remodelling, DNA methylation, and N⁶ -methyladenosine (m6A) RNA methylation underlying the evolutionary success of polyploids. Such research is expanding our understanding of the evolutionary adaptability of polyploids and the regulatory pathways that allow adaptive plasticity in a variety of plant species. Herein we review the roles of individual molecular regulatory mechanisms and their potential synergistic pathways underlying plant evolution and adaptation. Notably, increasing interest in m6A methylation has provided a new component in potential mechanistic coordination that is still predominantly unexplored. Future research should attempt to identify and functionally characterize the evolutionary impact of both individual and synergistic pathways in polyploid plant species.

Transcriptome sequencing and screening of genes related to sex determination of Trichosanthes kirilowii Maxim

Trichosanthes kirilowii Maxim. (TK) is a dioecious plant in the Cucurbitaceae for which different sexes have separate medicinal uses. In order to study the genes related to sex determination, transcriptome sequencing was performed on flower buds of male and female plants using the high-throughput sequencing technology. A total of 145,975 unigenes and 7110 DEGs were obtained. There were 6776 DEGs annotated to 1234 GO terms and enriched to 18 functional groups, including five biological processes related to sugar metabolism. KEGG pathway analysis indicated genes involved in hormone transduction, hormone synthesis and carbohydrate metabolism. Many DEGs of TK are involved in reproductive organ formation, hormone signal transduction and regulatory networks. Combining the results of GO, KEGG and qRT-PCR, 11 sex determining candidate genes of TK were selected, including MYB80, MYB108, CER1, CBL9, ABCB19, SERK1, HSP81-3, ACS9, SEP3, AUX1 and YUC6. The results provide a foundation for the study of sex differentiation in TK.

Sulfadiazine biodegradation by Phanerochaete chrysosporium: Mechanism and degradation product identification

Antibiotic contaminants have become a severe environmental problem in recent years and finding effective ways to deal with this issue is of great importance. In this study, Phanerochaete chrysosporium was used to degrade sulfadiazine (SDZ), which is frequently detected in the culture medium of isolates from soil and surface water systems. The results demonstrate that 10 mg L^-1 SDZ can be completely degraded by P. chrysosporium under conditions of pH 5.7 and 30 °C within 6 days. The Q-Exactive-MS/MS analysis identified and confirmed several different SDZ degradation intermediates, and four proposed degradation pathways of SDZ were deduced. Moreover, enzyme activity tests revealed that manganese peroxidase and ligninolytic peroxidase played important roles in SDZ degradation. Moreover, a transcriptome analysis method was performed to explore the mechanism and pathways of SDZ degradation by P. chrysosporium in greater detail. The results of GO and KEGG analysis strongly suggest that the metabolism pathway is significantly activated and plays an important role in antibiotic degradation. Further, this is the first study to identify SDZ degradation intermediates and two main intermediates were found to be involved in possible SDZ degradation pathways. This study is also the first report results from RNA sequencing to evaluate genome-wide changes of P. chrysosporium to further explore SDZ degradation mechanism.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4 bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3 bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B. carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.

Extensive changes in gene expression and alternative splicing due to homoeologous exchange in rice segmental allopolyploids

We report rampant homoeologous exchanges in progenies of a newly synthesized rice segmental allotetraploid and demonstrate their consequences to changes of gene expression and alternative splicing. Allopolyploidization is recurrent across the tree of angiosperms and known as a driving evolutionary force in both plants and animals. A salient feature of allopolyploidization is the induction of homoeologous exchange (HE) events between the constituent subgenomes, which may in turn cause changes in gene expression, transcript alternative splicing, and phenotypic novelty. However, this issue has been poorly studied, largely because lack of a system in which the exact parentage donating the subgenomes is known and the HE events are occurring in real time. Here, we employed whole-genome re-sequencing and RNA-seq-based transcriptome profiling in four randomly chosen progeny individuals (at the 10th-selfed generation) of segmental allotetraploids that were constructed by colchicine-mediated whole-genome doubling of F1 hybrids between the two subspecies (japonica and indica) of Asian cultivated Oryza sativa. We show that rampant HE events occurred in these tetraploid individuals, which converted most of the otherwise heterozygous genomic regions into a homogenized state of one parental subgenome. We demonstrate that genes within these homogenized genomic regions in the tetraploids showed high frequencies of altered expression and enhanced alternative splicing relative to their counterparts in the corresponding diploid parents in the embryo tissue. Intriguingly, limited overlaps between the differentially expressed genes and the differential alternative spliced genes were identified, which were partitioned to distinctly enriched gene ontology terms. Together, our results indicate that HE is a major mechanism to rapidly generate novelty in gene expression and transcriptome diversity, which may facilitate phenotypic innovation in nascent allopolyploids and relevant to allopolyploid crop breeding.

Screening and analysis on the differentially expression genes between diploid and autotetraploid watermelon by using of digital gene expression profile

Abstract Synthetic polyploids are key breeding materials for watermelon. Compared with diploid watermelon, the tetraploid watermelon often exhibit wide phenotypic differences and differential gene expression. Digital gene expression (DGE) profile technique was performed in this study to present gene expression patterns in an autotetraploid and its progenitor diploid watermelon, and deferentially expressed genes (DEGs) related to the abiotic and biotic stress were also addressed. Altogether, 4,985 DEGs were obtained in the autotetraploid against its progenitor diploid, and 66.02% DEGs is up-regulated. GO analysis shows that these DEGs mainly distributed in ‘metabolic process’, ‘cell’ and ‘catalytic activity’. KEGG analysis revealed that these DEGs mainly cover ‘metabolic pathways’, ‘secondary metabolites’ and ‘ribosome’. Moreover, 134 tolerance related DEGs were identified which cover osmotic adjustment substance, protective enzymes/protein, signaling proteins and pathogenesis-related proteins. This study present the differential expression of stress related genes and global gene expression patterns at background level in autotetraploid watermelons. These new evidences could supplement the molecular theoretical basis for the better resistance after the genome doubling in the gourd family.

Transcriptome-based gene expression profiling of diploid radish (Raphanus sativus L.) and the corresponding autotetraploid

Polyploidy is an important evolutionary factor in most land plant lineages which possess more than two complete sets of chromosomes. Radish (Raphanus sativus L.) is an economically annual/biennial root vegetable crop worldwide. However, the expression patterns of duplicated homologs involved in the autopolyploidization remains unclear. In present study, the autotetraploid radish plants (2n = 4x = 36) were produced with colchicine and exhibited an increase in the size of flowers, leaves, stomata and pollen grains. The differential gene expression (DGE) profiling was performed to investigate the differences in gene expression patterns between diploid and its corresponding autotetraploid by RNA-Sequencing (RNA-Seq). Totally, 483 up-regulated differentially expressed genes (DEGs) and 408 down-regulated DEGs were detected in diploid and autotetraploid radishes, which majorly involved in the pathways of hormones, photosynthesis and stress response. Moreover, the xyloglucan endotransglucosylase/hydrolase (XTH) and pectin methylesterases (PME) family members related to cell enlargement and cell wall construction were found to be enriched in GO enrichment analysis, of which XTH family members enriched in "apoplast" and "cell wall" terms, while PME family members enriched in "cell wall" term. Reverse-transcription quantitative PCR (RT-qPCR) analysis indicated that the expression profile of DEGs were consistent with results from the RNA-Seq analysis. The DEGs involved in cell wall construction and auxin metabolism were predicted to be associated with organs size increase of autotetraploid radishes in the present study. These results could provide valuable information for elucidating the molecular mechanism underlying polyploidization and facilitating further genetic improvements of important traits in radish breeding programs.

Which Genes in a Typical Intertidal Seagrass (Zostera japonica) Indicate Copper-, Lead-, and Cadmium Pollution?

Healthy seagrasses are considered a prime indicator of estuarine and coastal ecosystem function; however, as the only group of flowering plants recolonizing the sea, seagrasses are frequently exposed to anthropogenic heavy metal pollutants, which are associated with high levels of molecular damage. To determine whether biologically relevant concentrations of heavy metals cause systematic alterations in RNA expression patterns, we performed a gene expression study using transcriptome analyses (RNA-seq). We exposed the typical intertidal seagrass Zostera japonica to 0 and 50 μM of copper (Cu), lead (Pb), and cadmium (Cd) under laboratory conditions. A total of 18,266 differentially expressed genes (DEGs) were identified, of which 2001 co-expressed genes directly related by Cu, Pb, and Cd stress. We also examined the effects of short-term heavy metal Cu, Pb, and Cd pulses on the accumulation of metals in Z. japonica and showed metal concentrations were higher in the shoots than in roots. Twelve differentially expressed genes were further analyzed for expression differences using real-time quantitative polymerase chain reaction (RT-qPCR). Our data suggest that as coastal seawater pollution worsens, the sensitive genes identified in this study may be useful biomarkers of sublethal effects and provide fundamental information for Z. japonica resistant gene engineering.

Transcriptome analysis of alternative splicing in peanut (Arachis hypogaea L.)

BACKGROUND

Alternative splicing (AS) represents a mechanism widely used by eukaryotes for the post-transcriptional regulation of genes. The detailed exploration of AS in peanut has not been documented.

RESULTS

The strand-specific RNA-Seq technique was exploited to characterize the distribution of AS in the four samples of peanut (FH1-seed1, FH1-seed2, FH1-root and FH1-leaf). AS was detected as affecting around 37.2% of the full set of multi-exon genes. Some of these genes experienced AS throughout the plant, while in the case of others, the effect was organ-specific. Overall, AS was more frequent in the seed than in either the root or leaf. The predominant form of AS was intron retention, and AS in transcription start site and transcription terminal site were commonly identified in all the four samples. It is interesting that in genes affected by AS, the majority experienced only a single type of event. Not all of the in silico predicted transcripts appeared to be translated, implying that these are either degraded or sequestered away from the translation machinery. With respect to genes involved in fatty acid metabolism, about 61.6% were shown to experience AS.

CONCLUSION

Our report contributes significantly in AS analysis of peanut genes in general, and these results have not been mentioned before. The specific functions of different AS forms need further investigation.

Patterns of alternative splicing vary between species during heat stress

Plants have evolved a variety of mechanisms to respond and adapt to abiotic stress. High temperature stress induces the heat shock response. During the heat shock response a large number of genes are up-regulated, many of which code for chaperone proteins that prevent irreversible protein aggregation and cell death. However, it is clear that heat shock is not the only mechanism involved in the plant heat stress response. Alternative splicing (AS) is also important during heat stress since this post-transcriptional regulatory mechanism can produce significant transcriptome and proteome variation. In this study, we examine AS during heat stress in the model species Arabidopsis thaliana and in the highly thermotolerant native California mustard Boechera depauperata. Analyses of AS during heat stress revealed that while a significant number of genes undergo AS and are differentially expressed (DE) during heat stress, some undergo both AS and DE. Analysis of the functional categories of genes undergoing AS indicated that enrichment patterns are different in the two species. Categories enriched in B. depauperata included light response genes and numerous abiotic stress response genes. Categories enriched in A. thaliana, but not in B. depauperata, included RNA processing and nucleotide binding. We conclude that AS and DE are largely independent responses to heat stress. Furthermore, this study reveals significant differences in the AS response to heat stress in the two related mustard species. This indicates AS responses to heat stress are species-specific. Future studies will explore the role of AS of specific genes in organismal thermotolerance.

Resistance mechanisms and their difference between the root and leaf of Microsorum pteropus - A novel potential aquatic cadmium hyperaccumulator

Microsorum pteropus (M. pteropus), an aquatic Polypodiaceae fern, was identified as a novel potential cadmium (Cd) hyperaccumulator in our previous study. This study reveals the Cd-resistance mechanisms and their difference between the root and leaf of M. pteropus based on analyses of photosynthesis, antioxidant systems and gene expression. A high level of Cd at 500μM was used to treat the samples to test the effects of this compound. Superoxide dismutase (SOD), peroxidase (POD), malondialdehyde (MDA) and flavonoids were used as indicators for antioxidant system changes. Five chlorophyll fluorescent parameters including the maximal photochemical efficiency of photosystem II (F_v/F_m), effective quantum yield of photosystem II (Y(II)), photochemical quenching (qP), nonphotochemical quenching (qN) and electron transport rate (ETR) were measured to determine the photosynthetic changes. RNA-sequencing analysis was used to study the changes in gene expression. The results showed that after exposure to high levels of Cd, the concentrations of enzymatic oxidants (SOD and POD) were significantly increased, while the MDA levels were significantly decreased. There were no significant changes for the chlorophyll fluorescent parameters during Cd stress, which indicates that M. pteropus is highly effective at protecting itself. Certain functional genes, including photosystem genes and secondary metabolites, had significantly altered levels of expression. Different Cd-resistance mechanisms were found between the root and leaf tissues of M. pteropus. The root tissues of M. pteropus resist Cd damage using antioxidants, while its leaf tissues mainly protect themselves using photosystem self-protection.

De novo transcriptome assembly, development of EST-SSR markers and population genetic analyses for the desert biomass willow, Salix psammophila

Salix psammophila, a sandy shrub known as desert willow, is regarded as a potential biomass feedstock and plays an important role in maintaining local ecosystems. However, a lack of genomic data and efficient molecular markers limit the study of its population evolution and genetic breeding. In this study, chromosome counts, flow cytometry and SSR analyses indicated that S. psammophila is tetraploid. A total of 6,346 EST-SSRs were detected based on 71,458 de novo assembled unigenes from transcriptome data. Twenty-seven EST-SSR markers were developed to evaluate the genetic diversity and population structure of S. psammophila from eight natural populations in Northern China. High levels of genetic diversity (mean 10.63 alleles per locus; mean HE 0.689) were dectected in S. psammophila. The weak population structure and little genetic differentiation (pairwise FST = 0.006-0.016) were found among Population 1-Population 7 (Pop1-Pop7; Inner Mongolia and Shaanxi), but Pop8 (Ningxia) was clearly separated from Pop1-Pop7 and moderate differentiation (pairwise FST = 0.045-0.055) was detected between them, which may be influenced by local habitat conditions. Molecular variance analyses indicated that most of the genetic variation (94.27%) existed within populations. These results provide valuable genetic informations for natural resource conservation and breeding programme optimisation of S. psammophila.

Comparative Proteomic Analysis of Mature Pollen in Triploid and Diploid Populus deltoides

Ploidy affects plant growth vigor and cell size, but the relative effects of pollen fertility and allergenicity between triploid and diploid have not been systematically examined. Here we performed comparative analyses of fertility, proteome, and abundances of putative allergenic proteins of pollen in triploid poplar ‘ZhongHuai1’ (‘ZH1’, triploid) and ‘ZhongHuai2’ (‘ZH2’, diploid) generated from the same parents. The mature pollen was sterile in triploid poplar ‘ZH1’. By applying two-dimensional gel electrophoresis (2-DE), a total of 72 differentially expressed protein spots (DEPs) were detected in triploid poplar pollen. Among them, 24 upregulated and 43 downregulated proteins were identified in triploid poplar pollen using matrix-assisted laser desorption/ionisation coupled with time of-flight tandem mass spectrometer analysis (MALDI-TOF/TOF MS/MS). The main functions of these DEPs were related with “S-adenosylmethionine metabolism”, “actin cytoskeleton organization”, or “translational elongation”. The infertility of triploid poplar pollen might be related to its abnormal cytoskeletal system. In addition, the abundances of previously identified 28 putative allergenic proteins were compared among three poplar varieties (‘ZH1’, ‘ZH2’, and ‘2KEN8‘). Most putative allergenic proteins were downregulated in triploid poplar pollen. This work provides an insight into understanding the protein regulation mechanism of pollen infertility and low allergenicity in triploid poplar, and gives a clue to improving poplar polyploidy breeding and decreasing the pollen allergenicity.

Genome-wide identification of microRNAs in pomegranate (Punica granatum L.) by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs), a class of small non-coding endogenous RNAs that regulate gene expression post-transcriptionally, play multiple key roles in plant growth and development and in biotic and abiotic stress response. Knowledge and roles of miRNAs in pomegranate fruit development have not been explored.

RESULTS

Pomegranate, which accumulates a large amount of anthocyanins in skin and arils, is valuable to human health, mainly because of its antioxidant properties. In this study, we developed a small RNA library from pooled RNA samples from young seedlings to mature fruits and identified both conserved and pomegranate-specific miRNA from 29,948,480 high-quality reads. For the pool of 15- to 30-nt small RNAs, ~50 % were 24 nt. The miR157 family was the most abundant, followed by miR156, miR166, and miR168, with variants within each family. The base bias at the first position from the 5' end had a strong preference for U for most 18- to 26-nt sRNAs but a preference for A for 18-nt sRNAs. In addition, for all 24-nt sRNAs, the nucleotide U was preferred (97 %) in the first position. Stem-loop RT-qPCR was used to validate the expression of the predominant miRNAs and novel miRNAs in leaves, male and female flowers, and multiple fruit developmental stages; miR156, miR156a, miR159a, miR159b, and miR319b were upregulated during the later stages of fruit development. Higher expression of miR156 in later fruit developmental may positively regulate anthocyanin biosynthesis by reducing SPL transcription factor. Novel miRNAs showed variation in expression among different tissues. These novel miRNAs targeted different transcription factors and hormone related regulators. Gene ontology and KEGG pathway analyses revealed predominant metabolic processes and catalytic activities, important for fruit development. In addition, KEGG pathway analyses revealed the involvement of miRNAs in ascorbate and linolenic acid, starch and sucrose metabolism; RNA transport; plant hormone signaling pathways; and circadian clock.

CONCLUSION

Our first and preliminary report of miRNAs will provide information on the synthesis of biochemical compounds of pomegranate for future research. The functions of the targets of the novel miRNAs need further investigation.

Transcriptome profiling of differentially expressed genes in floral buds and flowers of male sterile and fertile lines in watermelon

Background

Male sterility is an important mechanism for the production of hybrid seeds in watermelon. Although fruit development has been studied extensively in watermelon, there are no reports on gene expression in floral organs. In this study, RNA-sequencing (RNA-seq) was performed in two near-isogenic watermelon lines (genic male sterile [GMS] line, DAH3615-MS and male fertile line, DAH3615) to identify the differentially expressed genes (DEGs) related to male sterility.

Results

DEG analysis showed that 1259 genes were significantly associated with male sterility at a FDR P-value of < 0.01. Most of these genes were only expressed in the male fertile line. In addition, 11 functional clusters were identified using DAVID functional classification analysis. Of detected genes in RNA-seq analysis, 19 were successfully validated by qRT-PCR.

Conclusions

In this study, we carried out a comprehensive floral transcriptome sequence comparison of a male fertile line and its near-isogenic male sterile line in watermelon. This analysis revealed essential genes responsible for stamen development, including pollen development and pollen tube elongation, and allowed their functional classification. These results provided new information on global mechanisms related to male sterility in watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2186-9) contains supplementary material, which is available to authorized users.

Transcriptome Analysis of Invasive Plants in Response to Mineral Toxicity of Reclaimed Coal-Mine Soil in the Appalachian Region

Efficient postmining reclamation requires successful revegetation. By using RNA sequencing, we evaluated the growth response of two invasive plants, goutweed (Aegopodium podagraria L.) and mugwort (Artemisia vulgaris), grown in two Appalachian acid-mine soils (MS-I and -II, pH ∼ 4.6). Although deficient in macronutrients, both soils contained high levels of plant-available Al, Fe and Mn. Both plant types showed toxicity tolerance, but metal accumulation differed by plant and site. With MS-I, Al accumulation was greater for mugwort than goutweed (385 ± 47 vs 2151 ± 251 μg g-1). Al concentration was similar between mine sites, but its accumulation in mugwort was greater with MS-I than MS-II, with no difference in accumulation by site for goutweed. An in situ approach revealed deregulation of multiple factors such as transporters, transcription factors, and metal chelators for metal uptake or exclusion. The two plant systems showed common gene expression patterns for different pathways. Both plant systems appeared to have few common heavy-metal pathway regulators addressing mineral toxicity/deficiency in both mine sites, which implies adaptability of invasive plants for efficient growth at mine sites with toxic waste. Functional genomics can be used to screen for plant adaptability, especially for reclamation and phytoremediation of contaminated soils and waters.

RNA-seq-based transcriptome profiling reveals differential gene expression in the lungs of Sprague-Dawley rats during early-phase acute hypobaric hypoxia

Individuals subjected to hypobaric hypoxia at high altitudes may exhibit differential physiological responses in terms of susceptibility and tolerance to the development of hypoxia-related disorders. We studied early-phase gene expression in the lungs of Sprague-Dawley rats exhibiting such differential physiological responses after exposure to acute hypobaric hypoxia for 1 h at a simulated altitude of 9144 m. RNA-seq transcriptome profiling of lung tissues revealed differential gene expression in tolerant and susceptible groups, subsequently validated by qRT-PCR for ten selected differentially expressed genes. The gene expression pattern indicated hypometabolism and negative regulation of vasoconstriction in all groups except susceptible rats, coupled with altered MAPK, p53 and JAK-STAT signaling. Upregulation of early-phase response genes including Dusp1 (dual specificity phosphatase), Cdkn1a (cyclin-dependent kinase inhibitor 1A), Txnip (thioredoxin-interacting protein), Rgs1 (regulator of G-protein signaling 1) and Rgs2 (regulator of G-protein signaling 2) in susceptible rats indicated a progression toward growth arrest and apoptosis. Enhanced expression of cell adhesion molecules, wound healing and repair bioprocesses was observed in tolerant males. Upregulated Kcnj15 (potassium inwardly rectifying channel subfamily j membrane 15) and Vsig4 (V-set and Ig domain containing 4) variants in tolerant females suggested adaptation to hypoxia possibly by fluid reabsorption to avoid edematous conditions and suppression of T cell proliferation to avoid acute lung inflammation. Our study might help in understanding the molecular-physiological mechanisms associated with progressive damage in the lung tissues of susceptible and tissue-protective measures in tolerant rats during acute hypobaric hypoxia.