Genome-wide identification of the COBRA-like gene family and expression pattern analysis under abiotic stresses of Sorghum bicolor (L.)

On the Expansion of “Dangerous” Gene Repertoires by Whole-Genome Duplications in Early Vertebrates

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Casparian strip membrane domain proteins like (CASPL), exhibit profound associations with root development, stress responsiveness and mineral element uptake in plants. Nonetheless, a comprehensive bioinformatics analysis of the ZmCASPL gene family in maize remains unreported. In the study, we have identified 47 ZmCASPL members at the whole-genome level, systematically classifying them into six distinct groups. Furthermore, our analysis revealed that the same group of ZmCASPL contains similar gene structures and conserved motifs. Duplication events showed whole genome duplication (WGD) and tandem duplication (TD) contribute to the generation of the ZmCASPL gene family together in maize, but the former plays a more prominent role. Furthermore, we observed that most ZmCASPL genes contain MYB-binding sites (CAACCA), which are associated with the Casparian strip. Utilizing RNA-seq data, we found that ZmCASPL21 and ZmCASPL47 are specifically highly expressed only in the roots. This finding implies that ZmCASPL21 and ZmCASPL47 may be involved in the Casparian strip development. Additionally, RNA-seq analysis illuminated that drought, salt, heat, cold stresses, low nitrogen and phosphorus conditions, as well as pathogen infection, significantly impact the expression patterns of ZmCASPL genes. RT-qPCR revealed that ZmCASPL 5/13/25/44 genes showed different expression patterns under PEG and NaCl treatments. Collectively, these findings provide a robust theoretical foundation for further investigations into the functional roles of the ZmCASPL gene family in maize.

Homeodomain-leucine zipper (HD-Zip) genes encode plant-specific transcription factors, which play important roles in plant growth, development, and response to environmental stress. These genes have not been fully studied in allopolyploid Brassica napus, an important kind of oil crop. In this study, 165 HD-Zip genes were identified in B. napus and classified into four subfamilies. If proteins belong to the same subfamily, they exhibit similarities in gene structure, motifs, and domain distribution patterns. BnHD-Zip genes were unevenly distributed in the An and Cn subgenomes. Whole genome triplication (WGT) events may be major mechanisms accounting for this gene family expansion. Orthologous gene analysis showed that the process of this gene family expansion was accompanied by domain loss. We further found three genes homologous to HB7 and three genes homologous to HB12, all induced by PEG, ABA, and NaCl treatment. HB7 could not form homodimers but could form heterodimers with HB12 based on yeast two-hybrid assays. The results of this study provide valuable information for further exploration of the HD-Zip gene family in B. napus.

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

BackgroundThe cytochrome P450 monooxygenase (P450) superfamily is involved in the biosynthesis of various primary and secondary metabolites. However, little is known about the effects of whole genome duplication (WGD) and tandem duplication (TD) events on the evolutionary history and functional divergence of P450s in Brassica after splitting from a common ancestor with Arabidopsis thaliana.ResultsUsing Hidden Markov Model search and manual curation, we detected that Brassica species have nearly 1.4-fold as many P450 members as A. thaliana. Most P450s in A. thaliana and Brassica species were located on pseudo-chromosomes. The inferred phylogeny indicated that all P450s were clustered into two different subgroups. Analysis of WGD event revealed that different P450 gene families had appeared after evolutionary events of species. For the TD event analyses, the P450s from TD events in Brassica species can be divided into ancient and recent parts. Our comparison of influence of WGD and TD events on the P450 gene superfamily between A. thaliana and Brassica species indicated that the family-specific evolution in the Brassica lineage can be attributed to both WGD and TD, whereas WGD was recognized as the major mechanism for the recent evolution of the P450 super gene family. Expression analysis of P450s from A. thaliana and Brassica species indicated that WGD-type P450s showed the same expression pattern but completely different expression with TD-type P450s across different tissues in Brassica species. Selection force analysis suggested that P450 orthologous gene pairs between A. thaliana and Brassica species underwent negative selection, but no significant differences were found between P450 orthologous gene pairs in A. thaliana–B. rapa and A. thaliana–B. oleracea lineages, as well as in different subgenomes in B. rapa or B. oleracea compared with A. thaliana.ConclusionsThis study is the first to investigate the effects of WGD and TD on the evolutionary history and functional divergence of P450 gene families in A. thaliana and Brassica species. This study provides a biology model to study the mechanism of gene family formation, particularly in the context of the evolutionary history of angiosperms, and offers novel insights for the study of angiosperm genomes.

The transcription factor known as TALE (three-amino acid loop extension) is essential for plant growth, cell differentiation and responses to environmental stresses. Although the TALE gene family has been identified in various plants, there has been a lack of comprehensive whole-genome identification and analysis in Gramineae species. In this study, 123 TALE family genes were identified in five Gramineae species, which can be categorized into two main subgroups: KONX and BELL. Most of the TALE genes in the same subgroup displayed analogous gene structures and conserved motifs. Furthermore, whole genome duplication (WGD) significantly contributes to the expansion of the TALE gene family in Gramineae. The promoter region of TALE genes in Gramineae contains a large number of cis-elements associated with abiotic stress and hormone response. Tissue-specific expression analysis indicated that most OsTALE, ZmTALE and AtTALE genes were highly expressed in stems and leaves. Additionally, RNA-seq data revealed that OsTALE, ZmTALE and AtTALE genes were found to respond to abiotic stress treatments. Furthermore, we found that the expression levels of SbTALE11/19 were up-regulated in response to PEG and NaCl treatment, respectively. This study provides a significant reference for further research on the biological function of TALE transcription factors in Gramineae plants.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae

Sugar transporter proteins (STPs), such as H+/sugar symporters, play essential roles in plants’ sugar transport, growth, and development, and possess an important potential to enhance plants’ performance of multiple agronomic traits, especially crop yield and stress tolerance. However, the evolutionary dynamics of this important gene family in Gramineae crops are still not well-documented and functional differentiation of rice STP genes remain unclear. To address this gap, we conducted a comparative genomic study of STP genes in seven representative Gramineae crops, which are Brachypodium distachyon (Bd), Hordeum vulgare (Hv), Setaria italica (Si), Sorghum bicolor (Sb), Zea mays (Zm), Oryza rufipogon (Or), and Oryza sativa ssp. japonica (Os). In this case, a total of 177 STP genes were identified and grouped into four clades. Of four clades, the Clade I, Clade III, and Clade IV showed an observable number expansion compared to Clade II. Our results of identified duplication events and divergence time of duplicate gene pairs indicated that tandem, Whole genome duplication (WGD)/segmental duplication events play crucial roles in the STP gene family expansion of some Gramineae crops (expect for Hv) during a long-term evolutionary process. However, expansion mechanisms of the STP gene family among the tested species were different. Further selective force studies revealed that the STP gene family in Gramineae crops was under purifying selective forces and different clades and orthologous groups with different selective forces. Furthermore, expression analysis showed that rice STP genes play important roles not only in flower organs development but also under various abiotic stresses (cold, high-temperature, and submergence stresses), blast infection, and wounding. The current study highlighted the expansion and evolutionary patterns of the STP gene family in Gramineae genomes and provided some important messages for the future functional analysis of Gramineae crop STP genes.

Elaeis guineensis is a tropical crop with high oil content, and the nutritional value of oil is very high. Amino acids not only affect the growth and development of plants but also act as intermediate metabolites of oils, which determine the final quality of oils. In this study, eleven amino acid permease genes (EgAAPs) in E. guineensis were identified from the amino acid transporter family. Real-time PCR results showed that seven EgAAPs (EgAAP1, EgAAP2, EgAAP3, EgAAP4, EgAAP6 EgAAP9, and EgAAP10) played an important role in vegetative growth because of their higher expression levels in roots and leaves, but EgAAP5, EgAAP7, EgAAP8 and EgAAP11 were also important for their higher expression levels in flowers or fruits. E. guineensis seedlings were treated with 20% PEG-6000 and 4 °C to induce drought stress and cold stress, respectively. The expression of six EgAAPs (EgAAP1, EgAAP2, EgAAP4, EgAAP6, EgAAP8 and EgAAP11) was decreased in both roots and leaves during cold treatment, and only the expression of EgAAP5 was increased in both roots and leaves after 6 h of cold treatment. The expression of six EgAAPs (EgAAP2, EgAAP4, EgAAP7, EgAAP8, EgAAP10 and EgAAP11) was decreased in roots but increased in leaves under PEG treatment, indicating an opposite pattern of expression levels of these EgAAPs in roots and leaves. However, only EgAAP5 had a similar pattern of expression levels between roots and leaves under PEG treatment. The findings provide information on how EgAAPs in E. guineensis are regulated during growth and development, and under various environmental stresses.

Genomic data have provided evidence of previously unknown ancient whole genome duplications (WGDs) and highlighted the role of WGDs in the evolution of many eukaryotic lineages. Ancient WGDs often are detected by examining distributions of synonymous substitutions per site (Ks) within a genome, or “Ks plots.” For example, WGDs can be detected from Ks plots by using univariate mixture models to identify peaks in Ks distributions. We performed gene family simulation experiments to evaluate the effects of different Ks estimation methods and mixture models on our ability to detect ancient WGDs from Ks plots. The simulation experiments, which accounted for variation in substitution rates and gene duplication and loss rates across gene families, tested the effects of WGD age and gene retention rates following WGD on inferring WGDs from Ks plots. Our simulations reveal limitations of Ks plot analyses. Strict interpretations of mixture model analyses often overestimate the number of WGD events, and Ks plot analyses typically fail to detect WGDs when ≤10% of the duplicated genes are retained following the WGD. However, WGDs can accurately be characterized over an intermediate range of Ks. The simulation results are supported by empirical analyses of transcriptomic data, which also suggest that biases in gene retention likely affect our ability to detect ancient WGDs. Although our results indicate mixture model results should be interpreted with great caution, using node-averaged Ks estimates and applying more appropriate mixture models can improve the accuracy of detecting WGDs.

High temperatures have a significant impact on plant growth and metabolism. In recent years, the fruit industry has faced a serious threat due to high-temperature stress on fruit plants caused by global warming. In the present study, we explored the molecular regulatory mechanisms that contribute to high-temperature tolerance in kiwifruit. A total of 36 Hsf genes were identified in the A. chinensis (Ac) genome, while 41 Hsf genes were found in the A. eriantha (Ae) genome. Phylogenetic analysis revealed the clustering of kiwifruit Hsfs into three distinct groups (groups A, B, and C). Synteny analysis indicated that the expansion of the Hsf gene family in the Ac and Ae genomes was primarily driven by whole genome duplication (WGD). Analysis of the gene expression profiles revealed a close relationship between the expression levels of Hsf genes and various plant tissues and stress treatments throughout fruit ripening. Subcellular localization analysis demonstrated that GFP-AcHsfA2a/AcHsfA7b and AcHsfA2a/AcHsfA7b -GFP were localized in the nucleus, while GFP-AcHsfA2a was also observed in the cytoplasm of Arabidopsis protoplasts. The results of real-time quantitative polymerase chain reaction (RT-qPCR) and dual-luciferase reporter assay revealed that the majority of Hsf genes, especially AcHsfA2a, were expressed under high-temperature conditions. In conclusion, our findings establish a theoretical foundation for analyzing the potential role of Hsfs in high-temperature stress tolerance in kiwifruit. This study also offers valuable information to aid plant breeders in the development of heat-stress-resistant plant materials.

In polyploids, whole genome duplication (WGD) played a significant role in genome expansion, evolution and diversification. Many gene families are expanded following polyploidization, with the duplicated genes functionally diversified by neofunctionalization or subfunctionalization. These mechanisms may support adaptation and have likely contributed plant survival during evolution. Flowering time is an important trait in plants, which affects critical features, such as crop yields. The flowering-time gene family is one of the largest expanded gene families in plants, with its members playing various roles in plant development. Here, we performed genome-wide identification and comparative analysis of flowering-time genes in three palnt families i.e., Malvaceae, Brassicaceae, and Solanaceae, which indicate these genes were expanded following the event/s of polyploidization. Duplicated genes have been retained during evolution, although genome reorganization occurred in their flanking regions. Further investigation of sequence conservation and similarity network analyses provide evidence for functional diversification of duplicated genes during evolution. These functionally diversified genes play important roles in plant development and provide advantages to plants for adaptation and survival in response to environmental changes encountered during evolution. Collectively, we show that flowering-time genes were expanded following polyploidization and retained as large gene family by providing advantages from functional diversification during evolution.

BackgroundThe circadian clock not only participates in regulating various stages of plant growth, development and metabolism, but confers plant environmental adaptability to stress such as drought. Pseudo-Response Regulators (PRRs) are important component of the central oscillator (the core of circadian clock) and play a significant role in plant photoperiod pathway. However, no systematical study about this gene family has been performed in cotton.MethodsPRR genes were identified in diploid and tetraploid cotton using bioinformatics methods to investigate their homology, duplication and evolution relationship. Differential gene expression, KEGG enrichment analysis and qRT-PCR were conducted to analyze PRR gene expression patterns under diurnal changes and their response to drought stress.ResultsA total of 44 PRR family members were identified in four Gossypium species, with 16 in G. hirsutum, 10 in G. raimondii, and nine in G. barbadense as well as in G. arboreum. Phylogenetic analysis indicated that PRR proteins were divided into five subfamilies and whole genome duplication or segmental duplication contributed to the expansion of Gossypium PRR gene family. Gene structure analysis revealed that members in the same clade are similar, and multiple cis-elements related to light and drought stress response were enriched in the promoters of GhPRR genes. qRT-PCR results showed that GhPRR genes transcripts presented four expression peaks (6 h, 9 h, 12 h, 15 h) during 24 h and form obvious rhythmic expression trend. Transcriptome data with PEG treatment, along with qRT-PCR verification suggested that members of clade III (GhPRR5a, b, d) and clade V (GhPRR3a and GhPRR3c) may be involved in drought response. This study provides an insight into understanding the function of PRR genes in circadian rhythm and in response to drought stress in cotton.

The LIM domain-containing proteins are dominantly found in plants and play a significant role in various biological processes such as gene transcription as well as actin cytoskeletal organization. Nevertheless, genome-wide identification as well as functional analysis of the LIM gene family have not yet been reported in the economically important plant sorghum (Sorghum bicolor L.). Therefore, we conducted an in silico identification and characterization of LIM genes in S. bicolor genome using integrated bioinformatics approaches. Based on phylogenetic tree analysis and conserved domain, we identified five LIM genes in S. bicolor (SbLIM) genome corresponding to Arabidopsis LIM (AtLIM) genes. The conserved domain, motif as well as gene structure analyses of the SbLIM gene family showed the similarity within the SbLIM and AtLIM members. The gene ontology (GO) enrichment study revealed that the candidate LIM genes are directly involved in cytoskeletal organization and various other important biological as well as molecular pathways. Some important families of regulating transcription factors such as ERF, MYB, WRKY, NAC, bZIP, C2H2, Dof, and G2-like were detected by analyzing their interaction network with identified SbLIM genes. The cis-acting regulatory elements related to predicted SbLIM genes were identified as responsive to light, hormones, stress, and other functions. The present study will provide valuable useful information about LIM genes in sorghum which would pave the way for the future study of functional pathways of candidate SbLIM genes as well as their regulatory factors in wet-lab experiments.