# Mapping of QTLs Associated with Biological Nitrogen Fixation Traits in Peanuts (Arachis hypogaea L.) Using an Interspecific Population Derived from the Cross between the Cultivated Species and Its Wild Ancestors

Darius T Nzepang, Djamel Gully, Joël R Nguepjop, Arlette Zaiya Zazou, Hodo-Abalo Tossim, Aissatou Sambou, Jean-François Rami, Valérie Hocher, Saliou Fall, Sergio Svistoonoff, et al.

► To cite this version:

Darius T Nzepang, Djamel Gully, Joël R Nguepjop, Arlette Zaiya Zazou, Hodo-Abalo Tossim, et al. Mapping of QTLs Associated with Biological Nitrogen Fixation Traits in Peanuts (Arachis hypogaea L.) Using an Interspecific Population Derived from the Cross between the Cultivated Species and Its Wild Ancestors. Genes, 2023, 14 (4), pp.797. 10.3390/genes14040797. [https://doi.org/10.3390/genes14040797](https://doi.org/10.3390/genes14040797)

HAL Id: hal-04119347

https://hal.inrae.fr/hal-04119347

Submitted on 6 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Mapping of QTLs Associated with Biological Nitrogen Fixation Traits in Peanuts (*Arachis hypogaea* L.) Using an Interspecific Population Derived from the Cross between the Cultivated Species and Its Wild Ancestors

Darius T. Nzepang<sup>1,2,3,4,5 ID</sup>, Djamel Gully<sup>2,4 ID</sup>, Joël R. Nguepjop<sup>1,5,6,7</sup>, Arlette Zaiya Zazou<sup>8</sup>, Hodo-Abalo Tossim<sup>1,5</sup>, Aissatou Sambou<sup>1,5</sup>, Jean-François Rami<sup>5,6,7</sup>, Valerie Hocher<sup>2,3,4</sup>, Saliou Fall<sup>3,4</sup>, Sergio Svistoonoff<sup>2,4 ID</sup> and Daniel Fonceka<sup>1,5,6,7,* ID</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Nzepang, D.T.; Gully, D.; Nguepjop, J.R.; Zaiya Zazou, A.; Tossim, H.-A.; Sambou, A.; Rami, J.-F.; Hocher, V.; Fall, S.; Svistoonoff, S.; et al. Mapping of QTLs Associated with Biological Nitrogen Fixation Traits in Peanuts (*Arachis hypogaea* L.) Using an Interspecific Population Derived from the Cross between the Cultivated Species and Its Wild Ancestors. *Genes* **2023**, *14*, 797.
https://doi.org/10.3390/genes14040797

Academic Editors: Liying Yan and Huaiyong Luo

Received: 20 February 2023

Revised: 16 March 2023

Accepted: 22 March 2023

Published: 26 March 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Peanuts (*Arachis hypogaea* L.) are an allotetraploid grain legume mainly cultivated by poor farmers in Africa, in degraded soil and with low input systems. Further understanding nodulation genetic mechanisms could be a relevant option to facilitate the improvement of yield and lift up soil without synthetic fertilizers. We used a subset of 83 chromosome segment substitution lines (CSSLs) derived from the cross between a wild synthetic tetraploid AiAd (*Arachis ipaensis* × *Arachis duranensis*)<sup>4×</sup> and the cultivated variety Fleur11, and evaluated them for traits related to BNF under shade-house conditions. Three treatments were tested: without nitrogen; with nitrogen; and without nitrogen, but with added0 *Bradyrhizobium vignae* strain ISRA400. The leaf chlorophyll content and total biomass were used as surrogate traits for BNF. We found significant variations for both traits specially linked to BNF, and four QTLs (quantitative trait loci) were consistently mapped. At all QTLs, the wild alleles decreased the value of the trait, indicating a negative effect on BNF. A detailed characterization of the lines carrying those QTLs in controlled conditions showed that the QTLs affected the nitrogen fixation efficiency, nodule colonization, and development. Our results provide new insights into peanut nodulation mechanisms and could be used to target BNF traits in peanut breeding programs.

**Keywords:** peanut; symbiosis; wild species; sustainable agriculture; plant nutrition

## 1. Introduction

Nitrogen is one of the most limiting nutrients in agriculture. Species of the legume family can associate with beneficial soil bacteria called “rhizobia”, which reduce atmospheric nitrogen ($N_2$) into nitrogen assimilable by the plant in exchange for a supply of carbon [1]. A legume-rhizobia symbiosis results in biological nitrogen fixation (BNF) inside a new organ called a nodule. This symbiosis provides approximately 40 million tons of nitrogen (N) annually to agricultural systems [2]. In farming systems using low amounts of synthetic fertilizers derived from fossil fuels (Haber-Bosh), BNF is key to the growth

<footnote n="1">¹ Centre d'Etudes Régional pour l'Amélioration de l'Adaptation à la Sécheresse, CERAAS-Route de Khombole, Thiès BP 3320, Senegal</footnote>

<footnote n="2">² PHIM Plant Health Institute, Univ Montpellier, IRD, CIRAD, INRAE, Institut Agro, Montpellier, France</footnote>

<footnote n="3">³ Laboratoire Commun de Microbiologie (LCM) (IRD/ISRA/UCAD), Centre de Recherche de Bel Air, Dakar BP 1386, Senegal</footnote>

<footnote n="4">⁴ Laboratoire Mixte International Adaptation des Plantes et Microorganismes Associés aux Stress Environnementaux (LAPSE), Centre de Recherche de Bel Air, Dakar CP 18524, Senegal</footnote>

<footnote n="5">⁵ Dispositif de Recherche et de Formation en Partenariat, Innovation et Amélioration Variétale en Afrique de l'Ouest (IAVAO), CERAAS Route de Khombole, Thiès BP 3320, Senegal</footnote>

<footnote n="6">⁶ CIRAD, UMR AGAP, F-34398 Montpellier, France</footnote>

<footnote n="7">⁷ AGAP, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, France</footnote>

<footnote n="8">⁸ Institute of Agricultural Research for Development (IRAD) (IRAD), Maroua, Cameroon</footnote>

<footnote>* Correspondence: daniel.fonceka@cirad.fr</footnote>
and production of legumes and also serves as a significant N source for associated or subsequent crops such as cereals [3]. Deciphering the genetic mechanisms that regulate this biological process is essential for its optimization and for the development of sustainable agriculture. However, although BNF is a well-studied process in model legumes and a few crops, breeding programs have rarely targeted this trait.

Domestication is the selection, modification, and adoption of wild plant species with traits of interest to humans [4]. Over 12,000 years ago, hundreds wild plant species have been transformed into improved crop plants [5]. Following domestication, modern agriculture practices have allowed yields to increase considerably, but have also led to a high reduction in the genetic diversity of cultivated plants as compared to wild relatives of crops [6,7]. A major consequence of the recent human selection pressure on crops is the alteration of the beneficial plant–microorganism associations in the soil and, consequently, an increase in the dependency on synthetic fertilizers, which are costly for farmers and harmful to the environment [8]. Mutch and Young [9] showed in peas (*Pisum sativum*) and broad beans (*Vicia faba* L.) that cultivated species are less promiscuous in their interaction with indigenous rhizobia compared to wild species. A reduction in host range was also reported between cultivated and wild species for chickpeas (*Cicer arietinum* L.) [10] and soybeans (*Glycine max* L.) [11], suggesting that the domestication process and human selection have affected the compatibility of the host–rhizobia interaction. Kiers et al. [12] showed that old soybean cultivars produce more seeds than newer cultivars when infected with a mixture of effective and ineffective rhizobia, suggesting that the older cultivars have an ability to select best-matching partners from a pool of compatible rhizobia. In peas, alfalfa (*Medicago sativa* L.), and common beans (*Phaseolus vulgaris* L.), the wild species and landraces exceed the improved varieties in nitrogen fixation efficiency [13]. These results suggest that wild legume species could be used as reservoirs of favorable alleles to improve the symbiotic potential of their cultivated counterparts [14,15]. In contrast, Muñoz et al. [16] reported that cultivated soybeans performed better than wild soybeans regarding BNF and identified quantitative trait loci (QTLs) for traits related to BNF, which showed a very low diversity among cultivated soybeans, possibly as a result of human selection during the domestication of soybeans. The overall impact of domestication on the symbioses of legume crops is difficult to assess, as very little information is available for other species.

Cultivated peanuts (*A. hypogaea* L.) are an oilseed legume crop that can meet up to 68% of its nitrogen needs via BNF [2], thanks to its interaction with rhizobia of the genus *Bradyrhizobium* [17,18]. Despite the low genetic diversity at the DNA level [19], several studies have reported phenotypic variations within the cultivated compartment of peanuts for key agronomic traits, including the ones linked to BNF such as the nodule weight and number, the biomass, nitrogenase activity, the total nitrogen content, and the amount of fixed nitrogen [20–23]. These studies suggest that the symbiotic potential of peanuts could be enhanced by exploiting the genetic variation in these traits. However, as BNF happens in nodules that are underground organs, it is not a straightforward trait to phenotype for; consequently, it remains difficult to target in breeding programs. In addition, BNF is subjected to environmental variation, meaning that multiple years and location trials are needed to better capture the trait variation. These constraints contribute to a limited knowledge of the genetic mechanisms linked to the nodulation of peanuts and other legume crops. The genetic and morphological characterization of non-nodulating and nodulating lines in peanuts suggests that the nodulation process is controlled by several genes [24–26]. Microsynteny and a transcriptome analysis showed that orthologs of almost all the symbiotic genes described in model legumes are present in *Arachis* sp. [27–30]. Peanuts belong to the dalbergioid clade, where rhizobia invade plant roots intercellularly through natural fissures at the lateral root base called “cracks” [31–33]. The molecular basis of this mode of rhizobia infection remains poorly understood.

Peanuts appeared about 9400 years ago from an interspecific hybridization between two wild diploid species, *A. duranensis* (genome A) and *A. ipaensis* (genome B), followed by a duplication of the chromosomes [34]. The combined effects of this single tetraploidization
event, domestication, and selection have greatly narrowed its genetic basis. Advances in
genetics and genomics offer opportunities to dive into feral peanut diversity to mine new
alleles for breeding [35,36]. Interspecific mapping populations such as advanced backcross
(AB-QTL) and chromosome segment substitution lines (CSSLs) were developed and used
in peanuts to identify several QTLs related to yield components, disease resistance, oil
quality, and plant morphology [37–40]. However, to our knowledge, interspecific mapping
populations have not been exploited to investigate the effect of wild species alleles on BNF.

In this study, we used a CSSL population derived from the cross between a cultivated
peanut variety and a synthetic tetraploid, created by gathering the genomes of the most
probable wild diploid ancestors of the cultivated species, to target traits related to BNF.
Genotype–phenotype associations allowed for the identification of QTLs for BNF varia-
tion. Our results provide new insights into the genetics of peanut nodulation and useful
information for pre-breeding programs.

**2. Materials and Methods**

2.1. *Plant Materials*

A subset of 83 CSSLs was selected from a population of 122 CSSLs (BC₄F₃) developed at CERAAS (Centre d'Etudes Régional de Recherche pour l'Amélioration de l'Adaptation à la Sécheresse). These CSSLs were derived from an interspecific cross between a cultivated Fleur11 variety, used as the female parent, and the synthetic tetraploid AiAd, which is a combination of the genomes of (A. ipaensis KG30076 (BB genome) × A. duranensis V14167 (AA genome))⁴×, used as the male parent [38]. Fleur11 is a Spanish-type variety with a short cycle (90 days) that is cultivated in West Africa for its agronomic performance and its moderate tolerance to drought. The selected CSSLs constitute the minimum set of lines that provide optimal genomic coverage. The wild introgression fragments have an average length of ~39.2 cM and cover ~88.7% of the genetic map. Information regarding the genotyping and markers used to develop this population is available in our previous work [38].

## 2.2. *Phenotyping of CSSL Population for Traits Associated with BNF*

The 83 CSSLs and the recurrent parent Fleur11, used as a control, were evaluated for their BNF capacity under shade-house conditions in December 2017 and March 2018 in Senegal. Seeds were sterilized with calcium hypochlorite and germinated in Petri dishes, and seedlings were transferred manually into pots containing a leached and sterilized substrate with a low nitrogen content, as described in [18]. Three treatments were applied: (i) without nitrogen (negative control or -N), where no external nitrogen supply was added to the substrate; (ii) with nitrogen (positive control or +N), where the equivalent of 100 kg/ha of urea was applied in two stages (at one and six weeks) to each pot; and (iii) without nitrogen + inoculation, where 5 mL of an inoculum of *B. vignae* strain ISRA400 (-N + ISRA400) was added to each pot. ISRA400 is an efficient bacterial strain that was isolated from the nodules of peanut plants grown in Senegal [18]. The inoculum was produced by cultivating the bacteria in yeast extract mannitol medium [41] in the dark for 4 days at 28 °C, shaken at 250 rpm, until obtaining cell suspensions with an absorbance of ~0.7 at 600 nm. The inoculum was applied one week after the seedlings were transferred to the pots. The experimental design used was an α-lattice with 3 replications and 6 blocks of 15 genotypes per replication. The studied factor was the "genotype" with 90 levels (83 CSSLs and the recurrent parent Fleur11 repeated 7 times), totaling 270 plants per treatment.

For each treatment, the leaf chlorophyll content and total dry biomass were evaluated.
The leaf chlorophyll content was indirectly estimated weekly on the third youngest leaf
using SPAD (soil plant analysis development; Konica-Minolta, Okasaka, Japan). After nine
weeks, the plants were harvested and dried at 65 °C for 72 h before weighing to obtain the
total dry biomass.
## 2.3. Statistical Analysis and QTL Mapping

Descriptive statistics (mean and standard deviation) were calculated for each trait and each treatment in both experiments. The normality and homogeneity of variance were checked with the Kolmogorov-Smirnov and Levene tests, respectively. An analysis of variance (ANOVA) was performed for each experiment to estimate the genotypes and replication effects on each trait using a standard linear model:

$$y_{ijk} = \mu + G_i + r_j + b_{jk} + \epsilon_{ijk}$$

where $y_{ijk}$ = observed value of genotype i in replication j and block k; $\mu$ = mean value of the population; $G_i$ = the effect of genotype i; $r_j$ = the effect of replication j; $b_{jk}$ = the effect of block k in replication j; and $\epsilon_{ijk}$ = residual error.

Estimates of broad-sense heritability were calculated for each treatment in each environment using the formula:

$$h^2 = 1 - (1/F)$$

where F is the F-value associated with the genotype effect [42].

When a significant genotype effect was found for a given trait, each CSSL was compared to the recurrent parent Fleur11 with Dunnett's multiple comparisons test at $p < 0.05$ using the package multcomp of the R software [43], as described in [38]. When a given CSSL carried one unique wild chromosome segment (i.e., target wild segment) and showed a significant difference compared to Fleur11, we assumed that the difference was linked to the presence of the wild segment, and consequently, that at least one QTL was located in the wild segment. When a given line harbored the target wild segment and other segments elsewhere in the genome, the QTL assignment was performed in two steps. We first compared the phenotype of the line with neighboring lines that had overlapping fragments. If the phenotypes were similar, the QTL was attributed to the target wild chromosome segment. If the phenotypes were different, the comparison was made with lines that carry similar regions elsewhere in the genome. In some particular cases, the genotype/phenotype comparison did not allow for the unambiguous assignment of a QTL in one wild fragment; therefore, the QTL was considered a putative QTL and was associated with all the wild chromosome fragments present in that line.

The QTL effect for a given trait was calculated according to the formula below [38]:

$$QTL \ effect(\%) = \frac{[Mean(CSSL) - Mean(Fleur11)] \times 100}{Mean(Fleur11)}$$

## 2.4. Cytological Analysis of Nodules

A cytological analysis was performed to provide information about the QTL effect on the shape, structure, and functioning of nodules. Fresh nodules obtained in controlled conditions were harvested, embedded in 4% agar, and sliced into 80 µm sections using a VT1000S vibratome (Leica Biosystems, Nussloch, Germany). The sections were successively incubated in live/dead staining solutions (5 µM SYTO 9 and 30 µM propidium iodide in 50 mM Tris (pH 7.0) buffer; Live/Dead BacLight, Invitrogen, Carlsbad, CA, USA) for 15 min, and calcofluor white M2R (0.01% wt/vol in 10 mM phosphate saline buffer, (Sigma, Germany) for 15 min in the dark [44]. The stained sections were mounted in saline-glycerol phosphate buffer (vol/vol) and observed using an LSM 700 confocal laser scanning microscope (Carl Zeiss, Germany). SYTO9 and propidium iodide were excited to 488 and 555 nm for an emission range of 490-522 nm and 555-700 nm, respectively. For calcofluor, the excitation wavelength was 405 nm and the emission range was 405-470 nm. Nodule sections were also observed under an AZ differential interference microscope (Nikon, Tokyo Japan) equipped with a DS-Ri2 camera to describe the nodule structure.
## 3. Results

### 3.1. Phenotypic Variation in BNF in the CSSL Population under Shade-House Conditions

In order to investigate the phenotypic variations in BNF, the CSSL population was grown in a shade house under three treatments (−N: without nitrogen, +N: with nitrogen, and −N+ISRA400: without nitrogen + *B. vignae* ISRA400), and two traits related to BNF, the leaf chlorophyll content and the biomass, were recorded. The experiment was conducted in two consecutive environments (December 2017 and March 2018). We observed a gradual decrease in the average leaf chlorophyll content in the negative control treatment (−N), while the chlorophyll content remained almost constant in the positive control treatment (+N) (Supplementary Figure S1). In the inoculated plants (−N+ISRA400), the leaf chlorophyll content decreased during the three first weeks, and then started to increase 4 weeks after inoculation, indicating the establishment of a functional symbiotic interaction that supplied nitrogen to the plants (Supplementary Figure S1). The fertilized and inoculated plants produced more biomass in comparison to the −N treatment (Table 1). The results of the ANOVAs showed a significant genotype effect for the chlorophyll content only in the +N and −N+ISRA400 treatments, while for the total biomass, a significant genotype effect was observed in all treatments (Table 1). The broad-sense heritability was estimated for each trait in each environment. We observed null, moderate, and high broad-sense heritability estimates for the chlorophyll content in the −N, +N, and −N+ISRA400 treatments, respectively, in both environments. For the total biomass, a moderate heritability was found for all the treatments in both environments (Table 1).

<table><caption>**Table 1.** Descriptive statistics and broad-sense heritability in 2017 and 2018 experiments for traits related to BNF: leaf chlorophyll content and total biomass.</caption><thead><tr><th rowspan="2">Trait</th><th rowspan="2">Treatment</th><th colspan="4">2017</th><th colspan="4">2018</th></tr><tr><th>Mean ± SD</th><th>F</th><th>Pr</th><th>h<sup>2</sup></th><th>Mean ± SD</th><th>F</th><th>Pr</th><th>h<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="3">SPAD</td><td>−N</td><td>26.2 ± 7.82</td><td>1.017</td><td>0.456</td><td>0</td><td>14.8 ± 5.10</td><td>1</td><td>0.583</td><td>0</td></tr><tr><td>+N</td><td>43.5 ± 3.79</td><td>1.742</td><td>0.001 **</td><td>0.43</td><td>40.6 ± 2.06</td><td>2</td><td>&lt;0.001 ***</td><td>0.50</td></tr><tr><td>−N+ISRA400</td><td>41.6 ± 6.42</td><td>3.843</td><td>&lt;0.001 ***</td><td>0.74</td><td>39.8 ± 6.48</td><td>11.358</td><td>&lt;0.001 ***</td><td>0.91</td></tr><tr><td rowspan="3">TB</td><td>−N</td><td>2.88 ± 1.26</td><td>1.611</td><td>0.005 **</td><td>0.38</td><td>2.69 ± 0.96</td><td>2.486</td><td>&lt;0.001 ***</td><td>0.60</td></tr><tr><td>+N</td><td>5.75 ± 1.62</td><td>1.310</td><td>0.073</td><td>0.24</td><td>7.98 ± 2.07</td><td>1.446</td><td>0.023 *</td><td>0.30</td></tr><tr><td>−N+ISRA400</td><td>4.36 ± 1.63</td><td>2.324</td><td>&lt;0.001 ***</td><td>0.57</td><td>3.90 ± 1.47</td><td>1.979</td><td>&lt;0.001 ***</td><td>0.49</td></tr></tbody></table>

SPAD: leaf chlorophyll content 49 days after inoculation; TB: total biomass (g); −N: without nitrogen; +N: with nitrogen; −N+ISRA400: without nitrogen and inoculated with *B. vignae* ISRA400; SD: standard deviation; F: value of genotype effect: Pr: *p*-value associated with the F-value of the fixed-effect genotype; *h*²: broad-sense heritability. *, **, and *** indicate values significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001, respectively.

### 3.2. QTL Mapping

To search for significant line × trait associations in each environment, the CSSLs were compared to Fleur11. We considered a significant line × trait association to be specifically linked to biological nitrogen fixation when it was not observed in the control treatments.

#### 3.2.1. Chlorophyll Content

In the 2017 environment, the Dunnett's test showed that six lines had a significantly lower leaf chlorophyll content compared to Fleur11 in the inoculated (−N+ISRA400) treatment, while in the negative (−N) and positive (+N) control treatments, no significant difference was observed (Supplementary Figure S2). The relative differences between these lines and Fleur11 ranged from −24.61 to −65.29% (Supplementary Table S1). When the experiment was repeated in 2018, a total of eight lines had a significantly lower leaf chlorophyll content compared to Fleur11 in the inoculated treatment (Supplementary Figure S3). The relative differences between these lines and Fleur11 ranged from −15.17 to −73.49% (Supplementary Table S2). One line (12CS_010) had a significantly higher chlorophyll content compared to Fleur11 in the +N treatment. All associations that were significant only under inoculated conditions were summarized and used to identify the
QTLs in each environment. In the 2017 environment, eight QTLs were linked to the leaf chlorophyll content (Table 2). One QTL was unambiguously located on each of the linkage groups (LGs) A02, A03, A04, A08, and B02 (Figure 1). In some particular cases, we were not able to unambiguously associate a QTL to one unique wild chromosome segment; therefore, the QTL was considered a putative QTL and was associated with all the fragments present in that line. This was, for example, the case for line 12CS_060, which had overlapping wild segments with several lines on the LGs A07, B06, and B08 (Figure 1). In 2018, eleven genomic regions distributed on nine LGs were identified (Table 2). One QTL was clearly mapped on each of the LGs A02, A03, A07, A08, and B02 (Figure 1). Six QTLs on the LGs A03, A04, A08, A09, B01, and B10 were considered putative QTLs. Finally, four QTLs located on the LGs A02, A03, A08, and B02 were consistently mapped in both environments.

### 3.2.2. Total Biomass (TB)

The comparison of the total biomass between the CSSLs and Fleur11 in 2017 showed that six lines had a significantly lower biomass than Fleur11 in the -N treatment, and 26 lines had a significantly lower biomass in the inoculated treatment (Supplementary Figure S4). However, among the lines found in the inoculated treatment, three lines (12CS_004, 12CS_042, and 12CS_098) were also found in the negative control treatment. Finally, a total of 23 lines were identified only in the inoculated treatment, with the relative differences ranging from -39.73 to -71.51% (Supplementary Table S1). In 2018, significant differences were also found between several lines and Fleur11: twenty-three lines in the -N treatment, one line in the +N treatment, and nine lines in the inoculated treatment (Supplementary Figure S5). Among the nine lines detected in the inoculated treatment, four (12CS_004, 12CS_034, 12CS_118, and 12CS_081) were also identified in the -N treatment. Finally, five lines were identified only in the inoculated treatment, with relative differences ranging from -55.10 to -60.58% (Supplementary Table S2). Regarding the mapping of QTLs in the 2017 environment, the line × trait associations allowed the identification of 20 QTLs located on fifteen linkage groups (Table 2). Three QTLs were mapped on LG A01, two QTLs on LG B07, and one QTL on each of the LGs A02, A05, A08, B02, B02, B04, B06, B08, B10, and B11 (Figure 1). Six QTLs were considered putative QTLs: two QTLs on LG A09 and one QTL on each of the LGs A04, A06, B03, and B10. In the 2018 environment, six QTLs located on five LGs were identified (Table 2). One QTL was unambiguously mapped on each of the LGs A02, A08, and B02 (Figure 1), whereas three QTLs located on the LGs A03, A04, and A08 were considered putative QTLs. Finally, three QTLs on the LGs A02, A08, and B02 were consistently identified in both environments.

In summary, 15 QTLs were identified for the leaf chlorophyll content and 25 QTLs were identified for the total biomass in 2017 and 2018 (Supplementary Table S3). Among these QTLs, four QTLs for the leaf chlorophyll content were stable (i.e., mapped in 2017 and 2018) and unambiguously mapped on the LGs A02, A03, A08, and B02, and three QTLs for the total biomass were stable and unambiguously mapped on the LGs A02, A08, and B02.

### 3.3. *Validation of QTLs in Controlled Conditions*

Five CSSLs that showed contrasting phenotypes compared to Fleur11 were selected to validate the QTLs. Among these lines, four (12CS_004, 12CS_044, 12CS_051, and 12CS_084) contained less chlorophyll and produced less biomass in the inoculated treatment, and one line (12CS_114) showed higher values for these traits, although these differences were not statistically significant (Supplementary Figures S2–S5).

#### Detailed Analysis of the Response to Inoculation

Differences in the plant growth between Fleur11 and the CSSLs were observed thirty-five days post-inoculation. For the leaf chlorophyll content, no difference was found between Fleur11 and the CSSLs under the -N treatment, while in the inoculated treatment, we observed a significant decrease in the chlorophyll content (~53–75%) for the lines 12CS_004, 12CS_044, 12CS_051, and 12CS_084 in the first and second experiments (Figure 2).
<table><caption>**Table 2.** Summary of the line × trait significant associations for the inoculated treatment in the 2017 and 2018 environments.</caption><thead><tr><th rowspan="2">Trait</th><th rowspan="2">CSSL</th><th rowspan="2">LG</th><th colspan="2">2017</th><th colspan="2">2018</th></tr><tr><th>QTL Effect (%)</th><th>Conf. Int (cM)</th><th>QTL Effect (%)</th><th>Conf. Int (cM)</th></tr></thead><tbody><tr><td rowspan="12">SPAD</td><td>12CS_051</td><td>A02</td><td>-65.29 ***</td><td>43.85-55.35</td><td>-43.26 ***</td><td>43.85-55.35</td></tr><tr><td>12CS_004</td><td>A03</td><td>-29.23 ***</td><td>0-16.6</td><td>-64.77 ***</td><td>0-16.6</td></tr><tr><td>12CS_103</td><td>A03</td><td>ns</td><td>0-48.975</td><td>-26.76 ***</td><td>57.2-74.2</td></tr><tr><td>12CS_027</td><td>A04</td><td>-24.61 *</td><td>0-48.975</td><td>ns</td><td></td></tr><tr><td>12CS_108</td><td>A04</td><td>ns</td><td></td><td>-22.79 ***</td><td>48.975-90.975</td></tr><tr><td>12CS_060</td><td>A07</td><td>-41.77 ***</td><td>12.1-28.525</td><td>ns</td><td></td></tr><tr><td>12CS_062</td><td>A07</td><td>ns</td><td></td><td>-26.32 ***</td><td>28.525-60.0</td></tr><tr><td>12CS_084</td><td>A08</td><td>-33.79 ***</td><td>75.525-85.375</td><td>-61.67 ***</td><td>75.525-85.375</td></tr><tr><td>12CS_103</td><td>A08</td><td>ns</td><td></td><td>-26.76 ***</td><td>46.15-75.525</td></tr><tr><td>12CS_028</td><td>A09</td><td>ns</td><td></td><td>-15.17 **</td><td>7.125-49.375</td></tr><tr><td>12CS_108</td><td>B01</td><td>ns</td><td></td><td>-22.79 ***</td><td>66.725-87.525</td></tr><tr><td>12CS_044</td><td>B02</td><td>-47.89 ***</td><td>24.55-43.25</td><td>-73.49 ***</td><td>24.55-43.25</td></tr><tr><td>12CS_060</td><td>B06</td><td>-41.77 ***</td><td>16.275-50.95</td><td>ns</td><td></td></tr><tr><td>12CS_060</td><td>B08</td><td>-41.77 ***</td><td>5.25-21.25</td><td>ns</td><td></td></tr><tr><td>12CS_028</td><td>B10</td><td>ns</td><td></td><td>-15.17 **</td><td>48.025-55.525</td></tr><tr><td rowspan="20">TB</td><td>12CS_023</td><td>A01</td><td>-48.76 **</td><td>0-9.675</td><td>ns</td><td></td></tr><tr><td>12CS_096</td><td>A01</td><td>-61.13 ***</td><td>42.325-47.7</td><td>ns</td><td></td></tr><tr><td>12CS_071</td><td>A01</td><td>-41.65 *</td><td>42.325-47.7</td><td>ns</td><td></td></tr><tr><td>12CS_016</td><td>A01</td><td>-58.32 **</td><td>111.125-119.975</td><td>ns</td><td></td></tr><tr><td>12CS_051</td><td>A02</td><td>-71.51 ***</td><td>55.35-63.125</td><td>-56.43 *</td><td>43.85-55.35</td></tr><tr><td>12CS_052</td><td>A02</td><td>-62.55 ***</td><td>55.35-63.125</td><td>ns</td><td></td></tr><tr><td>12CS_103</td><td>A03</td><td>ns</td><td></td><td>-60.58 **</td><td>57.2-74.2</td></tr><tr><td>12CS_027</td><td>A04</td><td>-47.38 **</td><td>90.975-113.15</td><td>ns</td><td></td></tr><tr><td>12CS_090</td><td>A04</td><td>-45.74 *</td><td>90.975-113.15</td><td>ns</td><td></td></tr><tr><td>12CS_098</td><td>A04</td><td>ns</td><td></td><td>-55.10 *</td><td>57.625-90.975</td></tr><tr><td>12CS_022</td><td>A05</td><td>-47.87 *</td><td>0-46.8</td><td>ns</td><td></td></tr><tr><td>12CS_118</td><td>A06</td><td>-52.02 *</td><td>29.825-43.925</td><td>ns</td><td></td></tr><tr><td>12CS_084</td><td>A08</td><td>-53.28 **</td><td>75.525-85.375</td><td>-58.35 *</td><td>75.525-85.375</td></tr><tr><td>12CS_103</td><td>A08</td><td>ns</td><td></td><td>-60.58 **</td><td>46.15-75.525</td></tr><tr><td>12CS_028</td><td>A09</td><td>-49.2 *</td><td>7.125-49.375</td><td>ns</td><td></td></tr><tr><td>12CS_118</td><td>A09</td><td>-52.02 *</td><td>72.1-96.1</td><td>ns</td><td></td></tr><tr><td>12CS_044</td><td>B02</td><td>-66.09 ***</td><td>24.55-43.25</td><td>-55.36 **</td><td>24.55-43.25</td></tr><tr><td>12CS_118</td><td>B03</td><td>-52.02 *</td><td>1.725-22.575</td><td>ns</td><td></td></tr><tr><td>12CS_059</td><td>B04</td><td>-49.12 **</td><td>40.875-57.225</td><td>ns</td><td></td></tr><tr><td>12CS_076</td><td>B06</td><td>-52.02 **</td><td>31.3-50.95</td><td>ns</td><td></td></tr><tr><td>12CS_060</td><td>B06</td><td>-69.63 ***</td><td>31.3-50.95</td><td>ns</td><td></td></tr><tr><td>12CS_047</td><td>B06</td><td>-57.13 ***</td><td>31.3-50.95</td><td>ns</td><td></td></tr><tr><td>12CS_081</td><td>B07</td><td>-60.27 **</td><td>0-2.05</td><td>ns</td><td></td></tr><tr><td>12CS_048</td><td>B07</td><td>-41.04 *</td><td>14.5-21.275</td><td>ns</td><td></td></tr><tr><td>12CS_006</td><td>B08</td><td>-50.57 *</td><td>26.95-29.8</td><td>ns</td><td></td></tr><tr><td>12CS_063</td><td>B10</td><td>-44.63 **</td><td>0-8.5</td><td>ns</td><td></td></tr><tr><td>12CS_064</td><td>B10</td><td>-47.1 **</td><td>0-8.5</td><td>ns</td><td></td></tr><tr><td>12CS_028</td><td>B10</td><td>-49.2 *</td><td>48.025-55.525</td><td>ns</td><td></td></tr><tr><td>12CS_039</td><td>B11</td><td>-39.73 *</td><td>0-15.1</td><td>ns</td><td></td></tr></tbody></table>

The QTL effect values for SPAD (leaf chlorophyll content) and TB (total biomass) represent the relative difference between the CSSLs and Fleur11 in each environment (2017 or 2018). LG indicates the linkage group, which carries the wild segment responsible for the phenotype. Conf. Int represents the confidence interval of the QTL on the chromosome. *, **, and *** indicate values significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001, respectively, using Dunnett's multiple comparisons test; ns: not significant.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Position of QTLs related to BNF in the genetic map of peanuts. Each column indicates a linkage group named from A01 to A10 (A genome) and B01 to B11 (B genome). Each row indicates a CSSL. The white background represents the genetic background of the recurrent parent Fleur11. The black areas represent the wild target chromosome segments, while the grey areas represent the wild supernumerary chromosome segments. CSSLs showing significant differences compared to Fleur11 are highlighted in orange. The green circles on the left of the figure represent the relative effect on the QTL. Their sizes are proportional to the QTL's relative effect. SPAD and TB correspond to the leaf chlorophyll content and the total biomass, respectively, evaluated in 2017 and 2018. All the QTLs identified had a negative effect on the measured traits. QTL locations are indicated by rectangles located at the bottom of the figure. Red = QTL with known location, orange = putative QTL.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Variation in leaf chlorophyll content 35 days after inoculation under controlled conditions. (a,b): First and second experiments, respectively. Without nitrogen and without nitrogen + ISRA400 indicate uninoculated plants (negative control) and plants inoculated with *B. vignae* ISRA400, respectively. *** indicate values significant at *p* < 0.001, respectively, using Dunnett's multiple comparisons test.</figcaption></figure>

Regarding the total biomass, except for line 12CS_004 in the first experiment, no difference was found between Fleur11 and the CSSLs under the –N treatment (Figure 3). In the inoculated treatment, we found that 12CS_114 produced more biomass compared to Fleur11, but this difference was statistically significant in only one experiment. Three lines (12CS_004, 12CS_044, and 12CS_051) exhibited a significantly lower biomass (23–43%) than Fleur11 in both experiments (Figure 3). Although the biomass decrease was less pronounced under controlled conditions, these results globally confirmed the observations made under the shade house, providing additional evidence that the corresponding QTLs are possibly involved in the symbiotic interaction with *B. vignae* ISRA400.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Variation in total biomass 35 days after inoculation under controlled conditions. (a,b): First and second experiments, respectively. Without nitrogen and without nitrogen + ISRA400 indicate uninoculated plants (negative control) and plants inoculated with *B. vignae* ISRA400, respectively. \* and \*\* indicate values significant at *p* < 0.05 and *p* < 0.01, respectively, using Dunnett's multiple comparisons test.</figcaption></figure>

To further investigate the phenotypes of these five lines, we evaluated the traits directly related to nodulation as indicators of symbiotic performance: the nitrogenase activity using an acetylene reduction assay, the number of nodules, the nodule dry weight, and the ratio between the nodule dry weight and the nodule number (Figure 4). The correlation matrix showed positive and significant Pearson coefficients between the nodulation traits, the leaf chlorophyll content, and the total biomass (Table 3), suggesting a strong relation between nodule development and general plant growth.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Significant line × trait of nodulation traits 35 days after inoculation between 5 CSSLs and the recurrent parent Fleur11. (a,b): First and second experiments, respectively. *, **, and *** indicate values significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001, respectively, using Dunnett's multiple comparisons test.</figcaption></figure>
<table><caption>**Table 3.** Relationships between traits related to BNF.</caption><thead><tr><th>Trait 1</th><th>Trait 2</th><th>Corr. Coef.</th></tr></thead><tbody><tr><td rowspan="4">Leaf chlorophyll content</td><td>Nitrogenase activity</td><td>0.92 ***</td></tr><tr><td>Nodule dry weight</td><td>0.89 ***</td></tr><tr><td>Nodule number</td><td>0.81 **</td></tr><tr><td>Nodule dry weight and number ratio</td><td>0.72 **</td></tr><tr><td rowspan="4">Total biomass</td><td>Nitrogenase activity</td><td>0.84 ***</td></tr><tr><td>Nodule dry weight</td><td>0.80 **</td></tr><tr><td>Nodule number</td><td>0.74 **</td></tr><tr><td>Nodule dry weight and number ratio</td><td>0.62 *</td></tr></tbody></table>

Corr. coef.: Pearson correlation coefficients obtained for 5 CSSLs and the recurrent parent Fleur11 using the average value per genotype. *, **, and *** indicate values significant at *p* < 0.05, *p* < 0.01, and *p* < 0.001.

Two independent experiments performed under controlled conditions showed that the line 12CS_004 exhibited a significant decrease in all the nodulation traits compared to Fleur11 (Figure 4). Similar results were also observed for the line 12CS_044, except for the nodule number. The line 12CS_084 produced as many nodules as Fleur11, but with a significant decrease in the nodule dry weight and nitrogenase activity. The line 12CS_051 showed a significantly lower nitrogenase activity compared to Fleur11 in both experiments. Moreover, in one experiment, we found a significant decrease in the nodule dry weight for line 12CS_051, and in the ratio between the nodule dry weight and the nodule number for line 12CS_084. No significant difference was found between Fleur11 and line 12CS_114 for all the nodulation traits.

During legume-rhizobia interactions, the symbiotic efficiency may depend on the nodule colonization by the rhizobia, and also on the development and functioning of the nodule. To better investigate the effect of QTLs, we examined histological sections of nodules at 35 days post-inoculation. In Fleur11, we observed spherical and determinate nodules of the aeschynomenoid type that exhibited central infected tissue corresponding to the fixation zone, surrounded by uninfected cells corresponding to the inner cortex or nodule parenchyma, the vascular bundles, and the outer cortex (Figure 5a). Similar anatomical features were observed for all CSSLs, including lines 12CS_004 and 12CS_044. However, the nodules of those lines were smaller compared to the ones formed by Fleur11 (Figure 5b-f). To further investigate the functional characteristics of these nodules, the bacteria were stained with a live/dead staining protocol and examined using confocal microscopy. For Fleur11, we observed numerous infected cells in the fixation zone. Those cells contained densely packed green-labelled spherical bacteroids, a landmark of functional nodules (Figure 5g,m) as reported by Sinharoy and DasGupta [45]. Similar features were observed for line 12CS_114, which was able to fix high amounts of nitrogen (Figure 5l,r). Lines with lower nitrogen fixation levels also exhibited numerous infected cells containing spherical bacteroids, but these were generally less dense and showed a higher proportion of red-labelled dead bacteroids compared to Fleur11 (Figure 5j,o-q). In line 12CS_051, we were able to detect zones where the nodule tissue collapsed, corresponding to a spontaneous degeneration of nodule tissue associated with defense responses as described by Gully et al. [46] and Songwattana et al. [47]. These observations indicate that the wild alleles at the QTLs can reduce the nodule number and/or size and also affect the maintenance of functional bacteroids in infected cells situated in the central tissue.

We screened the sequences of genomic regions corresponding to the QTLs using PeanutBase "https://peanutbase.org/gbrowse_peanut1.0" (accessed on 30 January 2023)" and Orthofinder "https://davidemms.github.io/" (accessed on 23 April 2021), and found 17 orthologs of the nodulation genes previously described in model legumes (Supplementary Table S4).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Longitudinal sections of the nodules of 5 CSSLs and the recurrent parent Fleur11 35 days after inoculation with *B. vignae* ISRA400. Differential interference contrast microscopic images (**a–f**) showing nodule structure and confocal microscopic images (**g–r**) showing live (green) bacteria, dead (red) bacteria, and plant cell walls (blue). Nodule sections of Fleur11 (**a,g,m**); line 12CS_004 (**b,h,n**); line 12CS_044 (**c,i,o**); line 12CS_051 (**d,j,p**); line 12CS_084 (**e,k,q**); and line 12CS_114 (**f,l,r**). FZ: fixation zone; NP: nodule parenchyma or inner cortex; VB: vascular bundles; NC: nucleus; NT: collapsed necrotic tissue; LB: live bacteroids; DB: dead bacteroids. Scale bar = 50 µm.</figcaption></figure>
**4. Discussion**

Biological nitrogen fixation (BNF) is a valuable component of sustainable agriculture. Understanding the genetics of the interaction between plants and rhizobia can facilitate the enhancement of BNF in grain legumes. In this study, 83 chromosome segment substitution lines (CSSLs) were used to map the QTLs of the leaf chlorophyll content and the total biomass under three treatments (−N, +N, and −N+ISRA400), and to explore the contribution of wild species to the improvement of BNF in cultivated peanuts.

**4.1. *Variability in Traits Related to BNF and QTL Mapping***

Genetic and phenotypic variations were found for the chlorophyll content and/or the total biomass in the CSSL population (Table 1). We hypothesize that these variations could be linked to the variability in nitrogen use efficiency (available or residual N uptake and transport) in the control treatments (−N and +N), and the nitrogen fixation efficiency in the inoculated treatment (−N + B. *vignae* ISRA400). Sahrawat et al. [48] showed that a non-nodulating peanut line is less efficient at acquiring N in comparison to nodulating genotypes, suggesting a relationship between nodulation and N-uptake abilities. Regarding the phenotyping for BNF, the leaf chlorophyll content has been successfully used as a surrogate trait to evaluate BNF in grain legumes, including peanuts [23,26,49–51].

We observed that some lines, such as 12CS_055 and 12CS_106, presented a total biomass that was significantly lower compared to Fleur11 in the -N treatment; however, in the treatment -N+ISRA400, these lines exhibited a total biomass similar to Fleur11 (Supplementary Tables S1 and S2). These findings could be explained by a strong sensitivity to nitrogen deficiency or by a positive effect of the wild fragments on the symbiotic performance. Thereby, for QTL detection, we assumed that significant associations were specifically linked to BNF when the phenotypic variation(s) they induced were not observed in the control treatments. Dunnett's test allowed for the identification of a total of twenty-four and nine lines that were significantly different from Fleur11 for the -N+ISRA400 treatment in the first and second shade-house experiments, respectively (Table 2). The variability in the response to inoculation observed between the two experiments was possibly due to temperature or light intensity differences in the shade-house trials (warmer days in 2017 compared to 2018). Several environmental factors, including the temperature, N concentration in the root zone, soil water content, and plant nutrient status (C and N), can affect biological nitrogen fixation in legumes [52]. Despite the variation in the experimental conditions between the 2017 and 2018 environments, four lines (12CS_051, 12CS_004, 12CS_084, and 12CS_044) had similar features for the chlorophyll content and/or the total biomass (Table 2).

In the 2017 and 2018 environments, the genotype × phenotype associations allowed for the identification of 15 and 25 QTLs for the leaf chlorophyll content and the total biomass, respectively (Supplementary Table S3). The QTLs were consistently mapped on the LGs A02, A03, A08, and B02 for the leaf chlorophyll content and/or the total biomass (Figure 1, Table 2). A co-localization for the two traits was found on the LGs A02, A08, and B02, in accordance with previous studies showing that a N deficiency reduces leaf photosynthesis and plant biomass accumulation in legumes [53] and other plants [54]. As with other economically and agronomically important traits, BNF is a complex trait regulated by several genes with different effects. Numerous QTLs involved in biological nitrogen fixation are described in other grain legumes, such as common beans [50,55], peas [56], and soybeans [16,57–60]. The simple-sequence repeat marker genotyping of two pairs of recombinant inbred lines and their parents allowed for the identification of 20 chromosome regions associated with the loss of the nodulation ability in peanuts [26]. Recently, map-based cloning and QTL-seq approaches using F₂ peanut populations have allowed for the identification of syntenic regions, including a pair of homoeologous genes coding for GRAS transcription factors involved in nodulation and nodulation signaling pathway 2 (NSP2) in the chromosomes A08 and B07 [61]. Due to the differences in the population and marker types used, it is difficult to precisely match our results with those
of previous genetic mapping experiments on peanut nodulation. Interestingly, we found that one QTL on chromosome A03 was linked to the total biomass in the inoculated treatment and in the negative control. This QTL is carried by line 12CS_004 (Supplementary Tables S1 and S2). Our findings suggest that two independent genes located in the same region are involved, for instance, in soil N uptake/transport and BNF, respectively, or that a single gene is involved in either downstream events such as nitrogen transport, metabolism, or remobilization or in low-nitrogen sensing. Examples of such genes have been reported in the literature. It is, for example, the case for the nitrate reductase gene, which plays a role in both soil N uptake [62] and in early nodulation steps [63]. The interactions among N fixation, nitrate reduction, and ammonium assimilation in legumes are well reviewed [64].

## 4.2. QTL Action on Nodulation Traits

The experiments under controlled conditions allowed for the confirmation of the shade-house results (Supplementary Figures S2 and S3) and pointed out a positive correlation between the surrogate traits, chlorophyll content and biomass, and the traits directly linked to nodulation (Table 3). We found that plants harboring wild alleles at the QTLs showed a reduced nodule number, a reduced nodule dry weight, and/or an alteration of nodule functioning, but all nodules, even small ones, exhibited similar anatomical features as those observed in the cultivated variety Fleur11. Moreover, all the genotypes exhibited determinate nodules of the aeschynomenoid type, as previously described in peanuts [65–67] and *Aeschynomene indica* [68]. Interestingly, we showed that these nodules contained numerous infected cells where we detected circular-shaped, well-differentiated bacteroids, indicating that the processes leading to bacterial recognition, nodule organogenesis, intracellular infection, and cell differentiation into bacteroids are functional. However, bacteroids were generally less dense in plant cells and often labelled in red, indicating the presence of dead bacteroids. It has been reported that in the legume-rhizobia symbiosis, the plant can sanction the rhizobia that fail to fix nitrogen inside the nodules by reducing the resource allocation (carbon), leading to decreased rhizobial fitness [69–71]. Furthermore, in line 12CS_051, we observed regions with collapsed necrotic tissue, which is often associated with incompatibility and can be assimilated to early nodule senescence triggered by plant defense mechanisms that target the invading bacteria. Necrotic zones and the presence of dead bacteria are associated with symbiotic incompatibilities and are characterized by low nitrogen fixation levels in *A. indica* [72] and in *Vigna unguiculata* [47]. These symbiotic incompatibilities were associated with mutations in the genes encoding the type III secretion system, which is a common and well-described mechanism used not only by bacterial pathogens, but also by bacterial symbionts to suppress plant immune functions [47]. Gully et al. [46] showed that necrotic zones appear in ineffective nodules of *Aeschynomene*, and they are associated with failed rhizobia differentiation.

Some QTL actions described here are close to those previously observed in functional studies performed using RNA interference to downregulate the orthologs of well-known nodulation genes in peanuts. Sinharoy and DasGupta [45] showed that the downregulation of Ca²⁺/calmodulin-dependent protein kinase (CCaMK) reduces nitrogenase activity by ~90% and affects rhizobial dissemination during nodule organogenesis. Likewise, the silencing of cytokinin receptor histidine kinase1 (HK1) led to a drastic reduction in the nodule number, nodule size, and nitrogenase activity [73]. In our study, we found 17 orthologs of nodulation genes located in the QTL regions, but at this stage, it is difficult to identify which of these genes are responsible for the phenotype(s) due to the large size of the wild introgressed chromosome segments in this CSSL population. A transcriptomic analysis and fine mapping of the major QTLs are underway to identify the underlying gene(s).

## 4.3. Wild Species, Domestication, and BNF

In our study, at all the QTLs, the wild alleles had a depressive effect on the plant growth and nodulation traits, indicating that when the efficient Senegalese strain ISRA400 is used as a symbiotic partner, *A. ipaensis* and *A. duranensis* are not valuable sources of
new alleles for increasing BNF. These results suggest either a potential positive impact of domestication and/or the breeding process in the improvement of peanut BNF or a disruption of the symbiotic mechanisms between the wild species and a foreign rhizobia strain as the result of peanut introduction and adaptation in Senegal. Bouznif et al. [74] suggested that cultivated peanuts acquired their current rhizobia by adapting to novel symbionts encountered within the area of introduction. Wild relative species of peanuts originating in South America exhibited a better nitrogen fixation ability when inoculated with native rhizobia strains from Argentina compared to the Senegalese strain *B. vignae* ISRA400 [75]. BNF is possibly the result of a coevolution and/or coadaptation process between symbiotic partners located in the same geographical area. Root nodule bacteria associated with cultivated chickpeas are phylogenetically different from those derived from the nodules of their wild relatives, indicating a coadaptation of rhizobia with legume domestication [76]. In common beans, plants were generally compatible with rhizobia isolated from their domestication centers [77], and specific polymorphisms were found both in bacterial and plant symbiotic genes [78]. However, no evidence for an adaptation to native rhizobia was found in *Medicago lupulina* [79]. In soybeans, a recombinant inbred line (RIL) population derived from a cross between a cultivated and wild germplasm showed that an improved BNF capacity could be associated with cultivated varieties [16]. Recently, the phenotyping of nodulation traits on 310 soybean accessions (10 wild soybeans, 71 landraces, and 229 improved cultivars) and 207 CSSLs allowed for the identification of the SNPs selected during the domestication and located in the regulated regions of the *GmCRP* gene encoding a rhizobial type III effector [60]. We hypothesize that the introgressions of the wild alleles of *A. duranensis* and *A. ipaensis* at the QTL region in the cultivated background of Fleur11 has partially disrupted the symbiotic association with ISRA400 in the absence of coevolution or coadaptation between the Senegalese strain and the wild species. To provide more insights into this assertion, it would be interesting to explore the response of the CSSL population when inoculated with strains coming from South America, where peanuts are native.

**5. Conclusions**

In this study, an interspecific population of 83 chromosome substitution segment lines was used to map genomic regions controlling biological nitrogen fixation traits in peanuts. We found four QTLs that were consistently associated with peanut nodulation. At these QTLs, the wild alleles had a negative effect on biological nitrogen traits. Then, screening these genomic regions allowed for the identification of 17 orthologs of the nodulation genes. However, more studies are needed to identify the underlying gene(s). Our results provide additional insights into the genetic aspects of peanut nodulation and will contribute to a better knowledge of the potential impact of domestication on the enhancement of the nitrogen fixation capacity of peanuts.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14040797/s1, Figure S1: Kinetics of average leaf chlorophyll content in the CSSL population during 49 days after inoculation in shade-house experiments; Figure S2: Significant line × leaf chlorophyll content associations observed between Fleur11 and CSSLs in the 2017 shade-house experiment; Figure S3: Significant line × leaf chlorophyll content associations observed between Fleur11 and CSSLs in the 2018 shade-house experiment; Figure S4: Significant line × total biomass associations observed between Fleur11 and CSSLs in the 2017 shade-house experiment; Figure S5: Significant line × total biomass associations observed between Fleur11 and CSSLs in the 2018 shade-house experiment; Table S1: Summary of significant line × trait associations with their QTL effects identified for all treatments in the 2017 environment; Table S2: Summary of significant line × trait associations with their QTL effects identified for all treatments in the 2018 environment; Table S3: Characteristics of the QTLs associated with leaf chlorophyll content and total biomass in the inoculated treatment; Table S4: Orthologs of nodulation genes located on the regions containing the QTLs.
**Author Contributions:** Conceptualization, D.T.N., S.F., S.S. and D.F.; methodology, D.T.N., A.Z.Z., V.H., D.G., A.S. and H.-A.T.; writing—original draft preparation, D.T.N.; writing—review and editing, D.T.N., J.R.N., D.G., A.Z.Z., V.H., J.-F.R., S.F., S.S. and D.F.; supervision, S.S., S.F. and D.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented and/or analyzed during the current study are available from the corresponding author upon reasonable request.

**Acknowledgments:** The authors thank the German Academic Exchange Service (DAAD), the Avril Foundation, and the Research Allowance for a Southern Thesis (ARTS) for providing a PhD scholarship.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Udvardi, M.; Poole, P.S. Transport and Metabolism in Legume-Rhizobia Symbioses. *Annu. Rev. Plant Biol.* **2013**, *64*, 781–805.
[CrossRef]

2. Herridge, D.F.; Peoples, M.B.; Boddey, R.M. Global Inputs of Biological Nitrogen Fixation in Agricultural Systems. *Plant Soil* **2008**, *311*, 1–18.
[CrossRef]

3. Brooker, R.W.; Bennett, A.E.; Cong, W.-F.; Daniell, T.J.; George, T.S.; Hallett, P.D.; Hawes, C.; Iannetta, P.P.M.; Jones, H.G.; Karley, A.J.; et al. Improving Intercropping: A Synthesis of Research in Agronomy, Plant Physiology and Ecology. *New Phytol.* **2015**, *206*, 107–117.
[CrossRef] [PubMed]

4. Gepts, P. Crop Domestication as a Long-Term Selection Experiment. *Plant Breed. Rev.* **2004**, *24*, 1–44.

5. Hufford, M.B.; Berny Mier y Teran, J.C.; Gepts, P. Crop Biodiversity: An Unfinished Magnum Opus of Nature. *Annu. Rev. Plant Biol.* **2019**, *70*, 727–751.
[CrossRef] [PubMed]

6. Muñoz, N.; Liu, A.; Kan, L.; Li, M.-W.; Lam, H.-M. Potential Uses of Wild Germplasms of Grain Legumes for Crop Improvement. *Int. J. Mol. Sci.* **2017**, *18*, 328.
[CrossRef] [PubMed]

7. Porter, S.S.; Sachs, J.L. Agriculture and the Disruption of Plant–Microbial Symbiosis. *Trends Ecol. Evol.* **2020**, *35*, 426–439.
[CrossRef] [PubMed]

8. Jaramillo, J.P.; Mendes, R.; Raaijmakers, J. Impact of Plant Domestication on Rhizosphere Microbiome Assembly and Functions. *Plant Mol. Biol.* **2015**, *90*, 635–644.
[CrossRef]

9. Mutch, L.A.; Young, J.P.W. Diversity and Specificity of Rhizobium Leguminosarum Biovar Viciae on Wild and Cultivated Legumes. *Mol. Ecol.* **2004**, *13*, 2435–2444.
[CrossRef]

10. Kim, D.H.; Kaashyap, M.; Rathore, A.; Das, R.R.; Parupalli, S.; Upadhyaya, H.D.; Gopalakrishnan, S.; Gaur, P.M.; Singh, S.; Kaur, J.; et al. Phylogenetic Diversity of Mesorhizobium in Chickpea. *J. Biosci.* **2014**, *39*, 513–517.
[CrossRef]

11. Chang, C.; Chen, W.; Luo, S.; Ma, L.; Li, X.; Tian, C. Rhizosphere Microbiota Assemblage Associated with Wild and Cultivated Soybeans Grown in Three Types of Soil Suspensions. *Arch. Agron. Soil Sci.* **2019**, *65*, 74–87.
[CrossRef]

12. Kiers, E.T.; Hutton, M.G.; Denison, R.F. Human Selection and the Relaxation of Legume Defences against Ineffective Rhizobia. *Proc. R. Soc. B Biol. Sci.* **2007**, *274*, 3119–3126.
[CrossRef]

13. Provorov, N.A.; Tikhonovich, I.A. Genetic Resources for Improving Nitrogen Fixation in Legume-Rhizobia Symbiosis. *Genet. Resour. Crop Evol.* **2003**, *50*, 89–99.
[CrossRef]

14. Hufford, M.B.; Lubinksy, P.; Pyhäjärvi, T.; Devengenzo, M.T.; Ellstrand, N.C.; Ross-Ibarra, J. The Genomic Signature of Crop-Wild Introgression in Maize. *PLoS Genet.* **2013**, *9*, e1003477.
[CrossRef]

15. Liu, J.; Yu, X.; Qin, Q.; Dinkins, R.D.; Zhu, H. The Impacts of Domestication and Breeding on Nitrogen Fixation Symbiosis in Legumes. *Front. Genet.* **2020**, *11*, 973.
[CrossRef]

16. Muñoz, N.; Qi, X.; Li, M.-W.; Xie, M.; Gao, Y.; Cheung, M.-Y.; Wong, F.-L.; Lam, H.-M. Improvement in Nitrogen Fixation Capacity Could Be Part of the Domestication Process in Soybean. *Heredity* **2016**, *117*, 84–93.
[CrossRef]

17. Jaiswal, S.K.; Msimbira, L.A.; Dakora, F.D. Phylogenetically Diverse Group of Native Bacterial Symbionts Isolated from Root Nodules of Groundnut (*Arachis hypogaea* L.) in South Africa. *Syst. Appl. Microbiol.* **2017**, *40*, 215–226.
[CrossRef] [PubMed]

18. Zaiya, Z.A.; Fonceka, D.; Fall, S.; Fabra, A.; Ibañez, F.; Pignoly, S.; Diouf, A.; Touré, O.; Faye, M.N.; Hocher, V.; et al. Genetic Diversity and Symbiotic Efficiency of Rhizobial Strains Isolated from Nodules of Peanut (*Arachis hypogaea* L.) in Senegal. *Agric. Ecosyst. Environ.* **2018**, *265*, 384–391.
[CrossRef]

19. Janila, P.; Nigam, S.N.; Pandey, M.K.; Nagesh, P.; Varshney, R.K. Groundnut Improvement: Use of Genetic and Genomic Tools. *Front. Plant Genet. Genom.* **2013**, *4*, 1–16.
[CrossRef]
20. Wynne, J.C.; Elkan, G.H.; Schneeweis, T.J. Increasing Nitrogen Fixation of the Groundnut by Strain and Host Selection. In Proceedings of the International Workshop on Groundnuts, Patancheru, India, 13–17 October 1980; ICRISAT Center: Hyderabad, India; pp. 95–109.

21. Arrenddell, S.; Wynne, J.C.; Elkan, G.H.; Isleib, T.G. Variation for Nitrogen Fixation Among Progenies of a Virginia × Spanish Peanut Cross1. Crop Sci. **1985**, *25*, 865–869. [CrossRef]

22. Nigam, S.N.; Dwivedi, S.L.; Nambiar, P.T.C.; Gibbons, R.W.; Dart, P.J. Combining Ability Analysis of N2-Fixation and Related Traits in Peanut. *Peanut Sci.* **1985**, *12*, 55–57. [CrossRef]

23. Pimratch, S.; Jogloy, S.; Toomsan, B.; Jaisil, P.; Sikhinarum, J.; Kesmala, T.; Patanothai, A. Evaluation of Seven Peanut Genotypes for Nitrogen Fixation and Agronomic Traits. *Songklanakarin J. Sci. Technol.* **2004**, *26*, 295–304.

24. Dutta, M.; Reddy, L.J. Further Studies on Genetics of Nonnodulation in Peanut. Crop Sci. **1988**, *28*, 60–62. [CrossRef]

25. Gallo-Meagher, M.; Dashiell, K.E.; Gorbet, D.W. Parental Effects in the Inheritance of Nonnodulation in Peanut. J. Hered. **2001**, *92*, 86–89. [CrossRef] [PubMed]

26. Peng, Z.; Tan, L.; Lopez, Y.; Maku, J.; Liu, F.; Zhou, H.; Tseng, Y.-C.; Yang, X.; Hsieh, Y.; Song, J.; et al. Morphological and Genetic Characterization of Non-Nodulating Peanut Recombinant Inbred Lines. Crop Sci. **2018**, *58*, 540–550. [CrossRef]

27. Griesmann, M.; Chang, Y.; Liu, X.; Song, Y.; Haberer, G.; Crook, M.B.; Billault-Penneteau, B.; Lauressergues, D.; Keller, J.; Imanishi, L.; et al. Phylogenomics reveals multiple losses of nitrogen-fixing root nodule symbiosis. Science **2018**, *361*, 144. [CrossRef] [PubMed]

28. Karmakar, K.; Kundu, A.; Rizvi, A.Z.; Dubois, E.; Severac, D.; Czernic, P.; Cartieaux, F.; DasGupta, M. Transcriptomic Analysis With the Progress of Symbiosis in 'Crack-Entry' Legume *Arachis hypogaea* Highlights Its Contrast With 'Infection Thread' Adapted Legumes. *MPMI* **2018**, *32*, 271–285. [CrossRef] [PubMed]

29. Roy, S.; Liu, W.; Nandety, R.S.; Crook, A.; Mysore, K.S.; Pislariu, C.I.; Frugoli, J.; Dickstein, R.; Udvardi, M.K. Celebrating 20 Years of Genetic Discoveries in Legume Nodulation and Symbiotic Nitrogen Fixation. Plant Cell **2020**, *32*, 15–41. [CrossRef]

30. Raul, B.; Bhattacharjee, O.; Ghosh, A.; Upadhyay, P.; Tembhare, K.; Singh, A.; Shaheen, T.; Ghosh, A.K.; Torres-Jerez, I.; Krom, N.; et al. Microscopic and Transcriptomic Analyses of Dalbergoid Legume Peanut Reveal a Divergent Evolution Leading to Nod-Factor-Dependent Epidermal Crack-Entry and Terminal Bacteroid Differentiation. *MPMI* **2022**, *35*, 131–145. [CrossRef]

31. Sprent, J.I.; James, E.K. Legume Evolution: Where Do Nodules and Mycorrhizas Fit In? Plant Physiol. **2007**, *144*, 575–581. [CrossRef]

32. Ibáñez, F.; Wall, L.; Fabra, A. Starting Points in Plant-Bacteria Nitrogen-Fixing Symbioses: Intercellular Invasion of the Roots. J. Exp. Bot. **2017**, *68*, 1905–1918. [CrossRef]

33. Quilbé, J.; Montiel, J.; Arrighi, J.-F.; Stougaard, J. Molecular Mechanisms of Intercellular Rhizobial Infection: Novel Findings of an Ancient Process. Front. Plant Sci. **2022**, *13*, 922982. [CrossRef]

34. Bertioli, D.J.; Cannon, S.B.; Froenicke, L.; Huang, G.; Farmer, A.D.; Cannon, E.K.S.; Liu, X.; Gao, D.; Clevenger, J.; Dash, S.; et al. The Genome Sequences of *Arachis duranensis* and *Arachis ipaensis*, the Diploid Ancestors of Cultivated Peanut. Nat. Genet. **2016**, *48*, 438–446. [CrossRef]

35. Fonceka, D.; Hodo-Abalo, T.; Rivallan, R.; Faye, I.; Sall, M.N.; Ndoye, O.; Fávero, A.P.; Bertioli, D.J.; Glaszmann, J.-C.; Courtois, B.; et al. Genetic Mapping of Wild Introgressions into Cultivated Peanut: A Way toward Enlarging the Genetic Basis of a Recent Allotetraploid. BMC Plant Biol. **2009**, *9*, 103. [CrossRef]

36. Mallikarjuna, N.; Senthilvel, S.; Hoisington, D. Development of new sources of tetraploid *Arachis* to broaden the genetic base of cultivated groundnut (*Arachis hypogaea* L.). Genet. Resour. Crop. Evol. **2010**, *58*, 889–907. [CrossRef]

37. Fonceka, D.; Tossim, H.-A.; Rivallan, R.; Vignes, H.; Faye, I.; Ndoye, O.; Moretzsohn, M.C.; Bertioli, D.J.; Glaszmann, J.-C.; Courtois, B.; et al. Fostered and left behind alleles in peanut: Interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding. BMC Plant Biol. **2012**, *12*, 26. [CrossRef] [PubMed]

38. Fonceka, D.; Tossim, H.-A.; Rivallan, R.; Vignes, H.; Lacut, E.; de Bellis, F.; Faye, I.; Ndoye, O.; Leal-Bertioli, S.C.M.; Valls, J.F.M.; et al. Construction of Chromosome Segment Substitution Lines in Peanut (*Arachis hypogaea* L.) Using a Wild Synthetic and QTL Mapping for Plant Morphology. PLoS ONE **2012**, *7*, e48642. [CrossRef] [PubMed]

39. Khera, P.; Pandey, M.K.; Mallikarjuna, N.; Sriswathi, M.; Roorkiwal, M.; Janila, P.; Sharma, S.; Shilpa, K.; Sudini, H.; Guo, B.; et al. Genetic imprints of domestication for disease resistance, oil quality, and yield component traits in groundnut (*Arachis hypogaea* L.). Mol. Genet. Genom. **2018**, *294*, 365–378. [CrossRef] [PubMed]

40. Tossim, H.-A.; Nguepjop, J.R.; Diatta, C.; Sambou, A.; Seye, M.; Sane, D.; Rami, J.-F.; Fonceka, D. Assessment of 16 Peanut (*Arachis hypogaea* L.) CSSLs Derived from an Interspecific Cross for Yield and Yield Component Traits: QTL Validation. *Agronomy* **2020**, *10*, 583. [CrossRef]

41. Vincent, J.M. A Manual for the Practical Study of the Root-Nodule Bacteria; Blackwell Scientific Publications: Hoboken, NJ, USA, 1970.

42. Kung, F.H.; Bey, C.F. Heritability Construction for Provenance and Family Selection; Department of Agriculture, Forest Service, North Central Forest Experiment Station: St. Paul, MN, USA, 1977.

43. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. *Biom. J.* **2008**, *50*, 346–363. [CrossRef]

44. Nagata, T.; Takebe, I. Cell Wall Regeneration and Cell Division in Isolated Tobacco Mesophyll Protoplasts. *Planta* **1970**, *92*, 301–308. [CrossRef] [PubMed]

45. Sinharoy, S.; DasGupta, M. RNA Interference Highlights the Role of CCaMK in Dissemination of Endosymbionts in the Aeschynomeneae Legume Arachis. *MPMI* **2009**, *22*, 1466–1475. [CrossRef] [PubMed]
46. Gully, D.; Gargani, D.; Bonaldi, K.; Grangeteau, C.; Chaintreuil, C.; Fardoux, J.; Nguyen, P.; Marchetti, R.; Nouwen, N.; Molinaro, A.; et al. A Peptidoglycan-Remodeling Enzyme Is Critical for Bacteroid Differentiation in *Bradyrhizobium* spp. During Legume Symbiosis. *MPMI* **2016**, *29*, 447–457. [CrossRef]

47. Songwattana, P.; Chaintreuil, C.; Wongdee, J.; Teulet, A.; Mbaye, M.; Piromyou, P.; Gully, D.; Fardoux, J.; Zoumman, A.M.A.; Camuel, A.; et al. Identification of Type III Effectors Modulating the Symbiotic Properties of *Bradyrhizobium vignae* Strain ORS3257 with Various Vigna Species. *Sci. Rep.* **2021**, *11*, 4874. [CrossRef]

48. Sahrawat, K.L.; Srinivasa Rao, B.; Nambiar, P.T.C. Macro- and Micronutrient Uptake by Nodulating and Non-Nodulating Peanut Lines. *Plant Soil* **1988**, *109*, 291–293. [CrossRef]

49. Vollmann, J.; Walter, H.; Sato, T.; Schweiger, P. Digital Image Analysis and Chlorophyll Metering for Phenotyping the Effects of Nodulation in Soybean. *Comput. Electron. Agric.* **2011**, *75*, 190–195. [CrossRef]

50. Ramaekers, L.; Galeano, C.H.; Garzón, N.; Vanderleyden, J.; Blair, M.W. Identifying Quantitative Trait Loci for Symbiotic Nitrogen Fixation Capacity and Related Traits in Common Bean. *Mol. Breed.* **2013**, *31*, 163–180. [CrossRef]

51. Hamawaki, R.L.; Kantartzi, S.K. Di-Nitrogen Fixation at the Early and Late Growth Stages of Soybean. *Acta Scientiarum. Agron.* **2018**, *40*, e36372. [CrossRef]

52. Liu, Y.; Wu, L.; Baddeley, J.A.; Watson, C.A. Models of biological nitrogen fixation of legumes. A review. *Agron. Sustain. Dev.* **2011**, *31*, 155–172. [CrossRef]

53. Herridge, D.F.; Robertson, M.J.; Cocks, B.; Peoples, M.B.; Holland, J.F.; Heuke, L. Low nodulation and nitrogen fixation of mungbean reduce biomass and grain yields. *Aust. J. Exp. Agric.* **2005**, *45*, 269–277. [CrossRef]

54. Nguyen, N.T.; Nakabayashi, K.; Mohapatra, P.K.; Thompson, J.; Fujita, K. Effect of Nitrogen Deficiency on Biomass Production, Photosynthesis, Carbon Partitioning, and Nitrogen Nutrition Status of *Melaleuca* and *Eucalyptus* Species. *Soil Sci. Plant Nutr.* **2003**, *49*, 99–109. [CrossRef]

55. Souza, A.A.; Boscariol, R.L.; Moon, D.H.; Camargo, L.E.A.; Tsai, S.M. Effects of Phaseolus Vulgaris QTL in Controlling Host-Bacteria Interactions under Two Levels of Nitrogen Fertilization. *Genet. Mol. Biol.* **2000**, *23*, 155–161. [CrossRef]

56. Bourion, V.; Rizvi, S.M.H.; Fournier, S.; de Larambergue, H.; Galmiche, F.; Marget, P.; Duc, G.; Burstin, J. Genetic dissection of nitrogen nutrition in pea through a QTL approach of root, nodule, and shoot variability. *Theor. Appl. Genet.* **2010**, *121*, 71–86. [CrossRef]

57. Nicolás, M.F.; Hungria, M.; Arias, C.A.A. Identification of Quantitative Trait Loci Controlling Nodulation and Shoot Mass in Progenies from Two Brazilian Soybean Cultivars. *Field Crops Res.* **2006**, *95*, 355–366. [CrossRef]

58. Santos, M.A.; Geraldi, I.O.; Garcia, A.A.F.; Bortolatto, N.; Schiavon, A.; Hungria, M. Mapping of QTLs Associated with Biological Nitrogen Fixation Traits in Soybean. *Hereditas* **2013**, *150*, 17–25. [CrossRef]

59. Ni, H.; Peng, Y.; Wang, J.; Wang, J.; Yuan, Y.; Fu, T.; Zhu, Z.; Zhang, J.; Pan, X.; Cui, Z.; et al. Mapping of Quantitative Trait Loci Underlying Nodule Traits in Soybean (*Glycine max* (L.) Merr.) and Identification of Genes Whose Expression Is Affected by the Sinorhizobium Fredii HH103 Effector Proteins NopL and NopT. *Agronomy* **2022**, *12*, 946. [CrossRef]

60. Wang, J.; Ma, C.; Ma, S.; Zheng, H.; Tian, H.; Wang, X.; Wang, Y.; Jiang, H.; Wang, J.; Zhang, Z.; et al. Genetic Variation in *GmCRP* Contributes to Nodulation in Soybean (*Glycine max* Merr.). *Crop J.* 2022; in press. [CrossRef]

61. Peng, Z.; Chen, H.; Tan, L.; Shu, H.; Varshney, R.K.; Zhou, Z.; Zhao, Z.; Luo, Z.; Chitikineni, A.; Wang, L.; et al. Natural Polymorphisms in a Pair of NSP2 Homoeologs Can Cause Loss of Nodulation in Peanut. *J. Exp. Bot.* **2021**, *72*, 1104–1118. [CrossRef] [PubMed]

62. Good, A.G.; Shrawat, A.K.; Muench, D.G. Can Less Yield More? Is Reducing Nutrient Input into the Environment Compatible with Maintaining Crop Production? *Trends Plant Sci.* **2004**, *9*, 597–605. [CrossRef]

63. Damiani, I.; Pauly, N.; Puppo, A.; Brouquisse, R.; Boscari, A. Reactive Oxygen Species and Nitric Oxide Control Early Steps of the Legume—Rhizobium Symbiotic Interaction. *Front. Plant Sci.* **2016**, *7*, 454. [CrossRef]

64. Murray, J.D.; Liu, C.-W.; Chen, Y.; Miller, A.J. Nitrogen Sensing in Legumes. *J. Exp. Bot.* **2017**, *68*, 1919–1926. [CrossRef] [PubMed]

65. Chandler, M.R. Some Observations on Infection of *Arachis hypogaea* L. by Rhizobium. *J. Exp. Bot.* **1978**, *29*, 749–755. [CrossRef]

66. Boogerd, F.C.; van Rossum, D. Nodulation of groundnut by *Bradyrhizobium*: A simple infection process by crack entry. *FEMS Microbiol. Rev.* **1997**, *21*, 5–27. [CrossRef]

67. Guinel, F.C. Getting around the Legume Nodule: I. The Structure of the Peripheral Zone in Four Nodule Types. *Botany* **2009**, *87*, 1117–1138. [CrossRef]

68. Bonaldi, K.; Gherbi, H.; Franche, C.; Bastien, G.; Fardoux, J.; Barker, D.; Giraud, E.; Cartieaux, F. The Nod Factor-Independent Symbiotic Signaling Pathway: Development of *Agrobacterium rhizogenes*-Mediated Transformation for the Legume *Aeschynomene indica*. *MPMI* **2010**, *23*, 1537–1544. [CrossRef]

69. Kiers, E.T.; Rousseau, R.A.; West, S.A.; Denison, R.F. Host Sanctions and the Legume-Rhizobium Mutualism. *Nature* **2003**, *425*, 78–81. [CrossRef]

70. Oono, R.; Denison, R.F.; Kiers, E.T. Controlling the Reproductive Fate of Rhizobia: How Universal Are Legume Sanctions? *New Phytol.* **2009**, *183*, 967–979. [CrossRef]

71. Oono, R.; Anderson, C.G.; Denison, R.F. Failure to Fix Nitrogen by Non-Reproductive Symbiotic Rhizobia Triggers Host Sanctions That Reduce Fitness of Their Reproductive Clonemates. *Proc. R. Soc. B Biol. Sci.* **2011**, *278*, 2698–2703. [CrossRef] [PubMed]
72. Teulet, A.; Busset, N.; Fardoux, J.; Gully, D.; Chaintreuil, C.; Cartieaux, F.; Jauneau, A.; Comorge, V.; Okazaki, S.; Kaneko, T.; et al. The rhizobial type III effector ErnA confers the ability to form nodules in legumes. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 21758–21768. [CrossRef] [PubMed]

73. Kundu, A.; DasGupta, M. Silencing of Putative Cytokinin Receptor Histidine Kinase1 Inhibits Both Inception and Differentiation of Root Nodules in *Arachis hypogaea*. *Mol. Plant Microbe Interact.* **2018**, *31*, 187–199. [CrossRef]

74. Bouznif, B.; Guefrachi, I.; Rodríguez de la Vega, R.C.; Hungria, M.; Mars, M.; Alunni, B.; Shykoff, J.A. Phylogeography of the *Bradyrhizobium* spp. Associated With Peanut, *Arachis hypogaea*: Fellow Travelers or New Associations? *Front. Microbiol.* **2019**, *10*, 2041. [CrossRef]

75. Zaiya, Z.A. Etude de La Composante Symbiotique Bactérienne Impliquée Dans l'adaptation de l'arachide Cultivée (*Arachis hypogaea* L.) Au Sénégal. Ph.D. Thesis, Cheikh Anta Diop de Dakar, Dakar, Senegal, 2018.

76. Greenlon, A.; Chang, P.L.; Damtew, Z.M.; Muleta, A.; Carrasquilla-Garcia, N.; Kim, D.; Nguyen, H.P.; Suryawanshi, V.; Krieg, C.P.; Yadav, S.K.; et al. Global-Level Population Genomics Reveals Differential Effects of Geography and Phylogeny on Horizontal Gene Transfer in Soil Bacteria. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 15200–15209. [CrossRef]

77. Aguilar, O.M.; Riva, O.; Peltzer, E. Analysis of *Rhizobium etli* and of its symbiosis with wild *Phaseolus vulgaris* supports coevolution in centers of host diversification. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 13548–13553. [CrossRef] [PubMed]

78. Aguilar, O.M.; Collavino, M.M.; Mancini, U. Nodulation competitiveness and diversification of symbiosis genes in common beans from the American centers of domestication. *Sci. Rep.* **2022**, *12*, 4591. [CrossRef] [PubMed]

79. Harrison, T.L.; Wood, C.W.; Borges, I.L.; Stinchcombe, J.R. No evidence for adaptation to local rhizobial mutualists in the legume *Medicago lupulina*. *Ecol. Evol.* **2017**, *7*, 4367–4376. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.