AbstractThe molecular mechanisms of symbiosis in cultivated peanut with a ‘crack entry’ infection process are largely understudied. In this study, we investigated the root transcriptional profiles of two pairs of non-nodulating (nod−) and nodulating (nod+) sister inbred peanut lines, E4/E5 and E7/E6, and their nod+ parents, F487A and PI262090 during rhizobial infection and nodule initiation by using RNA-seq technology. A total of 143, 101, 123, 215, 182, and 289 differentially expressed genes (DEGs) were identified in nod− E4, E7 and nod+ E5, E6, F487A, and PI262090 after inoculation with Bradyrhizobium sp. Different deficiencies at upstream of symbiotic signaling pathway were revealed in the two nod− genotypes. DEGs specific in nod+ genotypes included orthologs to some known symbiotic signaling pathway genes, such as NFR5, NSP2, NIN, ERN1, and many other novel and/or functionally unknown genes. Gene ontology (GO) enrichment analysis of nod+ specific DEGs revealed 54 significantly enriched GO terms, including oxidation-reduction process, metabolic process, and catalytic activity. Genes related with plant defense systems, hormone biosynthesis and response were particularly enriched. To our knowledge, this is the first report revealing symbiosis-related genes in a genome-wide manner in peanut representative of the ‘crack entry’ species.