PTGBase: an integrated database to study tandem duplicated genes in plants

Fine mapping and functional characterization of SiLAC3 in regulating brown seed coat pigmentation, integrated transcriptomic and metabolomic analyses in sesame (Sesamum indicum L.)

The color of the seed coat in sesame is a critical breeding trait, with significant implications for seed quality and consumer preferences. However, the genetic foundation and regulatory mechanisms underlying seed coat colors remain largely unexplored. Earlier, we identified a quantitative trait locus (QTL), qBSCchr6, associated with brown-colored seed coats in sesame. Currently, we utilized a recombinant inbred line (RIL) G3010, which harbors qBSCchr6 and crossed it with the white-seeded cultivar (YZ8) to develop an F2 population for QTL fine mapping. Through detailed recombinant analysis, qBSCchr6 was fine mapped to a narrow 83.4 kb genomic interval, co-segregating with the molecular marker BSC_SNV23. Within this interval, we identified SIN_1023239 for the brown-colored seed coat, which we have designated as SiLAC3. SiLAC3 encodes a laccase enzyme with three putative multicopper oxidase domains essential for its function. Sequence analysis revealed that a single nucleotide variant in the multicopper oxidase_3 domain is critical for the expression of brown seed coat phenotype. To further elucidate the role of SiLAC3, we employed a combination of chemical staining, transcriptomics, and metabolomics analyses. Our results indicate that the lignin biosynthesis pathway is crucial for brown seed coat pigmentation, with SiLAC3 facilitating oxidative polymerization of coniferyl alcohol in the seed coat. These findings advance our understanding of genetic and molecular mechanisms underlying seed coat color in sesame and offer a foundation to improve this trait in breeding programs.

Multiplex Editing of the Nucleoredoxin1 Tandem Array in Poplar: From Small Indels to Translocations and Complex Inversions

The CRISPR-Cas9 system has been deployed for precision mutagenesis in an ever-growing number of species, including agricultural crops and forest trees. Its application to closely linked genes with extremely high sequence similarities has been less explored. In this study, we used CRISPR-Cas9 to mutagenize a tandem array of seven Nucleoredoxin1 (NRX1) genes spanning ∼100 kb in Populus tremula × Populus alba. We demonstrated efficient multiplex editing with one single guide RNA in 42 transgenic lines. The mutation profiles ranged from small insertions and deletions and local deletions in individual genes to large genomic dropouts and rearrangements spanning tandem genes. We also detected complex rearrangements including translocations and inversions resulting from multiple cleavage and repair events. Target capture sequencing was instrumental for unbiased assessments of repair outcomes to reconstruct unusual mutant alleles. The work highlights the power of CRISPR-Cas9 for multiplex editing of tandemly duplicated genes to generate diverse mutants with structural and copy number variations to aid future functional characterization.

Combined genomic and transcriptomic analysis reveals the contribution of tandem duplication genes to low-temperature adaptation in perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is an agronomically important cool-season grass species that is widely used as forage for ruminant animal production and cultivated in temperate regions for the establishment of lawns. However, the underlying genetic mechanism of the response of L. perenne to low temperature is still unclear. In the present study, we performed a comprehensive study and identified 3,770 tandem duplication genes (TDGs) in L. perenne, and evolutionary analysis revealed that L. perenne might have undergone a duplication event approximately 7.69 Mya. GO and KEGG pathway functional analyses revealed that these TDGs were mainly enriched in photosynthesis, hormone-mediated signaling pathways and responses to various stresses, suggesting that TDGs contribute to the environmental adaptability of L. perenne. In addition, the expression profile analysis revealed that the expression levels of TDGs were highly conserved and significantly lower than those of all genes in different tissues, while the frequency of differentially expressed genes (DEGs) from TDGs was much higher than that of DEGs from all genes in response to low-temperature stress. Finally, in-depth analysis of the important and expanded gene family indicated that the members of the ELIP subfamily could rapidly respond to low temperature and persistently maintain higher expression levels during all low temperature stress time points, suggesting that ELIPs most likely mediate low temperature responses and help to facilitate adaptation to low temperature in L. perenne. Our results provide evidence for the genetic underpinning of low-temperature adaptation and valuable resources for practical application and genetic improvement for stress resistance in L. perenne.

Genome-wide identification and expression analysis of the Pisum sativum (L.) APETALA2/ethylene-responsive factor (AP2/ERF) gene family reveals functions in drought and cold stresses

AP2/ERF (APETALA2/Ethylene Response Factor) is a family of transcription factors that play essential roles in regulating gene expression in response to various environmental stimuli, including biotic and abiotic stresses, hormone signaling, and developmental processes. Pisum sativum (L.), commonly known as garden pea, is a winter crop sensitive to high temperatures and can also be affected by extreme cold and drought conditions. This study performed a genome-wide analysis of AP2/ERF genes and identified 153 AP2/ERF genes in P. sativum. Based on the conserved AP2/ERF domain and sequence homology, they were classified into AP2 (APETALA2), ERF (Ethylene Response Factor), DREB (Dehydration responsive element-binding), RAV (Related to Abscisic Acid Insensitive 3/ Viviparous 1) and Soloist subfamily. The DREB and ERF subfamily were further divided into groups A1-6 and B1-B6. Tandem and segmental duplication events were more frequent in the ERF subfamily, which can have important implications for their evolution and functional diversification. Under cold stress, the expression of DREB1A was highly induced in leaves, whereas DREB1B was suppressed. Similarly, the DREB2A, DREB2C, DREB2E, and DREB2F were induced in leaves under drought stress. The putative target genes of AP2/ERF transcription factors are highly diversified, suggesting that they play essential roles in various physiological responses in plants, including responses to biotic and abiotic stresses as well as developmental processes. Thus, this study of AP2/ERF genes and their functions provides valuable insight into how P. sativum responds to different environmental conditions, including cold and drought stresses.

Discovering favorable genes, QTLs, and genotypes as a genetic resource for sesame (Sesamum indicum L.) improvement

Sesame (Sesamum indicum L.) is an ancient diploid oilseed crop with high oil content, quality protein, and antioxidant characteristics that is produced in many countries worldwide. The genes, QTLs, and genetic resources of sesame are utilized by sesame researchers and growers. Researchers have identified the many useful traits of this crop, which are available on different platforms. The genes, genotypes, QTLs, and other genetic diversity data of sesame have been collected and stored in more than nine genomic resources, and five sesame crop marker databases are available online. However, data on phenotypic and genotypic variability, which would contribute to sesame improvements, are limited and not yet accessible. The present study comprehensively reviewed more than 110 original published research papers and scientifically incorporated the results. The candidate genes, genotypes, and QTLs of significantly important traits of sesame were identified. Genetic resources related to grain yield and yield component traits, oil content and quality, drought tolerance, salt tolerance, waterlogging resistance, disease resistance, mineral nutrient, capsule shattering resistance, and other agronomic important traits of sesame were studied. Numerous candidate genotypes, genes, QTLs, and alleles associated with those traits were summarized and discovered. The chromosome regions and linkage groups, maps associated with the best traits, and candidate genes were also included. The variability presented in this paper combined with sesame genetic information will help inform further sesame improvement.

An overview of online resources for intra-species detection of gene duplications

Gene duplication plays an important role in evolutionary mechanism, which can act as a new source of genetic material in genome evolution. However, detecting duplicate genes from genomic data can be challenging. Various bioinformatics resources have been developed to identify duplicate genes from single and/or multiple species. Here, we summarize the metrics used to measure sequence identity among gene duplicates within species, compare several computational approaches that have been used to predict gene duplicates, and review recent advancements of a Basic Local Alignment Search Tool (BLAST)-based web tool and database, allowing future researchers to easily identify intra-species gene duplications. This article is a quick reference guide for research tools used for detecting gene duplicates.

Genome-wide analysis of tandem duplicated genes and their expression under salt stress in seashore paspalum

Seashore paspalum (Paspalum vaginatum) is a halophytic, warm-season grass which is closely related to various grain crops. Gene duplication plays an important role in plant evolution, conferring significant plant adaptation at the genomic level. Here, we identified 2,542 tandem duplicated genes (TDGs) in the P. vaginatum genome and estimated the divergence time of pairs of TDGs based on synonymous substitution rates (Ks). Expression of P. vaginatum TDGs resulted in enrichment in many GO terms and KEGG pathways when compared to four other closely-related species. The GO terms included: "ion transmembrane transporter activity," "anion transmembrane transporter activity" and "cation transmembrane transport," and KEGG pathways included "ABC transport." RNA-seq analysis of TDGs showed tissue-specific expression under salt stress, and we speculated that P. vaginatum leaves became adapted to salt stress in the earlier whole-genome duplication (WGD; ~83.3 million years ago; Ma), whereas the entire P. vaginatum plant acquired a large number of TDGs related to salt stress in the second WGD (~23.3 Ma). These results can be used as a reference resource to accelerate salt-resistance research in other grasses and crops.

Comparative genomics analysis of bHLH genes in cucurbits identifies a novel gene regulating cucurbitacin biosynthesis

The basic helix-loop-helix (bHLH) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, and have undergone significant expansion during land plant evolution by gene duplications. In cucurbit crops, several bHLH genes have been found to be responsible for the agronomic traits such as bitterness. However, the characterization of bHLH genes across the genomes of cucurbit species has not been reported, and how they have evolved and diverged remains largely unanswered. Here we identified 1160 bHLH genes in seven cucurbit crops and performed a comprehensive comparative genomics analysis. We determined orthologous and paralogous bHLH genes across cucurbit crops by syntenic analysis between or within species. Orthology and phylogenetic analysis of the tandem-duplicated bHLH genes in the Bt cluster which regulate the biosynthesis of cucurbitacins suggest that this cluster is derived from three ancestral genes after the cucurbit-common tetraploidization event. Interestingly, we identified a new conserved cluster paralogous to the Bt cluster that includes two tandem bHLH genes, and the evolutionary history and expression profiles of these two genes in the new cluster suggest the involvement of one gene (Brp) in the regulation of cucurbitacin biosynthesis in roots. Further biochemical and transgenic assays in melon hairy roots support the function of Brp. This study provides useful information for further investigating the functions of bHLH TFs and novel insights into the regulation of cucurbitacin biosynthesis in cucurbit crops and other plants.

Global Identification and Characterization of C2 Domain-Containing Proteins Associated with Abiotic Stress Response in Rice (Oryza sativa L.)

C2 domain-containing proteins (C2DPs) have been identified in different genomes that contain single or multiple C2 domains in their C- or N-terminal. It possesses higher functional activity in the transmembrane regions. The identification of C2 domains were reported in a previous study, such as multiple C2 domains and transmembrane-region proteins (MCTPs) and N-terminal-TM-C2 domain proteins (NTMC2s) of rice, Arabidopsis thaliana, and cotton, whereas the C2DP gene family in rice has not been comprehensively studied, and the role of the C2DP gene in rice in response to abiotic stress is not yet fully understood. In this study, we identified 82 C2DPs in the rice genome and divided them into seven groups through phylogenetic analysis. The synteny analysis revealed that duplication events were either exhibited within the genome of rice or between the genomes of rice and other species. Through the analysis of cis-acting elements in promoters, expression profiles, and qRT-PCR results, the functions of OsC2DPs were found to be widely distributed in diverse tissues and were extensively involved in phytohormones-related and abiotic stresses response in rice. The prediction of the microRNA (miRNA) targets of OsC2DPs revealed the possibility of regulation by consistent miRNAs. Notably, OsC2DP50/51/52 as a co-tandem duplication exhibited similar expression variations and involved the coincident miRNA-regulation pathway. Moreover, the results of the genotypic variation and haplotype analysis revealed that OsC2DP17, OsC2DP29, and OsC2DP49 were associated with cold stress responses. These findings provided comprehensive insights for characterizations of OsC2DPs in rice as well as for their roles for abiotic stress.

Genome-Wide Identification of GDSL-Type Esterase/Lipase Gene Family in Dasypyrum villosum L. Reveals That DvGELP53 Is Related to BSMV Infection

GDSL-type esterase/lipase proteins (GELPs) characterized by a conserved GDSL motif at their N-terminus belong to the lipid hydrolysis enzyme superfamily. In plants, GELPs play an important role in plant growth, development and stress response. The studies of the identification and characterization of the GELP gene family in Triticeae have not been reported. In this study, 193 DvGELPs were identified in Dasypyrum villosum and classified into 11 groups (clade A–K) by means of phylogenetic analysis. Most DvGELPs contain only one GDSL domain, only four DvGELPs contain other domains besides the GDSL domain. Gene structure analysis indicated 35.2% DvGELP genes have four introns and five exons. In the promoter regions of the identified DvGELPs, we detected 4502 putative cis-elements, which were associated with plant hormones, plant growth, environmental stress and light responsiveness. Expression profiling revealed 36, 44 and 17 DvGELPs were highly expressed in the spike, the root and the grain, respectively. Further investigation of a root-specific expressing GELP, DvGELP53, indicated it was induced by a variety of biotic and abiotic stresses. The knockdown of DvGELP53 inhibited long-distance movement of BSMV in the tissue of D. villosum. This research provides a genome-wide glimpse of the D. villosum GELP genes and hints at the participation of DvGELP53 in the interaction between virus and plants.

Genome-Wide Comparative Analysis of Heat Shock Transcription Factors Provides Novel Insights for Evolutionary History and Expression Characterization in Cotton Diploid and Tetraploid Genomes

Heat shock transcription factors (HSFs) are involved in environmental stress response and plant development, such as heat stress and flowering development. According to the structural characteristics of the HSF gene family, HSF genes were classified into three major types (HSFA, HSFB, and HSFC) in plants. Using conserved domains of HSF genes, we identified 621 HSF genes among 13 cotton genomes, consisting of eight diploid and five tetraploid genomes. Phylogenetic analysis indicated that HSF genes among 13 cotton genomes were grouped into two different clusters: one cluster contained all HSF genes of HSFA and HSFC, and the other cluster contained all HSF genes of HSFB. Comparative analysis of HSF genes in Arabidopsis thaliana, Gossypium herbaceum (A1), Gossypium arboreum (A2), Gossypium raimondii (D5), and Gossypium hirsutum (AD1) genomes demonstrated that four HSF genes were inherited from a common ancestor, A0, of all existing cotton A genomes. Members of the HSF gene family in G. herbaceum (A1) genome indicated a significant loss compared with those in G. arboretum (A2) and G. hirsutum (AD1) A genomes. However, HSF genes in G. raimondii (D5) showed relative loss compared with those in G. hirsutum (AD1) D genome. Analysis of tandem duplication (TD) events of HSF genes revealed that protein-coding genes among different cotton genomes have experienced TD events, but only the two-gene tandem array was detected in Gossypium thurberi (D1) genome. The expression analysis of HSF genes in G. hirsutum (AD1) and Gossypium barbadense (AD2) genomes indicated that the expressed HSF genes were divided into two different groups, respectively, and the expressed HSF orthologous genes between the two genomes showed totally different expression patterns despite the implementation of the same abiotic stresses. This work will provide novel insights for the study of evolutionary history and expression characterization of HSF genes in different cotton genomes and a widespread application model for the study of HSF gene families in plants.

Genomic characterization and expression profiles of stress-associated proteins (SAPs) in castor bean (Ricinus communis)

Stress-associated proteins (SAPs) are known as response factors to multiple abiotic and biotic stresses in plants. However, the potential physiological and molecular functions of SAPs remain largely unclear. Castor bean (Ricinus communis L.) is one of the most economically valuable non-edible woody oilseed crops, able to be widely cultivated in marginal lands worldwide because of its broad adaptive capacity to soil and climate conditions. Whether SAPs in castor bean plays a key role in adapting diverse soil conditions and stresses remains unknown. In this study, we used the castor bean genome to identify and characterize nine castor bean SAP genes (RcSAP). Structural analysis showed that castor bean SAP gene structures and functional domain types vary greatly, differing in intron number, protein sequence, and functional domain type. Notably, the AN1-C2H2-C2H2 zinc finger domain within RcSAP9 has not been often observed in other plant families. High throughput RNA-seq data showed that castor bean SAP gene profiles varied among different tissues. In addition, castor bean SAP gene expression varied in response to different stresses, including salt, drought, heat, cold and ABA and MeJA, suggesting that the transcriptional regulation of castor bean SAP genes might operate independently of each other, and at least partially independent from ABA and MeJA signal pathways. Cis-element analyses for each castor bean SAP gene showed that no common cis-elements are shared across the nine castor bean SAP genes. Castor bean SAPs were localized to different regions of cells, including the cytoplasm, nucleus, and cytomembrane. This study provides a comprehensive profile of castor bean SAP genes that advances our understanding of their potential physiological and molecular functions in regulating growth and development and their responses to different abiotic stresses.

Influence of domestication on specialized metabolic pathways in fruit crops

During the process of plant domestication, the selection and traditional breeding for desired characters such as flavor, juiciness and nutritional value of fruits, probably have resulted in gain or loss of specialized metabolites contributing to these traits. Their appearance in fruits is likely due to the acquisition of novel and specialized metabolic pathways and their regulation, driven by systematic molecular evolutionary events facilitated by traditional breeding. Plants change their armory of specialized metabolism to adapt and survive in diverse ecosystems. This may occur through molecular evolutionary events, such as single nucleotide polymorphism, gene duplication and transposition, leading to convergent or divergent evolution of biosynthetic pathways producing such specialized metabolites. Breeding and selection for improved specific and desired traits (fruit size, color, taste, flavor, etc.) in fruit crops through conventional breeding approaches may further alter content and profile of specialized metabolites. Biosynthetic routes of these metabolites have been studied in various plants. Here, we explore the influence of plant domestication and breeding processes on the selection of biosynthetic pathways of favorable specialized metabolites in fruit crops. An orderly clustered arrangement of genes associated with their production is observed in many fruit crops. We further analyzed selection-based acquisition of specialized metabolic pathways comparing first the metabolic profiles and genes involved in their biosynthesis, followed by the genomic organization of such genes between wild and domesticated horticultural crops. Domestication of crop plants favored the acquisition and retention of metabolic pathways that enhanced the fruit value while eliminated those which produced toxic or unfavorable metabolites. Interestingly, unintentional reorganization of complex metabolic pathways by selection and traditional breeding processes has endowed us with flavorful, juicy and nutritionally rich fruits.

Genome-wide identification, and phylogenetic and expression profiling analyses, of XTH gene families in Brassica rapa L. and Brassica oleracea L

BACKGROUND

Xyloglucan endotransglucosylase/hydrolase genes (XTHs) are a multigene family and play key roles in regulating cell wall extensibility in plant growth and development. Brassica rapa and Brassica oleracea contain XTHs, but detailed identification and characterization of the XTH family in these species, and analysis of their tissue expression profiles, have not previously been carried out.

RESULTS

In this study, 53 and 38 XTH genes were identified in B. rapa and B. oleracea respectively, which contained some novel members not observed in previous studies. All XTHs of B. rapa, B. oleracea and Arabidopsis thaliana could be classified into three groups, Group I/II, III and the Early diverging group, based on phylogenetic relationships. Gene structures and motif patterns were similar within each group. All XTHs in this study contained two characteristic conserved domains (Glyco_hydro and XET_C). XTHs are located mainly in the cell wall but some are also located in the cytoplasm. Analyses of the mechanisms of gene family expansion revealed that whole-genome triplication (WGT) events and tandem duplication (TD) may have been the major mechanisms accounting for the expansion of the XTH gene family. Interestingly, TD genes all belonged to Group I/II, suggesting that TD was the main reason for the largest number of genes being in these groups. B. oleracea had lost more of the XTH genes, the conserved domain XET_C and the conserved active-site motif EXDXE compared with B. rapa, consistent with asymmetrical evolution between the two Brassica genomes. A majority of XTH genes exhibited different tissue-specific expression patterns based on RNA-seq data analyses. Moreover, there was differential expression of duplicated XTH genes in the two species, indicating that their functional differentiation occurred after B. rapa and B. oleracea diverged from a common ancestor.

CONCLUSIONS

We carried out the first systematic analysis of XTH gene families in B. rapa and B. oleracea. The results of this investigation can be used for reference in further studies on the functions of XTH genes and the evolution of this multigene family.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

Ab initio GO-based mining for non-tandem-duplicated functional clusters in three model plant diploid genomes

A functional Non-Tandem Duplicated Cluster (FNTDC) is a group of non-tandem-duplicated genes that are located closer than expected by mere chance and have a role in the same biological function. The identification of secondary-compounds–related FNTDC has gained increased interest in recent years, but little ab-initio attempts aiming to the identification of FNTDCs covering all biological functions, including primary metabolism compounds, have been carried out. We report an extensive FNTDC dataset accompanied by a detailed assessment on parameters used for genome scanning and their impact on FNTDC detection. We propose 70% identity and 70% alignment coverage as intermediate settings to exclude tandem duplicated genes and a dynamic scanning window of 24 genes. These settings were applied to rice, arabidopsis and grapevine genomes to call for FNTDCs. Besides the best-known secondary metabolism clusters, we identified many FNTDCs associated to primary metabolism ranging from macromolecules synthesis/editing, TOR signalling, ubiquitination, proton and electron transfer complexes. Using the intermediate FNTDC setting parameters (at P-value 1e^-6), 130, 70 and 140 candidate FNTDCs were called in rice, arabidopsis and grapevine, respectively, and 20 to 30% of GO tags associated to called FNTDC were common among the 3 genomes. The datasets developed along with this work provide a rich framework for pinpointing candidate FNTDCs reflecting all GO-BP tags covering both primary and secondary metabolism with large macromolecular complexes/metabolons as the most represented FNTDCs. Noteworthy, several FNTDCs are tagged with GOs referring to organelle-targeted multi-enzyme complex, a finding that suggest the migration of endosymbiont gene chunks towards nuclei could be at the basis of these class of candidate FNTDCs. Most FNTDC appear to have evolved prior of genome duplication events. More than one-third of genes interspersed/adjacent to called FNTDCs lacked any functional annotation; however, their co-localization may provide hints towards a candidate biological role.

Sustainable Sesame (Sesamum indicum L.) Production through Improved Technology: An Overview of Production, Challenges, and Opportunities in Myanmar

This paper aims to review the research achievements concerning sustainable sesame (Sesamum indicum L.) production and outlook on the production constraints and future perspectives for Myanmar sesame. Sesame is an economically and nutritionally important crop, and it is prized for oil. The global sesame market demand is rising with increasing health awareness. Meanwhile, there is high competition in the market among producing countries for an international trade. Smallholder farmers in developing countries cultivate sesame as a cash crop on marginal soils. The edible oilseed sectors currently face several challenges, including ones affecting sesame crops. For sustainable production of sesame, an integrated approach is needed to overcome these challenges and the critical limiting factors should be identified. In recent years, sesame genomic resources, including molecular markers, genetic maps, genome sequences, and online functional databases, are available for sesame genetic improvement programs. Since ancient times, sesame has been cultivated in Myanmar, but productivity is still lower than that of other sesame producing countries. Myanmar sesame production is limited by many factors, including production technology, research and development, etc. With integration of these genomic resources, crop production and protection techniques, postharvest practices, crop improvement programs, and capacity building will play a crucial role for improving sesame production in Myanmar.

Comprehensive phylogenomic analysis of ERF genes in sorghum provides clues to the evolution of gene functions and redundancy among gene family members

APETALA2/Ethylene-Responsive transcription factors (AP2/ERF), with their multifunctional roles in plant development, hormone signaling and stress tolerance, are important candidates for engineering crop plants. Here, we report identification and analysis of gene structure, phylogenetic distribution, expression, chromosomal localization and cis-acting promoter analysis of AP2/ERF genes in the C4 crop plant sorghum. We identified 158 ERF genes in sorghum with 52 of them encoding dehydration-responsive binding elements (DREB) while 106 code for ERF subfamily proteins. Phylogenetic analysis organized sorghum ERF proteins into 11 distinct groups exhibiting clade-specific expansion. About 68% ERF genes have paralogs indicating gene duplications as major cause of expansion of ERF family in sorghum. Analysis of spatiotemporal expression patterns using publicly available data revealed their tissue/genotype-preferential accumulation. In addition, 40 ERF genes exhibited differential accumulation in response to heat and/or drought stress. About 25% of the segmental gene pairs and eleven tandem duplicated genes exhibited high correlation (> 0.7) in their expression patterns indicating genetic redundancy. Comparative phylogenomic analysis of sorghum ERFs with 74 genetically characterized ERF genes from other plant species provided significant clues to sorghum ERF functions. Overall data generated here provides an overview of evolutionary relationship among ERF gene family members in sorghum and with respect to previously characterized ERF genes from other plant species. This information will be instrumental in initiating functional genomic studies of ERF candidates in sorghum.

Nanopore direct RNA sequencing maps the complexity of Arabidopsis mRNA processing and m6A modification

Understanding genome organization and gene regulation requires insight into RNA transcription, processing and modification. We adapted nanopore direct RNA sequencing to examine RNA from a wild-type accession of the model plant Arabidopsis thaliana and a mutant defective in mRNA methylation (m6A). Here we show that m6A can be mapped in full-length mRNAs transcriptome-wide and reveal the combinatorial diversity of cap-associated transcription start sites, splicing events, poly(A) site choice and poly(A) tail length. Loss of m6A from 3' untranslated regions is associated with decreased relative transcript abundance and defective RNA 3' end formation. A functional consequence of disrupted m6A is a lengthening of the circadian period. We conclude that nanopore direct RNA sequencing can reveal the complexity of mRNA processing and modification in full-length single molecule reads. These findings can refine Arabidopsis genome annotation. Further, applying this approach to less well-studied species could transform our understanding of what their genomes encode.

Characterization and Expression Analysis of ERF Genes in Fragaria vesca Suggest Different Divergences of Tandem ERF Duplicates

Ethylene-responsive factors (ERFs) play important roles in plant growth and development and in responses to abiotic stresses. However, little information was available about the ERF genes in woodland strawberry (Fragaria vesca), a genetic model plant for the Fragaria genus and Rosaceae family. In this study, 91 FveERF genes were identified, including 35 arrayed in tandem, indicating that tandem duplication is a major mechanism for the expansion of the FveERF family. According to their phylogenetic relationships with AtERFs from Arabidopsis thaliana, the tandem FveERF genes could be grouped into ancestral and lineage-specific tandem ones. The ancestral tandem FveERFs are likely derived from tandem duplications that occurred in the common ancestor of F. vesca and A. thaliana, whereas the lineage-specific ones are specifically present in the F. vesca lineage. The lineage-specific tandem FveERF duplicates are more conserved than the ancestral ones in sequence and structure. However, their expression in flowers and fruits is similarly diversified, indicating that tandem FveERFs have diverged rapidly after duplication in this respect. The lineage-specific tandem FveERFs display the same response patterns with only one exception under drought or cold, whereas the ancestral tandem ones are largely differentially expressed, suggesting that divergence of tandem FveERF expression under stress may have occurred later in the reproductive development. Our results provide evidence that the retention of tandem FveERF duplicates soon after their duplication may be related to their divergence in the regulation of reproductive development. In contrast, their further divergence in expression pattern likely contributes to plant response to abiotic stress.

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.

In Silico Genome-Wide Analysis of the ATP-Binding Cassette Transporter Gene Family in Soybean (Glycine max L.) and Their Expression Profiling

ATP-binding cassette (ABC) transporters constitute one of the largest gene families in all living organisms, most of which mediate transport across biological membranes by hydrolyzing ATP. However, detailed studies of ABC transporter genes in the important oil crop, soybean, are still lacking. In the present study, we carried out genome-wide identification and phylogenetic and transcriptional analyses of the ABC gene family in G. max. A total of 261 G. max ABC (GmABCs) genes were identified and unevenly localized onto 20 chromosomes. Referring to protein-domain orientation and phylogeny, the GmABC family could be classified into eight (ABCA-ABCG and ABCI) subfamilies and ABCG were the most abundantly present. Further, investigation of whole genome duplication (WGD) signifies the role of segmental duplication in the expansion of the ABC transporter gene family in soybean. The Ka/Ks ratio indicates that several duplicated genes are governed by intense purifying selection during evolution. In addition, in silico expression analysis based on RNA-sequence using publicly available database revealed that ABC transporters are differentially expressed in tissues and developmental stages and in dehydration. Overall, we provide an extensive overview of the GmABC transporter gene family and it promises the primary basis for the study in development and response to dehydration tolerance.

Demography and its effects on genomic variation in crop domestication

Over two thousand plant species have been modified morphologically through cultivation and human use. Here, we review three aspects of crop domestication that are currently undergoing marked revisions, due to analytical advancements and their application to whole genome resequencing (WGS) data. We begin by discussing the duration and demographic history of domestication. There has been debate as to whether domestication occurred quickly or slowly. The latter is tentatively supported both by fossil data and application of WGS data to sequentially Markovian coalescent methods that infer the history of effective population size. This history suggests the possibility of extended human impacts on domesticated lineages prior to their purposeful cultivation. We also make the point that demographic history matters, because it shapes patterns and levels of extant genetic diversity. We illustrate this point by discussing the evolutionary processes that contribute to the empirical observation that most crops examined to date have more putatively deleterious alleles than their wild relatives. These deleterious alleles may contribute to genetic load within crops and may be fitting targets for crop improvement. Finally, the same demographic factors are likely to shape the spectrum of structural variants (SVs) within crops. SVs are known to underlie many of the phenotypic changes associated with domestication and crop improvement, but we currently lack sufficient knowledge about the mechanisms that create SVs, their rates of origin, their population frequencies and their phenotypic effects.

GWAS Uncovers Differential Genetic Bases for Drought and Salt Tolerances in Sesame at the Germination Stage

Sesame has great potential as an industrial crop but its production is challenged by drought and salt stresses. To unravel the genetic variants leading to salinity and drought tolerances at the germination stage, genome-wide association studies of stress tolerance indexes related to NaCl-salt and polyethylene glycol-drought induced stresses were performed with a diversity panel of 490 sesame accessions. An extensive variation was observed for drought and salt responses in the population and most of the accessions were moderately tolerant to both stresses. A total of 132 and 120 significant Single Nucleotide Polymorphisms (SNPs) resolved to nine and 15 Quantitative trait loci (QTLs) were detected for drought and salt stresses, respectively. Only two common QTLs for drought and salt responses were found located on linkage groups 5 and 7, respectively. This indicates that the genetic bases for drought and salt responses in sesame are different. A total of 13 and 27 potential candidate genes were uncovered for drought and salt tolerance indexes, respectively, encoding transcription factors, antioxidative enzymes, osmoprotectants and involved in hormonal biosynthesis, signal transduction or ion sequestration. The identified SNPs and potential candidate genes represent valuable resources for future functional characterization towards the enhancement of sesame cultivars for drought and salt tolerances.

The Medicago truncatula Genome: Genomic Data Availability

Medicago truncatula emerged in 1990 as a model for legumes, comprising the third largest land plant family. Most legumes form symbiotic nitrogen-fixing root nodules with compatible soil bacteria and thus are important contributors to the global nitrogen cycle and sustainable agriculture. Legumes and legume products are important sources for human and animal protein as well as for edible and industrial oils. In the years since M. truncatula was chosen as a legume model, many genetic, genomic, and molecular resources have become available, including reference quality genome sequences for two widely used genotypes. Accessibility of genomic data is important for many different types of studies with M. truncatula as well as for research involving crop and forage legumes. In this chapter, we discuss strategies to obtain archived M. truncatula genomic data originally deposited into custom databases that are no longer maintained but are now accessible in general databases. We also review key current genomic databases that are specific to M. truncatula as well as those that contain M. truncatula data in addition to data from other plants.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Echinoderms, ubiquitous in the marine environment, are important from evolutionary, ecological, and socioeconomic perspectives. Together with chordates and hemichordates, they form the deuterostome clade, making them a crucial node in the study of chordate ancestry. Within echinoderms, class Holothuroidea is unique; its members (the sea cucumbers) display remarkable regenerative abilities and play key roles as sediment bioturbators and symbiotic hosts, and many are prized in the seafood and pharmaceutical industries. The sea cucumber genome therefore has the potential to significantly contribute to our understanding of important evolutionary and biological processes and help enhance aquaculture programs. Here we present a high-quality genome sequence for the economically important species Apostichopus japonicus. Through comparative analyses, we identified 763 echinoderm-specific gene families enriched in genes encoding membrane proteins, ion channels, and signal transduction proteins. Marker genes associated with the notochord and gill slits were also found, providing valuable insight into the origin of chordates. The reduced number and low expression levels of biomineralization genes reflect the skeletal degeneration seen in sea cucumbers. Importantly, 2 gene families appeared to be expanded in A. japonicus and may play crucial roles in its heightened regenerative potential. Together, findings from the sea cucumber genome provide important and novel insights into echinoderm and deuterostome biology.

Evolutionary history and functional divergence of the cytochrome P450 gene superfamily between Arabidopsis thaliana and Brassica species uncover effects of whole genome and tandem duplications

Background

The 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.

Results

Using 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.

Conclusions

This 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.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4094-7) contains supplementary material, which is available to authorized users.

Genome-wide analysis of WRKY gene family in the sesame genome and identification of the WRKY genes involved in responses to abiotic stresses

BACKGROUND

Sesame (Sesamum indicum L.) is one of the world's most important oil crops. However, it is susceptible to abiotic stresses in general, and to waterlogging and drought stresses in particular. The molecular mechanisms of abiotic stress tolerance in sesame have not yet been elucidated. The WRKY domain transcription factors play significant roles in plant growth, development, and responses to stresses. However, little is known about the number, location, structure, molecular phylogenetics, and expression of the WRKY genes in sesame.

RESULTS

We performed a comprehensive study of the WRKY gene family in sesame and identified 71 SiWRKYs. In total, 65 of these genes were mapped to 15 linkage groups within the sesame genome. A phylogenetic analysis was performed using a related species (Arabidopsis thaliana) to investigate the evolution of the sesame WRKY genes. Tissue expression profiles of the WRKY genes demonstrated that six SiWRKY genes were highly expressed in all organs, suggesting that these genes may be important for plant growth and organ development in sesame. Analysis of the SiWRKY gene expression patterns revealed that 33 and 26 SiWRKYs respond strongly to waterlogging and drought stresses, respectively. Changes in the expression of 12 SiWRKY genes were observed at different times after the waterlogging and drought treatments had begun, demonstrating that sesame gene expression patterns vary in response to abiotic stresses.

CONCLUSIONS

In this study, we analyzed the WRKY family of transcription factors encoded by the sesame genome. Insight was gained into the classification, evolution, and function of the SiWRKY genes, revealing their putative roles in a variety of tissues. Responses to abiotic stresses in different sesame cultivars were also investigated. The results of our study provide a better understanding of the structures and functions of sesame WRKY genes and suggest that manipulating these WRKYs could enhance resistance to waterlogging and drought.

The Emerging Oilseed Crop Sesamum indicum Enters the "Omics" Era

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops widely grown in Africa and Asia for its high-quality nutritional seeds. It is well adapted to harsh environments and constitutes an alternative cash crop for smallholders in developing countries. Despite its economic and nutritional importance, sesame is considered as an orphan crop because it has received very little attention from science. As a consequence, it lags behind the other major oil crops as far as genetic improvement is concerned. In recent years, the scenario has considerably changed with the decoding of the sesame nuclear genome leading to the development of various genomic resources including molecular markers, comprehensive genetic maps, high-quality transcriptome assemblies, web-based functional databases and diverse daft genome sequences. The availability of these tools in association with the discovery of candidate genes and quantitative trait locis for key agronomic traits including high oil content and quality, waterlogging and drought tolerance, disease resistance, cytoplasmic male sterility, high yield, pave the way to the development of some new strategies for sesame genetic improvement. As a result, sesame has graduated from an "orphan crop" to a "genomic resource-rich crop." With the limited research teams working on sesame worldwide, more synergic efforts are needed to integrate these resources in sesame breeding for productivity upsurge, ensuring food security and improved livelihood in developing countries. This review retraces the evolution of sesame research by highlighting the recent advances in the "Omics" area and also critically discusses the future prospects for a further genetic improvement and a better expansion of this crop.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

Background

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.

Results

We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.

Conclusion

We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

Electronic supplementary material

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

Genome evolutionary dynamics followed by diversifying selection explains the complexity of the Sesamum indicum genome

Background

Whole genome duplication (WGD) and tandem duplication (TD) provide two critical sources of raw genetic material for genome complexity and evolutionary novelty. Little is known about the complexity of the Sesamum indicum genome after it diverged from a common ancestor with the paleodiploid Vitis vinifera and further experienced WGD and TD events.

Results

Here, we analyzed the functional divergence of different classes of inter- and intra-genome gene pairs from ancestral events to uncover multiple-layers of evolutionary dynamics acting during the process of forming the modern S. indicum genome. Comprehensive inter-genome analyses revealed that 60% and 70% of syntenic orthologous gene pairs were retained among the two subgenomes in S. indicum compared to V. vinifera, although there was no evidence of significant differences under selection pressure. For the intra-genomic analyses, 5,932 duplicated gene pairs experienced fractionation, with the remaining 1,236 duplicated gene pairs having undergone functional divergence under diversifying selection. Analysis of the TD events indicated that 2,945 paralogous gene pairs, from 1,089 tandem arrays of 2–16 genes, experienced functional divergence under diversifying selection. Sequence diversification of different classes of gene pairs revealed that most of TD events occurred after the WGD event, with others following the ancestral gene order indicating ancient TD events at some time prior to the WGD event. Our comparison-of-function analyses for different classes of gene pairs indicated that the WGD and TD evolutionary events were both responsible for introducing genes that enabled exploration of novel and complementary functionalities, whilst maintaining individual plant ruggedness.

Conclusions

In this study, we first investigated functional divergence of different classes of gene pairs to characterize the dynamic processes associated with each evolutionary event in S. indicum. The data demonstrated massive and distinct functional divergence among different classes of gene pairs, and provided a genome-scale view of gene function diversification explaining the complexity of the S. indicum genome. We hope this provides a biological model to study the mechanism of plant species formation, particularly in the context of the evolutionary history of flowering plants, and offers novel insights for the study of angiosperm genomes.

Electronic supplementary material

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

Diversity, expansion, and evolutionary novelty of plant DNA-binding transcription factor families

Plant transcription factors (TFs) that interact with specific sequences via DNA-binding domains are crucial for regulating transcriptional initiation and are fundamental to plant development and environmental response. In addition, expansion of TF families has allowed functional divergence of duplicate copies, which has contributed to novel, and in some cases adaptive, traits in plants. Thus, TFs are central to the generation of the diverse plant species that we see today. Major plant agronomic traits, including those relevant to domestication, have also frequently arisen through changes in TF coding sequence or expression patterns. Here our goal is to provide an overview of plant TF evolution by first comparing the diversity of DNA-binding domains and the sizes of these domain families in plants and other eukaryotes. Because TFs are among the most highly expanded gene families in plants, the birth and death process of TFs as well as the mechanisms contributing to their retention are discussed. We also provide recent examples of how TFs have contributed to novel traits that are important in plant evolution and in agriculture.This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

Genome-Wide Identification, Evolutionary Analysis, and Stress Responses of the GRAS Gene Family in Castor Beans

Plant-specific GRAS transcription factors play important roles in regulating growth, development, and stress responses. Castor beans (Ricinus communis) are important non-edible oilseed plants, cultivated worldwide for its seed oils and its adaptability to growth conditions. In this study, we identified and characterized a total of 48 GRAS genes based on the castor bean genome. Combined with phylogenetic analysis, the castor bean GRAS members were divided into 13 distinct groups. Functional divergence analysis revealed the presence of mostly Type-I functional divergence. The gene structures and conserved motifs, both within and outside the GRAS domain, were characterized. Gene expression analysis, performed in various tissues and under a range of abiotic stress conditions, uncovered the potential functions of GRAS members in regulating plant growth development and stress responses. The results obtained from this study provide valuable information toward understanding the potential molecular mechanisms of GRAS proteins in castor beans. These findings also serve as a resource for identifying the genes that allow castor beans to grow in stressful conditions and to enable further breeding and genetic improvements in agriculture.