Placing Paleopolyploidy in Relation to Taxon Divergence: A Phylogenetic Analysis in Legumes Using 39 Gene Families

Phosphorus-Use-Efficiency Gene Identification in Fabaceae and RSL2 Expansion in Lupinus albus Is Associated with Low-Phosphorus Adaptation

Phosphorus is critical for plant growth but often becomes less accessible due to its precipitation with cations in soil. Fabaceae, a diverse plant family, exhibits robust adaptability and includes species like Lupinus albus, known for its efficient phosphorus utilization via cluster roots. Here, we systematically identified phosphorus-utilization-efficiency (PUE) gene families across 35 Fabaceae species, highlighting significant gene amplification in PUE pathways in Fabaceae. Different PUE pathways exhibited variable amplification, evolution, and retention patterns among various Fabaceae crops. Additionally, the number of homologous genes of the root hair development gene RSL2 in L. albus was far more than that in other Fabaceae species. Multiple copies of the RSL2 gene were amplified and retained in L. albus after whole genome triplication. The gene structure and motifs specifically retained in L. albus were different from homologous genes in other plants. Combining transcriptome analysis under low-phosphorus treatment, it was found that most of the homologous genes of RSL2 in L. albus showed high expression in the cluster roots, suggesting that the RSL2 gene family plays an important role in the adaptation process of L. albus to low-phosphorus environments and the formation of cluster roots.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

A chromosome-level genome assembly and evolutionary analysis of Coregonus ussuriensis Berg

Coregonus ussuriensis Berg, distributed widely in cold waters above 45° N latitude, is a savored freshwater whitefish that has been included in the list of endangered animals as a consequence of overfishing. Lack of genomic information seriously hampers evolutionary and genetic research on C. ussuriensis warranting the need to assemble a high-quality reference genome to promote its genetic breeding. We assembled and constructed a reference chromosome-level C. ussuriensis genome (sequence length, 2.51 Gb; contig N50 length, 4.27 Mb) using PacBio sequencing and Hi-C assembly technology, 3,109 contigs were assembled into scaffolds, resulting in a genome assembly with 40 chromosomes and a scaffold N50 length of 62.20 Mb. In addition, 43,320 protein-coding genes were annotated. The peak Ks position in the species comparison reflects the whole-genome replication event of C. ussuriensis. This chromosome-level genome provides reference data for further studies on the molecular breeding of C. ussuriensis.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Genome-wide analysis of the CML gene family and its response to melatonin in common bean (Phaseolus vulgaris L.)

Calmodulin-like proteins (CML) are important calcium signal transduction proteins in plants. CML genes have been analyzed in several plants. However, little information on CML in Phaseolus vulgare is available. In this study, we identified 111 PvCMLs distributed on eleven chromosomes. Phylogenetic analysis classified them into seven subfamilies. Cis-acting element prediction showed that PvCML contained elements related to growth and development, response to abiotic stress and hormones. Moreover, the majority of PvCMLs showed different expression patterns in most of the nine tissues and developmental stages which indicated the role of PvCML in the growth and development of common bean. Additionally, the common bean was treated with melatonin by seed soaking, and root transcriptome at the 5th day and qRT-PCR of different tissue at several stages were performed to reveal the response of PvCML to the hormone. Interestingly, 9 PvCML genes of subfamily VI were detected responsive to exogenous melatonin, and the expression dynamics of nine melatonin response PvCML genes after seed soaking with melatonin were revealed. Finally, the protein interaction network analysis of nine melatonin responsive PvCMLs was constructed. The systematic analysis of the PvCML gene family provides theoretical support for the further elucidation of their functions, and melatonin response molecular mechanism of the CML family in P. vulgaris.

Chromosome-scale genome assembly of Glycyrrhiza uralensis revealed metabolic gene cluster centred specialized metabolites biosynthesis

A high-quality genome assembly is imperative to explore the evolutionary basis of characteristic attributes that define chemotype and provide essential resources for a molecular breeding strategy for enhanced production of medicinal metabolites. Here, using single-molecule high-fidelity (HiFi) sequencing reads, we report chromosome-scale genome assembly for Chinese licorice (Glycyrrhiza uralensis), a widely used herbal and natural medicine. The entire genome assembly was achieved in eight chromosomes, with contig and scaffold N50 as 36.02 and 60.2 Mb, respectively. With only 17 assembly gaps and half of the chromosomes having no or one assembly gap, the presented genome assembly is among the best plant genomes to date. Our results showed an advantage of using highly accurate long-read HiFi sequencing data for assembling a highly heterozygous genome including its complexed repeat content. Additionally, our analysis revealed that G. uralensis experienced a recent whole-genome duplication at approximately 59.02 million years ago post a gamma (γ) whole-genome triplication event, which contributed to its present chemotype features. The metabolic gene cluster analysis identified 355 gene clusters, which included the entire biosynthesis pathway of glycyrrhizin. The genome assembly and its annotations provide an essential resource for licorice improvement through molecular breeding and the discovery of valuable genes for engineering bioactive components and understanding the evolution of specialized metabolites biosynthesis.

Genetic and Genomic Resources for Soybean Breeding Research

Soybean (Glycine max) is a legume species of significant economic and nutritional value. The yield of soybean continues to increase with the breeding of improved varieties, and this is likely to continue with the application of advanced genetic and genomic approaches for breeding. Genome technologies continue to advance rapidly, with an increasing number of high-quality genome assemblies becoming available. With accumulating data from marker arrays and whole-genome resequencing, studying variations between individuals and populations is becoming increasingly accessible. Furthermore, the recent development of soybean pangenomes has highlighted the significant structural variation between individuals, together with knowledge of what has been selected for or lost during domestication and breeding, information that can be applied for the breeding of improved cultivars. Because of this, resources such as genome assemblies, SNP datasets, pangenomes and associated databases are becoming increasingly important for research underlying soybean crop improvement.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

New insights on the function of plant acyl carrier proteins from comparative and evolutionary analysis

Acyl carrier proteins (ACPs) play a central role in both plastidial and mitochondrial Type II fatty acid synthesis in plant cells. However, a large proportion of plant ACPs remain functionally uncharacterized, and their evolutionary history remains elusive. In present study, 97 putative ACPs were identified from ten angiosperm species examined. Based on phylogenetic analysis, ACP genes were grouped into plastidial (cpACP: ACP1/2/3/4/5) and mitochondrial (mtACP: mtACP1/mtACP2/mtACP3) ACPs. Protein sequence (motifs and length), tertiary structure, and gene structure (exon number, average intron length, and intron phase) were highly conserved in different ACP subclades. The differentiation of ACPs into distinct types occurred 85-98 and 45-57 million years ago. A limited proportion of ACP genes experience tandem or segmental duplication, corresponding to two rounds of whole genome duplication. Ka/Ks ratios revealed that duplicated ACP genes underwent a purifying selection. Regarding expression patterns, most ACPs were expressed constitutively and tissue-specifically. Notably, the average expression levels of ACP1, mtACP3, and mtACP1 were positively correlated with those of ACP3, ACP4, and mtACP2, respectively. Analysis of cis-elements showed that seven motifs (CACTFTPPCA1, DOFCOREZM, GT1CONSENSUS, CAATBOX1, ARR1AT, POLLEN1LELAT52, and GATABOX) related to tissue-specific, ABA, and light-mediated gene regulation were ubiquitous in all ACPs investigated, which shed new light on the regulation patterns of these central enzymatic partners of the FAS system. This study presents a thorough overview of angiosperm ACP gene families and provides informative clues for the functional characterization of plant ACPs in the future.

Calmodulin-Like (CML) Gene Family in Medicago truncatula: Genome-Wide Identification, Characterization and Expression Analysis

Calcium is an important second messenger in mediating adaptation responses of plants to abiotic and biotic stresses. Calmodulin-like (CML) protein is an important calcium-signaling protein that can sense and decode Ca2+ signal in plants. Medicago truncatula is a model legume plant; however, investigations of MtCML proteins are limited. Using genome analysis and BLAST database searches, fifty MtCML proteins that possess EF-hand motifs were identified. Phylogenetic analysis showed that CML homologs between M. truncatula, Arabidopsis thaliana and Oryza sativa shared close relationships. Gene structure analysis revealed that these MtCML genes contained one to four conserved EF-hand motifs. All MtCMLs are localized to eight chromosomes and underwent gene duplication. In addition, MtCML genes were differentially expressed in different tissues of M. truncatula. Cis-acting elements in promoter region and expression analysis revealed the potential response of MtCML protein to abiotic stress and hormones. The results provide a basis of further functional research on the MtCML gene family and facilitate their potential use for applications in the genetic improvement on M. truncatula in drought, cold and salt stress environments.

Genome-wide study of C2H2 zinc finger gene family in Medicago truncatula

BACKGROUND

C2H2 zinc finger proteins (C2H2 ZFPs) play vital roles in shaping many aspects of plant growth and adaptation to the environment. Plant genomes harbor hundreds of C2H2 ZFPs, which compose one of the most important and largest transcription factor families in higher plants. Although the C2H2 ZFP gene family has been reported in several plant species, it has not been described in the model leguminous species Medicago truncatula.

RESULTS

In this study, we identified 218 C2H2 type ZFPs with 337 individual C2H2 motifs in M. truncatula. We showed that the high rate of local gene duplication has significantly contributed to the expansion of the C2H2 gene family in M. truncatula. The identified ZFPs exhibit high variation in motif arrangement and expression pattern, suggesting that the short C2H2 zinc finger motif has been adopted as a scaffold by numerous transcription factors with different functions to recognize cis-elements. By analyzing the public expression datasets and quantitative RT-PCR (qRT-PCR), we identified several C2H2 ZFPs that are specifically expressed in certain tissues, such as the nodule, seed, and flower.

CONCLUSION

Our genome-wide work revealed an expanded C2H2 ZFP gene family in an important legume M. truncatula, and provides new insights into the diversification and expansion of C2H2 ZFPs in higher plants.

Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa

Artificially improving traits of cultivated alfalfa (Medicago sativa L.), one of the most important forage crops, is challenging due to the lack of a reference genome and an efficient genome editing protocol, which mainly result from its autotetraploidy and self-incompatibility. Here, we generate an allele-aware chromosome-level genome assembly for the cultivated alfalfa consisting of 32 allelic chromosomes by integrating high-fidelity single-molecule sequencing and Hi-C data. We further establish an efficient CRISPR/Cas9-based genome editing protocol on the basis of this genome assembly and precisely introduce tetra-allelic mutations into null mutants that display obvious phenotype changes. The mutated alleles and phenotypes of null mutants can be stably inherited in generations in a transgene-free manner by cross pollination, which may help in bypassing the debate about transgenic plants. The presented genome and CRISPR/Cas9-based transgene-free genome editing protocol provide key foundations for accelerating research and molecular breeding of this important forage crop.

A Tale of Two Families: Whole Genome and Segmental Duplications Underlie Glutamine Synthetase and Phosphoenolpyruvate Carboxylase Diversity in Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been supplied with advanced genomic resources and, as such, has become a well-known model for molecular evolutionary studies within the legume family—a group of plants able to fix nitrogen from the atmosphere. The phylogenetic position of lupins in Papilionoideae and their evolutionary distance to other higher plants facilitates the use of this model species to improve our knowledge on genes involved in nitrogen assimilation and primary metabolism, providing novel contributions to our understanding of the evolutionary history of legumes. In this study, we present a complex characterization of two narrow-leafed lupin gene families—glutamine synthetase (GS) and phosphoenolpyruvate carboxylase (PEPC). We combine a comparative analysis of gene structures and a synteny-based approach with phylogenetic reconstruction and reconciliation of the gene family and species history in order to examine events underlying the extant diversity of both families. Employing the available evidence, we show the impact of duplications on the initial complement of the analyzed gene families within the genistoid clade and posit that the function of duplicates has been largely retained. In terms of a broader perspective, our results concerning GS and PEPC gene families corroborate earlier findings pointing to key whole genome duplication/triplication event(s) affecting the genistoid lineage.

The structure, functional evolution, and evolutionary trajectories of the H+-PPase gene family in plants

BACKGROUND

The H+-PPase (pyrophosphatase) gene family is an important class of proton transporters that play key roles in plant development and stress resistance. Although the physiological and biochemical functions of H+-PPases are well characterized, the structural evolution and functional differentiation of this gene family remain unclear.

RESULTS

We identified 124 H+-PPase members from 27 plant species using complete genomic data obtained from algae to angiosperms. We found that all analyzed plants carried H+-PPase genes, and members were not limited to the two main types (type I and II). Differentiation of this gene family occurred early in evolutionary history, probably prior to the emergence of algae. The type I and II H+-PPase genes were retained during the subsequent evolution of higher plants, and their copy numbers increased rapidly in some angiosperms following whole-genome duplication (WGD) events, with obvious expression pattern differentiation among the new copies. We found significant functional divergence between type I and II H+-PPase genes, with both showing evidence for positive selection pressure. We classified angiosperm type I H+-PPases into subtypes Ia and non-Ia, which probably differentiated at an early stage of angiosperm evolution. Compared with non-Ia subtype, the Ia subtype appears to confer some advantage in angiosperms, as it is highly conserved and abundantly expressed, but shows no evidence for positive selection.

CONCLUSIONS

We hypothesized that there were many types of H+-PPase genes in the plant ancestral genome, and that different plant groups retained different types of these genes. In the early stages of angiosperm evolution, the type I H+-PPase genes differentiated into various subtypes. In addition, the expression pattern varied not only among genes of different types or subtypes, but also among copies of the same subtype. Based on the expression patterns and copy numbers of H+-PPase genes in higher plants, we propose two possible evolutionary trajectories for this gene family.

The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection

Laccase is a widely used industrial oxidase for food processing, dye synthesis, paper making, and pollution remediation. At present, laccases used by industries come mainly from fungi. Plants contain numerous genes encoding laccase enzymes that show properties which are distinct from that of the fungal laccases. These plant-specific laccases may have better potential for industrial purposes. The aim of this work was to conduct a genome-wide search for the soybean laccase genes and analyze their characteristics and specific functions. A total of 93 putative laccase genes (GmLac) were identified from the soybean genome. All 93 GmLac enzymes contain three typical Cu-oxidase domains, and they were classified into five groups based on phylogenetic analysis. Although adjacent members on the tree showed highly similar exon/intron organization and motif composition, there were differences among the members within a class for both conserved and differentiated functions. Based on the expression patterns, some members of laccase were expressed in specific tissues/organs, while some exhibited a constitutive expression pattern. Analysis of the transcriptome revealed that some laccase genes might be involved in providing resistance to oomycetes. Analysis of the selective pressures acting on the laccase gene family in the process of soybean domestication revealed that 10 genes could have been under artificial selection during the domestication process. Four of these genes may have contributed to the transition of the soft and thin stem of wild soybean species into strong, thick, and erect stems of the cultivated soybean species. Our study provides a foundation for future functional studies of the soybean laccase gene family.

Reconstruction of ancestral genome reveals chromosome evolution history for selected legume species

Reconstruction of an ancestral genome for a set of plant species has been a challenging task because of complex histories that may include whole-genome duplications, segmental duplications, independent gene duplications or losses, diploidization and rearrangement events. Here, we describe the reconstruction a hypothetical ancestral genome for the papilionoid legumes (the largest subfamily within the third largest family in flowering plants), and evaluate the results relative to phylogenetic and chromosomal count data for this group of legumes, spanning 294 diverse papilionoid genera. To reconstruct the ancestral genomes for nine legume species with sequenced genomes, we used a maximum likelihood approach combined with a novel method for identifying informative markers for this purpose. Analyzing genomes from four species within the Phaseoleae, two in Dalbergieae, two in the 'inverted repeat loss' clade, and one in the Robinieae, we infer a common ancestral genome with nine chromosomes. The reconstructed genome structural histories are consistent with chromosomal and phylogenetic histories, but we also infer that a common ancestor with nine chromosomes was probably intermediate to an earlier state of 14 chromosomes following a whole-genome duplication that pre-dated the radiation of the papilionoid legumes, evidence for which is found in early-diverging papilionoid lineages.

Ancient Polyploidy and Genome Evolution in Palms

Mechanisms of genome evolution are fundamental to our understanding of adaptation and the generation and maintenance of biodiversity, yet genome dynamics are still poorly characterized in many clades. Strong correlations between variation in genomic attributes and species diversity across the plant tree of life suggest that polyploidy or other mechanisms of genome size change confer selective advantages due to the introduction of genomic novelty. Palms (order Arecales, family Arecaceae) are diverse, widespread, and dominant in tropical ecosystems, yet little is known about genome evolution in this ecologically and economically important clade. Here, we take a phylogenetic comparative approach to investigate palm genome dynamics using genomic and transcriptomic data in combination with a recent, densely sampled, phylogenetic tree. We find conclusive evidence of a paleopolyploid event shared by the ancestor of palms but not with the sister clade, Dasypogonales. We find evidence of incremental chromosome number change in the palms as opposed to one of recurrent polyploidy. We find strong phylogenetic signal in chromosome number, but no signal in genome size, and further no correlation between the two when correcting for phylogenetic relationships. Palms thus add to a growing number of diverse, ecologically successful clades with evidence of whole-genome duplication, sister to a species-poor clade with no evidence of such an event. Disentangling the causes of genome size variation in palms moves us closer to understanding the genomic conditions facilitating adaptive radiation and ecological dominance in an evolutionarily successful, emblematic tropical clade.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean’s agronomic traits.

Horizontal Gene Transfer in Five Parasite Plant Species in Orobanchaceae

We sequenced genomes of five parasite species in family Orobanchaceae to explore the evolutionary role of horizontal gene transfer in plants. Orobanche minor and Aeginetia indica are obligate parasites with no photosynthetic activity, whereas the other three (Pedicularis keiskei, Phtheirospermum japonicum, and Melampyrum roseum) are facultative parasites. By using reference genome sequences and/or transcriptomes of 14 species from Fabaceae and Poaceae, their major host families, we detected 106 horizontally transferred genes (HGT genes), only in the genomes of the two obligate parasites (22 and 84 for Oro. minor and Ae. indica, respectively), whereas none in the three facultative parasites. The HGT genes, respectively, account for roughly 0.1% and 0.2% of the coding genes in the two species. We found that almost all HGT genes retained introns at the same locations as their homologs in potential host species, indicating a crucial role of DNA-mediated gene transfer, rather than mRNA mediated retro transfer. Furthermore, some of the HGT genes might have transferred simultaneously because they located very closely in the host reference genome, indicating that the length of transferred DNA could exceed 100 kb. We confirmed that almost all introns are spliced in the current genome of the parasite species, and that about half HGT genes do not have any missense mutations or frameshift-causing indels, suggesting that some HGT genes may be still functional. Evolutionary analyses revealed that the nonsynonymous–synonymous substitution ratio is on average elevated on the lineage leading to HGT genes, due to either relaxation of selection or positive selection.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Genome-wide identification of the auxin response factor gene family in Cicer arietinum

Background

Auxin Response Factors act as critical components of the auxin-signaling pathway by regulating the transcription of auxin-responsive genes. The release of the chickpea reference genome provides an opportunity to identify and characterize the ARF gene family in this important legume by a data mining coupled by comparative genomics approaches.

Results

We performed a comprehensive characterization and analysis of 24 ARF genes in the chickpea reference genome. Comparative phylogenetic analysis of the ARF from chickpea, Medicago and Arabidopsis suggests that recent duplications have played a very limited role in the expansion of the ARF chickpea family. Gene structure analysis based on exon-intron organization provides additional evidence to support the evolutionary relationship among the ARF members. Conserved motif analysis shows that most of the proteins fit into the canonical ARF structure model, but 9 proteins lack or have a truncated dimerization domain. The mechanisms underlying the diversification of the ARF gene family are based on duplications, variations in domain organization and alternative splicing. Concerning duplications, segmental, but not tandem duplications, have contributed to the expansion of the gene family. Moreover, the duplicated pair genes have evolved mainly under the influence of purifying selection pressure with restricted functional divergence. Expression profiles responding to various environmental stimuli show a close relationship between tissue and expression patterns. Promoter sequence analysis reveals an enrichment of several cis-regulatory elements related to symbiosis, and modulation of plant gene expression during the interaction with microbes.

Conclusions

In conclusion, this study provides a comprehensive overview of the ARF gene family in chickpea. Globally, our data supports that auxin signaling pathway regulates a wide range of physiological processes and stress responses. Our findings could further provide new insights into the complexity of the regulation of ARF at the transcription level that may be useful to develop rational chickpea breeding strategies to improve development or stress responses. Our study also provides a foundation for comparative genomic analyses and a framework to trace the dynamic evolution of ARF genes on a large time-scale within the legume family.

Electronic supplementary material

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

Evolutionary history of the cobalamin-independent methionine synthase gene family across the land plants

Plants are successful paleopolyploids. The wide diversity of land plants is driven strongly by their gene duplicates undergoing distinct evolutionary fates after duplication. We used genomic resources from 35 model plant species to unravel the evolutionary fate of gene copies (paralogs) of the cobalamin-independent methionine synthase (metE) gene family across the land plants. To explore genealogical relationships and characterize positive selection as a driving force in the evolution of metE paralogs within a single species, we carried out complementary analyses on genomic data of 32 genotypes of soybean. The size of the metE gene family remained small across the land plants; most of the studied species possessed 1-6 paralogs. Gene products were either cytosolic or chloroplastic; this dual subcellular distribution arose early during the divergence of the land plants and reached all extant lineages. Biased gene loss and gene retention events took place multiple times; recurrent evolution remodeled redundant metE paralogs to recover and maintain the dual subcellular distribution of MetE. Shared whole-genome duplication events gave rise to the metE paralogs of both soybean and Medicago truncatula. In soybean, the ancestral paralog pair GlymaPP2A encoded a cytosolic isoform of MetE, was under strong purifying selection, and retained high levels of expression across eight RNA-seq expression libraries. The daughters GlymaPP1 and GlymaPP2B showed accelerated rates of evolution, accumulated many sites predicted to be under positive selection, and possessed low levels of expression. Our results suggest that the metE paralogs of soybean follow Ohno's neofunctionalization model of gene duplicate evolution.

Legume isoflavone synthase genes have evolved by whole-genome and local duplications yielding transcriptionally active paralogs

Isoflavone synthase (IFS) is the key enzyme of isoflavonoid biosynthesis. IFS genes were identified in numerous species, although their evolutionary patterns have not yet been reconstructed. To address this issue, we performed structural and functional genomic analysis. Narrow leafed lupin, Lupinus angustifolius L., was used as a reference species for the genus, because it has the most developed molecular tools available. Nuclear genome BAC library clones carrying IFS homologs were localized by linkage mapping and fluorescence in situ hybridization in three chromosome pairs. Annotation of BAC, scaffold and transcriptome sequences confirmed the presence of three full-length IFS genes in the genome. Microsynteny analysis and Bayesian inference provided clear evidence that IFS genes in legumes have evolved by lineage-specific whole-genome and tandem duplications. Gene expression profiling and RNA-seq data mining showed that the vast majority of legume IFS copies have maintained their transcriptional activity. L. angustifolius IFS homologs exhibited organ-specific expression patterns similar to those observed in other Papilionoideae. Duplicated lupin IFS homologs retained non-negligible levels of substitutions in conserved motifs, putatively due to positive selection acting during early evolution of the genus, before the whole-genome duplication. Strong purifying selection preserved newly arisen IFS duplicates from further nonsynonymous changes.

A tree of life based on ninety-eight expressed genes conserved across diverse eukaryotic species

Rapid advances in DNA sequencing technologies have resulted in the accumulation of large data sets in the public domain, facilitating comparative studies to provide novel insights into the evolution of life. Phylogenetic studies across the eukaryotic taxa have been reported but on the basis of a limited number of genes. Here we present a genome-wide analysis across different plant, fungal, protist, and animal species, with reference to the 36,002 expressed genes of the rice genome. Our analysis revealed 9831 genes unique to rice and 98 genes conserved across all 49 eukaryotic species analysed. The 98 genes conserved across diverse eukaryotes mostly exhibited binding and catalytic activities and shared common sequence motifs; and hence appeared to have a common origin. The 98 conserved genes belonged to 22 functional gene families including 26S protease, actin, ADP–ribosylation factor, ATP synthase, casein kinase, DEAD-box protein, DnaK, elongation factor 2, glyceraldehyde 3-phosphate, phosphatase 2A, ras-related protein, Ser/Thr protein phosphatase family protein, tubulin, ubiquitin and others. The consensus Bayesian eukaryotic tree of life developed in this study demonstrated widely separated clades of plants, fungi, and animals. Musa acuminata provided an evolutionary link between monocotyledons and dicotyledons, and Salpingoeca rosetta provided an evolutionary link between fungi and animals, which indicating that protozoan species are close relatives of fungi and animals. The divergence times for 1176 species pairs were estimated accurately by integrating fossil information with synonymous substitution rates in the comprehensive set of 98 genes. The present study provides valuable insight into the evolution of eukaryotes.

Genome-Wide Identification and Expression, Protein-Protein Interaction and Evolutionary Analysis of the Seed Plant-Specific BIG GRAIN and BIG GRAIN LIKE Gene Family

BIG GRAIN1 (BG1) is an auxin-regulated gene which functions in auxin pathway and positively regulates biomass, grain size and yield in rice. However, the evolutionary origin and divergence of these genes are still unknown. In this study, we found that BG genes are probably originated in seed plants. We also identified that seed plants evolved a class of BIG GRAIN LIKE (BGL) genes which share conserved middle and C-terminal motifs with BG. The BG genes were present in all monocot and eudicot species analyzed; however, the BGL genes were absent in few monocot lineages. Both BG and BGL were found to be serine-rich proteins; however, differences in expansion and rates of retention after whole genome duplication events were observed. Promoters of BG and BGL genes were found to be enriched with auxin-responsive elements and the Arabidopsis thaliana BG and BGL genes were found to be auxin-regulated. The auxin-induced expression of AthBG2 was found to be dependent on the conserved ARF17/19 module. Protein-protein interaction analysis identified that AthBG2 interact with regulators of splicing, transcription and chromatin remodeling. Taken together, this study provides interesting insights about BG and BGL genes and incentivizes future work in this gene family which has the potential to be used for crop manipulation.

Evolution in African tropical trees displaying ploidy-habitat association: The genus Afzelia (Leguminosae)

Polyploidy has rarely been documented in rain forest trees but it has recently been found in African species of the genus Afzelia (Leguminosae), which is composed of four tetraploid rain forest species and two diploid dry forest species. The genus Afzelia thus provides an opportunity to examine how and when polyploidy and habitat shift occurred in Africa, and whether they are associated. In this study, we combined three plastid markers (psbA, trnL, ndhF), two nuclear markers (ribosomal ITS and the single-copy PEPC E7 gene), plastomes (obtained by High Throughput Sequencing) and morphological traits, with an extensive taxonomic and geographic sampling to explore the evolutionary history of Afzelia. Both nuclear DNA and morphological vegetative characters separated diploid from tetraploid lineages. Although the two African diploid species were well differentiated genetically and morphologically, the relationships among the tetraploid species were not resolved. In contrast to the nuclear markers, plastid markers revealed that one of the diploid species forms a well-supported clade with the tetraploids, suggesting historical hybridisation, possibly in relation with genome duplication (polyploidization) and habitat shift from dry to rain forests. Molecular dating based on fossil-anchored gene phylogenies indicates that extant Afzelia started diverging c. 14.5 or 20Ma while extant tetraploid species started diverging c. 7.0 or 9.4Ma according to plastid and nuclear DNA, respectively. Additional studies of tropical polyploid plants are needed to assess whether the ploidy-habitat association observed in African Afzelia would reflect a role of polyploidization in niche divergence in the tropics.

Expansion of the phosphatidylethanolamine binding protein family in legumes: a case study of Lupinus angustifolius L. FLOWERING LOCUS T homologs, LanFTc1 and LanFTc2

BACKGROUND

The Arabidopsis FLOWERING LOCUS T (FT) gene, a member of the phosphatidylethanolamine binding protein (PEBP) family, is a major controller of flowering in response to photoperiod, vernalization and light quality. In legumes, FT evolved into three, functionally diversified clades, FTa, FTb and FTc. A milestone achievement in narrow-leafed lupin (Lupinus angustifolius L.) domestication was the loss of vernalization responsiveness at the Ku locus. Recently, one of two existing L. angustifolius homologs of FTc, LanFTc1, was revealed to be the gene underlying Ku. It is the first recorded involvement of an FTc homologue in vernalization. The evolutionary basis of this phenomenon in lupin has not yet been deciphered.

RESULTS

Bacterial artificial chromosome (BAC) clones carrying LanFTc1 and LanFTc2 genes were localized in different mitotic chromosomes and constituted sequence-specific landmarks for linkage groups NLL-10 and NLL-17. BAC-derived superscaffolds containing LanFTc genes revealed clear microsyntenic patterns to genome sequences of nine legume species. Superscaffold-1 carrying LanFTc1 aligned to regions encoding one or more FT-like genes whereas superscaffold-2 mapped to a region lacking such a homolog. Comparative mapping of the L. angustifolius genome assembly anchored to linkage map localized superscaffold-1 in the middle of a 15 cM conserved, collinear region. In contrast, superscaffold-2 was found at the edge of a 20 cM syntenic block containing highly disrupted collinearity at the LanFTc2 locus. 118 PEBP-family full-length homologs were identified in 10 legume genomes. Bayesian phylogenetic inference provided novel evidence supporting the hypothesis that whole-genome and tandem duplications contributed to expansion of PEBP-family genes in legumes. Duplicated genes were subjected to strong purifying selection. Promoter analysis of FT genes revealed no statistically significant sequence similarity between duplicated copies; only RE-alpha and CCAAT-box motifs were found at conserved positions and orientations.

CONCLUSIONS

Numerous lineage-specific duplications occurred during the evolution of legume PEBP-family genes. Whole-genome duplications resulted in the origin of subclades FTa, FTb and FTc and in the multiplication of FTa and FTb copy number. LanFTc1 is located in the region conserved among all main lineages of Papilionoideae. LanFTc1 is a direct descendant of ancestral FTc, whereas LanFTc2 appeared by subsequent duplication.

Homoeologs: What Are They and How Do We Infer Them?

The evolutionary history of nearly all flowering plants includes a polyploidization event. Homologous genes resulting from allopolyploidy are commonly referred to as 'homoeologs', although this term has not always been used precisely or consistently in the literature. With several allopolyploid genome sequencing projects under way, there is a pressing need for computational methods for homoeology inference. Here we review the definition of homoeology in historical and modern contexts and propose a precise and testable definition highlighting the connection between homoeologs and orthologs. In the second part, we survey experimental and computational methods of homoeolog inference, considering the strengths and limitations of each approach. Establishing a precise and evolutionarily meaningful definition of homoeology is essential for understanding the evolutionary consequences of polyploidization.

Integration of Lupinus angustifolius L. (narrow-leafed lupin) genome maps and comparative mapping within legumes

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been considered a reference genome for the Lupinus genus. In the present work, genetic and cytogenetic maps of L. angustifolius were supplemented with 30 new molecular markers representing lupin genome regions, harboring genes involved in nitrogen fixation during the symbiotic interaction of legumes and soil bacteria (Rhizobiaceae). Our studies resulted in the precise localization of bacterial artificial chromosomes (BACs) carrying sequence variants for early nodulin 40, nodulin 26, nodulin 45, aspartate aminotransferase P2, asparagine synthetase, cytosolic glutamine synthetase, and phosphoenolpyruvate carboxylase. Together with previously mapped chromosomes, the integrated L. angustifolius map encompasses 73 chromosome markers, including 5S ribosomal DNA (rDNA) and 45S rDNA, and anchors 20 L. angustifolius linkage groups to corresponding chromosomes. Chromosomal identification using BAC fluorescence in situ hybridization identified two BAC clones as narrow-leafed lupin centromere-specific markers, which served as templates for preliminary studies of centromere composition within the genus. Bioinformatic analysis of these two BACs revealed that centromeric/pericentromeric regions of narrow-leafed lupin chromosomes consisted of simple sequence repeats ordered into tandem repeats containing the trinucleotide and pentanucleotide simple sequence repeats AGG and GATAC, structured into long arrays. Moreover, cross-genus microsynteny analysis revealed syntenic patterns of 31 single-locus BAC clones among several legume species. The gene and chromosome level findings provide evidence of ancient duplication events that must have occurred very early in the divergence of papilionoid lineages. This work provides a strong foundation for future comparative mapping among legumes and may facilitate understanding of mechanisms involved in shaping legume chromosomes.

Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.)

Plant disease resistance genes (R-genes) play a critical role in the defense response to pathogens. Barley is one of the most important cereal crops, having a genome recently made available, for which the diversity and evolution of R-genes are not well understood. The main objectives of this research were to conduct a genome-wide identification of barley Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) genes and elucidate their evolutionary history. We employed a Hidden Markov Model using 52 Arabidopsis thaliana CNL reference sequences and analyzed for phylogenetic relationships, structural variation, and gene clustering. We identified 175 barley CNL genes nested into three clades, showing (a) evidence of an expansion of the CNL-C clade, primarily due to tandem duplications; (b) very few members of clade CNL-A and CNL-B; and (c) a complete absence of clade CNL-D. Our results also showed that several of the previously identified mildew locus A (MLA) genes may be allelic variants of two barley CNL genes, MLOC_66581 and MLOC_10425, which respond to powdery mildew. Approximately 23% of the barley CNL genes formed 15 gene clusters located in the extra-pericentromeric regions on six of the seven chromosomes; more than half of the clustered genes were located on chromosomes 1H and 7H. Higher average numbers of exons and multiple splice variants in barley relative to those in Arabidopsis and rice may have contributed to a diversification of the CNL-C members. These results will help us understand the evolution of R-genes with potential implications for developing durable resistance in barley cultivars.

A Conserved Cytochrome P450 Evolved in Seed Plants Regulates Flower Maturation

Global inspection of plant genomes identifies genes maintained in low copies across taxa and under strong purifying selection, which are likely to have essential functions. Based on this rationale, we investigated the function of the low-duplicated CYP715 cytochrome P450 gene family that appeared early in seed plants and evolved under strong negative selection. Arabidopsis CYP715A1 showed a restricted tissue-specific expression in the tapetum of flower buds and in the anther filaments upon anthesis. cyp715a1 insertion lines showed a strong defect in petal development, and transient alteration of pollen intine deposition. Comparative expression analysis revealed the downregulated expression of genes involved in pollen development, cell wall biogenesis, hormone homeostasis, and floral sesquiterpene biosynthesis, especially TPS21 and several key genes regulating floral development such as MYB21, MYB24, and MYC2. Accordingly, floral sesquiterpene emission was suppressed in the cyp715a1 mutants. Flower hormone profiling, in addition, indicated a modification of gibberellin homeostasis and a strong disturbance of the turnover of jasmonic acid derivatives. Petal growth was partially restored by the active gibberellin GA3 or the functional analog of jasmonoyl-isoleucine, coronatine. CYP715 appears to function as a key regulator of flower maturation, synchronizing petal expansion and volatile emission. It is thus expected to be an important determinant of flower-insect interaction.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

An advanced draft genome assembly of a desi type chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is an important pulse legume crop. We previously reported a draft genome assembly of the desi chickpea cultivar ICC 4958. Here we report an advanced version of the ICC 4958 genome assembly (version 2.0) generated using additional sequence data and an improved genetic map. This resulted in 2.7-fold increase in the length of the pseudomolecules and substantial reduction of sequence gaps. The genome assembly covered more than 94% of the estimated gene space and predicted the presence of 30,257 protein-coding genes including 2230 and 133 genes encoding potential transcription factors (TF) and resistance gene homologs, respectively. Gene expression analysis identified several TF and chickpea-specific genes with tissue-specific expression and displayed functional diversification of the paralogous genes. Pairwise comparison of pseudomolecules in the desi (ICC 4958) and the earlier reported kabuli (CDC Frontier) chickpea assemblies showed an extensive local collinearity with incongruity in the placement of large sequence blocks along the linkage groups, apparently due to use of different genetic maps. Single nucleotide polymorphism (SNP)-based mining of intra-specific polymorphism identified more than four thousand SNPs differentiating a desi group and a kabuli group of chickpea genotypes.

Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes

Unresolved questions about evolution of the large and diverse legume family include the timing of polyploidy (whole-genome duplication; WGDs) relative to the origin of the major lineages within the Fabaceae and to the origin of symbiotic nitrogen fixation. Previous work has established that a WGD affects most lineages in the Papilionoideae and occurred sometime after the divergence of the papilionoid and mimosoid clades, but the exact timing has been unknown. The history of WGD has also not been established for legume lineages outside the Papilionoideae. We investigated the presence and timing of WGDs in the legumes by querying thousands of phylogenetic trees constructed from transcriptome and genome data from 20 diverse legumes and 17 outgroup species. The timing of duplications in the gene trees indicates that the papilionoid WGD occurred in the common ancestor of all papilionoids. The earliest diverging lineages of the Papilionoideae include both nodulating taxa, such as the genistoids (e.g., lupin), dalbergioids (e.g., peanut), phaseoloids (e.g., beans), and galegoids (=Hologalegina, e.g., clovers), and clades with nonnodulating taxa including Xanthocercis and Cladrastis (evaluated in this study). We also found evidence for several independent WGDs near the base of other major legume lineages, including the Mimosoideae-Cassiinae-Caesalpinieae (MCC), Detarieae, and Cercideae clades. Nodulation is found in the MCC and papilionoid clades, both of which experienced ancestral WGDs. However, there are numerous nonnodulating lineages in both clades, making it unclear whether the phylogenetic distribution of nodulation is due to independent gains or a single origin followed by multiple losses.

Doubling down on genomes: polyploidy and crop plants

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

Evolutionary study of the isoflavonoid pathway based on multiple copies analysis in soybean

BACKGROUND

Previous studies suggest that the metabolic pathway structure influences the selection and evolution rates of involved genes. However, most of these studies have exclusively considered a single gene copy encoding each enzyme in the metabolic pathway. Considering multiple-copy encoding enzymes could provide direct evidence of gene evolution and duplication patterns in metabolic pathways. We conducted a detailed analysis of the phylogeny, synteny, evolutionary rate and selection pressure of the genes in the isoflavonoid metabolic pathway of soybeans.

RESULTS

The results revealed that 1) only the phenylalanine ammonia-lyase (PAL) gene family most upstream from the pathway preserved all of the ancient and recent segmental duplications and maintained a strongly conserved synteny among these duplicated segments; gene families encoding branch-point enzymes with higher pleiotropy tended to retain more types of duplication; and genes encoding chalcone reductase (CHR) and isoflavone synthase (IFS) specific for legumes retained only recent segmental duplications; 2) downstream genes evolved faster than upstream genes and were subject to positive selection or relaxed selection constraints; 3) gene members encoding enzymes with high pleiotropy at the branching points were more likely to have undergone evolutionary differentiation, which may correspond to their functional divergences.

CONCLUSIONS

We reconciled our results with existing controversies and proposed that gene copies at branch points with higher connectivity might be under stronger selective constraints and that the gene copies controlling metabolic flux allocation underwent positive selection. Our analyses demonstrated that the structure and function of a metabolic pathway shapes gene duplication and the evolutionary constraints of constituent enzymes.

Molecular and ultrastructural analysis of forisome subunits reveals the principles of forisome assembly

BACKGROUND AND AIMS

Forisomes are specialized structural phloem proteins that mediate sieve element occlusion after wounding exclusively in papilionoid legumes, but most studies of forisome structure and function have focused on the Old World clade rather than the early lineages. A comprehensive phylogenetic, molecular, structural and functional analysis of forisomes from species covering a broad spectrum of the papilionoid legumes was therefore carried out, including the first analysis of Dipteryx panamensis forisomes, representing the earliest branch of the Papilionoideae lineage. The aim was to study the molecular, structural and functional conservation among forisomes from different tribes and to establish the roles of individual forisome subunits.

METHODS

Sequence analysis and bioinformatics were combined with structural and functional analysis of native forisomes and artificial forisome-like protein bodies, the latter produced by expressing forisome genes from different legumes in a heterologous background. The structure of these bodies was analysed using a combination of confocal laser scanning microscopy (CLSM), scanning electron microscopy (SEM) and transmission electron microscopy (TEM), and the function of individual subunits was examined by combinatorial expression, micromanipulation and light microscopy.

KEY RESULTS

Dipteryx panamensis native forisomes and homomeric protein bodies assembled from the single sieve element occlusion by forisome (SEO-F) subunit identified in this species were structurally and functionally similar to forisomes from the Old World clade. In contrast, homomeric protein bodies assembled from individual SEO-F subunits from Old World species yielded artificial forisomes differing in proportion to their native counterparts, suggesting that multiple SEO-F proteins are required for forisome assembly in these plants. Structural differences between Medicago truncatula native forisomes, homomeric protein bodies and heteromeric bodies containing all possible subunit combinations suggested that combinations of SEO-F proteins may fine-tune the geometric proportions and reactivity of forisomes.

CONCLUSIONS

It is concluded that forisome structure and function have been strongly conserved during evolution and that species-dependent subsets of SEO-F proteins may have evolved to fine-tune the structure of native forisomes.

Genome-wide identification and characterization of RBR ubiquitin ligase genes in soybean

RBR (RING1-IBR-RING2) proteins play an important role in protein ubiquitination and are involved in many cellular processes. Recent studies showed plant RBR genes were induced by abiotic and biotic stresses. However, detailed studies on RBR genes in the important oil crop, soybean (Glycine max (L.) Merr.), is still lacking. Here we performed a genome-wide search and identified 24 RBR domain-containing genes from the soybean genome sequence and cloned 11 of them. Most soybean RBR proteins contain a highly conserved RBR supra-domain. Phylogenetic analyses indicated all 24 soybean RBR proteins are most related to the RBR proteins from Phaseolus vulgaris, and could be classified into seven groups including Ariadne A, Ariadne B, ARA54, Plant IIA, Plant IIB, Plant IIC, and Helicase. Tandem duplication and block duplication were found among the Ariadne B and Plant IIC group of soybean RBR genes. Despite the conserved RBR supra-domain, there are extensive variations in the additional protein motifs and exon-intron structures between different groups, which indicate they might have diverse functions. Most soybean RBR proteins are predicted to localize in nucleus, and four of them were experimentally confirmed by GFP fusion proteins. Soybean RBR genes are broadly expressed in many tissue types with a little more abundant in the roots and flowers than leaves, stems, and seeds. The expression of GmRTRTP3 (Plant IIB) and GmRTRTP5 (Plant IIC) are induced by NaCl treatment, which suggests these RBR genes might be involved in soybean response to abiotic stresses.

Nod factor perception protein carries weight in biotic interactions

Plant plasma membrane-bound receptors with extracellular lysin motif (LysM) domains participate in interactions with microorganisms. In Medicago truncatula, the LysM receptor-like kinase gene nodulation (Nod) factor perception (NFP) is a key gene that controls the perception of rhizobial lipochitooligosaccharide (LCO) Nod factors for the establishment of the Rhizobium-legume symbiosis. In this article, we review recent data that have refined our understanding of this function and that have revealed a role for NFP in the perception of arbuscular mycorrhizal (AM) symbiotic signals and plant pathogenic microorganisms. The dual role of NFP in symbiosis and immunity suggests that this receptor protein controls the perception of different signals and the activation of different downstream signalling pathways. These advances provide new insights into the evolution and functioning of this versatile plant protein.

Whole genome duplication events in plant evolution reconstructed and predicted using myosin motor proteins

Background

The evolution of land plants is characterized by whole genome duplications (WGD), which drove species diversification and evolutionary novelties. Detecting these events is especially difficult if they date back to the origin of the plant kingdom. Established methods for reconstructing WGDs include intra- and inter-genome comparisons, K_S age distribution analyses, and phylogenetic tree constructions.

Results

By analysing 67 completely sequenced plant genomes 775 myosins were identified and manually assembled. Phylogenetic trees of the myosin motor domains revealed orthologous and paralogous relationships and were consistent with recent species trees. Based on the myosin inventories and the phylogenetic trees, we have identified duplications of the entire myosin motor protein family at timings consistent with 23 WGDs, that had been reported before. We also predict 6 WGDs based on further protein family duplications. Notably, the myosin data support the two recently reported WGDs in the common ancestor of all extant angiosperms. We predict single WGDs in the Manihot esculenta and Nicotiana benthamiana lineages, two WGDs for Linum usitatissimum and Phoenix dactylifera, and a triplication or two WGDs for Gossypium raimondii. Our data show another myosin duplication in the ancestor of the angiosperms that could be either the result of a single gene duplication or a remnant of a WGD.

Conclusions

We have shown that the myosin inventories in angiosperms retain evidence of numerous WGDs that happened throughout plant evolution. In contrast to other protein families, many myosins are still present in extant species. They are closely related and have similar domain architectures, and their phylogenetic grouping follows the genome duplications. Because of its broad taxonomic sampling the dataset provides the basis for reliable future identification of further whole genome duplications.

A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.)

Cicer arietinum L. (chickpea) is the third most important food legume crop. We have generated the draft sequence of a desi-type chickpea genome using next-generation sequencing platforms, bacterial artificial chromosome end sequences and a genetic map. The 520-Mb assembly covers 70% of the predicted 740-Mb genome length, and more than 80% of the gene space. Genome analysis predicts the presence of 27,571 genes and 210 Mb as repeat elements. The gene expression analysis performed using 274 million RNA-Seq reads identified several tissue-specific and stress-responsive genes. Although segmental duplicated blocks are observed, the chickpea genome does not exhibit any indication of recent whole-genome duplication. Nucleotide diversity analysis provides an assessment of a narrow genetic base within the chickpea cultivars. We have developed a resource for genetic markers by comparing the genome sequences of one wild and three cultivated chickpea genotypes. The draft genome sequence is expected to facilitate genetic enhancement and breeding to develop improved chickpea varieties.

Comparative sequence analysis of nitrogen fixation-related genes in six legumes

Legumes play an important role as food and forage crops in international agriculture especially in developing countries. Legumes have a unique biological process called nitrogen fixation (NF) by which they convert atmospheric nitrogen to ammonia. Although legume genomes have undergone polyploidization, duplication and divergence, NF-related genes, because of their essential functional role for legumes, might have remained conserved. To understand the relationship of divergence and evolutionary processes in legumes, this study analyzes orthologs and paralogs for selected 20 NF-related genes by using comparative genomic approaches in six legumes i.e., Medicago truncatula (Mt), Cicer arietinum, Lotus japonicus, Cajanus cajan (Cc), Phaseolus vulgaris (Pv), and Glycine max (Gm). Subsequently, sequence distances, numbers of synonymous substitutions per synonymous site (Ks) and non-synonymous substitutions per non-synonymous site (Ka) between orthologs and paralogs were calculated and compared across legumes. These analyses suggest the closest relationship between Gm and Cc and the highest distance between Mt and Pv in six legumes. Ks proportional plots clearly showed ancient genome duplication in all legumes, whole genome duplication event in Gm and also speciation pattern in different legumes. This study also reports some interesting observations e.g., no peak at Ks 0.4 in Gm-Gm, location of two independent genes next to each other in Mt and low Ks values for outparalogs for three genes as compared to other 12 genes. In summary, this study underlines the importance of NF-related genes and provides important insights in genome organization and evolutionary aspects of six legume species analyzed.

Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

Three linkage maps of flax (Linum usitatissimum L.) were constructed from populations CDC Bethune/Macbeth, E1747/Viking and SP2047/UGG5-5 containing between 385 and 469 mapped markers each. The first consensus map of flax was constructed incorporating 770 markers based on 371 shared markers including 114 that were shared by all three populations and 257 shared between any two populations. The 15 linkage group map corresponds to the haploid number of chromosomes of this species. The marker order of the consensus map was largely collinear in all three individual maps but a few local inversions and marker rearrangements spanning short intervals were observed. Segregation distortion was present in all linkage groups which contained 1-52 markers displaying non-Mendelian segregation. The total length of the consensus genetic map is 1,551 cM with a mean marker density of 2.0 cM. A total of 670 markers were anchored to 204 of the 416 fingerprinted contigs of the physical map corresponding to ~274 Mb or 74 % of the estimated flax genome size of 370 Mb. This high resolution consensus map will be a resource for comparative genomics, genome organization, evolution studies and anchoring of the whole genome shotgun sequence.