The role of introns in the conservation of the metabolic genes of Arabidopsis thaliana

6Pgdh polymorphism in wild bulb mite populations: prevalence, environmental correlates and life history trade-offs

Genetic polymorphism in key metabolic genes plays a pivotal role in shaping phenotypes and adapting to varying environments. Polymorphism in the metabolic gene 6-phosphogluconate dehydrogenase (6Pgdh) in bulb mites, Rhizoglyphus robini is characterized by two alleles, S and F, that differ by a single amino acid substitution and correlate with male reproductive fitness. The S-bearing males demonstrate a reproductive advantage. Although the S allele rapidly fixes in laboratory settings, the persistence of polymorphic populations in the wild is noteworthy. This study examines the prevalence and stability of 6Pgdh polymorphism in natural populations across Poland, investigating potential environmental influences and seasonal variations. We found widespread 6Pgdh polymorphism in natural populations, with allele frequencies varying across locations and sampling dates but without clear geographical or seasonal clines. This widespread polymorphism and spatio-temporal variability may be attributed to population demography and gene flow between local populations. We found some correlation between soil properties, particularly cation content (Na, K, Ca, and Mg) and 6Pgdh allele frequencies, showcasing the connection between mite physiology and soil characteristics and highlighting the presence of environment-dependent balancing selection. We conducted experimental fitness assays to determine whether the allele providing the advantage in male-male competition has antagonistic effects on life-history traits and if these effects are temperature-dependent. We found that temperature does not differentially influence development time or juvenile survival in different 6Pgdh genotypes. This study reveals the relationship between genetic variation, environmental factors, and reproductive fitness in natural bulb mite populations, shedding light on the dynamic mechanisms governing 6Pgdh polymorphism.

Genome-wide identification and expression analysis of the glutamate receptor gene family in sweet potato and its two diploid relatives

Plant glutamate receptor (GLR) homologs are crucial calcium channels that play an important role in plant development, signal transduction, and response to biotic and abiotic stresses. However, the GLR gene family has not yet been thoroughly and systematically studied in sweet potato. In this study, a total of 37 GLR genes were identified in the cultivated hexaploid sweet potato (Ipomoea batatas), and 32 GLR genes were discovered in each of the two diploid relatives (Ipomoea trifida and Ipomoea triloba) for the first time. Based on their evolutionary relationships to those of Arabidopsis, these GLRs were split into five subgroups. We then conducted comprehensive analysis to explore their physiological properties, protein interaction networks, promoter cis-elements, chromosomal placement, gene structure, and expression patterns. The results indicate that the homologous GLRs of the cultivated hexaploid sweet potato and its two relatives are different. These variations are reflected in their functions related to plant growth, hormonal crosstalk, development of tuberous roots, resistance to root rot, and responses to abiotic stress factors, all of which are governed by specific individual GLR genes. This study offers a comprehensive analysis of GLR genes in sweet potato and its two diploid relatives. It also provides a theoretical basis for future research into their regulatory mechanisms, significantly influencing the field of molecular breeding in sweet potatoes.

Genome-Wide Identification of the CBF Gene Family and ICE Transcription Factors in Walnuts and Expression Profiles under Cold Conditions

Cold stress impacts woody tree growth and perennial production, especially when the temperature rapidly changes in late spring. To address this issue, we conducted the genome-wide identification of two important transcription factors (TFs), CBF (C-repeat binding factors) and ICE (inducers of CBF expression), in three walnut (Juglans) genomes. Although the CBF and ICE gene families have been identified in many crops, very little systematic analysis of these genes has been carried out in J. regia and J. sigillata. In this study, we identified a total of 16 CBF and 12 ICE genes in three Juglans genomes using bioinformatics analysis. Both CBF and ICE had conserved domains, motifs, and gene structures, which suggests that these two TFs were evolutionarily conserved. Most ICE genes are located at both ends of the chromosomes. The promoter cis-regulatory elements of CBF and ICE genes are largely involved in light and phytohormone responses. Based on 36 RNA sequencing of leaves from four walnut cultivars ('Zijing', 'Lvling', 'Hongren', and 'Liao1') under three temperature conditions (8 °C, 22 °C, and 5 °C) conditions in late spring, we found that the ICE genes were expressed more highly than CBFs. Both CBF and ICE proteins interacted with cold-related proteins, and many putative miRNAs had interactions with these two TFs. These results determined that CBF1 and ICE1 play important roles in the tolerance of walnut leaves to rapid temperature changes. Our results provide a useful resource on the function of the CBF and ICE genes related to cold tolerance in walnuts.

Genome wide mining of SNPs and INDELs through ddRAD sequencing in Sahiwal cattle

The study was conducted in Sahiwal cattle for genome wide identification and annotation of single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs) in Sahiwal cattle. The double digest restriction-site associated DNA (ddRAD) sequencing, a reduced representation method was used for the identification of variants at nucleotide level. A total of 1,615,211 variants were identified at RD10 and Q30 consisting of 1,480,930 SNPs and 134,281 INDELs with respect to the Bos taurus reference genome. The SNPs were annotated for their location, impact and functional class. The SNPs identified in Sahiwal cattle were found to be associated with a total of 26,229 genes. A total of 1819 SNPs were annotated for 209 candidate genes associated with different production and reproduction traits. The variants identified in the present study may be useful to strengthen the existing bovine SNP chips for reducing the biasness over the taurine cattle breeds. The diversity analysis provides the insight of the genetic architecture of the Sahiwal population Studied. The large genetic variations identified at the nucleotide level provide ample scope for implementing an effective and efficient breed improvement programme for increasing the productivity of Sahiwal cattle.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen

BACKGROUND

In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood.

RESULTS

We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture.

CONCLUSIONS

Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

Genome-Wide Identification and Expression Analysis of ACTIN Family Genes in the Sweet Potato and Its Two Diploid Relatives

ACTINs are structural proteins widely distributed in plants. They are the main components of microfilaments and participate in many crucial physiological activities, including the maintenance of cell shape and cytoplasmic streaming. Meanwhile, ACTIN, as a housekeeping gene, is widely used in qRT-PCR analyses of plants. However, ACTIN family genes have not been explored in the sweet potato. In this study, we identified 30, 39, and 44 ACTINs in the cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90) and its two diploid relatives, Ipomoea trifida (2n = 2x = 30) and Ipomoea triloba (2n = 2x = 30), respectively, via analysis of their genome structure and by phylogenetic characterization. These ACTINs were divided into six subgroups according to their phylogenetic relationships with Arabidopsis thaliana. The physiological properties of the protein, chromosome localization, phylogenetic relationship, gene structure, promoter cis-elements, protein interaction networks, and expression patterns of these 113 ACTINs were systematically investigated. The results suggested that homologous ACTINs are differentiated in the sweet potato and its two diploid relatives, and play various vital roles in plant growth, tuberous root development, hormone crosstalk, and abiotic stress responses. Some stable ACTINs that could be used as internal reference genes were found in the sweet potato and its two diploid relatives, e.g., IbACTIN18, -20, and -16.2; ItfACTIN2.2, -16, and -10; ItbACTIN18 and -19.1. This work provides a comprehensive comparison and furthers our understanding of the ACTIN genes in the sweet potato and its two diploid relatives, thereby supplying a theoretical foundation for their functional study and further facilitating the molecular breeding of sweet potatoes.

Integrated bioinformatic and physiological analyses reveal the pivotal role of hydrogen sulfide in enhancing low-temperature tolerance in alfalfa

Hydrogen sulfide (H₂ S) is an important gaseous signal molecule that regulates plant growth and stress resistance. However, research on the H₂ S synthase (HSase) genes is still limited in the model legume plant Medicago truncatula Gaertn. In the present study, a total of 40 HSase family members were first identified and analyzed in the M. truncatula genome, and these genes distributed across eight chromosomes and were clustered into five groups (I-V) based on their conserved gene structures and protein motifs. Expression analysis revealed that the MtHSase genes were expressed in all the tested abiotic stresses, albeit with expression level differences. This study also showed that H₂ S improves low temperature tolerance of alfalfa seedlings by regulating the antioxidant defense system and enhancing photosynthetic capacity. Thus, the study provides new insights into how the H₂ S signal regulates tolerance to low-temperature stress and provides the basis for further gene function and detection.

Genome-Wide Identification and Expression Analysis of SWEET Family Genes in Sweet Potato and Its Two Diploid Relatives

Sugar Will Eventually be Exported Transporter (SWEET) proteins are key transporters in sugar transportation. They are involved in the regulation of plant growth and development, hormone crosstalk, and biotic and abiotic stress responses. However, SWEET family genes have not been explored in the sweet potato. In this study, we identified 27, 27, and 25 SWEETs in cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90) and its two diploid relatives, Ipomoea trifida (2n = 2x = 30) and Ipomoea triloba (2n = 2x = 30), respectively. These SWEETs were divided into four subgroups according to their phylogenetic relationships with Arabidopsis. The protein physiological properties, chromosome localization, phylogenetic relationships, gene structures, promoter cis-elements, protein interaction networks, and expression patterns of these 79 SWEETs were systematically investigated. The results suggested that homologous SWEETs are differentiated in sweet potato and its two diploid relatives and play various vital roles in plant growth, tuberous root development, carotenoid accumulation, hormone crosstalk, and abiotic stress response. This work provides a comprehensive comparison and furthers our understanding of the SWEET genes in the sweet potato and its two diploid relatives, thereby supplying a theoretical foundation for their functional study and further facilitating the molecular breeding of sweet potato.

Interplay between gene expression and gene architecture as a consequence of gene and genome duplications: evidence from metabolic genes of Arabidopsis thaliana

Gene and genome duplications have been widespread during the evolution of flowering plant which resulted in the increment of biological complexity as well as creation of plasticity of a genome helping the species to adapt to changing environments. Duplicated genes with higher evolutionary rates can act as a mechanism of generating novel functions in secondary metabolism. In this study, we explored duplication as a potential factor governing the expression heterogeneity and gene architecture of Primary Metabolic Genes (PMGs) and Secondary Metabolic Genes (SMGs) of Arabidopsis thaliana. It is remarkable that different types of duplication processes controlled gene expression and tissue specificity differently in PMGs and SMGs. A complex relationship exists between gene architecture and expression patterns of primary and secondary metabolic genes. Our study reflects, expression heterogeneity and gene structure variation of primary and secondary metabolism in Arabidopsis thaliana are partly results of duplication events of different origins. Our study suggests that duplication has differential effect on PMGs and SMGs regarding expression pattern by controlling gene structure, epigenetic modifications, multifunctionality and subcellular compartmentalization. This study provides an insight into the evolution of metabolism in plants in the light of gene and genome scale duplication.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01188-2.

Comprehensive Evolutionary Analysis of CPP Genes in Brassica napus L. and Its Two Diploid Progenitors Revealing the Potential Molecular Basis of Allopolyploid Adaptive Advantage Under Salt Stress

Allopolyploids exist widely in nature and have strong environmental adaptability. The typical allopolyploid Brassica napus L. is a widely cultivated crop, but whether it is superior to its diploid progenitors in abiotic stress resistance and the key genes that may be involved are not fully understood. Cystein-rich polycomb-like protein (CPP) genes encode critical transcription factors involved in the response of abiotic stress, including salt stress. To explore the potential molecular basis of allopolyploid adaptation to salt stress, we comprehensively analyzed the characteristics and salt stress response of the CPP genes in B. napus and its two diploid progenitors in this study. We found some molecular basis that might be associated with the adaptability of B. napus, including the expansion of the CPP gene family, the acquisition of introns by some BnCPPs, and abundant cis-acting elements upstream of BnCPPs. We found two duplication modes (whole genome duplication and transposed duplication) might be the main reasons for the expansion of CPP gene family in B. napus during allopolyploidization. CPP gene expression levels and several physiological indexes were changed in B. napus and its diploid progenitors after salt stress, suggesting that CPP genes might play important roles in the response of salt stress. We found that some BnCPPs might undergo new functionalization or subfunctionalization, and some BnCPPs also show biased expression, which might contribute to the adaptation of B. napus under saline environment. Compared with diploid progenitors, B. napus showed stronger physiological responses, and BnCPP gene expression also showed higher changes after salt stress, indicating that the allopolyploid B. napus had an adaptive advantage under salt stress. This study could provide evidence for the adaptability of polyploid and provide important clues for the study of the molecular mechanism of salt stress resistance in B. napus.

The impacts of allopolyploidization on Methyl-CpG-Binding Domain (MBD) gene family in Brassica napus

BACKGROUND

Polyploidization promotes species formation and is widespread in angiosperms. Genome changes dramatically bring opportunities and challenges to plants after polyploidy. Methyl-CpG-Binding Domain (MBD) proteins can recognize and bind to methylation sites and they play an important role in the physiological process related to methylation in animals and plants. However, research on the influence of the allopolyploidization process on the MBD gene family is still lacking, so it is necessary to conduct a comprehensive analysis.

RESULTS

In this study, twenty-two, ten and eleven MBD genes were identified in the genome of allotetraploid B. napus and its diploid ancestors, B. rapa and B. oleracea, respectively. Based on the clades of the MBD gene in Arabidopsis, rice and maize, we divided the new phylogenetic tree into 8 clades. Among them, the true MBD genes in Brassica existed in only 5 clades. Clade IV and Clade VI were unique in term of MBD genes in dicotyledons. Ka/Ks calculations showed that MBD genes underwent purifying selection in Brassica and may retain genes through sequence or functional differentiation early in evolution. In the process of allopolyploidization, the number of MBD gene introns increased, and the protein motifs changed. The MBD proteins had their own special motifs in each clade, and the MBD domains were only conserved in their clades. At the same time, the MBD genes were expressed in flower, leaf, silique, and stem tissues, and the expression levels of the different genes were significantly different, while the tissue specificity was not obvious. The allopolyploidization process may increase the number of cis-acting elements and activate the transposable elements. During allopolyploidization, the expression pattern of the MBD gene changes, which may be regulated by cis-acting elements and transposable elements. The number imbalance of cis-acting elements and transposable elements in A_n and C_n subgenomes may also lead to biased A_n subgenome expression of the MBD gene in B. napus.

CONCLUSIONS

In this study, by evaluating the number, structure, phylogeny and expression of the MBD gene in B. napus and its diploid ancestors, we increased the understanding of MBD genes in allopolyploids and provided a reference for future analysis of allopolyploidization.

Identification and Characterization of Wall-Associated Kinase (WAK) and WAK-like (WAKL) Gene Family in Juglans regia and Its Wild Related Species Juglans mandshurica

Wall-associated kinase (WAK) and WAK-like kinase (WAKL) are receptor-like kinases (RLKs), which play important roles in signal transduction between the cell wall and the cytoplasm in plants. WAK/WAKLs have been studied in many plants, but were rarely studied in the important economic walnut tree. In this study, 27 and 14 WAK/WAKL genes were identified in Juglans regia and its wild related species Juglans mandshurica, respectively. We found tandem duplication might play a critical role in the expansion of WAK/WAKL gene family in J. regia, and most of the WAK/WAKL homologous pairs underwent purified selection during evolution. All WAK/WAKL proteins have the extracellular WAK domain and the cytoplasmic protein kinase domain, and the latter was more conserved than the former. Cis-acting elements analysis showed that WAK/WAKL might be involved in plant growth and development, plant response to abiotic stress and hormones. Gene expression pattern analysis further indicated that most WAK/WAKL genes in J. regia might play a role in the development of leaves and be involved in plant response to biotic stress. Our study provides a new perspective for the evolutionary analysis of gene families in tree species and also provides potential candidate genes for studying WAK/WAKL gene function in walnuts.

Genome-Wide Identification and Expression Analysis of JAZ Family Involved in Hormone and Abiotic Stress in Sweet Potato and Its Two Diploid Relatives

Jasmonate ZIM-domain (JAZ) proteins are key repressors of a jasmonic acid signaling pathway. They play essential roles in the regulation of plant growth and development, as well as environmental stress responses. However, this gene family has not been explored in sweet potato. In this study, we identified 14, 15, and 14 JAZs in cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90), and its two diploid relatives Ipomoea trifida (2n = 2x = 30) and Ipomoea triloba (2n = 2x = 30), respectively. These JAZs were divided into five subgroups according to their phylogenetic relationships with Arabidopsis. The protein physiological properties, chromosome localization, phylogenetic relationship, gene structure, promoter cis-elements, protein interaction network, and expression pattern of these 43 JAZs were systematically investigated. The results suggested that there was a differentiation between homologous JAZs, and each JAZ gene played different vital roles in growth and development, hormone crosstalk, and abiotic stress response between sweet potato and its two diploid relatives. Our work provided comprehensive comparison and understanding of the JAZ genes in sweet potato and its two diploid relatives, supplied a theoretical foundation for their functional study, and further facilitated the molecular breeding of sweet potato.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Cis-regulatory elements in conserved non-coding sequences of nuclear receptor genes indicate for crosstalk between endocrine systems

Nuclear receptors (NRs) are ligand-activated transcription factors that regulate gene expression when bound to specific DNA sequences. Crosstalk between steroid NR systems has been studied for understanding the development of hormone-driven cancers but not to an extent at a genetic level. This study aimed to investigate crosstalk between steroid NRs in conserved intron and exon sequences, with a focus on steroid NRs involved in prostate cancer etiology. For this purpose, we evaluated conserved intron and exon sequences among all 49 members of the NR Superfamily (NRS) and their relevance as regulatory sequences and NR-binding sequences. Sequence conservation was found to be higher in the first intron (35%), when compared with downstream introns. Seventy-nine percent of the conserved regions in the NRS contained putative transcription factor binding sites (TFBS) and a large fraction of these sequences contained splicing sites (SS). Analysis of transcription factors binding to putative intronic and exonic TFBS revealed that 5 and 16%, respectively, were NRs. The present study suggests crosstalk between steroid NRs, e.g., vitamin D, estrogen, progesterone, and retinoic acid endocrine systems, through cis-regulatory elements in conserved sequences of introns and exons. This investigation gives evidence for crosstalk between steroid hormones and contributes to novel targets for steroid NR regulation.

The functional analysis of ABCG transporters in the adaptation of pigeon pea (Cajanus cajan) to abiotic stresses

ATP-binding cassette (ABC) transporters are a class of proteins found in living organisms that mediate transmembrane transport by hydrolyzing ATP. They play a vital role in the physiological processes of growth and development in plants. The most numerous sub-type transporter in the ABC transporter family is the ABCG group and which have the most complex function in a plant’s response to abiotic stresses. Our study focused on the effect of ABCG transporters in the adaptation of the pigeon pea to adverse environments (such as drought, salt, temperature, etc.). We conducted a functional analysis of ABCG transporters in the pigeon pea and their role in response to abiotic stresses. A total of 51 ABCG genes (CcABCGs) were identified, and phylogenetic analysis was conducted. We also identified the physicochemical properties of the encoded proteins, predicted their subcellular localization, and identified of the conserved domains. Expression analysis showed that ABCG genes have different expression profiles with tissues and abiotic stresses. Our results showed that CcABCG28 was up-regulated at low temperatures, and CcABCG7 was up-regulated with drought and aluminum stress. The initial results revealed that ABCG transporters are more effective in the abiotic stress resistance of pigeon peas, which improves our understanding of their application in abiotic stress resistance.

Genome-wide and expression pattern analysis of JAZ family involved in stress responses and postharvest processing treatments in Camellia sinensis

The JASMONATE-ZIM DOMAIN (JAZ) family genes are key repressors in the jasmonic acid signal transduction pathway. Recently, the JAZ gene family has been systematically characterized in many plants. However, this gene family has not been explored in the tea plant. In this study, 13 CsJAZ genes were identified in the tea plant genome. Phylogenetic analysis showed that the JAZ proteins from tea and other plants clustered into 11 sub-groups. The CsJAZ gene transcriptional regulatory network predictive and expression pattern analyses suggest that these genes play vital roles in abiotic stress responses, phytohormone crosstalk and growth and development of the tea plant. In addition, the CsJAZ gene expression profiles were associated with tea postharvest processing. Our work provides a comprehensive understanding of the CsJAZ family and will help elucidate their contributions to tea quality during tea postharvest processing.

Association mapping of wheat Fusarium head blight resistance-related regions using a candidate-gene approach and their verification in a biparental population

Markers, located in Dicer1 and Ara6 genes, which are likely involved in cross-kingdom RNA trafficking, are associated with FHB resistance in GABI wheat population and were validated in biparental population. Association studies are a common approach to detect marker-trait associations for Fusarium head blight (FHB) resistance in wheat (Triticum aestivum), although verification of detected associations is exceptional. In the present study, candidate-gene association mapping (CG) of genes from silencing and secretory pathways, which may be involved in wheat resistance against FHB and cross-kingdom RNA trafficking, was performed. Fourteen markers, located in nine genes, were tested for association with FHB resistance in 356 lines from the GABI (genome analysis of the biological system of plants) wheat population. Three markers located in the genes Dicer1 and Ara6 were shown to be significantly associated with the studied trait. Verification of this finding was performed using the recombinant inbred lines (RILs) population 'Apache × Biscay', segregating for four of our 14 selected markers. We could show association of the Ara6 marker with plant height as well as association with FHB resistance for three markers located in Rab5-like GTPase gene Ara6 and Dicer1. These results confirmed the trait-marker associations detected also in the CG approach. Gene products of the associated genes are involved in response of the plant to pathogens, plant metabolism and may be involved in cross-kingdom RNA trafficking efficiency. The markers detected in the GABI wheat population, which were also validated in the biparental population, can potentially be used in wheat breeding.

Genome-wide identification and analysis of the EIN3/EIL gene family in allotetraploid Brassica napus reveal its potential advantages during polyploidization

BACKGROUND

Polyploidization is a common event in the evolutionary history of angiosperms, and there will be some changes in the genomes of plants other than a simple genomic doubling after polyploidization. Allotetraploid Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are a good group for studying the problems associated with polyploidization. On the other hand, the EIN3/EIL gene family is an important gene family in plants, all members of which are key genes in the ethylene signaling pathway. Until now, the EIN3/EIL gene family in B. napus and its diploid progenitors have been largely unknown, so it is necessary to comprehensively identify and analyze this gene family.

RESULTS

In this study, 13, 7 and 7 EIN3/EIL genes were identified in B. napus (2n = 4x = 38, A_nC_n), B. rapa (2n = 2x = 20, A_r) and B. oleracea (2n = 2x = 18, C_o). All of the identified EIN3/EIL proteins were divided into 3 clades and further divided into 8 sub-clades. Ka/Ks analysis showed that all identified EIN3/EIL genes underwent purifying selection after the duplication events. Moreover, gene structure analysis showed that some EIN3/EIL genes in B. napus acquired introns during polyploidization, and homolog expression bias analysis showed that B. napus was biased towards its diploid progenitor B. rapa. The promoters of the EIN3/EIL genes in B. napus contained more cis-acting elements, which were mainly involved in endosperm gene expression and light responsiveness, than its diploid progenitors. Thus, B. napus might have potential advantages in some biological aspects.

CONCLUSIONS

The results indicated allotetraploid B. napus might have potential advantages in some biological aspects. Moreover, our results can increase the understanding of the evolution of the EIN3/EIL gene family in B. napus, and provided more reference for future research about polyploidization.

Genome-wide analyses supported by RNA-Seq reveal non-canonical splice sites in plant genomes

BACKGROUND

Most eukaryotic genes comprise exons and introns thus requiring the precise removal of introns from pre-mRNAs to enable protein biosynthesis. U2 and U12 spliceosomes catalyze this step by recognizing motifs on the transcript in order to remove the introns. A process which is dependent on precise definition of exon-intron borders by splice sites, which are consequently highly conserved across species. Only very few combinations of terminal dinucleotides are frequently observed at intron ends, dominated by the canonical GT-AG splice sites on the DNA level.

RESULTS

Here we investigate the occurrence of diverse combinations of dinucleotides at predicted splice sites. Analyzing 121 plant genome sequences based on their annotation revealed strong splice site conservation across species, annotation errors, and true biological divergence from canonical splice sites. The frequency of non-canonical splice sites clearly correlates with their divergence from canonical ones indicating either an accumulation of probably neutral mutations, or evolution towards canonical splice sites. Strong conservation across multiple species and non-random accumulation of substitutions in splice sites indicate a functional relevance of non-canonical splice sites. The average composition of splice sites across all investigated species is 98.7% for GT-AG, 1.2% for GC-AG, 0.06% for AT-AC, and 0.09% for minor non-canonical splice sites. RNA-Seq data sets of 35 species were incorporated to validate non-canonical splice site predictions through gaps in sequencing reads alignments and to demonstrate the expression of affected genes.

CONCLUSION

We conclude that bona fide non-canonical splice sites are present and appear to be functionally relevant in most plant genomes, although at low abundance.