Genomics-based identification of a cold adapted clade in Deinococcus.

BACKGROUND: Microbes in the cold polar and alpine environments play a critical role in feedbacks that amplify the effects of climate change. Defining the cold adapted ecotype is one of the prerequisites for understanding the response of polar and alpine microbes to climate change. RESULTS: Here, we analysed 85 high-quality, de-duplicated genomes of Deinococcus, which can survive in a variety of harsh environments. By leveraging genomic and phenotypic traits with reverse ecology, we defined a cold adapted clade from eight Deinococcus strains isolated from Arctic, Antarctic and high alpine environments. Genome-wide optimization in amino acid composition and regulation and signalling enable the cold adapted clade to produce CO2 from organic matter and boost the bioavailability of mineral nitrogen. CONCLUSIONS: Based primarily on in silico genomic analysis, we defined a potential cold adapted clade in Deinococcus and provided an updated view of the genomic traits and metabolic potential of Deinococcus. Our study would facilitate the understanding of microbial processes in the cold polar and alpine environments.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera).

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species.

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

Pan-Genome Provides Insights into Vibrio Evolution and Adaptation to Deep-Sea Hydrothermal Vents.

This study delves into the genomic features of 10 Vibrio strains collected from deep-sea hydrothermal vents in the Pacific Ocean, providing insights into their evolutionary history and ecological adaptations. Through sequencing and pan-genome analysis involving 141 Vibrio species, we found that deep-sea strains exhibit larger genomes with unique gene distributions, suggesting adaptation to the vent environment. The phylogenomic reconstruction of the investigated isolates revealed the presence of 2 main clades: The first is monophyletic, consisting exclusively of Vibrio alginolyticus, while the second forms a monophyletic clade comprising both Vibrio antiquarius and Vibrio diabolicus species, which were previously isolated from deep-sea vents. All strains carry virulence and antibiotic resistance genes related to those found in human pathogenic Vibrio species which may play a wider ecological role other than host infection in these environments. In addition, functional genomic analysis identified genes potentially related to deep-sea survival and stress response, alongside candidate genes encoding for novel antimicrobial agents. Ultimately, the pan-genome we generated represents a valuable resource for future studies investigating the taxonomy, evolution, and ecology of Vibrio species.

The loci of environmental adaptation in a model eukaryote.

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Genetic admixture drives climate adaptation in the bank vole.

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Evaluation of different bermudagrass germplasm at physiological and molecular level under shade along longitudinal and latitudinal gradients.

Responses of turfgrass to shade vary in individual species, and the degree and quality of low light; therefore, the selection of low light tolerant cultivars of turfgrass is important and beneficial for turf management rather than other practices. The stolons of thirteen bermudagrass genotypes were planted with two treatments and three replications of each treatment to establish for one month in the Yangzhou University Jiangsu China greenhouse. The established plants were transferred outside of the greenhouse, and 50% shading was applied to them with a black net. After 30 days of stress treatment, the morpho-physiological and biochemical analyses were performed. The expression of genes such as HEMA, HY5, PIF4, and Cu/ZnSOD was assessed. Cynodon dactylon is a C4, and perennial that grows as lawn grass and is used as forage. Based on different indicator measurements, the most shade-tolerant germplasm was L01 and L06 along the longitudes and L09 and L10 along the latitudes. At the same time, L02 and L08 were more susceptible, respectively. However, germplasm showed greater tolerance in higher latitudes while longitudinal plants showed less stress response. The current study aimed (1) to screen out the most shade-tolerant Cynodon dactylon genotype among 13 along longitudinal and latitudinal gradients in China. (2) to examine morpho-physiological indicators of different bermudagrassgenotypes; (3) to evaluate if and how differences in various indicators of bermudagrass correlated with geographic region. This study will significantly advance the use of Cynodon germplasm in breeding, genomics, management, nomenclature, and phylogeographical study. It will decisively define whether natural selection and migration can drive evolutionary responses for populations to adapt to their new environments effectively.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen.

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Natural variation of domestication-related genes contributed to latitudinal expansion and adaptation in soybean.

Soybean is a major source of protein and edible oil worldwide. Originating from the Huang-Huai-Hai region, which has a temperate climate, soybean has adapted to a wide latitudinal gradient across China. However, the genetic mechanisms responsible for the widespread latitudinal adaptation in soybean, as well as the genetic basis, adaptive differentiation, and evolutionary implications of theses natural alleles, are currently lacking in comprehensive understanding. In this study, we examined the genetic variations of fourteen major gene loci controlling flowering and maturity in 103 wild species, 1048 landraces, and 1747 cultivated species. We found that E1, E3, FT2a, J, Tof11, Tof16, and Tof18 were favoured during soybean improvement and selection, which explained 75.5% of the flowering time phenotypic variation. These genetic variation was significantly associated with differences in latitude via the LFMM algorithm. Haplotype network and geographic distribution analysis suggested that gene combinations were associated with flowering time diversity contributed to the expansion of soybean, with more HapA clustering together when soybean moved to latitudes beyond 35°N. The geographical evolution model was developed to accurately predict the suitable planting zone for soybean varieties. Collectively, by integrating knowledge from genomics and haplotype classification, it was revealed that distinct gene combinations improve the adaptation of cultivated soybeans to different latitudes. This study provides insight into the genetic basis underlying the environmental adaptation of soybean accessions, which could contribute to a better understanding of the domestication history of soybean and facilitate soybean climate-smart molecular breeding for various environments.

Genome-wide identification and expression analysis revealed key transcription factors as potential regulators of high-temperature adaptation of Coriolopsis trogii.

Transcription factors (TFs) play a crucial role in gene expression, and studying them can lay the foundation for future research on the functional characterization of TFs involved in various biological processes. In this study, we conducted a genome-wide identification and analysis of TFs in the thermotolerant basidiomycete fungus, Coriolopsis trogii. The TF repertoire of C. trogii consisted of 350 TFs, with C2H2 and Zn2C6 being the largest TF families. When the mycelia of C. trogii were cultured on PDA and transferred from 25 to 35 °C, 14 TFs were up-regulated and 14 TFs were down-regulated. By analyzing RNA-seq data from mycelia cultured at different temperatures and under different carbon sources, we identified 22 TFs that were differentially expressed in more than three comparisons. Co-expression analysis revealed that seven differentially expressed TFs, including four Zn2C6s, one Hap4_Hap_bind, one HMG_box, and one Zinc_knuckle, showed significant correlation with 729 targeted genes. Overall, this study provides a comprehensive characterization of the TF family and systematically screens TFs involved in the high-temperature adaptation of C. trogii, laying the groundwork for further research into the specific roles of TFs in the heat tolerance mechanisms of filamentous fungi.

Deep metagenomic sequencing unveils novel SAR202 lineages and their vertical adaptation in the ocean.

SAR202 bacteria in the Chloroflexota phylum are abundant and widely distributed in the ocean. Their genome coding capacities indicate their potential roles in degrading complex and recalcitrant organic compounds in the ocean. However, our understanding of their genomic diversity, vertical distribution, and depth-related metabolisms is still limited by the number of assembled SAR202 genomes. In this study, we apply deep metagenomic sequencing (180 Gb per sample) to investigate microbial communities collected from six representative depths at the Bermuda Atlantic Time Series (BATS) station. We obtain 173 SAR202 metagenome-assembled genomes (MAGs). Intriguingly, 154 new species and 104 new genera are found based on these 173 SAR202 genomes. We add 12 new subgroups to the current SAR202 lineages. The vertical distribution of 20 SAR202 subgroups shows their niche partitioning in the euphotic, mesopelagic, and bathypelagic oceans, respectively. Deep-ocean SAR202 bacteria contain more genes and exhibit more metabolic potential for degrading complex organic substrates than those from the euphotic zone. With deep metagenomic sequencing, we uncover many new lineages of SAR202 bacteria and their potential functions which greatly deepen our understanding of their diversity, vertical profile, and contribution to the ocean's carbon cycling, especially in the deep ocean.

Discovering novel genomic regions explaining adaptation of bread wheat to conservation agriculture through GWAS.

To sustainably increase wheat yield to meet the growing world population's food demand in the face of climate change, Conservation Agriculture (CA) is a promising approach. Still, there is a lack of genomic studies investigating the genetic basis of crop adaptation to CA. To dissect the genetic architecture of 19 morpho-physiological traits that could be involved in the enhanced adaptation and performance of genotypes under CA, we performed GWAS to identify MTAs under four contrasting production regimes viz., conventional tillage timely sown (CTTS), conservation agriculture timely sown (CATS), conventional tillage late sown (CTLS) and conservation agriculture late sown (CALS) using an association panel of 183 advanced wheat breeding lines along with 5 checks. Traits like Phi2 (Quantum yield of photosystem II; CATS:0.37, CALS: 0.31), RC (Relative chlorophyll content; CATS:55.51, CALS: 54.47) and PS1 (Active photosystem I centers; CATS:2.45, CALS: 2.23) have higher mean values in CA compared to CT under both sowing times. GWAS identified 80 MTAs for the studied traits across four production environments. The phenotypic variation explained (PVE) by these QTNs ranged from 2.15 to 40.22%. Gene annotation provided highly informative SNPs associated with Phi2, NPQ (Quantum yield of non-photochemical quenching), PS1, and RC which were linked with genes that play crucial roles in the physiological adaptation under both CA and CT. A highly significant SNP AX94651261 (9.43% PVE) was identified to be associated with Phi2, while two SNP markers AX94730536 (30.90% PVE) and AX94683305 (16.99% PVE) were associated with NPQ. Identified QTNs upon validation can be used in marker-assisted breeding programs to develop CA adaptive genotypes.

Across two continents: The genomic basis of environmental adaptation in house mice (Mus musculus domesticus) from the Americas.

Replicated clines across environmental gradients can be strong evidence of adaptation. House mice (Mus musculus domesticus) were introduced to the Americas by European colonizers and are now widely distributed from Tierra del Fuego to Alaska. Multiple aspects of climate, such as temperature, vary predictably across latitude in the Americas. Past studies of North American populations across latitudinal gradients provided evidence of environmental adaptation in traits related to body size, metabolism, and behavior and identified candidate genes using selection scans. Here, we investigate genomic signals of environmental adaptation on a second continent, South America, and ask whether there is evidence of parallel adaptation across multiple latitudinal transects in the Americas. We first identified loci across the genome showing signatures of selection related to climatic variation in mice sampled across a latitudinal transect in South America, accounting for neutral population structure. Consistent with previous results, most candidate SNPs were in putatively regulatory regions. Genes that contained the most extreme outliers relate to traits such as body weight or size, metabolism, immunity, fat, eye function, and the cardiovascular system. We then compared these results with the results of analyses of published data from two transects in North America. While most candidate genes were unique to individual transects, we found significant overlap among candidate genes identified independently in the three transects. These genes are diverse, with functions relating to metabolism, immunity, cardiac function, and circadian rhythm, among others. We also found parallel shifts in allele frequency in candidate genes across latitudinal gradients. Finally, combining data from all three transects, we identified several genes associated with variation in body weight. Overall, our results provide strong evidence of shared responses to selection and identify genes that likely underlie recent environmental adaptation in house mice across North and South America.

Adaptive strategies to drought stress in grasses of the poaceae family under climate change: Physiological, genetic and molecular perspectives: A review.

Drought stress is one of the most critical abiotic factors which negatively impacts on growth, productivity, and survival of plants. Grass species have an important role in the sustainable intensification of cropping systems. This review focus on the specific drought tolerance characteristics in grass species and application of prevalent classical and molecular methods for genetic improvement of them to drought stress. Generally, grass species adapt to drought stress by utilizing more than one strategy including of changes in the root growth, photosynthetic pigments, activation of antioxidant enzymes, and accumulation of compatible osmolytes. They also have other specific characteristics consisted of summer dormancy, drought recovery, and persistence, which lead to drought adaptation after prolonged drought. Studies on different grasses, indicated that most of above mentioned traits usually have positive correlation with drought tolerance. Also, high heritability has been reported for most of them in different grasses. Therefore, an effective index might be considering in identification of drought tolerance genotypes. Recently, high-throughput imaging phenotyping and advanced molecular techniques such as genotyping-by-sequencing (GBS), RNA sequencing, genome-wide association study, and genome editing help conventional breeding methods to increase the accuracy, selection efficiency, genetic gains, and speed of breeding programs for developing drought tolerant cultivars.

Adaptation of the Invasive Plant Sphagneticola trilobata to Flooding Stress by Hybridization with Native Relatives.

Hybridization is common between invasive and native species and may produce more adaptive hybrids. The hybrid (Sphagneticola × guangdongensis) of Sphagneticola trilobata (an invasive species) and S. calendulacea (a native species) was found in South China. In this study, S. trilobata, S. calendulacea, and Sphagneticola × guangdongensis were used as research materials to explore their adaptability to flooding stress. Under flooding stress, the ethylene content and the expression of key enzyme genes related to ethylene synthesis in Sphagneticola × guangdongensis and S. calendulacea were significantly higher than those in S. trilobata. A large number of adventitious roots and aerenchyma were generated in Sphagneticola × guangdongensis and S. calendulacea. The contents of reactive oxygen species and malondialdehyde in Sphagneticola × guangdongensis and S. calendulacea were lower than those in S. trilobata, and the leaves of S. trilobata were the most severely damaged under flooding stress. The results indicate that hybridization catalyzed the tolerance of Sphagneticola × guangdongensis to flooding stress, and the responses of Sphagneticola × guangdongensis to flooding stress were more similar to that of its native parent. This suggests that hybridization with native relatives is an important way for invasive species to overcome environmental pressure and achieve invasion.

Patterns and effects of gene flow on adaptation across spatial scales: implications for management.

Gene flow can have rapid effects on adaptation and is an important evolutionary tool available when undertaking biological conservation and restoration. This tool is underused partly because of the perceived risk of outbreeding depression and loss of mean fitness when different populations are crossed. In this article, we briefly review some theory and empirical findings on how genetic variation is distributed across species ranges, describe known patterns of gene flow in nature with respect to environmental gradients, and highlight the effects of gene flow on adaptation in small or stressed populations in challenging environments (e.g., at species range limits). We then present a case study involving crosses at varying spatial scales among mountain populations of a trigger plant (Stylidium armeria: Stylidiaceae) in the Australian Alps to highlight how some issues around gene flow effects can be evaluated. We found evidence of outbreeding depression in seed production at greater geographic distances. Nevertheless, we found no evidence of maladaptive gene flow effects in likelihood of germination, plant performance (size), and performance variance, suggesting that gene flow at all spatial scales produces offspring with high adaptive potential. This case study demonstrates a path to evaluating how increasing sources of gene flow in managed wild and restored populations could identify some offspring with high fitness that could bolster the ability of populations to adapt to future environmental changes. We suggest further ways in which managers and researchers can act to understand and consider adaptive gene flow in natural and conservation contexts under rapidly changing conditions.

Genome evolution and transcriptome plasticity is associated with adaptation to monocot and dicot plants in Colletotrichum fungi.

BACKGROUND: Colletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes. RESULTS: Here we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation. CONCLUSIONS: These results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

Identifying critical growth stage and resilient genotypes in cowpea under drought stress contributes to enhancing crop tolerance for improvement and adaptation in Cameroon.

Drought stress following climate change is likely a scenario that will have to face crop growers in tropical regions. In mitigating this constraint, the best option should be the selection and use of resilient varieties that can withstand drought threats. Therefore, a pot experiment was conducted under greenhouse conditions at the Research and Teaching Farm of the Faculty of Agronomy and Agricultural Sciences of the University of Dschang. The objectives are to identify sensitive growth stage, to identify drought-tolerant genotypes with the help of yield-based selection indices and to identify suitable selection indices that are associated with yield under non-stress and stress circumstances. Eighty-eight cowpea genotypes from the sahelian and western regions of Cameroon were subjected to drought stress at vegetative (VDS) and flowering (FDS) stages by withholding water for 28 days, using a split plot design with two factors and three replications. Seed yields under stress (Ys) and non-stress (Yp) conditions were recorded. Fifteen drought indices were calculated for the two drought stress levels against the yield from non-stress plants. Drought Intensity Index (DII) under VDS and FDS were 0.71 and 0.84 respectively, indicating severe drought stress for both stages. However, flowering stage was significantly more sensitive to drought stress compared to vegetative stage. Based on PCA and correlation analysis, Stress Tolerance Index (STI), Relative Efficiency Index (REI), Geometric Mean Productivity (GMP), Mean Productivity (MP), Yield Index (YI) and Harmonic Mean (HM) correlated strongly with yield under stress and non-stress conditions and are therefore suitable to discriminate high-yielding and tolerant genotypes under both stress and non-stress conditions. Either under VDS and FDS, CP-016 exhibited an outstanding performance under drought stress and was revealed as the most drought tolerant genotype as shown by ranking, PCA and cluster analysis. Taking into account all indices, the top five genotypes namely CP-016, CP-021, MTA-22, CP-056 and CP-060 were identified as the most drought-tolerant genotypes under VDS. For stress activated at flowering stage (FDS), CP-016, CP-056, CP-021, CP-028 and MTA-22 were the top five most drought-tolerant genotypes. Several genotypes with insignificant Ys and irrelevant rank among which CP-037, NDT-001, CP-036, CP-034, NDT-002, CP-031, NDT-011 were identified as highly drought sensitive with low yield stability. This study identified the most sensitive stage and drought tolerant genotypes that are proposed for genetic improvement of cowpea.

Pseudomonas putida Dynamics of Adaptation under Prolonged Resource Exhaustion.

Many nonsporulating bacterial species survive prolonged resource exhaustion, by entering a state termed long-term stationary phase. Here, we performed long-term stationary phase evolutionary experiments on the bacterium Pseudomonas putida, followed by whole-genome sequencing of evolved clones. We show that P. putida is able to persist and adapt genetically under long-term stationary phase. We observed an accumulation of mutations within the evolving P. putida populations. Within each population, independently evolving lineages are established early on and persist throughout the 4-month-long experiment. Mutations accumulate in a highly convergent manner, with similar loci being mutated across independently evolving populations. Across populations, mutators emerge, that due to mutations within mismatch repair genes developed a much higher rate of mutation than other clones with which they coexisted within their respective populations. While these general dynamics of the adaptive process are quite similar to those we previously observed in the model bacterium Escherichia coli, the specific loci that are involved in adaptation only partially overlap between P. putida and E. coli.

Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions.

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.