Transcriptome analysis of the phosphate starvation response sheds light on strigolactone biosynthesis in rice.

Phosphorus (P) is a major element required for plant growth and development. To cope with P shortage, plants activate local and long-distance signaling pathways, such as an increase in the production and exudation of strigolactones (SLs). The role of the latter in mitigating P deficiency is, however, still largely unknown. To shed light on this, we studied the transcriptional response to P starvation and replenishment in wild-type rice and a SL mutant, dwarf10 (d10), and upon exogenous application of the synthetic SL GR24. P starvation resulted in major transcriptional alterations, such as the upregulation of P TRANSPORTER, SYG1/PHO81/XPR1 (SPX) and VACUOLAR PHOSPHATE EFFLUX TRANSPORTER. Gene Ontology (GO) analysis of the genes induced by P starvation showed enrichment in phospholipid catabolic process and phosphatase activity. In d10, P deficiency induced upregulation of genes enriched for sesquiterpenoid production, secondary shoot formation and metabolic processes, including lactone biosynthesis. Furthermore, several genes induced by GR24 treatment shared the same GO terms with P starvation-induced genes, such as oxidation reduction, heme binding and oxidoreductase activity, hinting at the role that SLs play in the transcriptional reprogramming upon P starvation. Gene co-expression network analysis uncovered a METHYL TRANSFERASE that displayed co-regulation with known rice SL biosynthetic genes. Functional characterization showed that this gene encodes an enzyme catalyzing the conversion of carlactonoic acid to methyl carlactonoate. Our work provides a valuable resource to further studies on the response of crops to P deficiency and reveals a tool for the discovery of SL biosynthetic genes.

Genome-Wide Identification, Characterization, and Expression Analysis of SPIRAL1 Family Genes in Legume Species.

The SPIRAL1 (SPR1) gene family encodes microtubule-associated proteins that are essential for the anisotropic growth of plant cells and abiotic stress resistance. Currently, little is known about the characteristics and roles of the gene family outside of Arabidopsis thaliana. This study intended to investigate the SPR1 gene family in legumes. In contrast to that of A. thaliana, the gene family has undergone shrinking in the model legume species Medicago truncatula and Glycine max. While the orthologues of SPR1 were lost, very few SPR1-Like (SP1L) genes were identified given the genome size of the two species. Specifically, the M. truncatula and G. max genomes only harbor two MtSP1L and eight GmSP1L genes, respectively. Multiple sequence alignment showed that all these members contain conserved N- and C-terminal regions. Phylogenetic analysis clustered the legume SP1L proteins into three clades. The SP1L genes showed similar exon-intron organizations and similar architectures in their conserved motifs. Many essential cis-elements are present in the promoter regions of the MtSP1L and GmSP1L genes associated with growth and development, plant hormones, light, and stress. The expression analysis revealed that clade 1 and clade 2 SP1L genes have relatively high expression in all tested tissues in Medicago and soybean, suggesting their function in plant growth and development. MtSP1L-2, as well as clade 1 and clade 2 GmSP1L genes, display a light-dependent expression pattern. The SP1L genes in clade 2 (MtSP1L-2, GmSP1L-3, and GmSP1L-4) were significantly induced by sodium chloride treatment, suggesting a potential role in the salt-stress response. Our research provides essential information for the functional studies of SP1L genes in legume species in the future.

Genome wide analysis of DWARF27 genes in soybean and functional characterization of GmD27c reveals eminent role of strigolactones in rhizobia interaction and nodulation in Glycine max.

BACKGROUND: Strigolactones (SLs) are newly identified hormones and their biosynthesis is stimulated under phosphate deprivation and accomplished by the action of several enzymes, including the beta-carotene isomerase DWARF27 (D27). Expression of D27 is well renowned to respond to phosphate insufficiency. However, the identification and functional analysis of the carotenoid isomerase D27 genes are not elucidated in soybean. METHODS AND RESULTS: A total of six D27 genes were identified in the soybean genome and designated on the basis of chromosomal localization. According to the findings, these genes were irregularly distributed on chromosomes, and segmental repetition led to the expansion of the soybean GmD27 gene family. Based on a neighbor-joining phylogenetic tree, the predicted D27 proteins of soybean were divided into three clades. Based on RNA seq data analysis, GmD27 genes were differently expressed in various tissues but GmD27c was the highest. Therefore, GmD27c was chosen for the additional functional study due to its rather obvious transcription in nodulation and roots. RT-qPCR results showed that GmD27c was highly expressed in different nodule stages and in response to rhizobia infection. Functional characterization of GmD27c revealed that overexpression of GmD27c led to higher nodule number, while GmD27c knockdown caused fewer nodules compared to GUS control. Furthermore, GmD27c overexpressed and knockdown lines oppositely regulated the expression of numerous nodulation genes, which are vital for the development of nodules. CONCLUSION: This study not only discovered that SL biosynthesis and signaling pathway genes are conserved, but it also revealed that SL biosynthesis gene GmD27c and legume rhizobia have close interactions in controlling plant nodule number.

Developing Genetic Engineering Techniques for Control of Seed Size and Yield.

Many signaling pathways regulate seed size through the development of endosperm and maternal tissues, which ultimately results in a range of variations in seed size or weight. Seed size can be determined through the development of zygotic tissues (endosperm and embryo) and maternal ovules. In addition, in some species such as rice, seed size is largely determined by husk growth. Transcription regulator factors are responsible for enhancing cell growth in the maternal ovule, resulting in seed growth. Phytohormones induce significant effects on entire features of growth and development of plants and also regulate seed size. Moreover, the vegetative parts are the major source of nutrients, including the majority of carbon and nitrogen-containing molecules for the reproductive part to control seed size. There is a need to increase the size of seeds without affecting the number of seeds in plants through conventional breeding programs to improve grain yield. In the past decades, many important genetic factors affecting seed size and yield have been identified and studied. These important factors constitute dynamic regulatory networks governing the seed size in response to environmental stimuli. In this review, we summarized recent advances regarding the molecular factors regulating seed size in Arabidopsis and other crops, followed by discussions on strategies to comprehend crops' genetic and molecular aspects in balancing seed size and yield.

Biochar supplementation regulates growth and heavy metal accumulation in tomato grown in contaminated soils.

Biochar application has recently gained increased attention to reclaim heavy metal degraded soils. In this context, the present study investigated the effects of biochar on the growth regulation and heavy metal accumulation in tomato grown on contaminated soils. A two-factorial design with factor A including three treatments with mine (contaminated soil) and garden soil in the following ratio viz., T1 = 1:2, T2 = 1:1, and T3 = 2:1, and garden soil only as control; whereas factor B consists of biochar amendments at three levels viz., B1 (3%), B2 (6%), and B3 (9%). Our results revealed significant negative effects of heavy metal-contaminated soil on plant growth, and besides resulted heavy metal accumulation in tomato fruit. Tomato plants showed maximum reduction of growth in T3 followed by T2, and lowest in T1, a similar pattern was found for accumulation of heavy metals in the fruit. However, the application of biochar reduced the bioavailability and accumulation of heavy metals in the tomato fruit, as well as improved plant growth in contaminated soils. Overall, among the three biochar treatments, B2 was determined as the optimum level for improved growth coupled with reduced heavy metal accumulation in the tomato fruit. Besides, biochar application decreased the daily intake of metals and human health risk index values, thus alleviating the health risk. Hence, the present study demonstrated a positive role of biochar in reclaiming heavy metal-contaminated soils and in increasing the plant growth.

Complete Genome Sequence of Cellulomonas sp. JZ18, a Root Endophytic Bacterium Isolated from the Perennial Desert Tussock-Grass Panicum turgidum.

Cellulomonas sp. JZ18 is a gram-positive, rod shaped bacterium that was previously isolated from the root endosphere of the perennial desert tussock-grass Panicum turgidum. Genome coverage of PacBio sequencing was approximately 199X. Genome assembly generated a single chromosome of 7,421,843 base pairs with a guanine-cytosine (GC) content of 75.60% with 3240 protein coding sequences, 361 pseudo genes, three ribosomal RNA operons, three non-coding RNAs and 45 transfer RNAs. Comparison of JZ18's genome with type strains from the same genus, using digital DNA-DNA hybridization and average nucleotide identity calculations, revealed that JZ18 might potentially belong to a new species. Functional analysis revealed the presence of genes that may complement previously observed biochemical and plant phenotypes. Furthermore, the presence of a number of enzymes could be of potential use in industrial processes as biocatalysts. Genome sequencing and analysis, coupled with comparative genomics, of endophytic bacteria for their potential plant growth promoting activities under different soil conditions will accelerate the knowledge and applications of biostimulants in sustainable agriculture.

Sequencing effort dictates gene discovery in marine microbial metagenomes.

Massive metagenomic sequencing combined with gene prediction methods were previously used to compile the gene catalogue of the ocean and host-associated microbes. Global expeditions conducted over the past 15 years have sampled the ocean to build a catalogue of genes from pelagic microbes. Here we undertook a large sequencing effort of a perturbed Red Sea plankton community to uncover that the rate of gene discovery increases continuously with sequencing effort, with no indication that the retrieved 2.83 million non-redundant (complete) genes predicted from the experiment represented a nearly complete inventory of the genes present in the sampled community (i.e., no evidence of saturation). The underlying reason is the Pareto-like distribution of the abundance of genes in the plankton community, resulting in a very long tail of millions of genes present at remarkably low abundances, which can only be retrieved through massive sequencing. Microbial metagenomic projects retrieve a variable number of unique genes per Tera base-pair (Tbp), with a median value of 14.7 million unique genes per Tbp sequenced across projects. The increase in the rate of gene discovery in microbial metagenomes with sequencing effort implies that there is ample room for new gene discovery in further ocean and holobiont sequencing studies.

Complete genome sequence of the endophytic bacterium Cellulosimicrobium sp. JZ28 isolated from the root endosphere of the perennial desert tussock grass Panicum turgidum.

Cellulosimicrobium sp. JZ28, a root endophytic bacterium from the desert plant Panicum turgidum, was previously identified as a plant growth-promoting bacterium. The genome of JZ28 consists of a 4378,193 bp circular chromosome and contains 3930 CDSs with an average GC content of 74.5%. Whole-genome sequencing analysis revealed that JZ28 was closely related to C. aquatile 3 bp. The genome harbors genes responsible for protection against oxidative, osmotic and salinity stresses, such as the production of osmoprotectants. It also contains genes with a role in the production of volatiles, such as hydrogen sulfide, which promote biotic and abiotic stress tolerance in plants. The presence of three copies of chitinase genes indicates a possible role of JZ28 as biocontrol agent against fungal pathogens, while a number of genes for the degradation of plant biopolymers indicates potential application in industrial processes. Genome sequencing and mining of culture-dependent collections of bacterial endophytes from desert plants provide new opportunities for biotechnological applications.

A promoter-level mammalian expression atlas.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

Complete Genome Sequence of Paenibacillus sp. JZ16, a Plant Growth Promoting Root Endophytic Bacterium of the Desert Halophyte Zygophyllum Simplex.

Paenibacillus sp. JZ16 is a gram-positive, rod-shaped, motile root endophytic bacterium of the pioneer desert halophytic plant Zygophyllum simplex. JZ16 was previously shown to promote salinity stress tolerance in Arabidopsis thaliana and possesses a highly motile phenotype on nutrient agar. JZ16 genome sequencing using PacBio generated 82,236 reads with a mean insert read length of 11,432 bp and an estimated genome coverage of 127X, resulting in a chromosome of 7,421,843 bp with a GC content of 49.25% encoding 6710 proteins, 8 rRNA operons, 117 ncRNAs and 73 tRNAs. Whole-genome sequencing analysis revealed a potentially new species for JZ16. Functional analysis revealed the presence of a number of enzymes involved in the breakdown of plant-based polymers. JZ16 could be of potential use in agricultural applications for promoting biotic and abiotic stress tolerance and for biotechnological processes (e.g., as biocatalysts for biofuel production). The culture-dependent collection of bacterial endophytes from desert plants combined with genome sequence mining provides new opportunities for industrial applications.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L.

BACKGROUND: Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana. RESULTS: We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments. CONCLUSIONS: Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Dynamic structure mediates halophilic adaptation of a DNA polymerase from the deep-sea brines of the Red Sea.

The deep-sea brines of the Red Sea are remote and unexplored environments characterized by high temperatures, anoxic water, and elevated concentrations of salt and heavy metals. This environment provides a rare system to study the interplay between halophilic and thermophilic adaptation in biologic macromolecules. The present article reports the first DNA polymerase with halophilic and thermophilic features. Biochemical and structural analysis by Raman and circular dichroism spectroscopy showed that the charge distribution on the protein's surface mediates the structural balance between stability for thermal adaptation and flexibility for counteracting the salt-induced rigid and nonfunctional hydrophobic packing. Salt bridge interactions via increased negative and positive charges contribute to structural stability. Salt tolerance, conversely, is mediated by a dynamic structure that becomes more fixed and functional with increasing salt concentration. We propose that repulsive forces among excess negative charges, in addition to a high percentage of negatively charged random coils, mediate this structural dynamism. This knowledge enabled us to engineer a halophilic version of Thermococcus kodakarensis DNA polymerase.-Takahashi, M., Takahashi, E., Joudeh, L. I., Marini, M., Das, G., Elshenawy, M. M., Akal, A., Sakashita, K., Alam, I., Tehseen, M., Sobhy, M. A., Stingl, U., Merzaban, J. S., Di Fabrizio, E., Hamdan, S. M. Dynamic structure mediates halophilic adaptation of a DNA polymerase from the deep-sea brines of the Red Sea.

Comprehensive Genomic Analyses of the OM43 Clade, Including a Novel Species from the Red Sea, Indicate Ecotype Differentiation among Marine Methylotrophs.

The OM43 clade within the family Methylophilaceae of Betaproteobacteria represents a group of methylotrophs that play important roles in the metabolism of C1 compounds in marine environments and other aquatic environments around the globe. Using dilution-to-extinction cultivation techniques, we successfully isolated a novel species of this clade (here designated MBRS-H7) from the ultraoligotrophic open ocean waters of the central Red Sea. Phylogenomic analyses indicate that MBRS-H7 is a novel species that forms a distinct cluster together with isolate KB13 from Hawaii (Hawaii-Red Sea [H-RS] cluster) that is separate from the cluster represented by strain HTCC2181 (from the Oregon coast). Phylogenetic analyses using the robust 16S-23S internal transcribed spacer revealed a potential ecotype separation of the marine OM43 clade members, which was further confirmed by metagenomic fragment recruitment analyses that showed trends of higher abundance in low-chlorophyll and/or high-temperature provinces for the H-RS cluster but a preference for colder, highly productive waters for the HTCC2181 cluster. This potential environmentally driven niche differentiation is also reflected in the metabolic gene inventories, which in the case of the H-RS cluster include those conferring resistance to high levels of UV irradiation, temperature, and salinity. Interestingly, we also found different energy conservation modules between these OM43 subclades, namely, the existence of the NADH:quinone oxidoreductase complex I (NUO) system in the H-RS cluster and the nonhomologous NADH:quinone oxidoreductase (NQR) system in the HTCC2181 cluster, which might have implications for their overall energetic yields.

BEACON: automated tool for Bacterial GEnome Annotation ComparisON.

BACKGROUND: Genome annotation is one way of summarizing the existing knowledge about genomic characteristics of an organism. There has been an increased interest during the last several decades in computer-based structural and functional genome annotation. Many methods for this purpose have been developed for eukaryotes and prokaryotes. Our study focuses on comparison of functional annotations of prokaryotic genomes. To the best of our knowledge there is no fully automated system for detailed comparison of functional genome annotations generated by different annotation methods (AMs). RESULTS: The presence of many AMs and development of new ones introduce needs to: a/ compare different annotations for a single genome, and b/ generate annotation by combining individual ones. To address these issues we developed an Automated Tool for Bacterial GEnome Annotation ComparisON (BEACON) that benefits both AM developers and annotation analysers. BEACON provides detailed comparison of gene function annotations of prokaryotic genomes obtained by different AMs and generates extended annotations through combination of individual ones. For the illustration of BEACON's utility, we provide a comparison analysis of multiple different annotations generated for four genomes and show on these examples that the extended annotation can increase the number of genes annotated by putative functions up to 27%, while the number of genes without any function assignment is reduced. CONCLUSIONS: We developed BEACON, a fast tool for an automated and a systematic comparison of different annotations of single genomes. The extended annotation assigns putative functions to many genes with unknown functions. BEACON is available under GNU General Public License version 3.0 and is accessible at: http://www.cbrc.kaust.edu.sa/BEACON/ .

Aroma Components in Horticultural Crops: Chemical Diversity and Usage of Metabolic Engineering for Industrial Applications.

Plants produce an incredible variety of volatile organic compounds (VOCs) that assist the interactions with their environment, such as attracting pollinating insects and seed dispersers and defense against herbivores, pathogens, and parasites. Furthermore, VOCs have a significant economic impact on crop quality, as well as the beverage, food, perfume, cosmetics and pharmaceuticals industries. These VOCs are mainly classified as terpenoids, benzenoids/phenylpropanes, and fatty acid derivates. Fruits and vegetables are rich in minerals, vitamins, antioxidants, and dietary fiber, while aroma compounds play a major role in flavor and quality management of these horticultural commodities. Subtle shifts in aroma compounds can dramatically alter the flavor and texture of fruits and vegetables, altering their consumer appeal. Rapid innovations in -omics techniques have led to the isolation of genes encoding enzymes involved in the biosynthesis of several volatiles, which has aided to our comprehension of the regulatory molecular pathways involved in VOC production. The present review focuses on the significance of aroma volatiles to the flavor and aroma profile of horticultural crops and addresses the industrial applications of plant-derived volatile terpenoids, particularly in food and beverages, pharmaceuticals, cosmetics, and biofuel industries. Additionally, the methodological constraints and complexities that limit the transition from gene selection to host organisms and from laboratories to practical implementation are discussed, along with metabolic engineering's potential for enhancing terpenoids volatile production at the industrial level.

Bioengineering Techniques to Improve Nitrogen Transformation and Utilization: Implications for Nitrogen Use Efficiency and Future Sustainable Crop Production.

Nitrogen (N) is crucial for plant growth and development, especially in physiological and biochemical processes such as component of different proteins, enzymes, nucleic acids, and plant growth regulators. Six categories, such as transporters, nitrate absorption, signal molecules, amino acid biosynthesis, transcription factors, and miscellaneous genes, broadly encompass the genes regulating NUE in various cereal crops. Herein, we outline detailed research on bioengineering modifications of N metabolism to improve the different crop yields and biomass. We emphasize effective and precise molecular approaches and technologies, including N transporters, transgenics, omics, etc., which are opening up fascinating opportunities for a complete analysis of the molecular elements that contribute to NUE. Moreover, the detection of various types of N compounds and associated signaling pathways within plant organs have been discussed. Finally, we highlight the broader impacts of increasing NUE in crops, crucial for better agricultural yield and in the greater context of global climate change.

Overexpression of GmWRKY172 enhances cadmium tolerance in plants and reduces cadmium accumulation in soybean seeds.

Introduction: Cadmium (Cd) stress is a significant threat to soybean production, and enhancing Cd tolerance in soybean is the focus of this study. The WRKY transcription factor family is associated with abiotic stress response processes. In this study, we aimed to identify a Cd-responsive WRKY transcription factor GmWRKY172 from soybean and investigate its potential for enhancing Cd tolerance in soybean. Methods: The characterization of GmWRKY172 involved analyzing its expression pattern, subcellular localization, and transcriptional activity. To assess the impact of GmWRKY172, transgenic Arabidopsis and soybean plants were generated and examined for their tolerance to Cd and Cd content in shoots. Additionally, transgenic soybean plants were evaluated for Cd translocation and various physiological stress indicators. RNA sequencing was performed to identify the potential biological pathways regulated by GmWRKY172. Results: GmWRKY172 was significantly upregulated by Cd stress, highly expressed in leaves and flowers, and localized to the nucleus with transcriptional activity. Transgenic plants overexpressing GmWRKY172 showed enhanced Cd tolerance and reduced Cd content in shoots compared to WT. Lower Cd translocation from roots to shoots and seeds was also observed in transgenic soybean. Under Cd stress, transgenic soybean accumulated less malondialdehyde (MDA) and hydrogen peroxide (H2O2) than WT plants, with higher flavonoid and lignin contents, and peroxidase (POD) activity. RNA sequencing analysis revealed that many stress-related pathways were regulated by GmWRKY172 in transgenic soybean, including flavonoid biosynthesis, cell wall synthesis, and peroxidase activity. Discussion: Our findings demonstrated that GmWRKY172 enhances Cd tolerance and reduces seed Cd accumulation in soybean by regulating multiple stress-related pathways, and could be a promising candidate for breeding Cd-tolerant and low Cd soybean varieties.

Mpox virus infection and drug treatment modelled in human skin organoids.

Mpox virus (MPXV) primarily infects human skin to cause lesions. Currently, robust models that recapitulate skin infection by MPXV are lacking. Here we demonstrate that human induced pluripotent stem cell-derived skin organoids are susceptible to MPXV infection and support infectious virus production. Keratinocytes, the predominant cell type of the skin epithelium, effectively support MPXV infection. Using transmission electron microscopy, we visualized the four stages of intracellular virus particle assembly: crescent formation, immature virions, mature virions and wrapped virions. Transcriptional analysis showed that MPXV infection rewires the host transcriptome and triggers abundant expression of viral transcripts. Early treatment with the antiviral drug tecovirimat effectively inhibits infectious virus production and prevents host transcriptome rewiring. Delayed treatment with tecovirimat also inhibits infectious MPXV particle production, albeit to a lesser extent. This study establishes human skin organoids as a robust experimental model for studying MPXV infection, mapping virus-host interactions and testing therapeutics.

Comprehensive Genomic Survey, Evolution, and Expression Analysis of GIF Gene Family during the Development and Metal Ion Stress Responses in Soybean.

The GIF gene family is one of the plant transcription factors specific to seed plants. The family members are expressed in all lateral organs produced by apical and floral meristems and contribute to the development of leaves, shoots, flowers, and seeds. This study identified eight GIF genes in the soybean genome and clustered them into three groups. Analyses of Ka/Ks ratios and divergence times indicated that they had undergone purifying selection during species evolution. RNA-sequence and relative expression patterns of these GmGIF genes tended to be conserved, while different expression patterns were also observed between the duplicated GIF members in soybean. Numerous cis-regulatory elements related to plant hormones, light, and stresses were found in the promoter regions of these GmGIF genes. Moreover, the expression patterns of GmGIF members were confirmed in soybean roots under cadmium (Cd) and copper (Cu) stress, indicating their potential functions in the heavy metal response in soybean. Our research provides valuable information for the functional characterization of each GmGIF gene in different legumes in the future.