Comparative genomic and transcriptomic analysis of tandemly and segmentally duplicated genes in rice

Genomic organization and expression profiles of nitrogen assimilation genes in Glycine max

Background

Glutamine synthetase (GS), glutamate synthase (GOGAT), and nitrate reductase (NR) are key enzymes involved in nitrogen assimilation and metabolism in plants. However, the systematic analysis of these gene families lacked reports in soybean (Glycine max (L.) Merr.), one of the most important crops worldwide.

Methods

In this study, we performed genome-wide identification and characterization of GS, GOGAT, and NR genes in soybean under abiotic and nitrogen stress conditions.

Results

We identified a total of 10 GS genes, six GOGAT genes, and four NR genes in the soybean genome. Phylogenetic analysis revealed the presence of multiple isoforms for each gene family, indicating their functional diversification. The distribution of these genes on soybean chromosomes was uneven, with segmental duplication events contributing to their expansion. Within the nitrogen assimilation genes (NAGs) group, there was uniformity in the exon-intron structure and the presence of conserved motifs in NAGs. Furthermore, analysis of cis-elements in NAG promoters indicated complex regulation of their expression. RT-qPCR analysis of seven soybean NAGs under various abiotic stresses, including nitrogen deficiency, drought-nitrogen, and salinity, revealed distinct regulatory patterns. Most NAGs exhibited up-regulation under nitrogen stress, while diverse expression patterns were observed under salt and drought-nitrogen stress, indicating their crucial role in nitrogen assimilation and abiotic stress tolerance. These findings offer valuable insights into the genomic organization and expression profiles of GS, GOGAT, and NR genes in soybean under nitrogen and abiotic stress conditions. The results have potential applications in the development of stress-resistant soybean varieties through genetic engineering and breeding.

Comparative genomics points to tandem duplications of SAD gene clusters as drivers of increased α-linolenic (ω-3) content in S. hispanica seeds

Salvia hispanica L. (chia) is a source of abundant ω-3 polyunsaturated fatty acids (ω-3-PUFAs) that are highly beneficial to human health. The genomic basis for this accrued ω-3-PUFA content in this emerging crop was investigated through the assembly and comparative analysis of a chromosome-level reference genome for S. hispanica. The highly contiguous 321.5-Mbp genome assembly covering all six chromosomes enabled the identification of 32,922 protein-coding genes. Two whole-genome duplications (WGD) events were identified in the S. hispanica lineage. However, these WGD events could not be linked to the high α-linolenic acid (ALA, ω-3) accumulation in S. hispanica seeds based on phylogenomics. Instead, our analysis supports the hypothesis that evolutionary expansion through tandem duplications of specific lipid gene families, particularly the stearoyl-acyl carrier protein desaturase (ShSAD) gene family, is the main driver of the abundance of ω-3-PUFAs in S. hispanica seeds. The insights gained from the genomic analysis of S. hispanica will help establish a molecular breeding target that can be leveraged through genome editing techniques to increase ω-3 content in oil crops.

Multi-algorithm cooperation research of WRKY genes under nitrogen stress in Panax notoginseng

WRKY transcription factors play an important role in the immune system and the innate defense response of plants. WRKY transcription factors have great feedback on nitrogen stress. In this study, bioinformatics was used to detect the WRKYs of Panax notoginseng (PnWRKYs). The response of PnWRKYs under nitrogen stress was also well studied. PnWRKYs were distributed on 11 chromosomes. According to PnWRKY and Arabidopsis thaliana WRKY (AtWRKY) domains, these PnWRKY proteins were divided into three groups by phylogenetic analysis. MEME analysis showed that almost every member contained motif 1 and motif 2. PlantCARE online predicted the cis-acting elements of the promoter. PnWRKY gene family members obtained 22 pairs of repeat fragments by collinearity analysis. The expression levels of PnWRKYs in different parts (roots, flowers, and leafs) were analyzed by the gene expression pattern. They reflected tissue-specific expressions. The qRT-PCR experiments were used to detect 74 PnWRKYs under nitrogen stress. The results showed that the expression levels of 8 PnWRKYs were significantly induced. The PnWRKY gene family may be involved in biotic/abiotic stresses and hormone induction. This study will not only lay the foundation to explore the functions of PnWRKYs but also provide candidate genes for the future improvement of P. notoginseng.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Genome-Wide Characterization of Nitrogenase Reductase (nifH) Genes in the Sweet Potato [Ipomoea batatas (L.) Lam] and Its Wild Ancestors

The sweet potato (Ipomoea batatas (L.) Lam.) is an important and widely grown crop, and the nitrogenase reductase (nifH) gene is the most widely sequenced marker gene used to identify nitrogen-fixing bacteria and archaea. There have been many examples of the isolation of the diazotrophic endophytes in sweet potatoes, and there has been no report on whether sweet potatoes and their wild ancestors harbored nifH genes. In this study, a comprehensive analysis of nifH genes has been conducted on these species by using bioinformatics and molecular biology methods. A total of 20, 19 and 17 nifH genes were identified for the first time in sweet potatoes, I. trifida and I. triloba, respectively. Based on a phylogenetic analysis, all of the nifH genes, except for g10233.t1, itf14g14040.t1 and itb14g15470.t1, were clustered into five independent clades: I, II, III, IV and V. The nifH genes clustered in the same phylogenetic branch showed a more similar distribution of conserved motifs and exons–introns than those of the other ones. All of the identified genes were further mapped on the 15 chromosomes of the sweet potato, I. trifida and I. triloba. No segmental duplication was detected in each genome of three Ipomoea species, and 0, 8 and 7 tandemly duplicated gene pairs were detected in the genome of the sweet potato, I. trifida and I. triloba, respectively. Synteny analysis between the three Ipomoea species revealed that there were 7, 7 and 8 syntenic gene pairs of nifH genes detected between the sweet potato and I. trifida, between the sweet potato and I. triloba and between I. trifida and I. triloba, respectively. All of the duplicated and syntenic nifH genes were subjected to purifying selection inside duplicated genomic elements during speciation, except for the tandemly duplicated gene pair itf11g07340.t2_itf11g07340.t3, which was subjected to positive selection. Different expression profiles were detected in the sweet potato, I. trifida and I. triloba. According to the above results, four nifH genes of the sweet potato (g950, g16683, g27094 and g33987) were selected for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in two sweet potato cultivars (Eshu 15 and Long 9) under nitrogen deficiency (N0) and normal (N1) conditions. All of them were upregulated in the N1 treatment and were consistent with the analysis of the RNA-seq data. We hope that these results will provide new insights into the nifH genes in the sweet potato and its wild ancestors and will contribute to the molecular breeding of sweet potatoes in the future.

The WRKY transcription factor family in cowpea: Genomic characterization and transcriptomic profiling under root dehydration

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most tolerant legume crops to drought and salt stresses. WRKY transcription factor (TF) family members stand out among plant transcriptional regulators related to abiotic stress tolerance. However, little information is currently available on the expression of the cowpea WRKY gene family (VuWRKY) in response to water deficit. Thus, we analyzed genomic and transcriptomic data from cowpea to identify VuWRKY members and characterize their structure and transcriptional response under root dehydration stress. Ninety-two complete VuWRKY genes were found in the cowpea genome based on their domain characteristics. They were clustered into three groups: I (15 members), II (58), and III (16), while three genes were unclassified. Domain analysis of the encoded proteins identified four major variants of the conserved heptapeptide motif WRKYGQK. In silico analysis of VuWRKY gene promoters identified eight candidate binding motifs of cis-regulatory elements, regulated mainly by six TF families associated with abiotic stress responses. Ninety-seven VuWRKY modulated splicing variants associated with 55 VuWRKY genes were identified via RNA-Seq analysis available at the Cowpea Genomics Consortium (CpGC) database. qPCR analyses showed that 22 genes are induced under root dehydration, with VuWRKY18, 21, and 75 exhibiting the most significant induction levels. Given their central role in activating signal transduction cascades in abiotic stress response, the data provide a foundation for the targeted modification of specific VuWRKY family members to improve drought tolerance in this important climate-resilient legume in the developing world and beyond.

SNORKEL Genes Relating to Flood Tolerance Were Pseudogenized in Normal Cultivated Rice

SNORKEL1 (SK1) and SNORKEL2 (SK2) are ethylene responsive factors that regulate the internode elongation of deepwater rice in response to submergence. We previously reported that normal cultivated rice lacks SK genes because the Chromosome 12 region containing SK genes was deleted from its genome. However, no study has analyzed how the genome defect occurred in that region by comparing normal cultivated rice and deepwater rice. In this study, comparison of the sequence of the end of Chromosome 12, which contains SK genes, between normal and deepwater rice showed that complicated genome changes such as insertions, deletions, inversions, substitutions, and translocation occurred frequently in this region. In addition to SK1 and SK2 of deepwater rice, gene prediction analysis identified four genes containing AP2/ERF domains in normal cultivated rice and six in deepwater rice; we called these genes SK-LIKE (SKL) genes. SKs and SKLs were present in close proximity to each other, and the SKLs in normal cultivated rice were in tandem. These predicted genes belong to the same AP2/ERF subfamily and were separated into four types: SK1, SK2, SKL3, and SKL4. Sequence comparison indicated that normal cultivated rice possesses a gene with high homology to SK2, which we named SKL1. However, none of the predicted SKLs except for SKL3s were expressed during submergence. Although SKL3s were expressed in both normal and deepwater rice, normal rice does not undergo internode elongation, suggesting that its expression does not contribute to internode elongation. Plants overexpressing SKL1, which showed the most homology to SK2, underwent internode elongation similar to plants overexpressing SK1 and SK2 under normal growth conditions. A yeast one-hybrid assay showed that the C-end of SKL1 has transcription activity, as do the C-ends of SK1 and SK2. Our results suggested that SKLs were derived via gene duplication, but were not expressed and pseudogenized in normal cultivated rice during sequence evolution.

Genetic Divergence of Lineage-Specific Tandemly Duplicated Gene Clusters in Four Diploid Potato Genotypes

Potato (Solanum tuberosum L.) is the most important non-grain food crop. Tandem duplication significantly contributes to genome evolution. The objectives of this study were to (i) identify tandemly duplicated genes and compare their genomic distributions across potato genotypes, (ii) investigate the bias in functional specificities, (iii) explore the relationships among coding sequence, promoter and expression divergences associated with tandemly duplicated genes, (iv) examine the role of tandem duplication in generating and expanding lineage-specific gene families, (v) investigate the evolutionary forces affecting tandemly duplicated genes, and (vi) assess the similarities and differences with respect to above mentioned aspects between cultivated genotypes and their wild-relative. In this study, we used well-annotated and chromosome-scale de novo genome assemblies of multiple potato genotypes. Our results showed that tandemly duplicated genes are abundant and dispersed through the genome. We found that several functional specificities, such as disease resistance, stress-tolerance, and biosynthetic pathways of tandemly duplicated genes were differentially enriched across multiple potato genomes. Our results indicated the existence of a significant correlation among expression, promoter, and protein divergences in tandemly duplicated genes. We found about one fourth of tandemly duplicated gene clusters as lineage-specific among multiple potato genomes, and these tended to localize toward centromeres and revealed distinct selection signatures and expression patterns. Furthermore, our results showed that a majority of duplicated genes were retained through sub-functionalization followed by genetic redundancy, while only a small fraction of duplicated genes was retained though neo-functionalization. The lineage-specific expansion of gene families by tandem duplication coupled with functional bias might have significantly contributed to potato's genotypic diversity, and, thus, to adaption to environmental stimuli.

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H₂O₂ caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (As^III and As^V) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of As^III and As^V treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Evolutionary history of the C-repeat binding factor/dehydration-responsive element-binding 1 (CBF/DREB1) protein family in 43 plant species and characterization of CBF/DREB1 proteins in Solanum tuberosum

Background

Plants are easily affected by temperature variations, and high temperature (heat stress) and low temperature (cold stress) will lead to poor plant development and reduce crop yields. Therefore, it is very important to identify resistance genes for improving the ability of plants to resist heat stress or cold stress by using modern biotechnology. Members of the C-repeat binding factor/Dehydration responsive element-binding 1 (CBF/DREB1) protein family are related to the stress resistance of many plant species. These proteins affect the growth and development of plants and play vital roles during environmental stress (cold, heat, drought, salt, etc.). In this study, we identified CBF/DREB1 genes from 43 plant species (including algae, moss, ferns, gymnosperms, angiosperms) by using bioinformatic methods to clarify the characteristics of the CBF/DREB1 protein family members and their functions in potato under heat and cold stresses.

Results

In this study, we identified 292 CBF/DREB1 proteins from 43 plant species. However, no CBF/DREB1 protein was found in algae, moss, ferns, or gymnosperms; members of this protein family exist only in angiosperms. Phylogenetic analysis of all the CBF/DREB1 proteins revealed five independent groups. Among them, the genes of group I do not exist in eudicots and are found only in monocots, indicating that these genes have a special effect on monocots. The analysis of motifs, gene duplication events, and the expression data from the PGSC website revealed the gene structures, evolutionary relationships, and expression patterns of the CBF/DREB1 proteins. In addition, analysis of the transcript levels of the 8 CBF/DREB1 genes in potato (Solanum tuberosum) under low-temperature and high-temperature stresses showed that these genes were related to temperature stresses. In particular, the expression levels of StCBF3 and StCBF4 in the leaves, stems, and roots significantly increased under high-temperature conditions, which suggested that StCBF3 and StCBF4 may be closely related to heat tolerance in potato.

Conclusion

Overall, members of the CBF/DREB1 protein family exist only in angiosperms and plays an important role in the growth and development of plants. In addition, the CBF/DREB1 protein family is related to the heat and cold resistance of potato. Our research revealed the evolution of the CBF/DREB1 family, and is useful for studying the precise functions of the CBF/DREB1 proteins when the plants are developing and are under temperature stress.

Evolutionary history of the heat shock protein 90 (Hsp90) family of 43 plants and characterization of Hsp90s in Solanum tuberosum

Heat shock protein 90 genes/proteins (Hsp90s) are related to the stress resistance found in various plant species. These proteins affect the growth and development of plants and have important effects on the plants under various stresses (cold, drought and salt) in the environment. In this study, we identified 334 Hsp90s from 43 plant species, and Hsp90s were found in all species. Phylogenetic tree and conserved domain database analysis of all Hsp90s showed three independent clades. The analysis of motifs, gene duplication events, and the expression data from PGSC website revealed the gene structures, evolution relationships, and expression patterns of the Hsp90s. In addition, analysis of the transcript levels of the 7 Hsp90s in potato (Solanum tuberosum) under low temperature and high temperature stresses showed that these genes were related to the temperature stresses. Especially StHsp90.2 and StHsp90.4, under high or low temperature conditions, the expression levels in leaves, stems, or roots were significantly up-regulated. Our findings revealed the evolution of the Hsp90s, which had guiding significance for further researching the precise functions of the Hsp90s.

A Comprehensive Survey on the Terpene Synthase Gene Family Provides New Insight into Its Evolutionary Patterns

Terpenes are organic compounds and play important roles in plant growth and development as well as in mediating interactions of plants with the environment. Terpene synthases (TPSs) are the key enzymes responsible for the biosynthesis of terpenes. Although some species were employed for the genome-wide identification and characterization of the TPS family, limited information is available regarding the evolution, expansion, and retention mechanisms occurring in this gene family. We performed a genome-wide identification of the TPS family members in 50 sequenced genomes. Additionally, we also characterized the TPS family from aromatic spearmint and basil plants using RNA-Seq data. No TPSs were identified in algae genomes but the remaining plant species encoded various numbers of the family members ranging from 2 to 79 full-length TPSs. Some species showed lineage-specific expansion of certain subfamilies, which might have contributed toward species or ecotype divergence or environmental adaptation. A large-scale family expansion was observed mainly in dicot and monocot plants, which was accompanied by frequent domain loss. Both tandem and segmental duplication significantly contributed toward family expansion and expression divergence and played important roles in the survival of these expanded genes. Our data provide new insight into the TPS family expansion and evolution and suggest that TPSs might have originated from isoprenyl diphosphate synthase genes.

Genome-wide identification, characterisation and functional evaluation of WRKY genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.) G. Don. under abiotic stresses

BACKGROUND

WRKY DNA-binding protein (WRKY) is a large gene family involved in plant responses and adaptation to salt, drought, cold and heat stresses. Sweet potato from the genus Ipomoea is a staple food crop, but the WRKY genes in Ipomoea species remain unknown to date. Hence, we carried out a genome-wide analysis of WRKYs in Ipomoea trifida (H.B.K.) G. Don., the wild ancestor of sweet potato.

RESULTS

A total of 83 WRKY genes encoding 96 proteins were identified in I. trifida, and their gene distribution, duplication, structure, phylogeny and expression patterns were studied. ItfWRKYs were distributed on 15 chromosomes of I. trifida. Gene duplication analysis showed that segmental duplication played an important role in the WRKY gene family expansion in I. trifida. Gene structure analysis showed that the intron-exon model of the ItfWRKY gene was highly conserved. Meanwhile, the ItfWRKYs were divided into five groups (I, IIa + IIb, IIc, IId + IIe and III) on the basis of the phylogenetic analysis on I. trifida and Arabidopsis thaliana WRKY proteins. In addition, gene expression profiles confirmed by quantitative polymerase chain reaction showed that ItfWRKYs were highly up-regulated or down-regulated under salt, drought, cold and heat stress conditions, implying that these genes play important roles in response and adaptation to abiotic stresses.

CONCLUSIONS

In summary, genome-wide identification, gene structure, phylogeny and expression analysis of WRKY gene in I. trifida provide basic information for further functional studies of ItfWRKYs and for the molecular breeding of sweet potato.

Genome Wide Identification, Molecular Characterization, and Gene Expression Analyses of Grapevine NHX Antiporters Suggest Their Involvement in Growth, Ripening, Seed Dormancy, and Stress Response

Plant NHX antiporters are critical for cellular pH, Na⁺, and K⁺ homeostasis and salt tolerance. Even though their genomic and functional studies have been conducted in many species, the grapevine NHX family has not been described yet. Our work highlights the presence of six VvNHX genes whose phylogenetic analysis revealed their classification in two distinct groups: group I vacuolar (VvNHX1-5) and group II endosomal (VvNHX6). Several cis-acting regulatory elements related to tissue-specific expression, transcription factor binding, abiotic/biotic stresses response, and light regulation elements were identified in their promoter. Expression profile analyses of VvNHX genes showed variable transcription within organs and tissues with diverse patterns according to biochemical, environmental, and biotic treatments. All VvNHXs are involved in berry growth, except VvNHX5 that seems to be rather implicated in seed maturation. VvNHX4 would be more involved in floral development, while VvNHX2 and 3 display redundant roles. QPCR expression analyses of VvNHX1 showed its induction by NaCl and KNO₃ treatments, whereas VvNHX6 was induced by ABA application and strongly repressed by PEG treatment. VvNHX1 plays a crucial role in a bunch of grape developmental steps and adaptation responses through mechanisms of phyto-hormonal signaling. Overall, VvNHX family members could be valuable candidate genes for grapevine improvement.

CLAVATA signaling pathway genes modulating flowering time and flower number in chickpea

A combinatorial genomic strategy delineated functionally relevant natural allele of a CLAVATA gene and its marker (haplotype)-assisted introgression led to development of the early-flowering chickpea cultivars with high flower number and enhanced yield/productivity. Unraveling the genetic components involved in CLAVATA (CLV) signaling is crucial for modulating important shoot apical meristem (SAM) characteristics and ultimately regulating diverse SAM-regulated agromorphological traits in crop plants. A genome-wide scan identified 142 CLV1-, 28 CLV2- and 6 CLV3-like genes, and their comprehensive genomic constitution and phylogenetic relationships were deciphered in chickpea. The QTL/fine mapping and map-based cloning integrated with high-resolution association analysis identified SNP loci from CaCLV3_01 gene within a major CaqDTF1.1/CaqFN1.1 QTL associated with DTF (days to 50% flowering) and FN (flower number) traits in chickpea, which was further ascertained by quantitative expression profiling. Molecular haplotyping of CaCLV3_01 gene, expressed specifically in SAM, constituted two major haplotypes that differentiated the early-DTF and high-FN chickpea accessions from late-DTF and low-FN. Enhanced accumulation of transcripts of superior CaCLV3_01 gene haplotype and known flowering promoting genes was observed in the corresponding haplotype-introgressed early-DTF and high-FN near-isogenic lines (NILs) with narrow SAM width. The superior haplotype-introgressed NILs exhibited early-flowering, high-FN and enhanced seed yield/productivity without compromising agronomic performance. These delineated molecular signatures can regulate DTF and FN traits through SAM proliferation and differentiation and thereby will be useful for translational genomic study to develop early-flowering cultivars with enhanced yield/productivity.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

Characterisation of pathogen-specific regions and novel effector candidates in Fusarium oxysporum f. sp. cepae

A reference-quality assembly of Fusarium oxysporum f. sp. cepae (Foc), the causative agent of onion basal rot has been generated along with genomes of additional pathogenic and non-pathogenic isolates of onion. Phylogenetic analysis confirmed a single origin of the Foc pathogenic lineage. Genome alignments with other F. oxysporum ff. spp. and non pathogens revealed high levels of syntenic conservation of core chromosomes but little synteny between lineage specific (LS) chromosomes. Four LS contigs in Foc totaling 3.9 Mb were designated as pathogen-specific (PS). A two-fold increase in segmental duplication events was observed between LS regions of the genome compared to within core regions or from LS regions to the core. RNA-seq expression studies identified candidate effectors expressed in planta, consisting of both known effector homologs and novel candidates. FTF1 and a subset of other transcription factors implicated in regulation of effector expression were found to be expressed in planta.

Plant transcriptomics and responses to environmental stress: an overview

Different stresses include nutrient deficiency, pathogen attack, exposure to toxic chemicals etc. Transcriptomic studies have been mainly applied to only a few plant species including the model plant, Arabidopsis thaliana. These studies have provided valuable insights into the genetic networks of plant stress responses. Transcriptomics applied to cash crops including barley, rice, sugarcane, wheat and maize have further helped in understanding physiological and molecular responses in terms of genome sequence, gene regulation, gene differentiation, posttranscriptional modifications and gene splicing. On the other hand, comparative transcriptomics has provided more information about plant's response to diverse stresses. Thus, transcriptomics, together with other biotechnological approaches helps in development of stress tolerance in crops against the climate change.

Comparative transcriptome resources of two Dysosma species (Berberidaceae) and molecular evolution of the CYP719A gene in Podophylloideae

Dysosma species (Berberidaceae, Podophylloideae) are of great medicinal pharmacogenetic importance and used as model systems to study the drivers and mechanisms of species diversification of temperate plants in East Asia. Recently, we have sequenced the transcriptome of the low-elevation D. versipellis. In this study, we sequenced the transcriptome of the high-elevation D. aurantiocaulis and used comparative genomic approaches to investigate the transcriptome evolution of the two species. We retrieved 53,929 unigenes from D. aurantiocaulis by de novo transcriptome assemblies using the Illumina HiSeq 2000 platform. Comparing the transcriptomes of both species, we identified 4593 orthologs. Estimation of Ka/Ks ratios for 3126 orthologs revealed that none had a Ka/Ks significantly greater than 1, whereas 1273 (Ka/Ks < 0.5, P < 0.05) were inferred to be under purifying selection. A total of 51 primer pairs were successfully designed from 461 EST-SSRs contained in 4593 orthologs. Marker validation assay revealed that 26 (51%) and 41 (80.4%) produced clear fragments with the expected sizes in all Podophylloideae species. Specifically, 19 different sequences of CYP719A were identified from PCR-amplified genomic DNA of all 12 species of Podophylloideae using primers designed from the assembled transcripts. The data further indicated that CYP719A was likely subject to strong selective constraints maintaining only one copy per genome. In Dysosma, there was relaxed purifying selection or more positive selection for high-elevation species. Overall, this study has generated a wealth of molecular resources potentially useful for pharmacogenetic and evolutionary studies in Dysosma and allied taxa.

Comprehensive Expression Profiling of Rice Tetraspanin Genes Reveals Diverse Roles During Development and Abiotic Stress

Tetraspanin family is comprised of evolutionarily conserved integral membrane proteins. The incredible ability of tetraspanins to form 'micro domain complexes' and their preferential targeting to membranes emphasizes their active association with signal recognition and communication with neighboring cells, thus acting as key modulators of signaling cascades. In animals, tetraspanins are associated with multitude of cellular processes. Unlike animals, the biological relevance of tetraspanins in plants has not been well investigated. In Arabidopsis tetraspanins are known to contribute in important plant development processes such as leaf morphogenesis, root, and floral organ formation. In the present study we investigated the genomic organization, chromosomal distribution, phylogeny and domain structure of 15 rice tetraspanin proteins (OsTETs). OsTET proteins had similar domain structure and signature 'GCCK/R' motif as reported in Arabidopsis. Comprehensive expression profiling of OsTET genes suggested their possible involvement during rice development. While OsTET9 and 10 accumulated predominantly in flowers, OsTET5, 8, and 12 were preferentially expressed in root tissues. Noticeably, seven OsTETs exhibited more than twofold up regulation at early stages of flag leaf senescence in rice. Furthermore, several OsTETs were differentially regulated in rice seedlings exposed to abiotic stresses, exogenous treatment of hormones and nutrient deprivation. Transient subcellular localization studies of eight OsTET proteins in tobacco epidermal cells showed that these proteins localized in plasma membrane. The present study provides valuable insights into the possible roles of tetraspanins in regulating development and defining response to abiotic stresses in rice. Targeted proteomic studies would be useful in identification of their interacting partners under different conditions and ultimately their biological function in plants.

JAZ Repressors: Potential Involvement in Nutrients Deficiency Response in Rice and Chickpea

Jasmonates (JA) are well-known phytohormones which play important roles in plant development and defense against pathogens. Jasmonate ZIM domain (JAZ) proteins are plant-specific proteins and act as transcriptional repressors of JA-responsive genes. JA regulates both biotic and abiotic stress responses in plants; however, its role in nutrient deficiency responses is very elusive. Although, JA is well-known for root growth inhibition, little is known about behavior of JAZ genes in response to nutrient deficiencies, under which root architectural alteration is an important adaptation. Using protein sequence homology and a conserved-domains approach, here we identify 10 novel JAZ genes from the recently sequenced Chickpea genome, which is one of the most nutrient efficient crops. Both rice and chickpea JAZ genes express in tissue- and stimuli-specific manners. Many of which are preferentially expressed in root. Our analysis further showed differential expression of JAZ genes under macro (NPK) and micronutrients (Zn, Fe) deficiency in rice and chickpea roots. While both rice and chickpea JAZ genes showed a certain level of specificity toward type of nutrient deficiency, generally majority of them showed induction under K deficiency. Generally, JAZ genes showed an induction at early stages of stress and expression declined at later stages of macro-nutrient deficiency. Our results suggest that JAZ genes might play a role in early nutrient deficiency response both in monocot and dicot roots, and information generated here can be further used for understanding the possible roles of JA in root architectural alterations for nutrient deficiency adaptations.

Evolution and expression analysis of the soybean glutamate decarboxylase gene family

Glutamate decarboxylase (GAD) is an enzyme that catalyses the conversion of L-glutamate into gamma-aminobutyric acid (GABA), which is a four-carbon non-protein amino acid present in all organisms. Although plant GAD plays important roles in GABA biosynthesis, our knowledge concerning GAD gene family members and their evolutionary relationship remains limited. Therefore, in this study, we have analysed the evolutionary mechanisms of soybean GAD genes and suggested that these genes expanded in the soybean genome partly due to segmental duplication events. The approximate dates of duplication events were calculated using the synonymous substitution rate, and we suggested that the segmental duplication of GAD genes in soybean originated 9.47 to 11.84 million years ago (Mya). In addition, all segmental duplication pairs (GmGAD1/3 and GmGAD2/4) are subject to purifying selection. Furthermore, GmGAD genes displayed differential expression either in their transcript abundance or in their expression patterns under abiotic stress conditions like salt, drought, and cold. The expression pattern of paralogous pairs suggested that they might have undergone neofunctionalization during the subsequent evolution process. Taken together, our results provide valuable information for the evolution of the GAD gene family and represent the basis for future research on the functional characterization of GAD genes in higher plants.

Evolution and expression analysis of the soybean glutamate decarboxylase gene family

Glutamate decarboxylase (GAD) is an enzyme that catalyses the conversion of L-glutamate into gamma-aminobutyric acid (GABA), which is a four-carbon non-protein amino acid present in all organisms. Although plant GAD plays important roles in GABA biosynthesis, our knowledge concerning GAD gene family members and their evolutionary relationship remains limited. Therefore, in this study, we have analysed the evolutionary mechanisms of soybean GAD genes and suggested that these genes expanded in the soybean genome partly due to segmental duplication events. The approximate dates of duplication events were calculated using the synonymous substitution rate, and we suggested that the segmental duplication of GAD genes in soybean originated 9.47 to 11.84 million years ago (Mya). In addition, all segmental duplication pairs (GmGAD1/3 and GmGAD2/4) are subject to purifying selection. Furthermore, GmGAD genes displayed differential expression either in their transcript abundance or in their expression patterns under abiotic stress conditions like salt, drought, and cold. The expression pattern of paralogous pairs suggested that they might have undergone neofunctionalization during the subsequent evolution process. Taken together, our results provide valuable information for the evolution of the GAD gene family and represent the basis for future research on the functional characterization of GAD genes in higher plants.

Genome-wide survey on genomic variation, expression divergence, and evolution in two contrasting rice genotypes under high salinity stress

Expression profiling is one of the most important tools for dissecting biological functions of genes and the upregulation or downregulation of gene expression is sufficient for recreating phenotypic differences. Expression divergence of genes significantly contributes to phenotypic variations. However, little is known on the molecular basis of expression divergence and evolution among rice genotypes with contrasting phenotypes. In this study, we have implemented an integrative approach using bioinformatics and experimental analyses to provide insights into genomic variation, expression divergence, and evolution between salinity-sensitive rice variety Nipponbare and tolerant rice line Pokkali under normal and high salinity stress conditions. We have detected thousands of differentially expressed genes between these two genotypes and thousands of up- or downregulated genes under high salinity stress. Many genes were first detected with expression evidence using custom microarray analysis. Some gene families were preferentially regulated by high salinity stress and might play key roles in stress-responsive biological processes. Genomic variations in promoter regions resulted from single nucleotide polymorphisms, indels (1–10 bp of insertion/deletion), and structural variations significantly contributed to the expression divergence and regulation. Our data also showed that tandem and segmental duplication, CACTA and hAT elements played roles in the evolution of gene expression divergence and regulation between these two contrasting genotypes under normal or high salinity stress conditions.

A strategy for genome-wide identification of gene based polymorphisms in rice reveals non-synonymous variation and functional genotypic markers

The genetic diversity of plants has traditionally been employed to improve crop plants to suit human needs, and in the future feed the increasing population and protect crops from environmental stresses and climate change. Genome-wide sequencing is a reality and can be used to make association to crop traits to be utilized by high-throughput marker based selection methods. This study describes a strategy of using next generation sequencing (NGS) data from the rice genome to make comparisons to the high-quality reference genome, identify functional polymorphisms within genes that might result in function changes and be used to study correlations to traits and employed in genetic mapping. We analyzed the NGS data of Oryza sativa ssp indica cv. G4 covering 241 Mb with ∼20X coverage and compared to the reference genome of Oryza sativa ssp. japonica to describe the genome-wide distribution of gene-based single nucleotide polymorphisms (SNPs). The analysis shows that the 63% covered genome consists of 1.6 million SNPs with 6.9 SNPs/Kb, and including 80,146 insertions and 92,655 deletions (INDELs) genome-wide. There are a total of 1,139,801 intergenic SNPs, 295,136 SNPs in intronic/non-coding regions, 195,098 in coding regions, 23,242 SNPs at the five-prime (5′) UTR regions and 22,686 SNPs at the three-prime (3′) UTR region. SNP variation was found in 40,761 gene loci, which include 75,262 synonymous and 119,836 non-synonymous changes, and functional reading frame changes through 3,886 inducing STOP-codon (isSNP) and 729 preventing STOP-codon (psSNP) variation. There are quickly evolving 194 high SNP hotspot genes (>100 SNPs/gene), and 1,513 out of 2,458 transcription factors displaying 2,294 non-synonymous SNPs that can be a major source of phenotypic diversity within the species. All data is searchable at https://plantstress-pereira.uark.edu/oryza2. We envision that this strategy will be useful for the identification of genes for crop traits and molecular breeding of rice cultivars.

Evolution of the tyrosinase gene family in bivalve molluscs: independent expansion of the mantle gene repertoire

Tyrosinase is a copper-containing enzyme that mediates the hydroxylation of monophenols and oxidation of o-diphenols to o-quinones. This enzyme is involved in a variety of biological processes, including pigment production, innate immunity, wound healing, and exoskeleton fabrication and hardening (e.g. arthropod skeleton and mollusc shell). Here we show that the tyrosinase gene family has undergone large expansions in pearl oysters (Pinctada spp.) and the Pacific oyster (Crassostrea gigas). Phylogenetic analysis reveals that pearl oysters possess at least four tyrosinase genes that are not present in the Pacific oyster. Likewise, C. gigas has multiple tyrosinase genes that are not orthologous to the Pinctada genes, indicating that this gene family has expanded independently in these bivalve lineages. Many of the tyrosinase genes in these bivalves are expressed at relatively high levels in the mantle, the organ responsible for shell fabrication. Detailed comparisons of tyrosinase gene expression in different regions of the mantle in two closely related pearl oysters, P. maxima and P. margaritifera, reveals that recently evolved orthologous tyrosinase genes can have markedly different expression profiles. The expansion of tyrosinase genes in these oysters and their co-option into the mantle's gene regulatory network is consistent with mollusc shell formation being underpinned by a rapidly evolving transcriptome.

Impact of recurrent gene duplication on adaptation of plant genomes

BACKGROUND

Recurrent gene duplication and retention played an important role in angiosperm genome evolution. It has been hypothesized that these processes contribute significantly to plant adaptation but so far this hypothesis has not been tested at the genome scale.

RESULTS

We studied available sequenced angiosperm genomes to assess the frequency of positive selection footprints in lineage specific expanded (LSE) gene families compared to single-copy genes using a dN/dS-based test in a phylogenetic framework. We found 5.38% of alignments in LSE genes with codons under positive selection. In contrast, we found no evidence for codons under positive selection in the single-copy reference set. An analysis at the branch level shows that purifying selection acted more strongly on single-copy genes than on LSE gene clusters. Moreover we detect significantly more branches indicating evolution under positive selection and/or relaxed constraint in LSE genes than in single-copy genes.

CONCLUSIONS

In this - to our knowledge -first genome-scale study we provide strong empirical support for the hypothesis that LSE genes fuel adaptation in angiosperms. Our conservative approach for detecting selection footprints as well as our results can be of interest for further studies on (plant) gene family evolution.

The large soybean (Glycine max) WRKY TF family expanded by segmental duplication events and subsequent divergent selection among subgroups

BACKGROUND

WRKY genes encode one of the most abundant groups of transcription factors in higher plants, and its members regulate important biological process such as growth, development, and responses to biotic and abiotic stresses. Although the soybean genome sequence has been published, functional studies on soybean genes still lag behind those of other species.

RESULTS

We identified a total of 133 WRKY members in the soybean genome. According to structural features of their encoded proteins and to the phylogenetic tree, the soybean WRKY family could be classified into three groups (groups I, II, and III). A majority of WRKY genes (76.7%; 102 of 133) were segmentally duplicated and 13.5% (18 of 133) of the genes were tandemly duplicated. This pattern was not apparent in Arabidopsis or rice. The transcriptome atlas revealed notable differential expression in either transcript abundance or in expression patterns under normal growth conditions, which indicated wide functional divergence in this family. Furthermore, some critical amino acids were detected using DIVERGE v2.0 in specific comparisons, suggesting that these sites have contributed to functional divergence among groups or subgroups. In addition, site model and branch-site model analyses of positive Darwinian selection (PDS) showed that different selection regimes could have affected the evolution of these groups. Sites with high probabilities of having been under PDS were found in groups I, II c, II e, and III. Together, these results contribute to a detailed understanding of the molecular evolution of the WRKY gene family in soybean.

CONCLUSIONS

In this work, all the WRKY genes, which were generated mainly through segmental duplication, were identified in the soybean genome. Moreover, differential expression and functional divergence of the duplicated WRKY genes were two major features of this family throughout their evolutionary history. Positive selection analysis revealed that the different groups have different evolutionary rates. Together, these results contribute to a detailed understanding of the molecular evolution of the WRKY gene family in soybean.