Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition

Systematic analysis of the Rboh gene family in seven gramineous plants and its roles in response to arbuscular mycorrhizal fungi in maize

BACKGROUND

Plant respiratory burst oxidase homolog (Rboh) gene family produces reactive oxygen species (ROS), and it plays key roles in plant-microbe interaction. Most Rboh gene family-related studies mainly focused on dicotyledonous plants; however, little is known about the roles of Rboh genes in gramineae.

RESULTS

A total of 106 Rboh genes were identified in seven gramineae species, including Zea mays, Sorghum bicolor, Brachypodium distachyon, Oryza sativa, Setaria italica, Hordeum vulgare, and Triticum aestivum. The Rboh protein sequences showed high similarities, suggesting that they may have conserved functions across different species. Duplication mode analysis detected whole-genome/segmental duplication (WGD)/(SD) and dispersed in the seven species. Interestingly, two local duplication (LD, including tandem and proximal duplication) modes were found in Z. mays, S. italica and H. vulgare, while four LD were detected in T. aestivum, indicating that these genes may have similar functions. Collinearity analysis indicated that Rboh genes are at a stable evolution state in all the seven species. Besides, Rboh genes from Z. mays were closely related to those from S. bicolor, consistent with the current understanding of plant evolutionary history. Phylogenetic analysis showed that the genes in the subgroups I and II may participate in plant-AM fungus symbiosis. Cis-element analysis showed that different numbers of elements are related to fungal induction in the promoter region. Expression profiles of Rboh genes in Z. mays suggested that Rboh genes had distinct spatial expression patterns. By inoculation with AM fungi, our transcriptome analysis showed that the expression of Rboh genes varies upon AM fungal inoculation. In particularly, ZmRbohF was significantly upregulated after inoculation with AM fungi. pZmRbohF::GUS expression analyses indicated that ZmRbohF was induced by arbuscular mycorrhizal fungi in maize. By comparing WT and ZmRbohF mutant, we found ZmRbohF had limited impact on the establishment of maize-AM fungi symbiosis, but play critical roles in regulating the proper development of arbuscules.

CONCLUSIONS

This study provides a comprehensive analysis of the evolution relationship of Rboh genes in seven gramineae species. Results showed that several Rboh genes regulate maize-AM fungal symbiosis process. This study provides valuable information for further studies of Rboh genes in gramineae.

Comprehensive characterization and expression analysis of enzymatic antioxidant gene families in passion fruit (Passiflora edulis)

Passion fruit, a valuable tropical fruit, faces climate-related growth challenges. Antioxidant enzymes are vital for both stress protection and growth regulation in plants. We first provided systemic analysis of enzymatic antioxidant gene families in passion fruit, identifying 90 members including 11 PeSODs, 45 PeAPXs, 8 PeCATs, 7 PeGPXs, 6 PeMDHARs, 8 PeDHARs, and 5 PeGRs. Gene members in each gene family with same subcellular localization showed closer phylogenetic relationship. Many antioxidant genes exhibited tissue- or developmental stage-specific expression patterns during floral and fruit development, with some widely expressed. Their co-expressed genes were linked to photosynthesis and energy metabolism, suggesting roles in protecting highly proliferating tissues from oxidative damage. Potential genes for enhancing temperature stress resistance were identified. The involvement of diverse regulatory factors including miRNAs, transcription factors, and CREs might contribute to the complex roles of antioxidant genes. This study informs future research on antioxidant genes and passion fruit breeding.

VDR promotes testosterone synthesis in mouse Leydig cells via regulation of cholesterol side chain cleavage cytochrome P450 (Cyp11a1) expression

BACKGROUND

The vitamin D receptor (VDR) mediates the pleiotropic biological actions that include osteoporosis, immune responses and androgen synthesis.VDR is widely expressed in testis cells such as Leydig cells, Sertoli cells, and sperm. The levels of steroids are critical for sexual development. In the early stage of steroidogenesis, cholesterol is converted to pregnenolone (precursor of most steroid hormones) by cholesterol side-chain lyase (CYP11A1), which eventually synthesizes the male hormone testosterone.

OBJECTIVE

This study aims to reveal how VDR regulates CYP11A1 expression and affects testosterone synthesis in murine Leydig cells.

METHODS

The levels of VDR, CYP11A1 were determined by quantitative real-time polymerase chain reaction (RT-qPCR) or western blot. Targeted relationship between VDR and Cyp11a1 was evaluated by dual-luciferase reporter assay. The levels of testosterone concentrations in cell culture media serum by enzyme-linked immunosorbent assay (ELISA).

RESULTS

Phylogenetic and motif analysis showed that the Cyp11a1 family had sequence loss, which may have special biological functions during evolution. The results of promoter prediction showed that vitamin D response element (VDRE) existed in the upstream promoter region of murine Cyp11a1. Dual-luciferase assay confirmed that VDR could bind candidate VDREs in upstream region of Cyp11a1, and enhance gene expression. Tissue distribution and localizatio analysis showed that Cyp11a1 was mainly expressed in testis, and dominantly existed in murine Leydig cells. Furthermore, over-expression VDR and CYP11A1 significantly increased testosterone synthesis in mice Leydig cells.

CONCLUSIONS

Active vitamin D3 (VD3) and Vdr interference treatment showed that VD3/VDR had a positive regulatory effect on Cyp11a1 expression and testosterone secretion. VDR promotes testosterone synthesis in male mice by up-regulating Cyp11a1 expression, which played an important role for male reproduction.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Malus domestica

The conserved BURP-containing proteins are specific to plants and play a crucial role in plant growth, development, and response to abiotic stresses. However, less is known about the systematic characterization of BURP-containing proteins in apple. This study aimed to identify and analyze all BURP-containing genes in the apple genome, as well as to examine their expression patterns through various bioinformatics methods. Eighteen members of BURP-containing genes were identified in apple, six members lacked signal peptides, and the secondary structure was mainly a Random coil of BURP-containing genes. Gene structure and Motif analysis showed that proteins have similar structures and are conserved at the C-terminal. Cis-acting element analysis revealed that the proteins contain phytohormone and stress response elements, and chromosomal localization revealed that the family is unevenly distributed across eight chromosomes, with duplication of fragments leading to the expansion of family proteins. Tissue expression showed that MdPG3 and MdPG4 were expressed in different tissues and different varieties, MdRD2 and MdRD7 were highly expressed in 'M74' fruits and MdRD7 in 'M49' leaves, while MdUSP1 was highly expressed in 'GD' roots. The quantitative real-time PCR analysis showed that the expressions of six and seven genes were significantly up-regulated under NaCl and PEG treatments, respectively, whereas MdRD7 was significantly up-regulated under NaCl and PEG treatment over time. This study offers a comprehensive identification and expression analysis of BURP-containing proteins in apple. The findings provide a theoretical foundation for further exploration of the functions of this protein family.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01393-7.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Superoxide dismutase gene family in cassava revealed their involvement in environmental stress via genome-wide analysis

Superoxide dismutase (SOD) is a crucial metal-containing enzyme that plays a vital role in catalyzing the dismutation of superoxide anions, converting them into molecular oxygen and hydrogen peroxide, essential for enhancing plant stress tolerance. We identified 8 SOD genes (4 CSODs, 2 FSODs, and 2 MSODs) in cassava. Bioinformatics analyses provided insights into chromosomal location, phylogenetic relationships, gene structure, conserved motifs, and gene ontology annotations. MeSOD genes were classified into two groups through phylogenetic analysis, revealing evolutionary connections. Promoters of these genes harbored stress-related cis-elements. Duplication analysis indicated the functional significance of MeCSOD2/MeCSOD4 and MeMSOD1/MeMSOD2. Through qRT-PCR, MeCSOD2 responded to salt stress, MeMSOD2 to drought, and cassava bacterial blight. Silencing MeMSOD2 increased XpmCHN11 virulence, indicating MeMSOD2 is essential for cassava's defense against XpmCHN11 infection. These findings enhance our understanding of the SOD gene family's role in cassava and contribute to strategies for stress tolerance improvement.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

Genome-wide identification of GH9 gene family and the assessment of its role during fruit abscission zone formation in Vaccinium ashei

KEY MESSAGE

The genomic location and stage-specific expression pattern of GH9 genes reveal their critical roles during fruit abscission zone formation in Vaccinium ashei. Glycosyl hydrolase family 9 (GH9) cellulases play a crucial role in both cellulose synthesis and hydrolysis during plant growth and development. Despite this importance, there is currently no study on the involvement of GH9-encoding genes, specifically VaGH9s, in abscission zone formation of rabbiteye blueberries (Vaccinium ashei). In this study, we identified a total of 61 VaGH9s in the genome, which can be classified into 3 subclasses based on conserved motifs and domains, gene structures, and phylogenetic analyses. Our synteny analysis revealed that VaGH9s are more closely related to the GH9s of Populus L. than to those of Arabidopsis, Vitis vinifera, and Citrus sinensis. In silico structural analysis predicted that most of VaGH9s are hydrophilic, and localized in cell membrane and/or cell wall, and the variable sets of cis-acting regulatory elements and functional diversity with four categories of stress response, hormone regulation, growth and development, and transcription factor-related elements are present in the promoter sequence of VaGH9s genes. Transcriptomic analysis showed that there were 22 differentially expressed VaGH9s in fruit abscission zone tissue at the veraison stage, and the expression of VaGH9B2 and VaGH9C10 was continuously increased during fruit maturation, which were in parallel with the increasing levels of cellulase activity and oxidative stress indicators, suggesting that they are involved in the separation stage of fruit abscission in Vaccinium ashei. Our work identified 22 VaGH9s potentially involved in different stages of fruit abscission and would aid further investigation into the molecular regulation of abscission in rabbiteye blueberries fruit.

Decoupling subgenomes within hybrid lavandin provide new insights into speciation and monoterpenoid diversification of Lavandula

Polyploidization and transposon elements contribute to shape plant genome diversity and secondary metabolic variation in some edible crops. However, the specific contribution of these variations to the chemo-diversity of Lamiaceae, particularly in economic shrubs, is still poorly documented. The rich essential oils (EOs) of Lavandula plants are distinguished by monoterpenoids among the main EO-producing species, L. angustifolia (LA), L. × intermedia (LX) and L. latifolia (LL). Herein, the first allele-aware chromosome-level genome was assembled using a lavandin cultivar 'Super' and its hybrid origin was verified by two complete subgenomes (LX-LA and LX-LL). Genome-wide phylogenetics confirmed that LL, like LA, underwent two lineage-specific WGDs after the γ triplication event, and their speciation occurred after the last WGD. Chloroplast phylogenetic analysis indicated LA was the maternal source of 'Super', which produced premium EO (higher linalyl/lavandulyl acetate and lower 1,8-cineole and camphor) close to LA. Gene expression, especially the monoterpenoid biosynthetic genes, showed bias to LX-LA alleles. Asymmetric transposon insertions in two decoupling 'Super' subgenomes were responsible for speciation and monoterpenoid divergence of the progenitors. Both hybrid and parental evolutionary analysis revealed that LTR (long terminal repeat) retrotransposon associated with AAT gene loss cause no linalyl/lavandulyl acetate production in LL, and multi-BDH copies retained by tandem duplication and DNA transposon resulted in higher camphor accumulation of LL. Advances in allelic variations of monoterpenoids have the potential to revolutionize future lavandin breeding and EO production.

Genome-wide identification and expression analysis of the trihelix transcription factor family in sesame (Sesamum indicum L.) under abiotic stress

BACKGROUND

The plant trihelix gene family is among the earliest discovered transcription factor families, and it is vital in modulating light, plant growth, and stress responses.

METHODS

The identification and characterization of trihelix family members in the sesame genome were analyzed by bioinformatics methods, and the expression patterns of sesame trihelix genes were assessed by quantitative real-time PCR.

RESULTS

There were 34 trihelix genes discovered in the genome of sesame, which were irregularly distributed among 10 linkage groups. Also, the genome contained 5 duplicate gene pairs. The 34 trihelix genes were divided into six sub-families through a phylogenetic study. A tissue-specific expression revealed that SiTH genes exhibited spatial expression patterns distinct from other trihelix genes in the same subfamily. The cis-element showed that the SiTHs gene promoter contained various elements associated with responses to hormones and multiple abiotic stresses. Additionally, the expression patterns of 8 SiTH genes in leaves under abiotic stresses demonstrated that all selected genes were significantly upregulated or downregulated at least once in the stress period. Furthermore, the SiTH4 gene was significantly induced in response to drought and salt stress, showing that SiTH genes may be engaged in the stress response mechanisms of sesame.

CONCLUSION

These findings establish a foundation for further investigation of the trihelix gene-mediated response to abiotic stress in sesame.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt

BACKGROUND

The receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense.

RESULTS

Seventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating "20B12" (cotton cultivar) with "V991" (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt.

CONCLUSIONS

These observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt.

Hsp70 Gene Family in Sebastiscus marmoratus: The Genome-Wide Identification and Transcriptome Analysis under Thermal Stress

Heat shock protein 70 kDa (Hsp70) is a highly conserved heat stress protein that is important in biotic processes and responses to abiotic stress. Hsp70 genes may be important in Sebastiscus marmoratus, for it is a kind of nearshore reef fish, and habitat temperature change is more drastic during development. However, genome-wide identification and expression analysis in the Hsp70 gene family of S. marmoratus are still lacking. Here, a total of 15 Hsp70 genes in the genome of S. marmoratus are identified, and their expression patterns were investigated using transcriptomic data from thermal stress experiments. The expansion and gene duplication events of Hsp70 genes from the Hspa4, Hspa8, and Hspa12a subfamilies in S. marmoratus are revealed by phylogenetic analysis. qRT-PCR expression patterns demonstrated that seven Hsp70 genes were significantly up-regulated and none were significantly down-regulated after heat treatment. Only the hsp70 gene was significantly up-regulated after cold treatment. The selection test further showed a purifying selection on the duplicated gene pairs, suggesting that these genes underwent subfunctionalization. Our results add novel insight to aquaculture and biological research on S. marmoratus, providing important information on how Hsp70 genes are regulated in Scorpaeniformes under thermal stress.

An Insight of Betula platyphylla SWEET Gene Family through Genome-Wide Identification, Expression Profiling and Function Analysis of BpSWEET1c under Cold Stress

SWEET proteins play important roles in plant growth and development, sugar loading in phloem and resistance to abiotic stress through sugar transport. In this study, 13 BpSWEET genes were identified from birch genome. Collinearity analysis showed that there were one tandem repeating gene pair (BpSWEET1b/BpSWEET1c) and two duplicative gene pairs (BpSWEET17a/BpSWEET17b) in the BpSWEET gene family. The BpSWEET gene promoter regions contained several cis-acting elements related to stress resistance, for example: hormone-responsive and low-temperature-responsive cis-elements. Analysis of transcriptome data showed that BpSWEET genes were highly expressed in several sink organs, and the most BpSWEET genes were rapidly up-regulated under cold stress. BpSWEET1c, which was highly expressed in cold stress, was selected for further analysis. It was found that BpSWEET1c was located on the cell membrane. After 6 h of 4 °C stress, sucrose content in the leaves and roots of transient overexpressed BpSWEET1c was significantly higher than that of the control. MDA content in roots was significantly lower than that of the control. These results indicate that BpSWEET1c may play a positive role in the response to cold stress by promoting the metabolism and transport of sucrose. In conclusion, 13 BpSWEET genes were identified from the whole genome level. Most of the SWEET genes of birch were expressed in the sink organs and could respond to cold stress. Transient overexpression of BpSWEET1c changed the soluble sugar content and improved the cold tolerance of birch.

Dosage-sensitivity shapes how genes transcriptionally respond to allopolyploidy and homoeologous exchange in resynthesized Brassica napus

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

RNA polymerase common subunit ZmRPABC5b is transcriptionally activated by Opaque2 and essential for endosperm development in maize

Maize (Zea mays) kernel size is an important factor determining grain yield; although numerous genes regulate kernel development, the roles of RNA polymerases in this process are largely unclear. Here, we characterized the defective kernel 701 (dek701) mutant that displays delayed endosperm development but normal vegetative growth and flowering transition, compared to its wild type. We cloned Dek701, which encoded ZmRPABC5b, a common subunit to RNA polymerases I, II and III. Loss-of-function mutation of Dek701 impaired the function of all three RNA polymerases and altered the transcription of genes related to RNA biosynthesis, phytohormone response and starch accumulation. Consistent with this observation, loss-of-function mutation of Dek701 affected cell proliferation and phytohormone homeostasis in maize endosperm. Dek701 was transcriptionally regulated in the endosperm by the transcription factor Opaque2 through binding to the GCN4 motif within the Dek701 promoter, which was subjected to strong artificial selection during maize domestication. Further investigation revealed that DEK701 interacts with the other common RNA polymerase subunit ZmRPABC2. The results of this study provide substantial insight into the Opaque2-ZmRPABC5b transcriptional regulatory network as a central hub for regulating endosperm development in maize.

A Comprehensive Genome-Wide Investigation of the Cytochrome 71 (OsCYP71) Gene Family: Revealing the Impact of Promoter and Gene Variants (Ser33Leu) of OsCYP71P6 on Yield-Related Traits in Indica Rice (Oryza sativa L.)

The cytochrome P450 (CYP450) gene family plays a critical role in plant growth and developmental processes, nutrition, and detoxification of xenobiotics in plants. In the present research, a comprehensive set of 105 OsCYP71 family genes was pinpointed within the genome of indica rice. These genes were categorized into twelve distinct subfamilies, where members within the same subgroup exhibited comparable gene structures and conserved motifs. In addition, 105 OsCYP71 genes were distributed across 11 chromosomes, and 36 pairs of OsCYP71 involved in gene duplication events. Within the promoter region of OsCYP71, there exists an extensive array of cis-elements that are associated with light responsiveness, hormonal regulation, and stress-related signaling. Further, transcriptome profiling revealed that a majority of the genes exhibited responsiveness to hormones and were activated across diverse tissues and developmental stages in rice. The OsCYP71P6 gene is involved in insect resistance, senescence, and yield-related traits in rice. Hence, understanding the association between OsCYP71P6 genetic variants and yield-related traits in rice varieties could provide novel insights for rice improvement. Through the utilization of linear regression models, a total of eight promoters were identified, and a specific gene variant (Ser33Leu) within OsCYP71P6 was found to be linked to spikelet fertility. Additionally, different alleles of the OsCYP71P6 gene identified through in/dels polymorphism in 131 rice varieties were validated for their allelic effects on yield-related traits. Furthermore, the single-plant yield, spikelet number, panicle length, panicle weight, and unfilled grain per panicle for the OsCYP71P6-1 promoter insertion variant were found to contribute 20.19%, 13.65%, 5.637%, 8.79%, and 36.86% more than the deletion variant, respectively. These findings establish a robust groundwork for delving deeper into the functions of OsCYP71-family genes across a range of biological processes. Moreover, these findings provide evidence that allelic variation in the promoter and amino acid substitution of Ser33Leu in the OsCYP71P6 gene could potentially impact traits related to rice yield. Therefore, the identified promoter variants in the OsCYP71P6 gene could be harnessed to amplify rice yields.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Exploring the Heat Shock Transcription Factor (HSF) Gene Family in Ginger: A Genome-Wide Investigation on Evolution, Expression Profiling, and Response to Developmental and Abiotic Stresses

Ginger is a valuable crop known for its nutritional, seasoning, and health benefits. However, abiotic stresses, such as high temperature and drought, can adversely affect its growth and development. Heat shock transcription factors (HSFs) have been recognized as crucial elements for enhancing heat and drought resistance in plants. Nevertheless, no previous study has investigated the HSF gene family in ginger. In this research, a total of 25 ZoHSF members were identified in the ginger genome, which were unevenly distributed across ten chromosomes. The ZoHSF members were divided into three groups (HSFA, HSFB, and HSFC) based on their gene structure, protein motifs, and phylogenetic relationships with Arabidopsis. Interestingly, we found more collinear gene pairs between ZoHSF and HSF genes from monocots, such as rice, wheat, and banana, than dicots like Arabidopsis thaliana. Additionally, we identified 12 ZoHSF genes that likely arose from duplication events. Promoter analysis revealed that the hormone response elements (MEJA-responsiveness and abscisic acid responsiveness) were dominant among the various cis-elements related to the abiotic stress response in ZoHSF promoters. Expression pattern analysis confirmed differential expression of ZoHSF members across different tissues, with most showing responsiveness to heat and drought stress. This study lays the foundation for further investigations into the functional role of ZoHSFs in regulating abiotic stress responses in ginger.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus

BACKGROUND

Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual.

RESULTS

Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture.

CONCLUSIONS

Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Conserved H3K27me3-associated chromatin looping mediates physical interactions of gene clusters in plants

Higher-order chromatin organization is essential for transcriptional regulation, genome stability maintenance, and other genome functions. Increasing evidence has revealed significant differences in 3D chromatin organization between plants and animals. However, the extent, pattern, and rules of chromatin organization in plants are still unclear. In this study, we systematically identified and characterized long-range chromatin loops in the Arabidopsis 3D genome. We identified hundreds of long-range cis chromatin loops and found their anchor regions are closely associated with H3K27me3 epigenetic modifications. Furthermore, we demonstrated that these chromatin loops are dependent on Polycomb group (PcG) proteins, suggesting that the Polycomb repressive complex 2 (PRC2) complex is essential for establishing and maintaining these novel loops. Although most of these PcG-medicated chromatin loops are stable, many of these loops are tissue-specific or dynamically regulated by different treatments. Interestingly, tandemly arrayed gene clusters and metabolic gene clusters are enriched in anchor regions. Long-range H3K27me3-marked chromatin interactions are associated with the coregulation of specific gene clusters. Finally, we also identified H3K27me3-associated chromatin loops associated with gene clusters in Oryza sativa and Glycine max, indicating that these long-range chromatin loops are conserved in plants. Our results provide novel insights into genome evolution and transcriptional coregulation in plants.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Homoeologous evolution of the allotetraploid genome of Poa annua L

BACKGROUND

Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty.

RESULTS

We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis.

CONCLUSIONS

The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Evolution of the 14-3-3 gene family in monocotyledons and dicotyledons and validation of MdGRF13 function in transgenic Arabidopsis thaliana

KEY MESSAGE

The 14-3-3 family is more highly conserved among monocotyledons, and overexpression of MdGRF13 improved drought and salt tolerance in transgenic Arabidopsis thaliana. The 14-3-3 are highly conserved regulatory proteins found in eukaryotes and play an essential role in plant growth, development and stress response. However, the 14-3-3 gene family evolution in monocotyledons and dicotyledons and the biological functions of the MdGRF13 under abiotic stress remain unknown. In our study, 195 members of the 14-3-3 family were identified from 12 species and divided into ε group and the Non-ε group. Synteny analysis within the 14-3-3 family indicated that segmental duplication events contributed to the expansion of the family. Selective pressure analysis indicated that purifying selection was a vital force in the 14-3-3 genes evolution, and monocotyledons had a lower million years ago (Mya) mean values than dicotyledons. Meanwhile, the codon adaptation index (CAI) and frequency of optical codons (FOP) are higher and the effective number of codons (Nc) is lower in monocotyledons 14-3-3 genes compared to dicotyledons. Moreover, the yeast two-hybrid (Y2H) demonstrated that MdGRF13 interacts with MdRD22, MdLHP1a and MdMORF1. Significantly, the malondialdehyde (MDA) content and relative conductivity were decreased, while the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities were increased in transgenic Arabidopsis compared to the wild type (WT) under drought and salt stress. These results suggest that overexpression of MdGRF13 significantly improved the tolerance to drought and salt stress in transgenic Arabidopsis. Thus, our results provide a theoretical basis for exploring the evolution and function of the 14-3-3 gene family in monocotyledons and dicotyledons.

The heavy metal-associated isoprenylated plant protein (HIPP) gene family plays a crucial role in cadmium resistance and accumulation in the tea plant (Camellia sinensis L.)

Heavy metal-associated isoprenylated plant proteins (HIPPs) are only distributed in vascular plants, and are essential for the detoxification and vascular transport of heavy metals in plants. However, the HIPP gene family has not been thoroughly explored in the tea plant (Camellia sinensis). In this study, we systematically identified 56C. sinensis CsHIPP genes from five groups and characterized their phylogeny, structures, and the features of the encoded proteins. The expression patterns of CsHIPP genes in various tissues of C. sinensis were investigated based on a previous RNA-seq data analysis. The expression patterns of CsHIPP genes were explored in cadmium (Cd)-treated C. sinensis roots using our RNA-seq data. Three CsHIPP genes (CsHIPP22, CsHIPP24, and CsHIPP36) with high expression levels in Cd-treated C. sinensis roots were selected as candidate genes associated with Cd tolerance. Overexpression of CsHIPP22, CsHIPP24, and CsHIPP36 in a yeast mutant (ycf1) rescued Cd-sensitive ycf1 yeast and increased the yeast resistance to Cd stress, implying that these three CsHIPPs might be involved in Cd tolerance. These findings will enable the roles of HIPPs in Cd absorption and detoxification to be better understood as well as improving our understanding of the Cd-resistance and Cd-accumulation mechanisms in tea plant.

Genome-Wide Identification of PEBP Gene Family in Solanum lycopersicum

The PEBP gene family is crucial for the growth and development of plants, the transition between vegetative and reproductive growth, the response to light, the production of florigen, and the reaction to several abiotic stressors. The PEBP gene family has been found in numerous species, but the SLPEBP gene family has not yet received a thorough bioinformatics investigation, and the members of this gene family are currently unknown. In this study, bioinformatics was used to identify 12 members of the SLPEBP gene family in tomato and localize them on the chromosomes. The physicochemical characteristics of the proteins encoded by members of the SLPEBP gene family were also examined, along with their intraspecific collinearity, gene structure, conserved motifs, and cis-acting elements. In parallel, a phylogenetic tree was built and the collinear relationships of the PEBP gene family among tomato, potato, pepper, and Arabidopsis were examined. The expression of 12 genes in different tissues and organs of tomato was analyzed using transcriptomic data. It was also hypothesized that SLPEBP3, SLPEBP5, SLPEBP6, SLPEBP8, SLPEBP9, and SLPEBP10 might be related to tomato flowering and that SLPEBP2, SLPEBP3, SLPEBP7, and SLPEBP11 might be related to ovary development based on the tissue-specific expression analysis of SLPEBP gene family members at five different stages during flower bud formation to fruit set. This article's goal is to offer suggestions and research directions for further study of tomato PEBP gene family members.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

Different evolutionary patterns of TIR1/AFBs and AUX/IAAs and their implications for the morphogenesis of land plants

BACKGROUND

The plant hormone auxin is widely involved in plant growth, development, and morphogenesis, and the TIR1/AFB and AUX/IAA proteins are closely linked to rapid auxin response and signal transmission. However, their evolutionary history, historical patterns of expansion and contraction, and changes in interaction relationships are still unknown.

RESULTS

Here, we analyzed the gene duplications, interactions, and expression patterns of TIR1/AFBs and AUX/IAAs to understand their underlying mechanisms of evolution. The ratios of TIR1/AFBs to AUX/IAAs range from 4:2 in Physcomitrium patens to 6:29 in Arabidopsis thaliana and 3:16 in Fragaria vesca. Whole-genome duplication (WGD) and tandem duplication have contributed to the expansion of the AUX/IAA gene family, but numerous TIR1/AFB gene duplicates were lost after WGD. We further analyzed the expression profiles of TIR1/AFBs and AUX/IAAs in different tissue parts of Physcomitrium patens, Selaginella moellendorffii, Arabidopsis thaliana and Fragaria vesca, and found that TIR1/AFBs and AUX/IAAs were highly expressed in all tissues in P. patens, S. moellendorffii. In A. thaliana and F. vesca, TIR1/AFBs maintained the same expression pattern as the ancient plants with high expression in all tissue parts, while AUX/IAAs appeared tissue-specific expression. In F. vesca, 11 AUX/IAAs interacted with TIR1/AFBs with different interaction strengths, and the functional specificity of AUX/IAAs was related to their ability to bind TIR1/AFBs, thus promoting the development of specific higher plant organs. Verification of the interactions among TIR1/AFBs and AUX/IAAs in Marchantia polymorpha and F. vesca also showed that the regulation of AUX/IAA members by TIR1/AFBs became more refined over the course of plant evolution.

CONCLUSIONS

Our results indicate that specific interactions and specific gene expression patterns both contributed to the functional diversification of TIR1/AFBs and AUX/IAAs.

Phylogenomic Analysis of micro-RNA Involved in Juvenile to Flowering-Stage Transition in Photophilic Rice and Its Sister Species

Vegetative to reproductive phase transition in phototropic plants is an important developmental process and is sequentially mediated by the expression of micro-RNA MIR172. To obtain insight into the evolution, adaptation, and function of MIR172 in photophilic rice and its wild relatives, we analyzed the genescape of a 100 kb segment harboring MIR172 homologs from 11 genomes. The expression analysis of MIR172 revealed its incremental accumulation from the 2-leaf to 10-leaf stage, with maximum expression coinciding with the flag-leaf stage in rice. Nonetheless, the microsynteny analysis of MIR172s revealed collinearity within the genus Oryza, but a loss of synteny was observed in (i) MIR172A in O. barthii (AA) and O. glaberima (AA); (ii) MIR172B in O. brachyantha (FF); and (iii) MIR172C in O. punctata (BB). Phylogenetic analysis of precursor sequences/region of MIR172 revealed a distinct tri-modal clade of evolution. The genomic information generated in this investigation through comparative analysis of MIRNA, suggests mature MIR172s to have evolved in a disruptive and conservative mode amongst all Oryza species with a common origin of descent. Further, the phylogenomic delineation provided an insight into the adaptation and molecular evolution of MIR172 to changing environmental conditions (biotic and abiotic) of phototropic rice through natural selection and the opportunity to harness untapped genomic regions from rice wild relatives (RWR).

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Chromosome-level reference genome of Tetrastigma hemsleyanum (Vitaceae) provides insights into genomic evolution and the biosynthesis of phenylpropanoids and flavonoids

Here, we present a high-quality chromosome-scale genome assembly (2.19 Gb) and annotation of Tetrastigma hemsleyanum, a perennial herbaceous liana native to subtropical China with diverse medicinal applications. Approximately 73% of the genome was comprised of transposable elements (TEs), of which long terminal repeat retrotransposons (LTR-RTs) were a predominant group (69% of the genome). The genome size increase of T. hemsleyanum (relative to Vitis species) was mostly due to the proliferation of LTR-RTs. Of the different modes of gene duplication identified, transposed duplication (TRD) and dispersed duplication (DSD) were the predominant ones. Genes, particularly those involved in the phenylpropanoid-flavonoid (PF) pathway and those associated with therapeutic properties and environmental stress resistance, were significantly amplified through recent tandem duplications. We dated the divergence of two intraspecific lineages in Southwest (SW) versus Central-South-East (CSE) China to the late Miocene (approximately 5.2 million years ago). Of those, the former showed more upregulated genes and metabolites. Based on resequencing data of 38 individuals representing both lineages, we identified various candidate genes related to 'response to stimulus' and 'biosynthetic process', including ThFLS11, which is putatively involved in flavonoid accumulation. Overall, this study provides abundant genomic resources for future evolutionary, ecological, and functional genomics studies in T. hemsleyanum and related species.

POInTbrowse: orthology prediction and synteny exploration for paleopolyploid genomes

We describe POInT_browse, a web portal that gives access to the orthology inferences made for polyploid genomes with POInT, the Polyploidy Orthology Inference Tool. Ancient, or paleo-, polyploidy events are widely distributed across the eukaryotic phylogeny, and the combination of duplicated and lost duplicated genes that these polyploidies produce can confound the identification of orthologous genes between genomes. POInT uses conserved synteny and phylogenetic models to infer orthologous genes between genomes with a shared polyploidy. It also gives confidence estimates for those orthology inferences. POInT_browse gives both graphical and query-based access to these inferences from 12 different polyploidy events, allowing users to visualize genomic regions produced by polyploidies and perform batch queries for each polyploidy event, downloading genes trees and coding sequences for orthologous genes meeting user-specified criteria. POInT_browse and the associated data are online at https://wgd.statgen.ncsu.edu .

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Evolution and expression analysis of the caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) gene family in jute (Corchorus L.)

BACKGROUND

Jute is considered one of the most important crops for fiber production and multipurpose usages. Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is a crucial enzyme involved in lignin biosynthesis in plants. The potential functions of CCoAOMT in lignin biosynthesis of jute have been reported in several studies. However, little is known about the evolution of the CCoAOMT gene family, and either their expression level at different developing stages in different jute cultivars, as well as under abiotic stresses including salt and drought stress.

RESULTS

In the present study, 66 CCoAOMT genes from 12 species including 12 and eight CCoAOMTs in Corchorus olitorius and C. capsularis were identified. Phylogenetic analysis revealed that CCoAOMTs could be divided into six groups, and gene expansion was observed in C. olitorius. Furthermore, gene expression analysis of developing jute fibers was conducted at different developmental stages (15, 30, 45, 60, and 90 days after sowing [DAS]) in six varieties (Jute-179 [J179], Lubinyuanguo [LB], and Qiongyueqing [QY] for C. capsularis; Funong No.5 [F5], Kuanyechangguo [KY], and Cvlv [CL] for C. olitorius). The results showed that CCoAOMT1 and CCoAOMT2 were the dominant genes in the CCoAOMT family. Of these two dominant CCoAOMTs, CCoAOMT2 showed a constitutive expression level during the entire growth stages, while CCoAOMT1 exhibited differential expression patterns. These two genes showed higher expression levels in C. olitorius than in C. capsularis. The correlation between lignin content and CCoAOMT gene expression levels indicated that this gene family influences the lignin content of jute. Using real-time quantitative reverse transcription PCR (qRT-PCR), a substantial up-regulation of CCoAOMTs was detected in stem tissues of jute 24 h after drought treatment, with an up to 17-fold increase in expression compared to that of untreated plants.

CONCLUSIONS

This study provides a basis for comprehensive genomic studies of the entire CCoAOMT gene family in C. capsularis and C. olitorius. Comparative genomics analysis among the CCoAOMT gene families of 12 species revealed the close evolutionary relationship among Corchorus, Theobroma cacao and Gossypium raimondii. This study also shows that CCoAOMTs are not only involved in lignin biosynthesis, but also are associated with the abiotic stress response in jute, and suggests the potential use of these lignin-related genes to genetically improve the fiber quality of jute.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Identification and analysis of lignin biosynthesis genes related to fruit ripening and stress response in banana (Musa acuminata L. AAA group, cv. Cavendish)

BACKGROUND

Lignin is a key component of the secondary cell wall of plants, providing mechanical support and facilitating water transport as well as having important impact effects in response to a variety of biological and abiotic stresses.

RESULTS

In this study, we identified 104 genes from ten enzyme gene families related to lignin biosynthesis in Musa acuminata genome and found the number of MaCOMT gene family was the largest, while MaC3Hs had only two members. MaPALs retained the original members, and the number of Ma4CLs in lignin biosynthesis was significantly less than that of flavonoids. Segmental duplication existed in most gene families, except for MaC3Hs, and tandem duplication was the main way to expand the number of MaCOMTs. Moreover, the expression profiles of lignin biosynthesis genes during fruit development, postharvest ripening stages and under various abiotic and biological stresses were investigated using available RNA-sequencing data to obtain fruit ripening and stress response candidate genes. Finally, a co-expression network of lignin biosynthesis genes was constructed by weighted gene co-expression network analysis to elucidate the lignin biosynthesis genes that might participate in lignin biosynthesis in banana during development and in response to stresses.

CONCLUSION

This study systematically identified the lignin biosynthesis genes in the Musa acuminata genome, providing important candidate genes for further functional analysis. The identification of the major genes involved in lignin biosynthesis in banana provides the basis for the development of strategies to improve new banana varieties tolerant to biological and abiotic stresses with high yield and high quality.

Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly

Background

The genus Rhododendron (Ericaceae), a species-rich and widely distributed genus of woody plants, is distinguished for the beautiful and diverse flowers. Rhododendron delavayi Franch. and Rhododendron irroratum Franch., are highly attractive species widely distributed in south-west China and abundant new varieties have been selected from their genetic resources.

Methods

We constructed chromosome-scale genome assemblies for Rhododendron delavayi and Rhododendron irroratum. Phylogenetic and whole-genome duplication analyses were performed to elucidate the evolutionary history of Rhododendron. Further, different types of gene duplications were identified and their contributions to gene family expansion were investigated. Finally, comprehensive characterization and evolutionary analysis of R2R3-MYB and NBS-encoding genes were conducted to explore their evolutionary patterns.

Results

The phylogenetic analysis classified Rhododendron species into two sister clades, 'rhododendrons' and 'azaleas'. Whole-genome duplication (WGD) analysis unveiled only one WGD event that occurred in Rhododendron after the ancestral γ triplication. Gene duplication and gene family expansion analyses suggested that the younger tandem and proximal duplications contributed greatly to the expansion of gene families involved in secondary metabolite biosynthesis and stress response. The candidate R2R3-MYB genes likely regulating anthocyanin biosynthesis and stress tolerance in Rhododendron will facilitate the breeding for ornamental use. NBS-encoding genes had undergone significant expansion and experienced species-specific gain and loss events in Rhododendron plants.

Conclusions

The reference genomes presented here will provide important genetic resources for molecular breeding and genetic improvement of plants in this economically important Rhododendron genus.

Genome-wide identification and characterization of members of the LEA gene family in Panax notoginseng and their transcriptional responses to dehydration of recalcitrant seeds

BACKGROUND

Late embryogenesis abundant (LEA) proteins play an important role in dehydration process of seed maturation. The seeds of Panax notoginseng (Burkill) F. H. Chen are typically characterized with the recalcitrance and are highly sensitive to dehydration. However, it is not very well known about the role of LEA proteins in response to dehydration stress in P. notoginseng seeds. We will perform a genome-wide analysis of the LEA gene family and their transcriptional responses to dehydration stress in recalcitrant P. notoginseng seeds.

RESULTS

In this study, 61 LEA genes were identified from the P. notoginseng genome, and they were renamed as PnoLEA. The PnoLEA genes were classified into seven subfamilies based on the phylogenetic relationships, gene structure and conserved domains. The PnoLEA genes family showed relatively few introns and was highly conserved. Unexpectedly, the LEA_6 subfamily was not found, and the LEA_2 subfamily contained 46 (75.4%) members. Within 19 pairs of fragment duplication events, among them 17 pairs were LEA_2 subfamily. In addition, the expression of the PnoLEA genes was obviously induced under dehydration stress, but the germination rate of P. notoginseng seeds decreased as the dehydration time prolonged.

CONCLUSIONS

We found that the lack of the LEA_6 subfamily, the expansion of the LEA_2 subfamily and low transcriptional levels of most PnoLEA genes might be implicated in the recalcitrant formation of P. notoginseng seeds. LEA proteins are essential in the response to dehydration stress in recalcitrant seeds, but the protective effect of LEA protein is not efficient. These results could improve our understanding of the function of LEA proteins in the response of dehydration stress and their contributions to the formation of seed recalcitrance.

Genome-Wide Identification of 2-Oxoglutarate and Fe (II)-Dependent Dioxygenase (2ODD-C) Family Genes and Expression Profiles under Different Abiotic Stresses in Camellia sinensis (L.)

The 2-oxoglutarate and Fe (II)-dependent dioxygenase (2ODD-C) family of 2-oxoglutarate-dependent dioxygenases potentially participates in the biosynthesis of various metabolites under various abiotic stresses. However, there is scarce information on the expression profiles and roles of 2ODD-C genes in Camellia sinensis. We identified 153 Cs2ODD-C genes from C. sinensis, and they were distributed unevenly on 15 chromosomes. According to the phylogenetic tree topology, these genes were divided into 21 groups distinguished by conserved motifs and an intron/exon structure. Gene-duplication analyses revealed that 75 Cs2ODD-C genes were expanded and retained after WGD/segmental and tandem duplications. The expression profiles of Cs2ODD-C genes were explored under methyl jasmonate (MeJA), polyethylene glycol (PEG), and salt (NaCl) stress treatments. The expression analysis showed that 14, 13, and 49 Cs2ODD-C genes displayed the same expression pattern under MeJA and PEG treatments, MeJA and NaCl treatments, and PEG and NaCl treatments, respectively. A further analysis showed that two genes, Cs2ODD-C36 and Cs2ODD-C21, were significantly upregulated and downregulated after MeJA, PEG, and NaCl treatments, indicating that these two genes played positive and negative roles in enhancing the multi-stress tolerance. These results provide candidate genes for the use of genetic engineering technology to modify plants by enhancing multi-stress tolerance to promote phytoremediation efficiency.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

Genome-Wide Identification and Expression Profiling of Glutathione S-Transferase Gene Family in Foxtail Millet (Setaria italica L.)

Glutathione S-transferases (GSTs) are a critical superfamily of multifunctional enzymes in plants. As a ligand or binding protein, GSTs regulate plant growth and development and detoxification. Foxtail millet (Setaria italica (L.) P. Beauv) could respond to abiotic stresses through a highly complex multi-gene regulatory network in which the GST family is also involved. However, GST genes have been scarcely studied in foxtail millet. Genome-wide identification and expression characteristics analysis of the foxtail millet GST gene family were conducted by biological information technology. The results showed that 73 GST genes (SiGSTs) were identified in the foxtail millet genome and were divided into seven classes. The chromosome localization results showed uneven distribution of GSTs on the seven chromosomes. There were 30 tandem duplication gene pairs belonging to 11 clusters. Only one pair of SiGSTU1 and SiGSTU23 were identified as fragment duplication genes. A total of ten conserved motifs were identified in the GST family of foxtail millet. The gene structure of SiGSTs is relatively conservative, but the number and length of exons of each gene are still different. The cis-acting elements in the promoter region of 73 SiGST genes showed that 94.5% of SiGST genes possessed defense and stress-responsive elements. The expression profiles of 37 SiGST genes covering 21 tissues suggested that most SiGST genes were expressed in multiple organs and were highly expressed in roots and leaves. By qPCR analysis, we found that 21 SiGST genes were responsive to abiotic stresses and abscisic acid (ABA). Taken together, this study provides a theoretical basis for identifying foxtail millet GST family information and improving their responses to different stresses.

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

Genome-wide identification and analysis of the PAL genes from the orchids Apostasia shenzhenica, Dendrobium catenatum and Phalaenopsis equestris

Phenylalanine ammonia-lyase (PAL) is a key gateway enzyme that connects the phenylpropanoid pathway to primary metabolism. The phenylpropanoid pathway plays a vital role in the growth and environmental adaptation of many plants leading to the production of valuable bioactive compounds with industrial and medical applications. In the present study, nine putative PAL genes from three orchids were identified; five in Apostasia shenzhenica and two each in Dendrobium catenatum and Phalaenopsis equestris. Eighteen motifs and four major conserved functional domains were identified as reported in PAL proteins of other species. All the nine PALs were stable based on their computed physicochemical properties and localized in the cytoplasm. The three-dimensional structures of PALs revealed a homo-tetrameric structure consisting of four identical subunits. A total of 21 cis-regulatory elements with known functions were identified from the promoter regions of all PALs which are responsible for various plant responses to light, stress and growth regulators like auxins, gibberellins and abscisic acid. Phylogenetic analysis showed that the studied PAL proteins clustered in two major clades (clade I and II), placing dicot and monocot PALs in two separate monophyletic clades. In silico gene expression of the identified PALs in different vegetative and reproductive tissues revealed the differential expressions based on tissue type and disclosed that the expression of PAL genes was upregulated in all the tissues examined with an exception of PePAL leaf samples where no expression was detected, however, the same being highly expressed in reproductive tissues (PePAL1-labellum; PePAL2-sepal). In case of AsPALs, the expression was found to be highest in reproductive tissues (AsPAL4-maximum in inflorescence). On the other hand, the expression of DcPALs was found to be highest in vegetative tissues (DcPAL2-maximum in root). Based on the medicinal importance of orchids and the significant role of PAL genes in synthesis of bioactive compounds, the functional characterization of PAL genes can be further exploited in genetic improvement of medicinal orchids.Communicated by Ramaswamy H. Sarma.

Zinc/iron-regulated transporter-like protein gene family in Theobroma cacao L: Characteristics, evolution, function and 3D structure analysis

The zinc/iron-regulated transporter-like protein (ZIP) gene family first identified in plants is highly distributed in the plant kingdom. This family has previously been reported to transport several essential and non-essential cationic elements, including those toxic to many economically important crops such as cacao (Theobroma cacao L.). In this article, we present a detailed study on physicochemical properties, evolution, duplication, gene structure, promoter region and TcZIP family three-dimensional protein structure. A total of 11 TcZIP genes have been identified to encode proteins from 309 to 435 aa, with localization in the plasma membrane and chloroplast, containing 6-9 putative domains (TM). Interspecies phylogenetic analysis subdivided the ZIP proteins into four groups. Segmental duplication events significantly contributed to the expansion of TcZIP genes. These genes underwent high pressure of purifying selection. The three-dimensional structure of the proteins showed that α helix conformations are predominant with several pocket sites, containing the metal binding site, with the residues leucine (LEU), alanine (ALA), glycine (GLY), serine (SER), lysine (LYS) and histidine (HIS) the most predicted. Regarding the analysis of the protein-protein interaction and enrichment of the gene ontology, four biological processes were assigned, the most important being the cation transport. These new discoveries expand the knowledge about the function, evolution, protein structures and interaction of ZIP family proteins in cacao and contribute to develop cacao genotypes enriched with important mineral nutrients as well as genotypes that bioaccumulate or exclude toxic metals.

Genomic Survey of Flavin Monooxygenases in Wild and Cultivated Rice Provides Insight into Evolution and Functional Diversities

The flavin monooxygenase (FMO) enzyme was discovered in mammalian liver cells that convert a carcinogenic compound, N-N'-dimethylaniline, into a non-carcinogenic compound, N-oxide. Since then, many FMOs have been reported in animal systems for their primary role in the detoxification of xenobiotic compounds. In plants, this family has diverged to perform varied functions like pathogen defense, auxin biosynthesis, and S-oxygenation of compounds. Only a few members of this family, primarily those involved in auxin biosynthesis, have been functionally characterized in plant species. Thus, the present study aims to identify all the members of the FMO family in 10 different wild and cultivated Oryza species. Genome-wide analysis of the FMO family in different Oryza species reveals that each species has multiple FMO members in its genome and that this family is conserved throughout evolution. Taking clues from its role in pathogen defense and its possible function in ROS scavenging, we have also assessed the involvement of this family in abiotic stresses. A detailed in silico expression analysis of the FMO family in Oryza sativa subsp. japonica revealed that only a subset of genes responds to different abiotic stresses. This is supported by the experimental validation of a few selected genes using qRT-PCR in stress-sensitive Oryza sativa subsp. indica and stress-sensitive wild rice Oryza nivara. The identification and comprehensive in silico analysis of FMO genes from different Oryza species carried out in this study will serve as the foundation for further structural and functional studies of FMO genes in rice as well as other crop types.

Acquisition of new function through gene duplication in the metallocarboxypeptidase family

Gene duplication is a key first step in the process of expanding the functionality of a multigene family. In order to better understand the process of gene duplication and its role in the formation of new enzymes, we investigated recent duplication events in the M14 family of proteolytic enzymes. Within vertebrates, four of 23 M14 genes were frequently found in duplicate form. While AEBP1, CPXM1, and CPZ genes were duplicated once through a large-scale, likely whole-genome duplication event, the CPO gene underwent many duplication events within fish and Xenopus lineages. Bioinformatic analyses of enzyme specificity and conservation suggested a greater amount of neofunctionalization and purifying selection in CPO paralogs compared with other CPA/B enzymes. To examine the functional consequences of evolutionary changes on CPO paralogs, the four CPO paralogs from Xenopus tropicalis were expressed in Sf9 and HEK293T cells. Immunocytochemistry showed subcellular distribution of Xenopus CPO paralogs to be similar to that of human CPO. Upon activation with trypsin, the enzymes demonstrated differential activity against three substrates, suggesting an acquisition of new function following duplication and subsequent mutagenesis. Characteristics such as gene size and enzyme activation mechanisms are possible contributors to the evolutionary capacity of the CPO gene.