A novel QTL <Emphasis Type="Italic">qSPP2.2</Emphasis> controlling spikelet per panicle identified from <Emphasis Type="Italic">Oryza longistaminata</Emphasis> (A. Chev. et Roehr.), mapped and transferred to <Emphasis Type="Italic">Oryza sativa</Emphasis> (L.)

Increasing the rice productivity from the current 10 to 12 tons/ha to meet the demand of estimated 8.8 billion people in 2035 is posing a major challenge. Wild relatives of rice contain some novel genes which can help in improving rice yield. Spikelet per panicle (SPP) is a valuable trait for determining yield potential in rice. In this study, a major QTL for increasing SPP has been identified, mapped, and transferred from African wild rice O. longistaminata to O. sativa (L.). The QTL was mapped on the long arm of chromosome 2 in a 167.1 kb region flanked by SSR markers RM13743 and RM13750, which are 1.0 cM apart, and is designated as qSPP2.2. The QTL explained up to 30% of phenotypic variance in different generations/seasons and showed positive additive effect of allele contributed by O. longistaminata. In addition, O. longistaminata allele in qSPP2.2 contributed to increase in grains per panicle, but decrease in the tillers per plant. The 167.1 kb region contains 23 predicted genes. Based on the functional annotation, three genes, LOC_Os02g44860, LOC_Os02g44990, and LOC_Os02g45010, were selected as putative candidates for characterization. Sequence analysis of the three genes revealed functional variations between the parental lines for LOC_Os02g44990 and a variation in 5′UTR for LOC_Os02g45010 which will help further to identify putative candidate gene(s). This is the first yield component QTL to be identified, mapped, and transferred from O. longistaminata.     

Fine mapping of QTL <Emphasis Type="Italic">qCTB10</Emphasis>-<Emphasis Type="Italic">2</Emphasis> that confers cold tolerance at the booting stage in rice

Key message

The QTL qCTB10 - 2 controlling cold tolerance at the booting stage in rice was delimited to a 132.5 kb region containing 17 candidate genes and 4 genes were cold-inducible.

Abstract

Low temperature at the booting stage is a major abiotic stress-limiting rice production. Although some QTL for cold tolerance in rice have been reported, fine mapping of those QTL effective at the booting stage is few. Here, the near-isogenic line ZL31-2, selected from a BC₇F₂ population derived from a cross between cold-tolerant variety Kunmingxiaobaigu (KMXBG) and the cold-sensitive variety Towada, was used to map a QTL on chromosome 10 for cold tolerance at the booting stage. Using BC₇F₃ and BC₇F₄ populations, we firstly confirmed qCTB10-2 and gained confidence that it could be fine mapped. QTL qCTB10-2 explained 13.9 and 15.9% of the phenotypic variances in those two generations, respectively. Using homozygous recombinants screened from larger BC₇F₄ and BC₇F₅ populations, qCTB10-2 was delimited to a 132.5 kb region between markers RM25121 and MM0568. 17 putative predicted genes were located in the region and only 5 were predicted to encode expressed proteins. Expression patterns of these five genes demonstrated that, except for constant expression of LOC_Os10g11820, LOC_Os10g11730, LOC_Os10g11770, and LOC_Os10g11810 were highly induced by cold stress in ZL31-2 compared to Towada, while LOC_Os10g11750 showed little difference. Our results provide a basis for identifying the genes underlying qCTB10-2 and indicate that markers linked to the qCTB10-2 locus can be used to improve the cold tolerance of rice at the booting stage by marker-assisted selection.

     

Tandem duplications of receptor-like kinase genes reshaped the homologous chromosome 1 regions in <Emphasis Type="Italic">Oryza sativa</Emphasis> ssp. <Emphasis Type="Italic">indica</Emphasis> and <Emphasis Type="Italic">japonica</Emphasis> and their possible ancestral genomes

Plant receptor-like kinase (Rlk) genes form a large family, each encoding a protein with a signal motif, a single transmembrane region, and a cytoplasmic kinase domain. Various gene duplications have contributed to the establishment and expansion of the family. Here, we characterized the formation and evolution of the Rlk gene family in cultivated rice and their possible progenitors. Using wheat Rlk gene sequences, we identified orthologs from the genomes of domesticated rice subspecies Oryza sativa ssp. japonica and ssp. indica and their putative progenitors O. glaberrima and O. rufipogon. The four chromosome 1 orthologous regions ranged from 103 to 281 kb comprising 181 syntenic blocks with 75 to 100% sequence identity. These regions contained 11–19 Triticum aestivum kinases (Taks) and 10–15 Lr10 receptor-like kinases (Lrks) organized in clusters and 3–12 transposable elements (TEs). Dot plot analyses showed that the 4 regions had 21–37 conserved catalytic domains, mainly in protein kinases (PKs) and tyrosine kinases (TyrKs) in coupling state. Over 50% of the sequences of glaberrima/rufipogon and japonica/indica pairs were colinear, while japonica/indica displayed a marked sequence expansion with duplicated genes and TEs. A total of 2312 single nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) were identified between japonica and indica. Duplication of the Rlk genes in O. glaberrima and O. rufipogon occurred after the grass species radiation and before the divergence of O. rufipogon from O. glaberrima; the orthologous Rlk genes from O. japonica and O. indica duplicated after O. sativa separated from O. rufipogon; paralogs, obtained through extensive duplication, happened after the separation of rice from maize. Tandem duplication was the major factor contributing to the gene copy number variation and genome size expansion.     

Fine mapping and candidate gene analysis of <Emphasis Type="Italic">qFL-chr1</Emphasis>, a fiber length QTL in cotton

Key message

A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length.

Abstract

Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC₄F₂ plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC₄F₃ plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F₂ population, supporting these genes as candidates for qFL-chr1.

     

<Emphasis Type="Italic">Trueperella pyogenes</Emphasis> isolated from a brain abscess of an adult roebuck (<Emphasis Type="Italic">Capreolus capreolus</Emphasis>)

The present study was designed to characterize phenotypically and genotypically a Trueperella pyogenes strain isolated from a brain abscess of an adult roebuck (Capreolus capreolus). The species identity could be confirmed by phenotypical investigations, by MALDI-TOF MS analysis, and by sequencing the 16S ribosomal RNA (rRNA) gene, the 16S–23S rRNA intergenic spacer region (ISR); by sequencing the target genes rpoB, gap, and tuf; and by detection of T. pyogenes chaperonin-encoding gene cpn60 with a previously developed loop-mediated isothermal amplification (LAMP) assay. The T. pyogenes strain could additionally be characterized by PCR-mediated amplification of several known and putative virulence factor-encoding genes which revealed the presence of the genes plo encoding pyolysin and nanH and nanP encoding neuraminidases; the genes fimA, fimC, and fimE encoding the fimbrial subunits FimA, FimC, and FimE; and the gene cbpA encoding collagen-binding protein CbpA. The present data give a detailed characterization of a T. pyogenes strain isolated from a brain abscess of a roebuck. However, the route of infection of the roebuck remains unclear.     

Complete genome sequence and pathogenic genes analysis of <Emphasis Type="Italic">Pectobacterium atroseptica</Emphasis> JG10-08

Pectobacterium atroseptica is known as a rod-shaped gram-negative bacterial pathogen associated with the blackleg of potato. P. atroseptica has been widely identified as the predominant agent causing tuber rot in temperate regions, a disease that leads to severe economic losses to potato industry. In this study, we provide the complete genome sequence of P. atroseptica JG10-08, which revealed that P. atroseptica strain JG10-08 carries a single 5,004,926 bp chromosome with 51.15% G+C content and harbors 4252 predicted coding genes. Phylogenetic analysis based on the genome sequences showed a close evolutionary relationship between P. atroseptica and Pectobacterium wasabiae. We discovered total 168 genes were potentially related to pathogenesis including 9 strain-specific genes encoding toxins on the genome of JG10-08. Further comparison with other species in Pectobacterium revealed a better understanding of pathogenic factors, especially secretion systems in P. atroseptica JG10-08. Collectively, the results of this research provide a solid foundation for discovering the underlying pathogenic mechanisms of P. atroseptica and offer the information to develop more effective strategies against blackleg of potatoes.     

A rice <Emphasis Type="Italic">White-stripe leaf3</Emphasis> (<Emphasis Type="Italic">wsl3</Emphasis>) mutant lacking an HD domain-containing protein affects chlorophyll biosynthesis and chloroplast development

Leaf-color mutants are ideal genetic materials for understanding the mechanism of chloroplast development and chlorophyll (Chl) biosynthesis. Here we isolated and identified a new leaf-color mutant of rice, named white-stripe leaf3 (wsl3), from a ⁶⁰Co-irradiated mutant pool. The wsl3 mutant displayed a visible white-stripe leaf in both young seedlings and flag leaves of mature plant. Chl content in homozygous wsl3 mutant was approximately 47% of that in the wild type. Besides, chloroplast development in the mutant was severely arrested. By a map-based cloning strategy, the wsl3 gene was finely confined to a 50.8 kb region on chromosome 1. Moreover, a 9-bp deletion was identified in the genomic region of LOC_Os01g01920, which encodes an HD (histidine and aspartic acid) domaincontaining protein. Genetic complementation confirmed that LOC_Os01g01920 could recover the lesion of wsl3 mutation. Real-time PCR analyses showed that the expression levels of WSL3 were the highest in young and flag leaves among various tissues, and most of the genes associated with Chl biosynthesis were significantly down-regulated in the wsl3 mutant. Meanwhile, in contrast to many nuclear gene-encoded phage-type RNA polymerase(s) (NEP) transcribed genes were up-regulated, most of plastid-encoded bacterialtype RNA polymerase (PEP) transcribed genes were downregulated. These results demonstrated that the WSL3 gene, as an HD domain-containing protein, is involved in chl biosynthesis and chloroplast development in rice.     

The causal deletions in the second exon of <Emphasis Type="Italic">An-3</Emphasis> closely associated with awn development and rice yield

Awn is one of important traits during rice domestication. To understand the development of rice awn and the roles it played in rice domestication, we preliminary mapped a major QTL An-3 for awn development using chromosome segment substitution line CSSL138 developed by introgressed genomic fragments of long-awned Guangxi common wild rice (GXCWR, Oryza rufipogon Griff.) into genetic background of short-awned indica cultivar 93–11. An-3 was then fine mapped to a 7-kb region of chromosome 8. An epidermal patterning factor-like protein gene was identified as the single candidate gene corresponding to this QTL. An-3 was showed to be an allele of RAE2 and GAD1, and negatively regulated 1000-grains weight, grain length, and length–width ratio. Comparing with the coding sequences of An-3 from CSSL138, a 2- and 4-bp frame-shift deletions in the second exon were identified in 93–11 and Nipponbare, respectively. Taken together, our results provide valuable natural variation in the alleles of An-3 between common wild rice and cultivated rice, which will be helpful in clarifying the mechanism of awn development and promoting the application of an-3 in genetic improvement of rice yield traits.     

Mapping of days to flower and seed yield in spring oilseed <Emphasis Type="Italic">Brassica napus</Emphasis> carrying genome content introgressed from <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Earliness of flowering and maturity and high seed yield are important objectives of breeding spring Brassica napus canola. Previously, we have introgressed earliness of flowering from Brassica oleracea into spring B. napus canola through interspecific crossing between these two species. In this paper, we report quantitative trait locus (QTL) mapping of days to flower and seed yield by use of publicly available markers and markers designed based on flowering time genes and a doubled haploid population, derived from crossing of the spring canola parent and an early flowering line developed from a B. napus × B. oleracea cross, tested in nine field trials for over 5 years. Five genomic regions associated with days to flower were identified on C1, C2, C3, and C6 of which the single QTL of C1 was detected in all trials; in all cases, the allele introgressed from B. oleracea reduced the number of days to flower. BLASTn search in the Brassica genomes located the physical position of the QTL markers and identified putative flowering time genes in these regions. In the case of seed yield, ten QTL from eight linkage groups were detected; however, none could be consistently detected in all trials. The QTL region of C1 associated with days to flower did not show significant association with seed yield in more than 80% of the field trials; however, in a single trial, the allele introgressed from B. oleracea exerted a negative effect on seed yield. Thus, the genomic regions and molecular markers identified in this research could potentially be used in breeding for the development of early flowering B. napus canola cultivars without affecting seed yield in a majority of the environments.     

Genetic dissection of seed weight by QTL analysis and detection of allelic variation in Indian and east European gene pool lines of <Emphasis Type="Italic">Brassica juncea</Emphasis>

Key message

Seed weight QTL identified in different populations were synthesized into consensus QTL which were shown to harbor candidate genes by in silico mapping. Allelic variation inferred would be useful in breeding B. juncea lines with high seed weight.

Abstract

Seed weight is an important yield influencing trait in oilseed Brassicas and is a multigenic trait. Among the oilseed Brassicas, Brassica juncea harbors the maximum phenotypic variation wherein thousand seed weight varies from around 2.0 g to more than 7.0 g. In this study, we have undertaken quantitative trait locus/quantitative trait loci (QTL) analysis of seed weight in B. juncea using four bi-parental doubled-haploid populations. These four populations were derived from six lines (three Indian and three east European lines) with parental phenotypic values for thousand seed weight ranging from 2.0 to 7.6 g in different environments. Multi-environment QTL analysis of the four populations identified a total of 65 QTL ranging from 10 to 25 in each population. Meta-analysis of these component QTL of the four populations identified six ‘consensus’ QTL (C-QTL) in A3, A7, A10 and B3 by merging 33 of the 65 component Tsw QTL from different bi-parental populations. Allelic diversity analysis of these six C-QTL showed that Indian lines, Pusajaikisan and Varuna, hold the most positive allele in all the six C-QTL. In silico mapping of candidate genes with the consensus QTL localized 11 genes known to influence seed weight in Arabidopsis thaliana and also showed conserved crucifer blocks harboring seed weight QTL between the A subgenomes of B. juncea and B. rapa. These findings pave the way for a better understanding of the genetics of seed weight in the oilseed crop B. juncea and reveal the scope available for improvement of seed weight through marker-assisted breeding.

     

Construction of a sequence-based bin map and mapping of QTLs for downy mildew resistance at four developmental stages in Chinese cabbage (<Emphasis Type="Italic">Brassica rapa</Emphasis> L. ssp. <Emphasis Type="Italic">pekinensis</Emphasis>)

Specific-locus amplified fragment sequencing is a high-resolution method for genetic mapping, genotyping, and single nucleotide polymorphism (SNP) marker discovery. Previously, a major QTL for downy mildew resistance, BraDM, was mapped to linkage group A08 in a doubled-haploid population derived from Chinese cabbage lines 91–112 and T12–19. The aim of the present study was to improve the linkage map and identify the genetic factors involved in downy mildew resistance. We detected 53,692 high quality SLAFs, of which 7230 were polymorphic, and 3482 of the polymorphic markers were used in genetic map construction. The final map included 1064 bins on ten linkage groups and was 858.98 cM in length, with an average inter-locus distance of 0.81 cM. We identified six QTLs that are involved in downy mildew resistance. The four major QTLs, sBrDM8, yBrDM8, rBrDM8, and hBrDM8, for resistance at the seedling, young plant, rosette, and heading stages were mapped to A08, and are identical to BraDM. The two minor resistance QTLs, rBrDM6 (A06) and hBrDM4 (A04), were active at the rosette and heading stages. The major QTL sBrDM8 defined a physical interval of ~228 Kb on A08, and a serine/threonine kinase family gene, Bra016457, was identified as the possible candidate gene. We report here the first high-density bin map for Chinese cabbage, which will facilitate mapping QTLs for economically important traits and SNP marker development. Our results also expand knowledge of downy mildew resistance in Chinese cabbage and provide three SNP markers (A08-709, A08-028, and A08-018) that we showed to be effective when used in MAS to breed for downy mildew resistance in B. rapa.     

Acid phosphatase gene <Emphasis Type="Italic">GmHAD1</Emphasis> linked to low phosphorus tolerance in soybean,through fine mapping

Key message

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress.

Abstract

Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23–13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80–13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4–16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

     

A single nucleotide mutation of <Emphasis Type="Italic">IspF</Emphasis> gene involved in the MEP pathway for isoprenoid biosynthesis causes yellow-green leaf phenotype in rice

Key message

We identified IspF gene through yellow-green leaf mutant 505ys in rice. OsIspF was expressed in all tissues detected, and its encoded protein was targeted to the chloroplast. On expression levels of genes in this mutant, OsIspF itself and the genes encoding other enzymes of the MEP pathway and chlorophyll synthase were all up-regulated, however, among eight genes associated with photosynthesis, only psaA, psaN and psbA genes for three reaction center subunits of photosystem obviously changed.

Abstract

Isoprenoids are the most abundant natural compounds in all organisms, which originate from the basic five-carbon units isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). In plants, IPP and DMAPP are synthesized through two independent pathways, the mevalonic acid pathway in cytoplasm and the 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway in plastids. The MEP pathway comprises seven enzymatic steps, in which IspF is the fifth enzyme. So far, no IspF gene has been identified in monocotyledonous plants. In this study, we isolated a leaf-color mutant, 505ys, in rice (Oryza sativa). The mutant displayed yellow-green leaf phenotype, reduced level of photosynthetic pigments, and arrested development of chloroplasts. By map-based cloning of this mutant, we identified OsIspF gene (LOC_Os02g45660) showing significant similarity to IspF gene of Arabidopsis, in which a missense mutation occurred in the mutant, resulting in an amino acid change in the encoded protein. OsIspF gene was expressed in all tissues detected, and its encoded protein was targeted to the chloroplast. Further, the mutant phenotype of 505ys was complemented by transformation with the wild-type OsIspF gene. Therefore, we successfully identified an IspF gene in monocotyledonous plants. In addition, real-time quantitative RT-PCR implied that a positive regulation could exist between the OsIspF gene and the genes encoding other enzymes of the MEP pathway and chlorophyll synthase. At the same time, it also implied that the individual genes involved in the MEP pathway might differentially regulated expression levels of the genes associated with photosynthesis.

     

Dissection of a major QTL <Emphasis Type="Italic">qhir1</Emphasis> conferring maternal haploid induction ability in maize

Key message

Among the qhir11 and qhir12 sub-regions of a major QTL qhir1, only qhir11 has significant effect on maternal haploid induction, segregation distortion and kernel abortion.

Abstract

In vivo haploid induction in maize can be triggered in high frequencies by pollination with special genetic stocks called haploid inducers. Several genetic studies with segregating populations from non-inducer x inducer crosses identified a major QTL, qhir1, on chromosome 1.04 contributing to in vivo haploid induction. A recent Genome Wide Association Study using 51 inducers and 1482 non-inducers also identified two sub-regions within the qhir1 QTL region, named qhir11 and qhir12; qhir12 was proposed to be mandatory for haploid induction because the haplotype of qhir11 was also present in some non-inducers and putative candidate genes coding for DNA and amino acid binding proteins were identified in the qhir12 region. To characterize the effects of each sub-region of qhir1 on haploid induction rate, F₂ recombinants segregating for one of the sub-regions and fixed for the other were identified in a cross between CML269 (non-inducer) and a tropicalized haploid inducer TAIL8. To quantify the haploid induction effects of qhir11 and qhir12, selfed progenies of recombinants between these sub-regions were genotyped. F₃ plants homozygous for qhir11 and/or qhir12 were identified, and crossed to a ligueless tester to determine their haploid induction rates. The study revealed that only the qhir11 sub-region has a significant effect on haploid induction ability, besides causing significant segregation distortion and kernel abortion, traits that are strongly associated with maternal haploid induction. The results presented in this study can guide fine mapping efforts of qhir1 and in developing new inducers efficiently using marker assisted selection.

     

Evolution of <Emphasis Type="Italic">fixNOQP</Emphasis> genes encoding cytochrome oxidase with high affinity to oxygen in rhizobia and related bacteria

Many bacteria belonging to the order Rhizobiales have fixNOQP genes which encode cytochrome oxidase with high affinity to oxygen required for oxidative phosphorylation in microaerophilic conditions. There is one copy of the identified fixNOQP operon in ancestral forms of rhizobia (Bradyrhizobium), as well as in their putative evolutionary predecessors (bacteria related to Rhodopseudomonas). At the same time, forms deeply specialized in symbiosis (Rhizobium leguminosarum, Sinorhizobium meliloti) have multiple (2–3) copies, some of them have a high similarity (>90%) to fixNOQP genes of Bradyrhizobium and Rhodopseudomonas, and others have only 30–50% similarity. Two divergent copies fixNOQP are detected in Tardiphaga, which is a representative of the Bradyrhizobiaceae family, lacking the ability to fix N₂ (lack of nif genes encoding the synthesis of nitrogenase) and to induce the formation of nodules on legumes roots (lack of nod genes encoding the synthesis of signal Nod factors activating symbiosis development). The presence of Tardiphaga in nodule bacterial communities from a range of legumes, including Vavilovia formosa (relic representative of the tribe Fabeae, for which R. leguminosarum bv. viciae is the main microsymbiont), suggests that the ancestral gene duplication and subsequent divergence of fixNOQP operon in bacteria related to Tardiphaga opened the possibility of wide dissemination of functionally different copies of this cluster among symbiotically active forms of Rhizobiales. It is possible that the acquisition of fixNOQP genes determines adaptation of bacteria to microaerophilic niches not only in plants nodules but also in their environment (the rhizosphere, rhizoplane, internal portions of soil aggregates).     

Genome-wide identification and expression analysis of sulphate transporter (<Emphasis Type="Italic">SULTR</Emphasis>) genes under sulfur deficiency in <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

Sulphur is an important mineral element for plant growth and development. It involves in a number of metabolic processes with crucial functions. This study has performed a genome-wide analysis of sulfate transporter (SULTR) genes in Brachypodium distachyon. Ten putative SULTR genes were identified in Brachypodium genome. BdSULTR genes included 6–17 exons encoding a protein of 647–693 residues with basic nature. BdSULTR proteins included both sulfate_transp (PF00916) and STAS (PF01740) domains. BdSULTRs were classified into 4 groups based on the phylogenetic distribution. Promoter regions of all BdSULTR genes, except for BdSULTR3;3 and 3;5 included the SURECOREATSULTR11 elements. A considerable structural overlap was identified between superimposed SULTR1;3 and 3;1 proteins, indicating that SULTR1 members may also involve in plant stress response/tolerance like SULTR3 members. Microarray and RNA-Seq analyses also revealed the differential expression of SULTR 1 and 3 genes under different biotic/abiotic stresses. Protein–protein interaction partners of BdSULTRs were mainly related with adenylyl-sulfate kinases, 5′-adenylylsulfate reductases, ATP sulfurylases, and acyl carrier proteins. Moreover, expression profiles of identified BdSULTR genes under S-deficiency were analyzed using RT-qPCR. It was revealed that BdSULTR1;1 and 3;1 are highly expressed in plant roots as ~tenfold and ~fivefold, respectively, while BdSULTR2 (~15-fold) and 3;1 (~twofold) are abundantly expressed in leaf tissues.