A systematic search for positive selection in higher plants (Embryophytes)

Identification and characterization of FBA genes in moso bamboo reveals PeFBA8 related to photosynthetic carbon metabolism

Fructose 1,6-bisphosphate aldolase (FBA) is a pivotal enzyme, which plays a critical role in fixing CO2 through the process of in the Calvin cycle. In this study, a comprehensive exploration of the FBA family genes in moso bamboo (Phyllostachys edulis) was conducted by the bioinformatics and biological analyses. A total of nine FBA genes (PeFBA1-PeFBA9) were identified in the moso bamboo genome. The expression patterns of PeFBAs across diverse tissues of moso bamboo suggested that they have multifaceted functionality. Notably, PeFBA8 might play an important role in regulating photosynthetic carbon metabolism. Co-expression and cis-element analyses demonstrated that PeFBA8 was regulated by a photosynthetic regulatory transcription factor (PeGLK1), which was confirmed by yeast one-hybrid and dual-luciferase assays. In-planta gene editing analysis revealed that the edited PeFBA8 mutants displayed compromised photosynthetic functionality, characterized by reduced electron transport rate and impaired photosystem I, leading to decreased photosynthesis rate overall, compared to the unedited control. The recombinant protein of PeFBA8 from prokaryotic expression exhibited enzymatic catalytic function. The findings suggest that the expression of PeFBA8 can affect photosynthetic efficiency of moso bamboo leaves, which underlines the potential of leveraging PeFBA8's regulatory mechanism to breed bamboo varieties with enhanced carbon fixation capability.

Genome-Wide Identification and Expression Analysis of the Class III Peroxidase Gene Family under Abiotic Stresses in Litchi (Litchi chinensis Sonn.)

Class III peroxidases (CIII PRXs) are plant-specific enzymes with high activity that play key roles in the catalysis of oxidation-reduction reactions. In plants, CIII PRXs can reduce hydrogen peroxide to catalyze oxidation-reduction reactions, thereby affecting plant growth, development, and stress responses. To date, no systematic analysis of the CIII PRX gene family in litchi (Litchi chinensis Sonn.) has been documented, although the genome has been reported. In this study, a total of 77 CIII PRX (designated LcPRX) gene family members were predicted in the litchi genome to provide a reference for candidate genes in the responses to abiotic stresses during litchi growth and development. All of these LcPRX genes had different numbers of highly conserved PRX domains and were unevenly distributed across fifteen chromosomes. They were further clustered into eight clades using a phylogenetic tree, and almost every clade had its own unique gene structure and motif distribution. Collinearity analysis confirmed that there were eleven pairs of duplicate genes among the LcPRX members, and segmental duplication (SD) was the main driving force behind the LcPRX gene expansion. Tissue-specific expression profiles indicated that the expression levels of all the LcPRX family members in different tissues of the litchi tree were significantly divergent. After different abiotic stress treatments, quantitative real-time PCR (qRT-PCR) analysis revealed that the LcPRX genes responded to various stresses and displayed differential expression patterns. Physicochemical properties, transmembrane domains, subcellular localization, secondary structures, and cis-acting elements were also analyzed. These findings provide insights into the characteristics of the LcPRX gene family and give valuable information for further elucidating its molecular function and then enhancing abiotic stress tolerance in litchi through molecular breeding.

Molecular genetic variability of Cryphonectria hypovirus 1 associated with Cryphonectria parasitica in South Tyrol (northern Italy)

Cryphonectria hypovirus 1 (CHV-1) has been widely studied and used as a biocontrol agent because of its ability to infect the chestnut blight fungus, Cryphonectria parasitica, and to reduce its virulence. Knowledge about the hypovirus, its presence, and diversity is completely lacking in South Tyrol (northern Italy), which may obstruct biocontrol measures for chestnut blight based on CHV-1. This work aimed to study the occurrence of CHV-1 infecting C. parasitica in South Tyrol and to perform a genetic characterization of the hypovirus. In South Tyrol, CHV-1 was found to occur in 29.2% of the fungal isolates investigated, varying in frequency between different regions and chestnut stands. Twenty-three haplotypes based on partial cDNA (complementary DNA) sequences of open reading frame (ORF)-A and 30 haplotypes based on partial cDNA sequences of ORF-B were identified among 47 and 56 hypovirulent fungal isolates, respectively. Phylogenetic analysis showed that all the haplotypes belonged to the Italian subtype of CHV-1 and that they were closely related to the populations of Italy, Switzerland, Croatia and Slovenia. Evidence of recombination was not found in the sequences and point mutations were the main source of diversity. Overall, this study indicated that the prevalence of CHV-1 in South Tyrol is low compared to many other central and western European populations and determined a need to actively impose biocontrol measures. Using sequence analysis, we identified some variants of interest of CHV-1 that should be studied in detail for their potential use in biocontrol.

Genome-wide identification, evolution, and expression analysis of the NAC gene family in chestnut (Castanea mollissima)

The NAC gene family is one of the most important transcription factor families specific to plants, responsible for regulating many biological processes, including development, stress response, and signal transduction. However, it has not yet been characterized in chestnut, an important nut tree species. Here, we identified 115 CmNAC genes in the chestnut genome, which were divided into 16 subgroups based on the phylogenetic analysis. Numerous cis-acting elements related to auxin, gibberellin, and abscisic acid were identified in the promoter region of CmNACs, suggesting that they play an important role in the growth and development of chestnut. The results of the collinear analysis indicated that dispersed duplication and whole-genome-duplication were the main drivers of CmNAC gene expansion. RNA-seq data of developmental stages of chestnut nut, bud, and ovule revealed the expression patterns of CmNAC genes. Additionally, qRT-PCR experiments were used to verify the expression levels of some CmNAC genes. The comprehensive analysis of the above results revealed that some CmNAC members may be related to chestnut bud and nut development, as well as ovule fertility. The systematic analysis of this study will help to increase understanding of the potential functions of the CmNAC genes in chestnut growth and development.

WRKY transcription factors in passion fruit analysis reveals key PeWRKYs involved in abiotic stress and flavonoid biosynthesis

WRKY transcription factors (TFs) are a superfamily of regulators involved in plant responses to pathogens and abiotic stress. Passion fruit is famous for its unique flavor and nutrient-rich juice, but its growth is limited by environmental factors and pathogens. In this study, 55 WRKY genes were identified from the Passiflora edulis genome. The structure and evolutionary characteristics of PeWRKYs were analyzed using a bioinformatics approach. PeWRKYs were classified into seven subgroups (I, IIa, IIb, IIc, IId, IIe, III) according to their homologs in Arabidopsis thaliana. Group IIa PeWRKY48 gene was highly up-regulated under cold stress by RNA expression analysis, and transgenic PeWRKY48 in yeast and Arabidopsis showed resistance exposure to cold, salt, and drought stress. Metabolome and transcriptome co-expression analysis of two different disease resistance genotypes of P. edulis identified PeWRKY30 as a key TF co-expressed with flavonoid accumulation in yellow fruit P. edulis, which may contribute to biotic or abiotic resistance. The qRT-PCR verified the expression of key genes in different tissues of P. edulis and in different species of Passiflora. This study provides a set of WRKY candidate genes that will facilitate the genetic improvement of disease and abiotic tolerance in passion fruit.

Genome-wide identification and expression analysis of xyloglucan endotransglucosylase/hydrolase genes family in Salicaceae during grafting

BACKGROUND

Poplar (Populus cathayana)and willow (Salix rehderiana) are important fast-growing trees in China. Grafting plays an important role in improving plant stress resistance and construction of ornamental plants. It is found that willow scions grafted onto poplar rootstocks can form ornamental plants. However, this grafted combination has a low survival rate. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in the healing process of grafts.

RESULTS

A total of 38 PtrXTHs and 32 SpuXTHs were identified in poplar and willow respectively, and were classified into three subfamilies. Tandem duplication was the main reason for the expansion of the PtrXTHs. Grafting treatment and Quantitative real time PCR (RT-qPCR) analysis revealed that five XTH genes differentially expressed between self-grafted and reciprocal grafted combinations. Specifically, the high expression levels of SrXTH16, SrXTH17, SrXTH25, PcXTH22 and PcXTH17 may contribute to the high survival rate of the grafted combination with willow scion and poplar rootstock. Subcellular localization identified that the SrXTH16, SrXTH17, SrXTH25, PcXTH17 and PcXTH22 proteins were located on the cell walls. Transcription factors (NAC, MYB and DOF) may regulate the five XTH genes.

CONCLUSIONS

This study provides a new understanding of the roles of PcXTH and SrXTH genes and their roles in grafting. Our results will give some hints to explore the molecular mechanisms of PcXTH and SrXTH genes involved in grafting in the future.

Genome-Wide Analysis of Q-Type C2H2 ZFP Genes in Response to Biotic and Abiotic Stresses in Sugar Beet

A plant's Q-type C2H2-type ZFP plays key roles in plant growth and development and responses to biotic and abiotic stresses. Sugar beet (Beta vulgaris L.) is an important crop for sugar production. Salt stress and viral infection significantly reduce the root yield and sugar content of sugar beet. However, there is a lack of comprehensive genome-wide analyses of Q-type C2H2 ZFPs and their expression patterns in sugar beet under stress. In this study, 35 sugar beet Q-type C2H2 ZFPs (BvZFPs) containing at least one conserved "QALGGH" motif were identified via bioinformatics techniques using TBtools software. According to their evolutionary relationship, the BvZFPs were classified into five subclasses. Within each subclass, the physicochemical properties and motif compositions showed strong similarities. A Ka/Ks analysis indicated that the BvZFPs were conserved during evolution. Promoter cis-element analysis revealed that most BvZFPs are associated with elements related to phytohormone, biotic or abiotic stress, and plant development. The expression data showed that the BvZFPs in sugar beet are predominantly expressed in the root. In addition, BvZFPs are involved in the response to abiotic and biotic stresses, including salt stress and viral infection. Overall, these results will extend our understanding of the Q-type C2H2 gene family and provide valuable information for the biological breeding of sugar beet against abiotic and biotic stresses in the future.

Genome-wide identification and expression analysis of the Auxin-Response factor (ARF) gene family in Medicago sativa under abiotic stress

BACKGROUND

Alfalfa (Medicago sativa) is the most widely planted legume forage and one of the most economically valuable crops in the world. The periodic changes in its growth and development and abiotic stress determine its yield and economic benefits. Auxin controls many aspects of alfalfa growth by regulating gene expression, including organ differentiation and stress response. Auxin response factors (ARF) are transcription factors that play an essential role in auxin signal transduction and regulate the expression of auxin-responsive genes. However, the function of ARF transcription factors is unclear in autotetraploid-cultivated alfalfa.

RESULT

A total of 81 ARF were identified in the alfalfa genome in this study. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, identifying that ARF genes are mainly involved in transcriptional regulation and plant hormone signal transduction pathways. Phylogenetic analysis showed that MsARF was divided into four clades: I, II, III, and IV, each containing 52, 13, 7, and 9 genes, respectively. The promoter region of the MsARF gene contained stress-related elements, such as ABRE, TC-rich repeats, MBS, LTR. Proteins encoded by 50 ARF genes were localized in the nucleus without guide peptides, signal peptides, or transmembrane structures, indicating that most MsARF genes are not secreted or transported but only function in the nucleus. Protein structure analysis revealed that the secondary and tertiary structures of the 81 MsARF genes varied. Chromosomal localization analysis showed 81 MsARF genes were unevenly distributed on 25 chromosomes, with the highest distribution on chromosome 5. Furthermore, 14 segmental duplications and two sets of tandem repeats were identified. Expression analysis indicated that the MsARF was differentially expressed in different tissues and under various abiotic stressors. The quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that the expression profiles of 23 MsARF genes were specific to abiotic stresses such as drought, salt, high temperature, and low temperature, as well as tissue-specific and closely related to the duration of stress.

CONCLUSION

This study identified MsARF in the cultivated alfalfa genome based on the autotetraploid level, which GO, KEGG analysis, phylogenetic analysis, sequence characteristics, and expression pattern analysis further confirmed. Together, these findings provide clues for further investigation of MsARF functional verification and molecular breeding of alfalfa. This study provides a novel approach to systematically identify and characterize ARF transcription factors in autotetraploid cultivated alfalfa, revealing 23 MsARF genes significantly involved in response to various stresses.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Genome-Wide Identification of Expansin Gene Family and Their Response under Hormone Exposure in Ginkgo biloba L

Expansins are pH-dependent enzymatic proteins that irreversibly and continuously facilitate cell-wall loosening and extension. The identification and comprehensive analysis of Ginkgo biloba expansins (GbEXPs) are still lacking. Here, we identified and investigated 46 GbEXPs in Ginkgo biloba. All GbEXPs were grouped into four subgroups based on phylogeny. GbEXPA31 was cloned and subjected to a subcellular localization assay to verify our identification. The conserved motifs, gene organization, cis-elements, and Gene Ontology (GO) annotation were predicted to better understand the functional characteristics of GbEXPs. The collinearity test indicated segmental duplication dominated the expansion of the GbEXPA subgroup, and seven paralogous pairs underwent strong positive selection during expansion. A majority of GbEXPAs were mainly expressed in developing Ginkgo kernels or fruits in transcriptome and real-time quantitative PCR (qRT-PCR). Furthermore, GbEXLA4, GbEXLA5, GbEXPA5, GbEXPA6, GbEXPA8, and GbEXPA24 were inhibited under the exposure of abiotic stresses (UV-B and drought) and plant hormones (ABA, SA, and BR). In general, this study expanded our understanding for expansins in Ginkgo tissues' growth and development and provided a new basis for studying GbEXPs in response to exogenous phytohormones.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

Genome-wide identification of 14-3-3 gene family and characterization of their expression in developmental stages of Solanum tuberosum under multiple biotic and abiotic stress conditions

GF14 proteins are a family of conserved proteins involved in many cellular processes including transport, growth, metabolism, and stress response. However, only few reports are available regarding the 14-3-3 genes in potato. In this study, twelve 14-3-3 genes were detected in the potato genome. Based on their phylogenetic relationships, the StGF14 family members were categorized into two classes. Gene expression analysis demonstrated that StGF14h, StGF14a, and StGF14k had the highest gene expression, induced by abiotic and biotic stresses in all three tissues. The number of exons in 14-3-3 genes ranged from four to seven and most of these genes in the same subfamily had similar exon-intron patterns. The results of our study showed that the conserved motifs are similar in most of the proteins in each group. The intron-exon patterns and the composition of conserved motifs validated the 14-3-3 gene phylogenetic classification. According to the genome distribution results, 14-3-3 genes were located unevenly on the 12 Solanum tuberosum chromosomes. We find out 97 orthologous gene pairs between potato and Arabidopsis as well as 15 paralogous genes among potato genomes. Our results showed that GF-14 genes have an effective role in functional and molecular mechanisms in response to environmental stresses.

Systematic Identification and Expression Analysis of the Auxin Response Factor (ARF) Gene Family in Ginkgo biloba L

Auxin participates in various physiological and molecular response-related developmental processes and is a pivotal hormone that regulates phenotypic formation in plants. Auxin response factors (ARFs) are vital transcription factors that mediate downstream auxin signaling by explicitly binding to auxin-responsive genes' promoters. Here, to investigate the possible developmental regulatory functions of ARFs in Ginkgo biloba, through employing comprehensive bioinformatics, we recognized 15 putative GbARF members. Conserved domains and motifs, gene and protein structure, gene duplication, GO enrichment, transcriptome expression profiles, and qRT-PCR all showed that Group I and III members were highly conserved. Among them, GbARF10b and GbARF10a were revealed as transcriptional activators in the auxin response for the development of Ginkgo male flowers through sequences alignment, cis-elements analysis and GO annotation; the results were corroborated for the treatment of exogenous SA. Moreover, the GbARFs expansion occurred predominantly by segmental duplication, and most GbARFs have undergone purifying selection. The Ka/Ks ratio test identified the functional consistence of GbARF2a and GbARF2c, GbARF10b, and GbARF10a in tissue expression profiles and male flower development. In summary, our study established a new research basis for exploring Ginkgo GbARF members' roles in floral organ development and hormone response.

Genome-Wide Characterization of Superoxide Dismutase (SOD) Genes in Daucus carota: Novel Insights Into Structure, Expression, and Binding Interaction With Hydrogen Peroxide (H2O2) Under Abiotic Stress Condition

Superoxide dismutase (SOD) proteins are important antioxidant enzymes that help plants to grow, develop, and respond to a variety of abiotic stressors. SOD gene family has been identified in a number of plant species but not yet in Daucus carota. A total of 9 DcSOD genes, comprising 2 FeSODs, 2 MnSODs, and 5 Cu/ZnSODs, are identified in the complete genome of D. carota, which are dispersed in five out of nine chromosomes. Based on phylogenetic analysis, SOD proteins from D. carota were categorized into two main classes (Cu/ZnSODs and MnFeSODs). It was predicted that members of the same subgroups have the same subcellular location. The phylogenetic analysis was further validated by sequence motifs, exon-intron structure, and 3D protein structures, with each subgroup having a similar gene and protein structure. Cis-regulatory elements responsive to abiotic stresses were identified in the promoter region, which may contribute to their differential expression. Based on RNA-seq data, tissue-specific expression revealed that DcCSD2 had higher expression in both xylem and phloem. Moreover, DcCSD2 was differentially expressed in dark stress. All SOD genes were subjected to qPCR analysis after cold, heat, salt, or drought stress imposition. SODs are antioxidants and play a critical role in removing reactive oxygen species (ROS), including hydrogen peroxide (H2O2). DcSODs were docked with H2O2 to evaluate their binding. The findings of this study will serve as a basis for further functional insights into the DcSOD gene family.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Implications of divergence of methionine adenosyltransferase in archaea

Methionine adenosyltransferase (MAT) catalyzes the biosynthesis of S-adenosyl methionine from l-methionine and ATP. MAT enzymes are ancient, believed to share a common ancestor, and are highly conserved in all three domains of life. However, the sequences of archaeal MATs show considerable divergence compared with their bacterial and eukaryotic counterparts. Furthermore, the structural significance and functional significance of this sequence divergence are not well understood. In the present study, we employed structural analysis and ancestral sequence reconstruction to investigate archaeal MAT divergence. We observed that the dimer interface containing the active site (which is usually well conserved) diverged considerably between the bacterial/eukaryotic MATs and archaeal MAT. A detailed investigation of the available structures supports the sequence analysis outcome: The protein domains and subdomains of bacterial and eukaryotic MAT are more similar than those of archaea. Finally, we resurrected archaeal MAT ancestors. Interestingly, archaeal MAT ancestors show substrate specificity, which is lost during evolution. This observation supports the hypothesis of a common MAT ancestor for the three domains of life. In conclusion, we have demonstrated that archaeal MAT is an ideal system for studying an enzyme family that evolved differently in one domain compared with others while maintaining the same catalytic activity.

Characterisation of the class III peroxidase gene family in carrot taproots and its role in anthocyanin and lignin accumulation

Plant class III peroxidases (CIII Prxs) are involved in numerous essential plant life processes, such as plant development and differentiation, lignification and seed germination, and defence against pathogens. However, there is limited information about the structure-function relationships of Prxs in carrots. This study identified 75 carrot peroxidases (DcPrxs) and classified them into seven subgroups based on phylogenetic analysis. Gene structure analysis revealed that these DcPrxs had between one and eight introns, while conserved motif analysis showed a typical motif composition and arrangement for CIII Prx. In addition, eighteen tandem duplication events, but only eight segmental duplications, were identified among these DcPrxs, indicating that tandem duplication was the main contributor to the expansion of this gene family. Histochemical analyses showed that lignin was mainly localised in the cell walls of xylem, and Prx activity was determined in the epidermal region of taproots. The xylem always showed higher lignin concentration and lower Prx activity compared to the phloem in the taproots of both carrot cultivars. Combining these observations with RNA sequencing, some Prx genes were identified as candidate genes related to lignification and pigmentation. Three peroxidases (DcPrx30, DcPrx32, DcPrx62) were upregulated in the phloem of both genotypes. Carrot taproots are an attractive resource for natural food colourants and this study elucidated genome-wide insights of Prx for the first time, developing hypotheses concerning their involvement with lignin and anthocyanin in purple carrots. The findings provide an essential foundation for further studies of Prx genes in carrot, especially on pigmentation and lignification mechanisms.

Genome Characterization, Comparison and Phylogenetic Analysis of Complete Mitochondrial Genome of Evolvulus alsinoides Reveals Highly Rearranged Gene Order in Solanales

Mitogenome sequencing provides an understanding of the evolutionary mechanism of mitogenome formation, mechanisms driving plant gene order, genome structure, and migration sequences. Data on the mitochondrial genome for family Convolvulaceae members is lacking. E. alsinoides, also known as shankhpushpi, is an important medicinal plant under the family Convolvulaceae, widely used in the Ayurvedic system of medicine. We identified the mitogenome of E. alsinoides using the Illumina mate-pair sequencing platform, and annotated using bioinformatics approaches in the present study. The mitogenome of E. alsinoides was 344184 bp in length and comprised 46 unique coding genes, including 31 protein-coding genes (PCGs), 12 tRNA genes, and 3 rRNA genes. The secondary structure of tRNAs shows that all the tRNAs can be folded into canonical clover-leaf secondary structures, except three trnW, trnG, and trnC. Measurement of the skewness of the nucleotide composition showed that the AT and GC skew is positive, indicating higher A's and G's in the mitogenome of E. alsinoides. The Ka/Ks ratios of 11 protein-coding genes (atp1, ccmC, cob, cox1, rps19, rps12, nad3, nad9, atp9, rpl5, nad4L) were <1, indicating that these genes were under purifying selection. Synteny and gene order analysis were performed to identify homologous genes among the related species. Synteny blocks representing nine genes (nad9, nad2, ccmFc, nad1, nad4, nad5, matR, cox1, nad7) were observed in all the species of Solanales. Gene order comparison showed that a high level of gene rearrangement has occurred among all the species of Solanales. The mitogenome data obtained in the present study could be used as the Convolvulaceae family representative for future studies, as there is no complex taxonomic history associated with this plant.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Genome-Wide Identification and Comparative Analysis of MYB Transcription Factor Family in Musa acuminata and Musa balbisiana

MYB transcription factors (TFs) make up one of the most important TF families in plants. These proteins play crucial roles in processes related to development, metabolism, and stimulus-response; however, very few studies have been reported for the characterization of MYB TFs from banana. The current study identified 305 and 251 MYB genes from Musa acuminata and Musa balbisiana, respectively. Comprehensive details of MYBs are reported in terms of gene structure, protein domain, chromosomal localization, phylogeny, and expression patterns. Based on the exon-intron arrangement, these genes were classified into 12 gene models. Phylogenetic analysis of MYBs involving both species of banana, Oryza sativa, and Arabidopsis thaliana distributed these genes into 27 subfamilies. This highlighted not only the conservation, but also the gain/loss of MYBs in banana. Such genes are important candidates for future functional investigations. The MYB genes in both species exhibited a random distribution on chromosomes with variable densities. Estimation of gene duplication events revealed that segmental duplications represented the major factor behind MYB gene family expansion in banana. Expression profiles of MYB genes were also explored for their potential involvement in acetylene response or development. Collectively, the current comprehensive analysis of MYB genes in both species of banana will facilitate future functional studies.

Thermospermine Synthase (ACL5) and Diamine Oxidase (DAO) Expression Is Needed for Zygotic Embryogenesis and Vascular Development in Scots Pine

Unlike in flowering plants, the detailed roles of the enzymes in the polyamine (PA) pathway in conifers are poorly known. We explored the sequence conservation of the PA biosynthetic genes and diamine oxidase (DAO) in conifers and flowering plants to reveal the potential functional diversification of the enzymes between the plant lineages. The expression of the genes showing different selective constraints was studied in Scots pine zygotic embryogenesis and early seedling development. We found that the arginine decarboxylase pathway is strongly preferred in putrescine production in the Scots pine as well as generally in conifers and that the reduced use of ornithine decarboxylase (ODC) has led to relaxed purifying selection in ODC genes. Thermospermine synthase (ACL5) genes evolve under strong purifying selection in conifers and the DAO gene is also highly conserved in pines. In developing Scots pine seeds, the expression of both ACL5 and DAO increased as embryogenesis proceeded. Strong ACL5 expression was present in the procambial cells of the embryo and in the megagametophyte cells destined to die via morphologically necrotic cell death. Thus, the high sequence conservation of ACL5 genes in conifers may indicate the necessity of ACL5 for both embryogenesis and vascular development. Moreover, the result suggests the involvement of ACL5 in morphologically necrotic cell death and supports the view of the genetic regulation of necrosis in Scots pine embryogenesis and in plant development. DAO transcripts were located close to the cell walls and between the walls of adjacent cells in Scots pine zygotic embryos and in the roots of young seedlings. We propose that DAO, in addition to the role in Put oxidation for providing H2O2 during the cell-wall structural processes, may also participate in cell-to-cell communication at the mRNA level. To conclude, our findings indicate that the PA pathway of Scots pines possesses several special functional characteristics which differ from those of flowering plants.

The expression pattern of the Pho1a genes encoding plastidic starch phosphorylase correlates with the degradation of starch during fruit ripening in green-fruited and red-fruited tomato species

Genes encoding plastidic starch phosphorylase Pho1a were identified in 10 tomato species (Solanum section Lycopersicon). Pho1a genes showed higher variability in green-fruited than in red-fruited tomato species, but had an extremely low polymorphism level compared with other carbohydrate metabolism genes and an unusually low ratio of intron to exon single nucleotide polymorphisms (SNPs). In red-fruited species, Pho1a was expressed in all analysed tissues, including fruit at different developmental stages, with the highest level in mature green fruit, which is strong sink organ importing sucrose and accumulating starch. In green-fruited species Solanum peruvianum and Solanum arcanum, the Pho1a expression level was similar in mature green and ripe fruit, whereas in Solanum chmielewskii, it was higher in ripe fruit, and in Solanum habrochaites, the dynamics of fruit-specific Pho1a expression was similar to that in red-fruited tomatoes. During fruit development, in red-fruited Solanum lycopersicum, sucrose level was low, the monosaccharide content increased; in green-fruited S. peruvianum, the sucrose concentration increased and those of monosaccharides decreased. In both species, the starch content and Pho1a expression were downregulated. The evolutionary topology based on Pho1a sequences was consistent with the current division of tomatoes into red-fruited and green-fruited species, except for S. habrochaites.

Structure of proteins: Evolution with unsolved mysteries

Evolution of macromolecules could be considered as a milestone in the history of life. Nucleic acids are the long stretches of nucleotides that contain all the possible codes and information of life. On the other hand, proteins are their actual translated outcomes, or reflections of modifications in their structure that have occurred at a slow, but steady rate over a very long period of evolution. Over the years of research, biophysicists, biochemists, molecular and structural biologists have unfurled several layers of the structural convolutions in these chemical molecules; however evolutionists look over their structures through a different prism, which may or may not coincide with others. There remains a need to outline several well-known, but less discussed features of protein structures, like intrinsically disordered states, degron signals and different types of ubiquitin chains providing degradation signals, which help the cellular proteolytic machinery to identify and target the proteins towards degradation pathways. There are several important factors, which are critical for folding of proteins into their native three-dimensional conformations by the cytoplasmic chaperones; but in real time how the chaperones fold the newly synthesized polypeptide sequences into a particular three-dimensional shape within a fraction of second is still a mystery for biologists as well as mathematicians. Multiple similar unsolved or unaddressed questions need to be addressed in detail so that future line of research can dig deeper into the finer details of these structures of the proteins.

A Spotlight on Rad52 in Cyanidiophytina (Rhodophyta): A Relic in Algal Heritage

The RADiation sensitive52 (RAD52) protein catalyzes the pairing between two homologous DNA sequences' double-strand break repair and meiotic recombination, mediating RAD51 loading onto single-stranded DNA ends, and initiating homologous recombination and catalyzing DNA annealing. This article reports the characterization of RAD52 homologs in the thermo-acidophilic Cyanidiophyceae whose genomes have undergone extensive sequencing. Database mining, phylogenetic inference, prediction of protein structure and evaluation of gene expression were performed in order to determine the functionality of the RAD52 protein in Cyanidiophyceae. Its current function in Cyanidiophytina could be related to stress damage response for thriving in hot and acidic environments as well as to the genetic variability of these algae, in which, conversely to extant Rhodophyta, sexual mating was never observed.

Structural and Evolutionary Insights within the Polysaccharide Deacetylase Gene Family of Bacillus anthracis and Bacillus cereus

Functional and folding constraints impose interdependence between interacting sites along the protein chain that are envisaged through protein sequence evolution. Studying the influence of structure in phylogenetic models requires detailed and reliable structural models. Polysaccharide deacetylases (PDAs), members of the carbohydrate esterase family 4, perform mainly metal-dependent deacetylation of O- or N-acetylated polysaccharides such as peptidoglycan, chitin and acetylxylan through a conserved catalytic core termed the NodB homology domain. Genomes of Bacillus anthracis and its relative Bacillus cereus contain multiple genes of putative or known PDAs. A comparison of the functional domains of the recently determined PDAs from B. anthracis and B. cereus and multiple amino acid and nucleotide sequence alignments and phylogenetic analysis performed on these closely related species showed that there were distinct differences in binding site formation, despite the high conservation on the protein sequence, the folding level and the active site assembly. This may indicate that, subject to biochemical verification, the binding site-forming sequence fragments are under functionally driven evolutionary pressure to accommodate and recognize distinct polysaccharide residues according to cell location, use, or environment. Finally, we discuss the suggestion of the paralogous nature of at least two genes of B. anthracis, ba0330 and ba0331, via specific differences in gene sequence, protein structure, selection pressure and available localization patterns. This study may contribute to understanding the mechanisms under which sequences evolve in their structures and how evolutionary processes enable structural variations.

The Adaptive Evolution Database (TAED): A New Release of a Database of Phylogenetically Indexed Gene Families from Chordates

With the large collections of gene and genome sequences, there is a need to generate curated comparative genomic databases that enable interpretation of results in an evolutionary context. Such resources can facilitate an understanding of the co-evolution of genes in the context of a genome mapped onto a phylogeny, of a protein structure, and of interactions within a pathway. A phylogenetically indexed gene family database, the adaptive evolution database (TAED), is presented that organizes gene families and their evolutionary histories in a species tree context. Gene families include alignments, phylogenetic trees, lineage-specific dN/dS ratios, reconciliation with the species tree to enable both the mapping and the identification of duplication events, mapping of gene families onto pathways, and mapping of amino acid substitutions onto protein structures. In addition to organization of the data, new phylogenetic visualization tools have been developed to aid in interpreting the data that are also available, including TreeThrasher and TAED Tree Viewer. A new resource of gene families organized by species and taxonomic lineage promises to be a valuable comparative genomics database for molecular biologists, evolutionary biologists, and ecologists. The new visualization tools and database framework will be of interest to both evolutionary biologists and bioinformaticians.

Adaptive sequence evolution is driven by biotic stress in a pair of orchid species (Dactylorhiza) with distinct ecological optima

The orchid family is the largest in the angiosperms, but little is known about the molecular basis of the significant variation they exhibit. We investigate here the transcriptomic divergence between two European terrestrial orchids, Dactylorhiza incarnata and Dactylorhiza fuchsii, and integrate these results in the context of their distinct ecologies that we also document. Clear signals of lineage-specific adaptive evolution of protein-coding sequences are identified, notably targeting elements of biotic defence, including both physical and chemical adaptations in the context of divergent pools of pathogens and herbivores. In turn, a substantial regulatory divergence between the two species appears linked to adaptation/acclimation to abiotic conditions. Several of the pathways affected by differential expression are also targeted by deviating post-transcriptional regulation via sRNAs. Finally, D. incarnata appears to suffer from insufficient sRNA control over the activity of RNA-dependent DNA polymerase, resulting in increased activity of class I transposable elements and, over time, in larger genome size than that of D. fuchsii. The extensive molecular divergence between the two species suggests significant genomic and transcriptomic shock in their hybrids and offers insights into the difficulty of coexistence at the homoploid level. Altogether, biological response to selection, accumulated during the history of these orchids, appears governed by their microenvironmental context, in which biotic and abiotic pressures act synergistically to shape transcriptome structure, expression and regulation.

Evolution of jasmonate biosynthesis and signaling mechanisms

Jasmonates are phytohormones that modulate a wide spectrum of plant physiological processes, especially defense against herbivores and necrotrophs. The molecular mechanisms of jasmonate biosynthesis and signaling have been well characterized in model plants. In this review, we provide an in-depth analysis and overview of the origin and evolution of the jasmonate biosynthesis and signaling pathways. Furthermore, we discuss the striking parallels between jasmonate and auxin signaling mechanisms, which reveals a common ancestry of these signaling mechanisms. Finally, we highlight the importance of studying jasmonate biosynthesis and signaling in lower plants.

Characterization of parasite-specific indels and their proposed relevance for selective anthelminthic drug targeting

Insertions and deletions (indels) are important sequence variants that are considered as phylogenetic markers that reflect evolutionary adaptations in different species. In an effort to systematically study indels specific to the phylum Nematoda and their structural impact on the proteins bearing them, we examined over 340,000 polypeptides from 21 nematode species spanning the phylum, compared them to non-nematodes and identified indels unique to nematode proteins in more than 3000 protein families. Examination of the amino acid composition revealed uneven usage of amino acids for insertions and deletions. The amino acid composition and cost, along with the secondary structure constitution of the indels, were analyzed in the context of their biological pathway associations. Species-specific indels could enable indel-based targeting for drug design in pathogens/parasites. Therefore, we screened the spatial locations of the indels in the parasite's protein 3D structures, determined the location of the indel and identified potential unique drug targeting sites. These indels could be confirmed by RNA-Seq data. Examples are presented illustrating the close proximity of some indels to established small-molecule binding pockets that can potentially facilitate selective targeting to the parasites and bypassing their host, thus reducing or eliminating the toxicity of the potential drugs. This study presents an approach for understanding the adaptation of pathogens/parasites at a molecular level, and outlines a strategy to identify such nematode-selective targets that remain essential to the organism. With further experimental characterization and validation, it opens a possible channel for the development of novel treatments with high target specificity, addressing both host toxicity and resistance concerns.

Positive selection in the leucine-rich repeat domain of Gro1 genes in Solanum species

In pathogen resistant plants, solvent-exposed residues in the leucine-rich repeat (LRR) proteins are thought to mediate resistance by recognizing plant pathogen elicitors. In potato, the gene Gro1-4 confers resistance to Globodera rostochiensis. The investigation of variability in different copies of this gene represents a good model for the verification of positive selection mechanisms. Two datasets of Gro1 LRR sequences were constructed, one derived from the Gro1-4 gene, belonging to different cultivated and wild Solanum species, and the other belonging to paralogues of a resistant genotype. Analysis of nonsynonymous to synonymous substitution rates (K(a)/K(s)) highlighted 14 and six amino acids with K(a)/K(s) >1 in orthologue and paralogue datasets, respectively. Selection analysis revealed that the leucine-rich regions accumulate variability in a very specific way, and we found that some combinations of amino acids in these sites might be involved in pathogen recognition. The results confirm previous studies on positive selection in the LRR domain of R protein in Arabidopsis and other model plants and extend these to wild Solanum species. Moreover, positively selected sites in the Gro1 LRR domain show that coevolution mainly occurred in two regions on the internal surface of the three-dimensional horseshoe structure of the domain, albeit with different evolutionary forces between paralogues and orthologues.

Transcriptome Sequencing of Lima Bean (Phaseolus lunatus) to Identify Putative Positive Selection in Phaseolus and Legumes

The identification of genes under positive selection is a central goal of evolutionary biology. Many legume species, including Phaseolus vulgaris (common bean) and Phaseolus lunatus (lima bean), have important ecological and economic value. In this study, we sequenced and assembled the transcriptome of one Phaseolus species, lima bean. A comparison with the genomes of six other legume species, including the common bean, Medicago, lotus, soybean, chickpea, and pigeonpea, revealed 15 and 4 orthologous groups with signatures of positive selection among the two Phaseolus species and among the seven legume species, respectively. Characterization of these positively selected genes using Non redundant (nr) annotation, gene ontology (GO) classification, GO term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes are mostly involved in thylakoids, photosynthesis and metabolism. This study identified genes that may be related to the divergence of the Phaseolus and legume species. These detected genes are particularly good candidates for subsequent functional studies.

Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index

Swarna is a popular cultivated indica rice variety with low glycemic index (GI) but its genetic basis is not known. The whole genome of Swarna was sequenced using Illumina’s paired-end technology, and the reads were mapped to the Nipponbare reference genome. Overall, 65,984 non-synonymous SNPs were identified in 20,350 genes, and in silico analysis predicted that 4,847 of them in 2,214 genes may have deleterious effect on protein functions. Polymorphisms were found in all the starch biosynthesis genes, except the gene for branching enzyme IIa. It was found that T/G SNP at position 246, ‘A’ at position 2,386, and ‘C’ at position 3,378 in the granule bound starch synthase I gene, and C/T SNP at position 1,188 in the glucose-6-phosphate translocator gene may contribute to the low GI phenotype in Swarna. All these variants were also found in the genome of another low GI indica rice variety from Columbia, Fedearroz 50. The whole genome analysis of Swarna helped to understand the genetic basis of GI in rice, which is a complex trait involving multiple factors.

Darwin and Fisher meet at biotech: on the potential of computational molecular evolution in industry

Background

Today computational molecular evolution is a vibrant research field that benefits from the availability of large and complex new generation sequencing data – ranging from full genomes and proteomes to microbiomes, metabolomes and epigenomes. The grounds for this progress were established long before the discovery of the DNA structure. Specifically, Darwin’s theory of evolution by means of natural selection not only remains relevant today, but also provides a solid basis for computational research with a variety of applications. But a long-term progress in biology was ensured by the mathematical sciences, as exemplified by Sir R. Fisher in early 20th century. Now this is true more than ever: The data size and its complexity require biologists to work in close collaboration with experts in computational sciences, modeling and statistics.

Results

Natural selection drives function conservation and adaptation to emerging pathogens or new environments; selection plays key role in immune and resistance systems. Here I focus on computational methods for evaluating selection in molecular sequences, and argue that they have a high potential for applications. Pharma and biotech industries can successfully use this potential, and should take the initiative to enhance their research and development with state of the art bioinformatics approaches.

Conclusions

This review provides a quick guide to the current computational approaches that apply the evolutionary principles of natural selection to real life problems – from drug target validation, vaccine design and protein engineering to applications in agriculture, ecology and conservation.

Comparative transcriptome analysis within the Lolium/Festuca species complex reveals high sequence conservation

Background

The Lolium-Festuca complex incorporates species from the Lolium genera and the broad leaf fescues, both belonging to the subfamily Pooideae. This subfamily also includes wheat, barley, oat and rye, making it extremely important to world agriculture. Species within the Lolium-Festuca complex show very diverse phenotypes, and many of them are related to agronomically important traits. Analysis of sequenced transcriptomes of these non-model species may shed light on the molecular mechanisms underlying this phenotypic diversity.

Results

We have generated de novo transcriptome assemblies for four species from the Lolium-Festuca complex, ranging from 52,166 to 72,133 transcripts per assembly. We have also predicted a set of proteins and validated it with a high-confidence protein database from three closely related species (H. vulgare, B. distachyon and O. sativa). We have obtained gene family clusters for the four species using OrthoMCL and analyzed their inferred phylogenetic relationships. Our results indicate that VRN2 is a candidate gene for differentiating vernalization and non-vernalization types in the Lolium-Festuca complex. Grouping of the gene families based on their BLAST identity enabled us to divide ortholog groups into those that are very conserved and those that are more evolutionarily relaxed. The ratio of the non-synonumous to synonymous substitutions enabled us to pinpoint protein sequences evolving in response to positive selection. These proteins may explain some of the differences between the more stress tolerant Festuca, and the less stress tolerant Lolium species.

Conclusions

Our data presents a comprehensive transcriptome sequence comparison between species from the Lolium-Festuca complex, with the identification of potential candidate genes underlying some important phenotypical differences within the complex (such as VRN2). The orthologous genes between the species have a very high %id (91,61%) and the majority of gene families were shared for all of them. It is likely that the knowledge of the genomes will be largely transferable between species within the complex.

Electronic supplementary material

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

Phylogenomic analyses of nuclear genes reveal the evolutionary relationships within the BEP clade and the evidence of positive selection in Poaceae

BEP clade of the grass family (Poaceae) is composed of three subfamilies, i.e. Bambusoideae, Ehrhartoideae, and Pooideae. Controversies on the phylogenetic relationships among three subfamilies still persist in spite of great efforts. However, previous evidence was mainly provided from plastid genes with only a few nuclear genes utilized. Given different evolutionary histories recorded by plastid and nuclear genes, it is indispensable to uncover their relationships based on nuclear genes. Here, eleven species with whole-sequenced genome and six species with transcriptomic data were included in this study. A total of 121 one-to-one orthologous groups (OGs) were identified and phylogenetic trees were reconstructed by different tree-building methods. Genes which might have undergone positive selection and played important roles in adaptive evolution were also investigated from 314 and 173 one-to-one OGs in two bamboo species and 14 grass species, respectively. Our results support the ((B, P) E) topology with high supporting values. Besides, our findings also indicate that 24 and nine orthologs with statistically significant evidence of positive selection are mainly involved in abiotic and biotic stress response, reproduction and development, plant metabolism and enzyme etc. from two bamboo species and 14 grass species, respectively. In summary, this study demonstrates the power of phylogenomic approach to shed lights on the evolutionary relationships within the BEP clade, and offers valuable insights into adaptive evolution of the grass family.

Identification of positive selection in disease response genes within members of the Poaceae

Millions of years of coevolution between plants and pathogens can leave footprints on their genomes and genes involved on this interaction are expected to show patterns of positive selection in which novel, beneficial alleles are rapidly fixed within the population. Using information about upregulated genes in maize during Colletotrichum graminicola infection and resources available in the Phytozome database, we looked for evidence of positive selection in the Poaceae lineage, acting on protein coding sequences related with plant defense. We found six genes with evidence of positive selection and another eight with sites showing episodic selection. Some of them have already been described as evolving under positive selection, but others are reported here for the first time including genes encoding isocitrate lyase, dehydrogenases, a multidrug transporter, a protein containing a putative leucine-rich repeat and other proteins with unknown functions. Mapping positively selected residues onto the predicted 3-D structure of proteins showed that most of them are located on the surface, where proteins are in contact with other molecules. We present here a set of Poaceae genes that are likely to be involved in plant defense mechanisms and have evidence of positive selection. These genes are excellent candidates for future functional validation.

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

The molecular diversity and evolution of Rice tungro bacilliform virus from Indian perspective

Rice tungro disease is caused by a combination of two viruses: Rice tungro spherical virus and Rice tungro bacilliform virus (RTBV). This study was performed with the objective to decipher the molecular variability and evolution of RTBV isolates present in the tungro-affected states of Indian subcontinent. Phylogenetic analysis based on ORF-I, ORF-II, and ORF-IV sequences showed distinct divergence of Indian RTBV isolates into two groups; one consisted isolates from Hyderabad (Andhra Pradesh), Cuttack (Orissa), and Puducherry and another from West Bengal, Chinsura West Bengal, and Kanyakumari (Tamil Nadu). The results obtained from phylogenetic analysis were further supported with the single nucleotide polymorphisms (SNPs), insertion and deletion (INDELs) and evolutionary distance analysis. In addition, sequence difference count matrix revealed a maximum of 56 (ORF-I), 13 (ORF-II) and 73 (ORF-IV) nucleotides differences among all the Indian RTBV isolates taken in this study. However, at the protein level these differences were not significant as revealed by K (a)/K (s) ratio calculation. Sequence identity at nucleotide and amino acid level was 92-100 % (ORF-I), 96-100 % (ORF-II), 94-100 % (ORF-IV) and 86-100 % (ORF-I), 98-100 % (ORF-II) and 95-100 % (ORF-IV), respectively, among Indian isolates of RTBV. The divergence of RTBV isolates into two independent clusters of Indian and non-Indian was shown with the help of the data obtained from phylogeny, SNPs, and INDELs, evolutionary distance analysis, and conserved motifs analysis. The important role of ORF-I and ORF-IV in RTBV diversification and adaptation to different rice growing regions is also discussed.

The plant proteome folding project: structure and positive selection in plant protein families

Despite its importance, relatively little is known about the relationship between the structure, function, and evolution of proteins, particularly in land plant species. We have developed a database with predicted protein domains for five plant proteomes (http://pfp.bio.nyu.edu) and used both protein structural fold recognition and de novo Rosetta-based protein structure prediction to predict protein structure for Arabidopsis and rice proteins. Based on sequence similarity, we have identified ∼15,000 orthologous/paralogous protein family clusters among these species and used codon-based models to predict positive selection in protein evolution within 175 of these sequence clusters. Our results show that codons that display positive selection appear to be less frequent in helical and strand regions and are overrepresented in amino acid residues that are associated with a change in protein secondary structure. Like in other organisms, disordered protein regions also appear to have more selected sites. Structural information provides new functional insights into specific plant proteins and allows us to map positively selected amino acid sites onto protein structures and view these sites in a structural and functional context.

Rates of evolution in stress-related genes are associated with habitat preference in two Cardamine lineages

Background

Elucidating the selective and neutral forces underlying molecular evolution is fundamental to understanding the genetic basis of adaptation. Plants have evolved a suite of adaptive responses to cope with variable environmental conditions, but relatively little is known about which genes are involved in such responses. Here we studied molecular evolution on a genome-wide scale in two species of Cardamine with distinct habitat preferences: C. resedifolia, found at high altitudes, and C. impatiens, found at low altitudes. Our analyses focussed on genes that are involved in stress responses to two factors that differentiate the high- and low-altitude habitats, namely temperature and irradiation.

Results

High-throughput sequencing was used to obtain gene sequences from C. resedifolia and C. impatiens. Using the available A. thaliana gene sequences and annotation, we identified nearly 3,000 triplets of putative orthologues, including genes involved in cold response, photosynthesis or in general stress responses. By comparing estimated rates of molecular substitution, codon usage, and gene expression in these species with those of Arabidopsis, we were able to evaluate the role of positive and relaxed selection in driving the evolution of Cardamine genes. Our analyses revealed a statistically significant higher rate of molecular substitution in C. resedifolia than in C. impatiens, compatible with more efficient positive selection in the former. Conversely, the genome-wide level of selective pressure is compatible with more relaxed selection in C. impatiens. Moreover, levels of selective pressure were heterogeneous between functional classes and between species, with cold responsive genes evolving particularly fast in C. resedifolia, but not in C. impatiens.

Conclusions

Overall, our comparative genomic analyses revealed that differences in effective population size might contribute to the differences in the rate of protein evolution and in the levels of selective pressure between the C. impatiens and C. resedifolia lineages. The within-species analyses also revealed evolutionary patterns associated with habitat preference of two Cardamine species. We conclude that the selective pressures associated with the habitats typical of C. resedifolia may have caused the rapid evolution of genes involved in cold response.

Genetic variation of coat protein gene among the isolates of Rice tungro spherical virus from tungro-endemic states of the India

Rice tungro disease, one of the major constraints to rice production in South and Southeast Asia, is caused by a combination of two viruses: Rice tungro spherical virus (RTSV) and Rice tungro bacilliform virus (RTBV). The present study was undertaken to determine the genetic variation of RTSV population present in tungro endemic states of Indian subcontinent. Phylogenetic analysis based on coat protein sequences showed distinct divergence of Indian RTSV isolates into two groups; one consisted isolates from Hyderabad (Andhra Pradesh), Cuttack (Orissa), and Puducherry and another from West Bengal, Coimbatore (Tamil Nadu), and Kanyakumari (Tamil Nadu). The results obtained from phylogenetic study were further supported with the SNPs (single nucleotide polymorphism), INDELs (insertion and deletion) and evolutionary distance analysis. In addition, sequence difference count matrix revealed 2-68 nucleotides differences among all the Indian RTSV isolates taken in this study. However, at the protein level these differences were not significant as revealed by Ka/Ks ratio calculation. Sequence identity at nucleotide and amino acid level was 92-100% and 97-100%, respectively, among Indian isolates of RTSV. Understanding of the population structure of RTSV from tungro endemic regions of India would potentially provide insights into the molecular diversification of this virus.

Discovery of polymorphisms in starch-related genes in rice germplasm by amplification of pooled DNA and deeply parallel sequencing

High-throughput sequencing of pooled DNA was applied to polymorphism discovery in candidate genes involved in starch synthesis. This approach employed semi- to long-range PCR (LR-PCR) followed by next-generation sequencing technology. A total of 17 rice starch synthesis genes encoding seven classes of enzymes, including ADP-glucose pyrophosphorylase (AGPase), granule starch synthase (GBSS), soluble starch synthase (SS), starch branching enzyme (BE), starch debranching enzyme (DBE) and starch phosphorylase (SPHOL) and phosphate translocator (GPT1) from 233 genotypes were PCR amplified using semi- to long-range PCR. The amplification products were equimolarly pooled and sequenced using massively parallel sequencing technology (MPS). By detecting single nucleotide polymorphism (SNP)/Indels in both coding and noncoding areas of the genes, we identified genetic differences and characterized the SNP/Indel variation and distribution patterns among individual starch candidate genes. Approximately, 60.9 million reads were generated, of which 54.8 million (90%) mapped to the reference sequences. The average coverage rate ranged from 12,708 to 38,300 times for SSIIa and SSIIIb, respectively. SNPs and single/multiple-base Indels were analysed in a total assembled length of 116,403 bp. In total, 501 SNPs and 113 Indels were detected across the 17 starch-related loci. The ratio of synonymous to nonsynonymous SNPs (Ka/Ks) test indicated GBSSI and isoamylase 1 (ISA1) as the least diversified (most purified) and conservative genes as the studied populations have been through cycles of selection. This report demonstrates a useful strategy for screening germplasm by MPS to discover variants in a specific target group of genes.

The evolution of protein structures and structural ensembles under functional constraint

Protein sequence, structure, and function are inherently linked through evolution and population genetics. Our knowledge of protein structure comes from solved structures in the Protein Data Bank (PDB), our knowledge of sequence through sequences found in the NCBI sequence databases (http://www.ncbi.nlm.nih.gov/), and our knowledge of function through a limited set of in-vitro biochemical studies. How these intersect through evolution is described in the first part of the review. In the second part, our understanding of a series of questions is addressed. This includes how sequences evolve within structures, how evolutionary processes enable structural transitions, how the folding process can change through evolution and what the fitness impacts of this might be. Moving beyond static structures, the evolution of protein kinetics (including normal modes) is discussed, as is the evolution of conformational ensembles and structurally disordered proteins. This ties back to a question of the role of neostructuralization and how it relates to selection on sequences for functions. The relationship between metastability, the fitness landscape, sequence divergence, and organismal effective population size is explored. Lastly, a brief discussion of modeling the evolution of sequences of ordered and disordered proteins is entertained.

An EST-based analysis identifies new genes and reveals distinctive gene expression features of Coffea arabica and Coffea canephora

BACKGROUND

Coffee is one of the world's most important crops; it is consumed worldwide and plays a significant role in the economy of producing countries. Coffea arabica and C. canephora are responsible for 70 and 30% of commercial production, respectively. C. arabica is an allotetraploid from a recent hybridization of the diploid species, C. canephora and C. eugenioides. C. arabica has lower genetic diversity and results in a higher quality beverage than C. canephora. Research initiatives have been launched to produce genomic and transcriptomic data about Coffea spp. as a strategy to improve breeding efficiency.

RESULTS

Assembling the expressed sequence tags (ESTs) of C. arabica and C. canephora produced by the Brazilian Coffee Genome Project and the Nestlé-Cornell Consortium revealed 32,007 clusters of C. arabica and 16,665 clusters of C. canephora. We detected different GC3 profiles between these species that are related to their genome structure and mating system. BLAST analysis revealed similarities between coffee and grape (Vitis vinifera) genes. Using KA/KS analysis, we identified coffee genes under purifying and positive selection. Protein domain and gene ontology analyses suggested differences between Coffea spp. data, mainly in relation to complex sugar synthases and nucleotide binding proteins. OrthoMCL was used to identify specific and prevalent coffee protein families when compared to five other plant species. Among the interesting families annotated are new cystatins, glycine-rich proteins and RALF-like peptides. Hierarchical clustering was used to independently group C. arabica and C. canephora expression clusters according to expression data extracted from EST libraries, resulting in the identification of differentially expressed genes. Based on these results, we emphasize gene annotation and discuss plant defenses, abiotic stress and cup quality-related functional categories.

CONCLUSION

We present the first comprehensive genome-wide transcript profile study of C. arabica and C. canephora, which can be freely assessed by the scientific community at http://www.lge.ibi.unicamp.br/coffea. Our data reveal the presence of species-specific/prevalent genes in coffee that may help to explain particular characteristics of these two crops. The identification of differentially expressed transcripts offers a starting point for the correlation between gene expression profiles and Coffea spp. developmental traits, providing valuable insights for coffee breeding and biotechnology, especially concerning sugar metabolism and stress tolerance.

Ancient and recent adaptive evolution of primate non-homologous end joining genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.

Because all cells experience DNA damage, they must also have mechanisms for repairing DNA. When the proteins that repair DNA malfunction, mutation and disease often result. Based on their fundamental importance, DNA repair proteins would be expected to be well preserved over evolutionary time in order to ensure optimal DNA repair function. However, a previous genome-wide study of molecular evolution in Saccharomyces yeast identified the non-homologous end joining (NHEJ) DNA repair pathway as one of the two most rapidly evolving pathways in the yeast genome. In order to analyze the evolution of this pathway in humans, we have generated large evolutionary sequence sets of NHEJ genes from our primate relatives. Similar to the scenario in yeast, several genes in this pathway are evolving rapidly in primate genomes and in modern human populations. Thus, complex and seemingly opposite selective forces are shaping the evolution of these important DNA repair genes. The finding that NHEJ genes are rapidly evolving in species groups as diverse as yeasts and primates indicates a systematic perturbation of the NHEJ pathway, one that is potentially important to human health.

Lineage-specific differences in the amino acid substitution process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size. To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.

Comparative genomic analysis of C4 photosynthetic pathway evolution in grasses

BACKGROUND

Sorghum is the first C4 plant and the second grass with a full genome sequence available. This makes it possible to perform a whole-genome-level exploration of C4 pathway evolution by comparing key photosynthetic enzyme genes in sorghum, maize (C4) and rice (C3), and to investigate a long-standing hypothesis that a reservoir of duplicated genes is a prerequisite for the evolution of C4 photosynthesis from a C3 progenitor.

RESULTS

We show that both whole-genome and individual gene duplication have contributed to the evolution of C4 photosynthesis. The C4 gene isoforms show differential duplicability, with some C4 genes being recruited from whole genome duplication duplicates by multiple modes of functional innovation. The sorghum and maize carbonic anhydrase genes display a novel mode of new gene formation, with recursive tandem duplication and gene fusion accompanied by adaptive evolution to produce C4 genes with one to three functional units. Other C4 enzymes in sorghum and maize also show evidence of adaptive evolution, though differing in level and mode. Intriguingly, a phosphoenolpyruvate carboxylase gene in the C3 plant rice has also been evolving rapidly and shows evidence of adaptive evolution, although lacking key mutations that are characteristic of C4 metabolism. We also found evidence that both gene redundancy and alternative splicing may have sheltered the evolution of new function.

CONCLUSIONS

Gene duplication followed by functional innovation is common to evolution of most but not all C4 genes. The apparently long time-lag between the availability of duplicates for recruitment into C4 and the appearance of C4 grasses, together with the heterogeneity of origins of C4 genes, suggests that there may have been a long transition process before the establishment of C4 photosynthesis.