Genome-Wide In Silico Analysis and Expression Profiling of Phosphoenolpyruvate Carboxylase Genes in Loquat, Apple, Peach, Strawberry and Pear

Phosphoenolpyruvate carboxylase (PEPC) genes have multiple potential roles in plant metabolism such as regulation and accumulation of organic acids in fruits, movement of guard cells and stress tolerance, etc. However, the systematic identification and characterization of PEPC genes in Rosaceae species i.e., loquat, apple, peach, strawberry, and pear are yet to be performed. In present study, 27 putative PEPC genes (loquat 4, apple 6, peach 3, strawberry 9, and pear 5) were identified. To further investigate the role of those PEPC genes, comprehensive bioinformatics and expression analysis were performed. In bioinformatic analysis, the physiochemical properties, conserved domains, gene structure, conserved motif, phylogenetic and syntenic analysis of PEPC genes were performed. The result revealed that the PEPcase superfamily domain was conserved in all examined PEPC proteins. Most of the PEPC proteins were predicted to be localized in cytonuclear. Genomic structural and motif analysis showed that the exon and motif number of each PEPC gene ranged dramatically, from 8 to 20, and 7 to 10, respectively. Syntenic analysis indicated that the segmental or whole-genome duplication played a vital role in extension of PEPC gene family in Rosacea species. The Ka and Ks values of duplicated genes depicted that PEPC genes have undergone a strong purifying selection. Furthermore, the expression analysis of PEPC genes in root, mature leaf, stem, full-bloom flower, and ripened fruit of loquat, apple, peach, strawberry, and pear was performed. Some genes were differentially expressed in aforementioned plant tissues, signifying their role in plant metabolism. This study provides the first genome-wide identification, characterization, and expression profiling of PEPC gene family in Rosaceae species, and provides the foundation for further functional analysis.

Genome-Wide Identification, Characterization and Expression Profiling of Aluminum-Activated Malate Transporters in Eriobotrya japonica Lindl.

Aluminum-activated malate transporters (ALMTs) have multiple potential roles in plant metabolism such as regulation of organic acids in fruits, movement of guard cells and inducing tolerance against aluminum stress. However, the systematic characterization of ALMT genes in loquat is yet to be performed. In the current study, 24 putative ALMT genes were identified in the genome of Eriobotrya japonica Lindl. To further investigate the role of those ALMT genes, comprehensive bioinformatics and expression analysis were performed. In bioinformatics analysis, the physiochemical properties, conserved domains, gene structure, conserved motif, phylogenetic and syntenic analysis of EjALMT genes were conducted. The result revealed that the ALMT superfamily domain was conserved in all EjALMT proteins. EjALMT proteins were predicted to be localized in the plasma membrane. Genomic structural and motif analysis showed that the exon and motif number of each EjALMT gene ranged dramatically, from 5 to 7, and 6 to 10, respectively. Syntenic analysis indicated that the segmental or whole-genome duplication played a vital role in extension of the EjALMT gene family. The Ka and Ks values of duplicated genes depicted that EjALMT genes have undergone a strong purifying selection. Furthermore, the expression analysis of EjALMT genes was performed in the root, mature leaf, stem, full-bloom flower and ripened fruit of loquat. Some genes were expressed differentially in examined loquat tissues, signifying their differential role in plant growth and development. This study provides the first genome-wide identification, characterization, and relative expression of the ALMT gene family in loquat and provides the foundation for further functional analysis.

Dissecting the subtropical adaptation traits and cuticle synthesis pathways via the genome of the subtropical blueberry Vaccinium darrowii

Vaccinium darrowii is a subtropical wild blueberry species, which was used to breed economically important southern highbush cultivars. The adaptation traits of V. darrowii to subtropical climate would provide valuable information for breeding blueberry and perhaps other plants, especially against the background of global warming. Here, we assembled the V. darrowii genome into 12 pseudochoromosomes using Oxford Nanopore long reads complemented with Hi-C scaffolding technologies, and predicted 41 815 genes using RNAseq evidence. Syntenic analysis across three Vaccinium species revealed a highly conserved genome structure, with the highest collinearity between V. darrowii and V. corymbosum. This conserved genome structure may explain the high fertilization during crossbreeding between V. darrowii and other blueberry cultivars. Gene expansion and tandem duplication analysis indicated possible roles of defense and flowering associated genes in adaptation of V. darrowii to the subtropics. The possible SOC1 genes in V. darrowii were identified with phylogeny and expression analysis. Blueberries are covered in a thick cuticle layer and contain anthocyanins, which confer their powdery blue color. Using RNA-sequencing, the cuticle biosynthesis pathways of Vaccinium species were delineated here in V. darrowii. This result could serve as a reference for breeding berries with customer-desired colors. The V. darrowii reference genome, together with the unique traits of this species, including diploid genome, short vegetative phase, and high compatibility in hybridization with other blueberries, make V. darrowii a potential research model for blueberry species.

The Mitragyna speciosa (Kratom) Genome: a resource for data-mining potent pharmaceuticals that impact human health

Mitragyna speciosa (kratom) produces numerous compounds with pharmaceutical properties including the production of bioactive monoterpene indole and oxindole alkaloids. Using a linked-read approach, a 1,122,519,462 bp draft assembly of M. speciosa “Rifat” was generated with an N50 scaffold size of 1,020,971 bp and an N50 contig size of 70,448 bp that encodes 55,746 genes. Chromosome counting revealed that “Rifat” is a tetraploid with a base chromosome number of 11, which was further corroborated by orthology and syntenic analysis of the genome. Analysis of genes and clusters involved in specialized metabolism revealed genes putatively involved in alkaloid biosynthesis. Access to the genome of M. speciosa will facilitate an improved understanding of alkaloid biosynthesis and accelerate the production of bioactive alkaloids in heterologous hosts.

Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family Under Abiotic Stresses in Medicago truncatula

The trihelix transcription factor (GT) family is widely involved in regulating plant growth and development, and most importantly, responding to various abiotic stresses. Our study first reported the genome-wide identification and analysis of GT family genes in Medicago truncatula. Overall, 38 trihelix genes were identified in the M. truncatula genome and were classified into five subfamilies (GT-1, GT-2, SH4, GTγ and SIP1). We systematically analyzed the phylogenetic relationship, chromosomal distribution, tandem and segmental duplication events, gene structures and conserved motifs of MtGTs. Syntenic analysis revealed that trihelix family genes in M. truncatula had the most collinearity relationship with those in soybean followed by alfalfa, but very little collinearity with those in the maize and rice. Additionally, tissue-specific expression analysis of trihelix family genes suggested that they played various roles in the growth and development of specific tissues in M. truncatula. Moreover, the expression of some MtGT genes, such as MtGT19, MtGT20, MtGT22, and MtGT33, was dramatically induced by drought, salt, and ABA treatments, illustrating their vital roles in response to abiotic stresses. These findings are helpful for improving the comprehensive understanding of trihelix family; additionally, the study provides candidate genes for achieving the genetic improvement of stress resistance in legumes.

The Genome Sequence of Five Highly Pathogenic Isolates of Fusarium oxysporum f. sp. lini

Fusarium wilt is the most destructive fungal disease in flax, limiting flax cultivation in all the main flax and linseed growing countries. The causative agent is seedborne and soilborne fungus F. oxysporum f. sp. lini. Here, we report, for the first time, genome assemblies of five highly pathogenic isolates of Fusarium oxysporum f. sp. lini, namely monoisolate 39 and strains F329, F324, F282, F287. In addition, syntenic analysis provided a powerful approach to distinguish between core and lineage-specific parts of the genome. These results lay a solid foundation for comparative genomics studies of plant fungal pathogens, evolution of pathogenicity, and virulence factors underlying the dynamics of host-pathogen interactions, thus eventually offering solutions to Fusarium disease control.

Coevolution of the Spexin/Galanin/Kisspeptin Family: Spexin Activates Galanin Receptor Type II and III

The novel neuropeptide spexin (SPX) was discovered using bioinformatics. The function of this peptide is currently under investigation. Here, we identified SPX along with a second SPX gene (SPX2) in vertebrate genomes. Syntenic analysis and relocating SPXs and their neighbor genes on reconstructed vertebrate ancestral chromosomes revealed that SPXs reside in the near vicinity of the kisspeptin (KISS) and galanin (GAL) family genes on the chromosomes. Alignment of mature peptide sequences showed some extent of sequence similarity among the 3 peptide groups. Gene structure analysis indicated that SPX is more closely related to GAL than KISS. These results suggest that the SPX, GAL, and KISS genes arose through local duplications before 2 rounds (2R) of whole-genome duplication. Receptors of KISS and GAL (GAL receptor [GALR]) are phylogenetically closest among rhodopsin-like G protein-coupled receptors, and synteny revealed the presence of 3 distinct receptor families KISS receptor, GALR1, and GALR2/3 before 2R. A ligand-receptor interaction study showed that SPXs activate human, Xenopus, and zebrafish GALR2/3 family receptors but not GALR1, suggesting that SPXs are natural ligands for GALR2/3. Particularly, SPXs exhibited much higher potency toward GALR3 than GAL. Together, these results identify the coevolution of SPX/GAL/KISS ligand genes with their receptor genes. This study demonstrates the advantage of evolutionary genomics to explore the evolutionary relationship of a peptide gene family that arose before 2R by local duplications.

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis

Myosin heavy chain genes ( MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu ( Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH M743-2 and MYH M5 of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.