Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.

     

Exploiting the <Emphasis Type="Italic">Rht</Emphasis> portfolio for hybrid wheat breeding

Key message

The portfolio of available Reduced height loci (Rht-B1, Rht-D1, and Rht24) can be exploited for hybrid wheat breeding to achieve the desired heights in the female and male parents, as well as in the hybrids, without adverse effects on other traits relevant for hybrid seed production.

Abstract

Plant height is an important trait in wheat line breeding, but is of even greater importance in hybrid wheat breeding. Here, the height of the female and male parental lines must be controlled and adjusted relative to each other to maximize hybrid seed production. In addition, the height of the resulting hybrids must be fine-tuned to meet the specific requirements of the farmers in the target regions. Moreover, this must be achieved without adversely impacting traits relevant for hybrid seed production. In this study, we explored Reduced height (Rht) loci effective in elite wheat and exploited their utilization for hybrid wheat breeding. We performed association mapping in a panel of 1705 wheat hybrids and their 225 parental lines, which besides the Rht-B1 and Rht-D1 loci revealed Rht24 as a major QTL for plant height. Furthermore, we found that the Rht-1 loci also reduce anther extrusion and thus cross-pollination ability, whereas Rht24 appeared to have no adverse effect on this trait. Our results suggest different haplotypes of the three Rht loci to be used in the female or male pool of a hybrid breeding program, but also show that in general, plant height is a quantitative trait controlled by numerous small-effect QTL. Consequently, marker-assisted selection for the major Rht loci must be complemented by phenotypic selection to achieve the desired height in the female and male parents as well as in the wheat hybrids.

     

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.     

Characterisation of a novel quantitative trait locus, <Emphasis Type="Italic">GN4</Emphasis>-<Emphasis Type="Italic">1</Emphasis>, for grain number and yield in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Key message

A novel QTL for grain number, GN4-1, was identified and fine-mapped to an ~ 190-kb region on the long arm of rice chromosome 4.

Abstract

Rice grain yield is primarily determined by three components: number of panicles per plant, grain number per panicle and grain weight. Among these traits, grain number per panicle is the major contributor to grain yield formation and is a crucial trait for yield improvement. In this study, we identified a major quantitative trait locus (QTL) responsible for rice grain number on chromosome 4, designated GN4-1 (a QTL for Grain Number on chromosome 4), using advanced segregating populations derived from the crosses between an elite indica cultivar ‘Zhonghui 8006’ (ZH8006) and a japonica rice ‘Wuyunjing 8’ (WYJ8). GN4-1 was delimited to an ~ 190-kb region on chromosome 4. The genetic effect of GN4-1 was estimated using a pair of near-isogenic lines. The GN4-1 gene from WYJ8 promoted accumulation of cytokinins in the inflorescence and increased grain number per panicle by ~ 17%. More importantly, introduction of the WYJ8 GN4-1 gene into ZH8006 increased grain yield by ~ 14.3 and ~ 11.5% in the experimental plots in 2014 and 2015, respectively. In addition, GN4-1 promoted thickening of the culm and may enhance resistance to lodging. These results demonstrate that GN4-1 is a potentially valuable gene for improvement of yield and lodging resistance in rice breeding.

     

Construction of a genome-anchored,high-density genetic map for melon (<Emphasis Type="Italic">Cucumis melo</Emphasis> L.) and identification of <Emphasis Type="Italic">Fusarium oxysporum</Emphasis> f. sp. <Emphasis Type="Italic">melonis</Emphasis> race 1 resistance QTL

Key message

Four QTLs and an epistatic interaction were associated with disease severity in response to inoculation with Fusarium oxysporum f. sp. melonis race 1 in a recombinant inbred line population of melon.

Abstract

The USDA Cucumis melo inbred line, MR-1, harbors a wealth of alleles associated with resistance to several major diseases of melon, including powdery mildew, downy mildew, Alternaria leaf blight, and Fusarium wilt. MR-1 was crossed to an Israeli cultivar, Ananas Yok’neam, which is susceptible to all of these diseases, to generate a recombinant inbred line (RIL) population of 172 lines. In this study, the RIL population was genotyped to construct an ultra-dense genetic linkage map with 5663 binned SNPs anchored to the C. melo genome and exhibits the overall high quality of the assembly. The utility of the densely genotyped population was demonstrated through QTL mapping of a well-studied trait, resistance to Fusarium wilt caused by Fusarium oxysporum f. sp. melonis (Fom) race 1. A major QTL co-located with the previously validated resistance gene Fom-2. In addition, three minor QTLs and an epistatic interaction contributing to Fom race 1 resistance were identified. The MR-1 × AY RIL population provides a valuable resource for future QTL mapping studies and marker-assisted selection of disease resistance in melon.

     

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.     

QTL analysis and the development of closely linked markers for days to flowering in spring oilseed rape (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Days to flowering (DTF) is an important trait impacting cultivar performance in oilseed rape (Brassica napus L.), but the interaction of all loci controlling this trait in spring-type oilseed rape is not fully understood. We identified quantitative trait loci (QTL) for variation in DTF in a doubled haploid (DH) population from the Qinghai–Tibet Plateau that includes 217 lines derived from a cross between spring-type oilseed rape (B. napus L.) line No. 5246 and line No. 4512, the latter of which is responsive to the effective accumulated temperature (EAT). A linkage map was constructed for the DH population, using 202 SSR and 293 AFLP markers. At least 22 DTF QTL were found in multiple environments. Four major QTL were located on linkage groups A7, C2, C8 and C8. Among these QTL, cqDTFA7a and cqDTFC2a were identified in five environments and individually explained 10.4 and 23.0 % of the trait variation, respectively. cqDTFC8, a major QTL observed in spring environments, and a unique winter environment QTL, qDTFC8-3, were identified; these QTL explained 10.0 and 46.5 % of the phenotypic variation, respectively. Minor QTL (for example, cqDTFC2c) and epistatic interactions seemed evident in this population. Two closely linked SSR markers for cqDTFA7a and cqDTFC8 were developed (G1803 and S034). BnAP1, a B. napus gene with homology to Arabidopsis thaliana that was identified as a cqDTFA7a candidate gene, played a major role in this study. The allelic effects of the major and minor QTL on DTF were further validated in the DH population and in 93 breeding genotypes.     

Integrated analysis of leaf morphological and color traits in different populations of Chinese cabbage (<Emphasis Type="Italic">Brassica rapa</Emphasis> ssp. <Emphasis Type="Italic">pekinensis</Emphasis>)

Key message

QTLs and candidate gene markers associated with leaf morphological and color traits were identified in two immortalized populations of Brassica rapa, which will provide genetic information for marker-assisted breeding.

Abstract

Brassica rapa is an important leafy vegetable consumed worldwide and morphology is a key character for its breeding. To enhance genetic control, quantitative trait loci (QTLs) for leaf color and plant architecture were identified using two immortalized populations with replications of 2 and 4 years. Overall, 158 and 80 QTLs associated with 23 and 14 traits were detected in the DH and RIL populations, respectively. Among them, 23 common robust-QTLs belonging to 12 traits were detected in common loci over the replications. Through comparative analysis, five crucifer genetic blocks corresponding to morphology trait (R, J&U, F and E) and color trait (F, E) were identified in three major linkage groups (A2, A3 and A7). These might be key conserved genomic regions involved with the respective traits. Through synteny analysis with Arabidopsis, 64 candidate genes involved in chlorophyll biosynthesis, cell proliferation and elongation were co-localized within QTL intervals. Among them, SCO3, ABI3, FLU, HCF153, HEMB1, CAB3 were mapped within QTLs for leaf color; and CYCD3;1, CYCB2;4, AN3, ULT1 and ANT were co-localized in QTL regions for leaf size. These robust QTLs and their candidate genes provide useful information for further research into leaf architecture with crop breeding.

     

QTL-seq and marker development for resistance to <Emphasis Type="Italic">Fusarium oxysporum</Emphasis> f. sp. <Emphasis Type="Italic">niveum</Emphasis> race 1 in cultivated watermelon

Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. niveum (Fon), is one of the predominant diseases of watermelon. Resistance to Fon race 1 is conferred by a single major quantitative trait locus (QTL), Fo-1.1, but resolution of this region has been poor due to low marker density. In this study, a combination of whole genome resequencing of bulked segregants (QTL-seq analysis) followed by QTL mapping with kompetitive allele specific PCR (KASP) markers developed across Fo-1.1 successfully increased the resolution from 2.03 to 1.56 Mb and 315 kb, respectively. The linkage of the KASP markers to Fon race 1 resistance across a wide range of watermelon germplasm was validated in a set of elite watermelon cultivars. The linked markers described here provide a breeder-friendly toolkit immediately available for high-throughput genotyping in large-scale breeding programs for fine mapping and incorporation of Fon race 1 resistance in watermelon.     

Fine mapping a major QTL <Emphasis Type="Italic">qFCC7</Emphasis><Subscript><Emphasis Type="Italic">L</Emphasis></Subscript> for chlorophyll content in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) cv. PA64s

Chlorophyll (Chl) content is an important agronomic trait directly affecting the photosynthetic rate. Using a high-density genetic map of 132 recombinant inbred lines (RILs) derived from the cross between 93-11 and PA64s, we detected the quantitative trait loci (QTLs) for Chl content of the top three leaves under two nitrogen (N) conditions at two developmental stages. A total of 32 main-effect QTLs located on chromosomes 1, 4, 5, 6, 7, 8, and 12 were identified, and these QTLs individually accounted for 6.0–20.8?% of the total phenotypic variation. A major QTL qFCC7 _L affecting the Chl content under low N condition was identified, and its positive allele came from PA64s. This QTL might be associated with the ability to tolerate low-N stress in rice. The chromosomal segment substitution line (CSSL) with the corresponding segment from PA64s had a higher SPAD value and photosynthetic rate than 93-11 and showed a lower specific leaf area (SLA). We performed a fine-mapping using a BC₄F₂ population via marker-assisted backcross and finally mapped this QTL to a 124.5 kb interval on the long arm of chromosome 7. Candidate gene analysis showed that there were sequence variations and expression differences in the predicted candidate gene between the two parents. These results suggest that the QTL qFCC7 _L may be useful for breeding the rice varieties with higher photosynthetic rate and grain yield.     

Rapid and efficient <Emphasis Type="Italic">Agrobacterium</Emphasis>-mediated transformation of sorghum (<Emphasis Type="Italic">Sorghum bicolor</Emphasis>) employing standard binary vectors and <Emphasis Type="Italic">bar</Emphasis> gene as a selectable marker

Key message

A rapid and efficient Agrobacterium -mediated transformation system in sorghum has been developed employing standard binary vectors and bar gene as a selectable marker.

Abstract

Sorghum (Sorghum bicolor) is an important food and biofuel crop worldwide, for which improvements in genetic transformation are needed to study its biology and facilitate agronomic and commercial improvement. Here, we report optimization of regeneration and transformation of public sorghum genotype P898012 using standard binary vectors and bar gene as a selectable marker. The tissue culture regeneration time frame has been reduced to 7–12 weeks with a yield of over 18 plants per callus, and the optimized transformation system employing Agrobacterium tumefaciens strain AGL1 and the bar with a MAS promoter achieved an average frequency over 14 %. Of randomly analyzed independent transgenic events, 40–50 % carry single copy of integrated T-DNA. Some independent transgenic events were derived from the same embryogenic callus lines, but a 3:1 Mendelian segregation ratio was found in all transgenic events with single copy as estimated by Southern blots. The system described here should facilitate studies of sorghum biology and agronomic improvement.

     

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.

     

<Emphasis Type="Italic">STAYGREEN</Emphasis> (<Emphasis Type="Italic">CsSGR</Emphasis>) is a candidate for the anthracnose (<Emphasis Type="Italic">Colletotrichum orbiculare</Emphasis>) resistance locus <Emphasis Type="Italic">cla</Emphasis> in Gy14 cucumber

Key message

Map-based cloning identified a candidate gene for resistance to the anthracnose fungal pathogen Colletotrichum orbiculare in cucumber, which reveals a novel function for the highly conserved STAYGREEN family genes for host disease resistance in plants.

Abstract

Colletotrichum orbiculare is a hemibiotrophic fungal pathogen that causes anthracnose disease in cucumber and other cucurbit crops. No host resistance genes against the anthracnose pathogens have been cloned in crop plants. Here, we reported fine mapping and cloning of a resistance gene to the race 1 anthracnose pathogen in cucumber inbred lines Gy14 and WI 2757. Phenotypic and QTL analysis in multiple populations revealed that a single recessive gene, cla, was underlying anthracnose resistance in both lines, but WI2757 carried an additional minor-effect QTL. Fine mapping using 150 Gy14?×?9930 recombinant inbred lines and 1043 F₂ individuals delimited the cla locus into a 32 kb region in cucumber Chromosome 5 with three predicted genes. Multiple lines of evidence suggested that the cucumber STAYGREEN (CsSGR) gene is a candidate for the anthracnose resistance locus. A single nucleotide mutation in the third exon of CsSGR resulted in the substitution of Glutamine in 9930 to Arginine in Gy14 in CsSGR protein which seems responsible for the differential anthracnose inoculation responses between Gy14 and 9930. Quantitative real-time PCR analysis indicated that CsSGR was significantly upregulated upon anthracnose pathogen inoculation in the susceptible 9930, while its expression was much lower in the resistant Gy14. Investigation of allelic diversities in natural cucumber populations revealed that the resistance allele in almost all improved cultivars or breeding lines of the U.S. origin was derived from PI 197087. This work reveals an unknown function for the highly conserved STAYGREEN (SGR) family genes for host disease resistance in plants.

     

Identification and fine mapping of a stay-green gene (<Emphasis Type="Italic">Brnye1</Emphasis>) in pakchoi (<Emphasis Type="Italic">Brassica campestris</Emphasis> L. ssp. <Emphasis Type="Italic">chinensis</Emphasis>)

Key message

Using bulked segregant analysis combined with next-generation sequencing, we delimited the Brnye1 gene responsible for the stay-green trait of nye in pakchoi. Sequence analysis identified Bra019346 as the candidate gene.

Abstract

“Stay-green” refers to a plant trait whereby leaves remain green during senescence. This trait is useful in the cultivation of pakchoi (Brassica campestris L. ssp. chinensis), which is marketed as a green leaf product. This study aimed to identify the gene responsible for the stay-green trait in pakchoi. We identified a stay-green mutant in pakchoi, which we termed “nye”. Genetic analysis revealed that the stay-green trait is controlled by a single recessive gene, Brnye1. Using the BSA-seq method, a 3.0-Mb candidate region was mapped on chromosome A03, which helped us localize Brnye1 to an 81.01-kb interval between SSR markers SSRWN27 and SSRWN30 via linkage analysis in an F₂ population. We identified 12 genes in this region, 11 of which were annotated based on the Brassica rapa annotation database, and one was a functionally unknown gene. An orthologous gene of the Arabidopsis gene AtNYE1, Bra019346, was identified as the potential candidate for Brnye1. Sequence analysis revealed a 40-bp insertion in the second exon of Bra019346 in nye, which generated the TAA stop codon. A candidate gene-specific Indel marker in 1561 F₂ individuals showed perfect cosegregation with Brnye1 in the nye mutant. These results provide a foundation for uncovering the molecular mechanism of the stay-green trait in pakchoi.

     

Genetic mapping of <Emphasis Type="Italic">psl</Emphasis> locus and quantitative trait loci for angular leaf spot resistance in cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.)

One of the most important cucumber diseases is bacterial angular leaf spot (ALS), whose increased occurrence in open-field production has been observed over the last years. To map ALS resistance genes, a recombinant inbred line (RIL) mapping population was developed from a narrow cross of cucumber line Gy14 carrying psl resistance gene and susceptible B10 line. Parental lines and RILs were tested under growth chamber conditions as well as in the field for angular leaf spot symptoms. Based on simple sequence repeat and DArTseq, genotyping a genetic map was constructed, which contained 717 loci in seven linkage groups, spanning 599.7 cM with 0.84 cM on average between markers. Monogenic inheritance of the lack of chlorotic halo around the lesions, which is typical for ALS resistance and related with the presence of recessive psl resistance gene, was confirmed. The psl locus was mapped on cucumber chromosome 5. Two major quantitative trait loci (QTL) psl5.1 and psl5.2 related to disease severity were found and located next to each other on chromosome 5; moreover, psl5.1 was co-located with psl locus. Identified QTL were validated in the field experiment. Constructed genetic map and markers linked to ALS resistance loci are novel resources that can contribute to cucumber breeding programs.     

<Emphasis Type="Italic">qRfg3</Emphasis>, a novel quantitative resistance locus against <Emphasis Type="Italic">Gibberella</Emphasis> stalk rot in maize

Key message

A quantitative trait locus  qRfg3 imparts recessive resistance to maize Gibberella stalk rot. qRfg3 has been mapped into a 350-kb interval and could reduce the disease severity index by ~26.6%.

Abstract

Gibberella stalk rot, caused by the fungal pathogen Fusarium graminearum, severely affects maize yield and grain quality worldwide. To identify more resistance quantitative trait loci (QTLs) against this disease, we analyzed a recombinant inbred line (RIL) population derived from a cross between resistant H127R and susceptible C7-2 inbred lines. Within this population, maize resistance to Gibberella stalk rot had high broad-sense heritability. A major QTL, qRfg3, on chromosome 3 was consistently detected across three field trials, accounting for 10.7–19.4% of the total phenotypic variation. Using a progeny-based sequential fine-mapping strategy, we narrowed qRfg3 down to an interval of ~350 kb. We further demonstrated that qRfg3 is a recessive resistance locus to Gibberella stalk rot that reduced the disease severity index by ~26.6%. Both the gene location and recessive genetic mode distinguish qRfg3 from other stalk rot resistance loci. Hence, qRfg3 is valuable as a complement to existing resistance QTLs to improve maize resistance to Gibberella stalk rot.

     

Genome-wide identification and cadmium induced expression profiling of sulfate transporter (<Emphasis Type="Italic">SULTR</Emphasis>) genes in sorghum (<Emphasis Type="Italic">Sorghum bicolor</Emphasis> L.)

Sulfur is an essential element for all living organisms. Plants can convert inorganic sulfur into organic sulfur compounds by complex enzymatic steps. In this study, we conducted a genome-wide analysis of sulfate transporter genes (SULTRs) in the sorghum (Sorghum bicolor) genome and examined expression profiles of SbSULTR genes under 200 µM cadmium (Cd) exposure. As a result of sorghum genome analysis, 11 SULTR genes were identified, including SbSULTR1;1, SbSULTR1;2, SbSULTR1;3, SbSULTR2;1, SbSULTR2;2, SbSULTR3;1, SbSULTR3;2, SbSULTR3;3, SbSULTR3;4, SbSULTR3;5, and SbSULTR4. Given names are based on phylogeny and chromosomal locations. Except SbSULTR4, all SbSULTR proteins contained Sulfate_transp (PF00916), STAS (PF01740) domains and 12 trans-membrane domains. Phylogenetic analysis revealed that four major groups were identified such as SULTR1, 2, 3, and 4 groups and SULTR4 group was separated to other SULTR groups. In promotor sequences of SbSULTR genes, many diverse cis-acting elements were found mainly related with physiological processes such as light, stress and hormone responsiveness. The expression profiles of SbSULTR genes showed that SULTR1;2, 1;3, 3;3, and 3;5 genes up-regulated in root, while expression level of SULTR4 decreased under 200 µM Cd exposure. The predicted 3D structures of SULTR proteins showed some conformational changes, suggesting functional diversities of SbSULTRs. Finally, results of this study may contribute towards understanding SbSULTR genes and their regulations and roles in Cd stress in sorghum.