Development of microsatellite markers and characterization of simple sequence length polymorphism (SSLP) in rice (Oryza sativa L.)

Genotypic and phenotypic characterization of thermo-sensitive genic male sterile (TGMS) rice lines using simple sequence repeat (SSR) markers and population structure analysis

Thermo-sensitive genic male sterile (TGMS) lines in rice are crucial for hybrid breeding, enhancing genetic diversity by eliminating the need for manual emasculation and restorer genes. These lines induce sterility at high temperatures and restore fertility at low temperatures, in contrast to cytoplasmic male sterility (CMS) systems that require specific restorative genes. This temperature-sensitive mechanism allows for greater flexibility in pairing parent lines, increasing genetic diversity and enabling recombination of beneficial traits in hybrids. A randomized block design (RBD) with three replications was employed for the evaluation of these TGMS rice lines. This study investigates the molecular diversity and genetic variability among TGMS rice lines. Traits such as single plant yield, grains per panicle, glume angle, and pollen fertility showed significant phenotypic and genotypic variation, indicated by high coefficients of variation (PCV and GCV), heritability estimates, and genetic advance as a percentage of mean (GAM). These results highlight substantial genetic variation and selection potential. Euclidean distance matrix analysis of morphological data revealed notable genetic differences. TNAU 137S 1 and TNAU 137S 2 were the most genetically similar, while TNAU 112S and TNAU 114S showed the greatest divergence. Principal component analysis (PCA) revealed distinct genetic profiles among lines such as TNAU 136S, TNAU 113S, TNAU 142S, and TNAU 126S, important for hybrid development. Molecular diversity analysis using simple sequence repeat (SSR) markers identified 90 alleles and eight genetic clusters. Bayesian analysis further confirmed two major subpopulations with significant genetic divergence. These findings support the selective use of parent lines for hybrid rice breeding.

Revisiting the emerging pathosystem of rice sheath blight: deciphering the Rhizoctonia solani virulence, host range, and rice genotype-based resistance

Sheath blight, caused by Rhizoctonia solani AG1 IA, is a challenging disease of rice worldwide. In the current study, nine R. solani isolates, within the anastomosis group AG-1 IA, were isolated, characterized based on their macroscopic and microscopic features, as well as their ability to produce cell wall degrading enzymes (CWDEs), and further molecularly identified via ITS sequencing. Although all isolates were pathogenic and produced typical sheath blight symptoms the susceptible rice cultivar, Sakha 101, R. solani AG1 IA -isolate SHBP9 was the most aggressive isolate. The virulence of isolate SHBP9 was correlated with its overproduction of CWDEs, where it had the highest pectinase, amylase, and cellulase activity in vitro. R. solani AG1 IA -isolate SHBP9 was able to infect 12 common rice-associated weeds from the family Poaceae, as well as over 25 economic crops from different families, except chickpea (Cicer arietinum) from Fabaceae, Rocket (Eruca sativa) from Brassicaceae, and the four crops from Solanaceae. Additionally, rice genotype-based resistance was evaluated using 11 rice genotypes for their response to R. solani isolates, morphological traits, yield components, and using 12 SSR markers linked to sheath blight resistance. Briefly, the tested 11 rice genotypes were divided into three groups; Cluster "I" included only two resistant genotypes (Egyptian Yasmine and Giza 182), Cluster "II" included four moderately resistant genotypes (Egyptian hybrid 1, Giza 178, 181, and 183), whereas Cluster "III" included five susceptible (Sakha 104, 101, 108, Super 300 and Giza 177). Correspondingly, only surface-mycelium growth was microscopically noticed on the resistant cultivar Egyptian Yasmine, as well as the moderately resistant Egyptian hybrid 1, however, on the susceptible Sakha 104, the observed mycelium was branched, shrunk, and formed sclerotia. Accordingly, Indica and Indica/Japonica rice genotypes showed more resistance to R. solani than Japonica genotypes. These findings provide insights into its pathogenicity mechanisms and identify potential targets for disease control which ultimately contributes to the development of sustainable eco-friendly disease management strategies. Moreover, our findings might pave the way for developing resistant rice varieties by using more reliable resistance sources of non-host plants, as well as, rice genotype-based resistance as a genetic resource.

QTL detection for grain shape and fine mapping of two novel locus qGL4 and qGL6 in rice

Rice grain size and grain weight, which have a great influence on rice quality and yield, are complex quantitative traits that are mediated by grain length (GL), grain width (GW), length-to-width ratio (LWR), and grain thickness (GT). In this study, the BC1F2 and BC1F2:3 populations derived from a cross between two indica rice varieties, Guangzhan 63-4S (GZ63-4S) and Dodda, were used to locate quantitative trait loci (QTL) related to grain size. A total of 30 QTL associated with GL, GW and LWR were detected, of which six QTL were scanned repeatedly in both populations. Two QTL, qGL4 and qGL6, were selected for genetic effect validation and were subsequently fine mapped to 2.359 kb and 176 kb, respectively. LOC_Os04g52240 (known as OsKS2/OsKSL2), which encoding an ent-beyerene synthase and as the only gene found in 2.359 kb interval, was proposed to be the candidate for qGL4. Moreover, the grains of qGL4 homozygous mutant plants generated by the CRISPR-Cas9 system became shorter and wider. In addition, the qGL4 allele from GZ63-4S contributes to the increase of yield per plant. Our study not only laid the foundation for further functional study of qGL4 and map-based cloning of qGL6, but also provided genetic resources for the development of high yield and good quality rice varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01502-8.

Marker-assisted enhancement of bacterial blight (Xanthomonas oryzae pv. oryzae) resistance in a salt-tolerant rice variety for sustaining rice production of tropical islands

Introduction

Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae is a major disease of rice, specially in the tropical regions of the world. Developing rice varieties with host resistance against the disease is the most effective and economical solution for managing the disease.

Methods

Pyramiding resistance genes (Xa4, xa5, xa13,and Xa21) in popular rice varieties using marker-assisted backcross breeding (MABB) has been demonstrated as a cost-effective and sustainable approach for establishing durable BB resistance. Here, we report our successful efforts in introgressing four resistance genes (Xa4, xa5, xa13, and Xa21) from IRBB60 to CARI Dhan 5, a popular salt-tolerant variety developed from a somaclonal variant of Pokkali rice, through functional MABB.

Results and discussion

Both BB and coastal salinity are among the major challenges for rice production in tropical island and coastal ecosystems. Plants with four, three, and two gene pyramids were generated, which displayed high levels of resistance to the BB pathogen at the BC3F2 stage. Under controlled salinity microplot environments, the line 131-2-175-1223 identified with the presence of three gene pyramid (Xa21+xa13+xa5) displayed notable resistance across locations and years as well as exhibited a salinity tolerance comparable to the recurrent parent, CARI Dhan 5. Among two BB gene combinations (Xa21+xa13), two lines, 17-1-69-334 and 46-3-95-659, demonstrated resistance across locations and years, as well as salt tolerance and grain production comparable to CARI Dhan 5. Besides salinity tolerance, five lines, 17-1-69-179, 46-3-95-655, 131-2-190-1197, 131-2-175-1209, and 131-2-175-1239, exhibited complete resistance to BB disease. Following multilocation testing, potential lines have been identified that can serve as a prospective candidate for producing varieties for the tropical Andaman and Nicobar Islands and other coastal locations, which are prone to BB and coastal salinity stresses.

Wild and cultivated allele effects on rice phenotypic traits in reciprocal backcross populations between Oryza rufipogon and two cultivars, O. sativa Nipponbare and IR36

A total of four populations of reciprocal backcross recombinant inbred lines were produced from a cross between a wild accession of Oryza rufipogon W630 and two major cultivars, O. sativa Japonica Nipponbare and Indica IR36. Using these populations, quantitative trait locus (QTL) analysis for eight morphological traits (culm length, panicle length, days to heading, panicle shape, pericarp color, hull color, seed shattering and seed awning) was carried out, and the putative QTL regions were compared among the populations. The QTLs with strong allele effects were commonly detected for culm length, panicle shape, pericarp color and hull color in all four populations, and their peak locations were close to the major genes of sd1, Spr3, Rc and Bh4, respectively. For panicle length and days to heading, some QTL regions overlapped between two or three populations. In the case of seed shattering and seed awning, strong wild allele effects at major loci were observed only in the populations with cultivated backgrounds. Since the wild and cultivated alleles have never been evaluated in the reciprocal genetic backgrounds, the present results provide new information on gene effects in breeding and domestication studies.

Chromosomal Location of xa19, a Broad-Spectrum Rice Bacterial Blight Resistant Gene from XM5, a Mutant Line from IR24

Bacterial blight is an important rice disease caused by bacteria named Xanthomonas oryzae pv. oryzae (Xoo). XM5 is an Xoo resistant mutant line with the genetic background of IR24, an Indica Xoo susceptible cultivar, induced by a chemical mutagen N-methyl-N-nitrosourea (MNU). XM5 carries a recessive Xoo resistant gene, xa19. Trisomic analysis was conducted using the cross between XM5 and the trisomic series under the genetic background of IR24, showing that xa19 was located on chromosome 7. The approximate chromosomal location was found using 37 surely resistant plants in the F₂ population from XM5 × Kinmaze, which was susceptible to most Japanese Xoo races. The IAS44 line carries a Japonica cultivar Asominori chromosomal segment covering the xa19 locus under the IR24 genetic background. Linkage analysis using the F₂ population from the cross between XM5 and IAS44 revealed that xa19 was located within the 0.8 cM region between RM8262 and RM6728. xa19 is not allelic to the known Xoo resistant genes. However, its location suggests that it might be allelic to a lesion-mimic mutant gene spl5, some alleles of which are resistant to several Xoo races. Together with xa20 and xa42, three Xoo resistant genes were induced from IR24 by MNU. The significance of chemical mutagen as a source of Xoo resistance was discussed.

Variations in Grain Traits among Local Rice Varieties Collected More Than Half-Century Ago in Indochinese Countries

More than half-century ago, local rice varieties were collected from Indochinese countries (Cambodia, Thailand, Laos, and Vietnam). Of these, 162 local varieties were examined for 7 grain-size traits: seed length/width/thickness, brown rice length/width/thickness, and 100-seed weight. Since these traits varied considerably, a survey of functional mutations was performed in the genes related to these traits. In total, 19 markers (12 InDel and 7 dCAPS markers) were used to investigate the mutations at 14 grain-size loci of GW2, GS2, qLGY3, GS3, GL3.1, TGW3, GS5, GW5, GS6, TGW6, GW6a, GLW7, GL7, and GW8. Significant allele effects were observed with six markers detecting base substitution mutations at GW2 and GS3 and insertion/deletion mutations at GS5, GW5, and GW6a, suggesting that these mutations might have affected the grain trait and caused variation among local varieties in the Indochinese countries. In addition to grain size, the hull color, grain color, and glutinosity were also examined using a survey of loss-of-function mutations at major responsible loci. Most phenotypes were reflected based on functional mutations at these loci. Since the local varieties have wide genetic variation, they are important genetic resources for future rice breeding.

Fine mapping of qTGW2b and qGL9, two minor QTL conferring grain size and weight in rice

Rice grain size is a key determinant of both grain yield and quality. In this study, we conducted QTL mapping on grain size using a recombinant inbred line (RIL) population derived from a cross between japonica variety Beilu130 (BL130) and indica variety Jin23B (J23B). A total of twenty-two QTL related to grain length (GL), grain width (GW), grain length-to-width ratio (LWR), grain thickness (GT), and thousand grain weight (TGW) were detected under two environments, and 14 of them were repeatedly detected. Two minor QTL, qTGW2b and qGL9, were validated and further delimited to regions of 631 kb and 272 kb, respectively. Parental sequence comparison of genes expressed in inflorescence in corresponding candidate regions identified frameshifts in the exons of LOC_Os02g38690 and LOC_Os02g38780, both of which encode protein phosphatase 2C-containing protein, and LOC_Os09g29930, which encodes a BIM2 protein. Scanning electron microscopy (SEM) analysis revealed that the increase of cell size rather than cell number caused the differences in grain size between NILs of qTGW2b and qGL9. Quantitative RT-PCR analysis showed that the expression levels of EXPA4, EXPA5, EXPA6, EXPB3, EXPB4, and EXPB7 were significantly different in both qTGW2b NILs and qGL9 NILs. Our results lay the foundation for the cloning of qTGW2b and qGL9, and provide genetic materials for the improvement of rice yield and quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01328-2.

Fine Mapping and Cloning of a Major QTL qph12, Which Simultaneously Affects the Plant Height, Panicle Length, Spikelet Number and Yield in Rice (Oryza sativa L.)

Plant height is one of the most important agronomical traits in rice (Oryza sativa L.). Introducing the semidwarf rice in the 1960s significantly enhanced the rice yield potential in Asia. Implementing near-isogenic lines (NILs) is the most powerful tool for the identification and fine mapping of quantitative trait loci (QTLs). In this study, 176 NILs were produced from the crossing and back-crossing of two rice cultivars. Specifically, the indica rice cultivar Jiafuzhan served as a recipient, and the restorer japonica cultivar Hui1586 served as a donor. Using the 176 NILs, we identified a novel major QTL for reduced plant height in the NIL36 line. The qph12 QTL was mapped to a 31 kb genomic region between the indel markers Indel12-29 and Indel12-31. The rice genome annotation indicated the presence of three candidate genes in this genomic region. Through gene prediction and cDNA sequencing, we confirmed that LOC_Os12g40890 (qPH12) is the target gene in the NIL36 line. Further analysis showed that the qph12 QTL is caused by a 1 bp deletion in the first exon that resulted in premature termination of the qPH12. Knockout experiments showed that the qph12 QTL is responsible for the reduced plant height phenotype of the NIL36 line. Although the qph12 gene from the NIL36 line showed a shorter panicle length, fewer spikelets per panicle and a lower plant grain yield, the plant also exhibited a lower plant height. Taken together, our results revealed that the qph12 have good specific application prospects in future rice breeding.

Development of Wide-Compatible Indica Lines by Pyramiding Multiple Neutral Alleles of Indica-Japonica Hybrid Sterility Loci

Since the development of indica hybrid rice in the 1970s, great success has been achieved in hybrid rice production in China and around the world. The utilization of inter-subspecific indica-japonica hybrid rice has always been considered due to its stronger heterosis characteristics. However, indica-japonica hybrids face a serious problem of sterility, which hinders the exploitation of their heterosis. In the past decades, the genetic basis of indica-japonica hybrid sterility has been well studied. It was found that in sterile indica-japonica hybrids, female sterility was mainly controlled by the S5 locus and male sterility by the Sa, Sb, Sc, Sd, and Se loci. In this study, we developed wide-compatible indica lines (WCILs) by pyramiding multiple neutral (n) alleles of the hybrid sterility loci. First, we identified Sⁿ alleles of the loci in single-segment substitution lines (SSSLs) in the genetic background of indica Huajingxian 74 (HJX74). Then, the Sⁿ alleles of S5, Sb, Sc, Sd, and Se loci in SSSLs were pyramided in the HJX74 genetic background. The WCILs carrying Sⁿ alleles at the S5, Sb, Sc, Sd, and Se loci showed wide compatibility with indica and japonica rice varieties. Therefore, the WCILs will be used to develop inter-subspecific indica-japonica hybrid rice with normal fertility.

Fine Mapping of Two Interacting Loci for Transmission Ratio Distortion in Rice (Oryza sativa L.)

Transmission ratio distortion (TRD) denotes the observed allelic or genotypic frequency deviation from the expected Mendelian segregation ratios in the offspring of a heterozygote. TRD can severely hamper gene flow between and within rice species. Here, we report the fine mapping and characterization of two loci (TRD4.1 and TRD4.2) for TRD using large F₂ segregating populations, which are derived from rice chromosome segment substitution lines, each containing a particular genomic segment introduced from the japonica cultivar Nipponbare (NIP) into the indica cultivar Zhenshan (ZS97). The two loci exhibited a preferential transmission of ZS97 alleles in the derived progeny. Reciprocal crossing experiments using near-isogenic lines harboring three different alleles at TRD4.1 suggest that the gene causes male gametic selection. Moreover, the transmission bias of TRD4.2 was diminished in heterozygotes when they carried homozygous TRD4.1 ^ZS97. This indicates an epistatic interaction between these two loci. TRD4.2 was mapped into a 35-kb region encompassing one candidate gene that is specifically expressed in the reproductive organs in rice. These findings broaden the understanding of the genetic mechanisms of TRD and offer an approach to overcome the barrier of gene flow between the subspecies in rice, thus facilitating rice improvement by introgression breeding.

QTL Analysis Revealed One Major Genetic Factor Inhibiting Lesion Elongation by Bacterial Blight (Xanthomonas oryzae pv. oryzae) from a japonica Cultivar Koshihikari in Rice

Xanthomonas oryzae pv. oryzae (Xoo) is a pathogen that has ravaged the rice industry as the causal agent of bacterial blight (BB) diseases in rice. Koshihikari (KO), an elite japonica cultivar, and ARC7013 (AR), an indica cultivar, are both susceptible to Xoo. Their phenotypic characteristics reveal that KO has shorter lesion length than that of AR. The F₂ population from KO × AR results in continuous distribution of lesion length by inoculation of an Xoo race (T7147). Consequently, quantitative trait loci (QTL) mapping of the F₂ population is conducted, covering 12 chromosomes with 107 simple sequence repeat (SSR) and insertion/deletion (InDel) genetic markers. Three QTLs are identified on chromosomes 2, 5, and 10. Of them, qXAR5 has the strongest resistance variance effect of 20.5%, whereas qXAR2 and qXAR10 have minor QTL effects on resistance variance, with 3.9% and 2.3%, respectively, for a total resistance variance of 26.7%. The QTLs we examine for this study differ from the loci of BB resistance genes from earlier studies. Our results can help to facilitate understanding of genetic and morphological fundamentals for use in rice breeding programs that are more durable against evolving Xoo pathogens and uncertain climatic temperature.

NILs of Cold Tolerant Japonica Cultivar Exhibited New QTLs for Mineral Elements in Rice

Chilling stress at booting stage can cause floret deterioration and sterility by limiting the supply of food chain and the accumulation of essential mineral elements resulting in reduction of yield and grain quality attributes in rice. Genomic selection of chilling tolerant rice with reference to the accumulation of mineral elements will have great potential to cope with malnutrition and food security in times of climate change. Therefore, a study was conducted to explore the genomic determinants of cold tolerance and mineral elements content in near-isogenic lines (NILs) of japonica rice subjected to chilling stress at flowering stage. Detailed morphological analysis followed by quantitative analysis of 17 mineral elements revealed that the content of phosphorus (P, 3,253 mg/kg) and potassium (K, 2,485 mg/kg) were highest while strontium (Sr, 0.26 mg/kg) and boron (B, 0.34 mg/kg) were lowest among the mineral elements. The correlation analysis revealed extremely positive correlation of phosphorus (P) and copper (Cu) with most of the cold tolerance traits. Among all the effective ear and the second leaf length correlation was significant with half of the mineral elements. As a result of comparative analysis, some QTLs (qBRCC-1, qBRCIC-2, qBRZC-6, qBRCHC-6, qBRMC-6, qBRCIC-6a, qBRCIC-6b, qBRCHC-6, and qBRMC-6) identified for calcium (Ca), zinc (Zn), chromium (Cr) and magnesium (Mg) on chromosome number 1, 2, and 6 while, a novel QTL (qBCPC-1) was identified on chromosome number 1 for P element only. These findings provided bases for the identification of candidate genes involved in mineral accumulation and cold tolerance in rice at booting stage.

Identifying Wild Versus Cultivated Gene-Alleles Conferring Seed Coat Color and Days to Flowering in Soybean

Annual wild soybean (G. soja) is the ancestor of the cultivated soybean (G. max). To reveal the genetic changes from soja to max, an improved wild soybean chromosome segment substitution line (CSSL) population, SojaCSSLP5, composed of 177 CSSLs with 182 SSR markers (SSR-map), was developed based on SojaCSSLP1 generated from NN1138-2(max)×N24852(soja). The SojaCSSLP5 was genotyped further through whole-genome resequencing, resulting in a physical map with 1366 SNPLDBs (SNP linkage-disequilibrium blocks), which are composed of more markers/segments, shorter marker length and more recombination breakpoints than the SSR-map and caused 721 new wild substituted segments. Using the SNPLDB-map, two loci co-segregating with seed-coat color (SCC) and six loci for days to flowering (DTF) with 88.02% phenotypic contribution were identified. Integrated with parental RNA-seq and DNA-resequencing, two SCC and six DTF candidate genes, including three previously cloned (G, E2 and GmPRR3B) and five newly detected ones, were predicted and verified at nucleotide mutant level, and then demonstrated with the consistency between gene-alleles and their phenotypes in SojaCSSLP5. In total, six of the eight genes were identified with the parental allele-pairs coincided to those in 303 germplasm accessions, then were further demonstrated by the consistency between gene-alleles and germplasm phenotypes. Accordingly, the CSSL population integrated with parental DNA and RNA sequencing data was demonstrated to be an efficient platform in identifying candidate wild vs. cultivated gene-alleles.

Fine Mapping of a Novel Major Quantitative Trait Locus, qPAA7, That Controls Panicle Apical Abortion in Rice

The panicle apical abortion (PAA) causes severe yield losses in rice production, but details about its development and molecular basis remain elusive. Here, we detected PAA quantitative trait loci (QTLs) in three environments using a set of chromosome segment substitution lines (CSSLs) that was constructed with indica Changhui121 as the recurrent parent and japonica Koshihikari as the donor parent. First, we identified a novel major effector quantitative trait locus, qPAA7, and selected a severe PAA line, CSSL176, which had the highest PAA rate among CSSLs having Koshihikari segments at this locus. Next, an F₂ population was constructed from a cross between CSS176 and CH121. Using F₂ to make recombinantion analysis, qPAA7 was mapped to an 73.8-kb interval in chromosome 7. Among nine candidate genes within this interval, there isn't any known genes affecting PAA. According to the gene annotation, gene expression profile and alignment of genomic DNA, LOC_Os07g41220 and LOC_Os07g41280 were predicted as putative candidate genes of qPAA7. Our study provides a foundation for cloning and functional characterization of the target gene from this locus.

Characterization of Quantitative Trait Loci for Germination and Coleoptile Length under Low-Temperature Condition Using Introgression Lines Derived from an Interspecific Cross in Rice

Previously, five putative quantitative trait loci (QTLs) for low-temperature germination (LTG) have been detected using 96 BC₃F₈ lines derived from an interspecific cross between the Korean japonica cultivar “Hwaseong” and Oryza rufipogon. In the present study, two introgression lines, CR1517 and CR1518, were used as parents to detect additional QTLs and analyze interactions among QTLs for LTG. The F₂ population (154 plants) along with parental lines, Hwaseong and O. rufipogon, were evaluated for LTG and coleoptile length under low-temperature conditions (13 °C). Among five QTLs for LTG, two major QTLs, qLTG1 and qLTG3, were consistently detected at 6 and 7 days after incubation. Three minor QTLs were detected on chromosomes 8 and 10. Two QTLs, qLTG10.1 and qLTG10.2, showing linkage on chromosome 10, exerted opposite effects with the Hwaseong allele at qLTG10.2 and the O. rufipogon allele at qLTG10.1 respectively, in turn, increasing LTG. Interactions among QTLs were not significant, implying that the QTLs act in an additive manner. Near-isogenic line plants with the combination of favorable alleles from O. rufipogon and Hwaseong exhibited higher LTG than two introgression lines. With regard to coleoptile length, three QTLs observed on chromosomes 1, 3, and 8 were colocalized with QTLs for LTG, suggesting the pleiotropy of the single gene at each locus. According to the results, the introgression of favorable O. rufipogon alleles could hasten the development of rice with high LTG and high coleoptile elongation in japonica cultivars.

Development of Chromosome Segment Substitution Lines (CSSLs) Derived from Guangxi Wild Rice (Oryza rufipogon Griff.) under Rice (Oryza sativa L.) Background and the Identification of QTLs for Plant Architecture, Agronomic Traits and Cold Tolerance

Common wild rice contains valuable resources of novel alleles for rice improvement. It is well known that genetic populations provide the basis for a wide range of genetic and genomic studies. In particular, chromosome segment substitution lines (CSSLs) ais a powerful tool for fine mapping of quantitative traits, new gene discovery and marker-assisted breeding. In this study, 132 CSSLs were developed from a cultivated rice (Oryza sativa) cultivar (93-11) and common wild rice (Oryza rufipogon Griff. DP30) by selfing-crossing, backcrossing and marker-assisted selection (MAS). Based on the high-throughput sequencing of the 93-11 and DP30, 285 pairs of Insertion-deletions (InDel) markers were selected with an average distance of 1.23 Mb. The length of this DP30-CSSLs library was 536.4 cM. The coverage rate of substitution lines cumulatively overlapping the whole genome of DP30 was about 91.55%. DP30-CSSLs were used to analyze the variation for 17 traits leading to the detection of 36 quantitative trait loci (QTLs) with significant phenotypic effects. A cold-tolerant line (RZ) was selected to construct a secondary mapping F2 population, which revealed that qCT2.1 is in the 1.7 Mb region of chromosome 2. These CSSLs may, therefore, provide powerful tools for genome wide large-scale gene discovery in wild rice. This research will also facilitate fine mapping and cloning of QTLs and genome-wide study of wild rice. Moreover, these CSSLs will provide a foundation for rice variety improvement.

InDel Marker Based Estimation of Multi-Gene Allele Contribution and Genetic Variations for Grain Size and Weight in Rice (Oryza sativa L.)

The market success of any rice cultivar is exceedingly dependent on its grain appearance, as well as its grain yield, which define its demand by consumers as well as growers. The present study was undertaken to explore the contribution of nine major genes, qPE9~1, GW2, SLG7, GW5, GS3, GS7, GW8, GS5, and GS2, in regulating four size and weight related traits, i.e., grain length (GL), grain width (GW), grain thickness (GT), and thousand grain weight (TGW) in 204 diverse rice germplasms using Insertion/Deletion (InDel) markers. The studied germplasm displayed wide-ranging variability in the four studied traits. Except for three genes, all six genes showed considerable association with these traits with varying strengths. Whole germplasm of 204 genotypes could be categorized into three major clusters with different grain sizes and weights that could be utilized in rice breeding programs where grain appearance and weight are under consideration. The study revealed that TGW was 24.9% influenced by GL, 37.4% influenced by GW, and 49.1% influenced by GT. Hence, assuming the trend of trait selection, i.e., GT > GW > GL, for improving TGW in the rice yield enhancement programs. The InDel markers successfully identified a total of 38 alleles, out of which 27 alleles were major and were found in more than 20 genotypes. GL was associated with four genes (GS3, GS7, GW8, and GS2). GT was also found to be regulated by four different genes (GS3, GS7, GW8, and GS2) out of the nine studied genes. GW was found to be under the control of three studied genes (GW5, GW8, and GS2), whereas TGW was found to be under the influence of four genes (SLG7, GW5, GW8, and GS5) in the germplasm under study. The Unweighted Pair Group Method with Arithmetic means (UPGMA) tree based on the studied InDel marker loci segregated the whole germplasm into three distinct clusters with dissimilar grain sizes and weights. A two-dimensional scatter plot constructed using Principal Coordinate Analysis (PCoA) based on InDel markers further separated the 204 rice germplasms into four sub-populations with prominent demarcations of extra-long, long, medium, and short grain type germplasms that can be utilized in breeding programs accordingly. The present study could help rice breeders to select a suitable InDel marker and in formulation of breeding strategies for improving grain appearance, as well as weight, to develop rice varieties to compete international market demands with higher yield returns. This study also confirms the efficient application of InDel markers in studying diverse types of rice germplasm, allelic frequencies, multiple-gene allele contributions, marker-trait associations, and genetic variations that can be explored further.

Data on genetic potentiality of folk rice (Oryza sativa L.) genotypes from Koraput, India in reference to drought tolerance traits

Precise physiological and molecular marker-based assessment provides information about the extent of genetic diversity, which helps for effective breeding programmes. We have conducted detailed physiological and molecular marker-based assessment of selected eight indigenous rice landraces from Koraput, India along with tolerant (N22) and susceptible (IR64) check varieties under control and simulated drought stress using polyethylene glycol (PEG) 6000. After exposure to different levels of drought stress, relative germination performance (RGP), seedling vigour index (SVI) and relative growth index (RGI) were significantly declined in all the rice landraces compared to the control plants and significant varietal differences were observed. Genetic relationship among the studied rice landraces was assessed with 24 previously reported drought tolerance linked Simple Sequence Repeat (SSR) markers. A total of 53 alleles were detected at the loci of the 24 markers across the 10 rice accessions. The Nei's gene diversity (He) and the polymorphism information content (PIC) ranged from 0 to 0.665 and 0 to 0.687, respectively. Six SSR loci, RM276, RM411, RM3, RM263, RM216 and RM28199, provided the highest PIC values and are potential for exploring the genetic diversity of studied rice lines for drought tolerance. Four rice genotypes (Butkichudi, Haldichudi, Machakanta and Kalajeera) showed the highest genetic distance with tolerant check variety (N22) and can be considered as valuable genetic resources for drought breeding program.

Identification of differentially expressed transcripts at critical developmental stages in sorghum [Sorghum bicolor (L.) Moench] in relation to grain yield heterosis

Evaluation of a set of 10 F₁ hybrids along with their female (27A and 7A) and male parents (C 43, RS 673, RS 627, CB 26, and CB 29) for grain yield and its component traits revealed that grain yield/plant followed by panicle weight, primary branches/panicle, and 100-seed weight exhibited high levels of heterosis. Eight hybrids exhibited 50% or more mid-parent heterosis for grain yield/plant, of which, one hybrid (27A × RS673) recorded heterobeltiosis above 50% (73.61%). Differential display analysis generated about 2995 reproducible transcripts, which were categorized as UPF₁-expressed in any one of the parents and F₁ (10.53-14.76%), BPnF₁-expressed in both parents but not in F₁ (4.56-11.44%), UPnF₁-expressed in either of the parents and not in F₁ (17.95-27.40%), F₁nBP-expressed only in F₁ but not in either of the parents (14.39-20.54%), and UET-expressed in both parents and F₁ (34.52-42.43%). A comparison between high and low heterotic hybrids revealed that the proportions of UPF₁ and F₁nBP transcript patterns were much higher in the former (21.31% and 45.24%) as compared to the latter (16.67% and 32.14%) at the booting and flowering stage, respectively, indicating the role of over-dominance and dominance in the manifestation of grain yield heterosis. Significant positive correlations were observed for differential transcript patterns with mid-parent and better-parent heterosis for the components of grain yield such as primary branches (0.63 and 0.61 at p < 0.01) and 100-seed weight (0.64 and 0.52 at p < 0.01). Cloning and sequence analysis of 16 transcripts that were differentially expressed in hybrids and their parental lines revealed that they code for genes involved in basic cellular processes, cellulose biosynthesis, and assimilate partitioning between various organs and allocation between various pathways, pyrimidine, and polyamine biosynthesis, enhancing ATP production and regulation of plant growth and development.

The potentiality of rice microsatellite markers in assessment of cross-species transferability and genetic diversity of rice and its wild relatives

The main aim of this study is to assess the potentiality of SSR markers for the identification of the cross-species transferability frequency in a large set of the diverse genome types of wild relative rice along with cultivated rice. Here, we used 18 different rice genotypes representing nine different genome types with 70 SSR markers to investigate the potentiality of cross-species transferability rate. The overall cross-species transferability of SSR markers across the 18 rice genotypes ranged from 38.9% (RM280 and RM447) to 100% (RM490, RM318, RM279, RM18877 and RM20033, RM19303) with an average of 76.58%. Also, cross-species transferability across chromosome ranged from 54.4% (chromosome 4) to 86.5% (chromosome 2) with an average of 74.35%. The polymorphism information content of the markers varied from 0.198 (RM263) to 0.868 (RM510) with a mean of 0.549 ± 0.153, showing high discriminatory power. The highest rate of cross-transferability was observed in O. rufipogon (97%), The highest rate of cross-species transferability was in O. rufipogon (97.00%), followed by O. glaberrima (94.20%), O. nivara (92.80%), Swarna (92.80%), O. longistaminata (91.40%), O. eichingeri (90%), O. barthii (88.50%), O. alta (82.80%), O. australiensis (77.10%), O. grandiglumis (74.20%), O. officinalis (74.20%), Zizania latifolia (70.00%), O. latifolia (68.50%), O. brachyantha (62.80%), Leersia perrieri (57.10%) and O. ridleyi (41.40%) with least in O. coarctata (28.50%). A total of 341 alleles from 70 loci were detected with the number of alleles per locus ranged from 2 to 12. Based on dendrogram analysis, the AA genome groups was separated as distinct group from the rest of the genome types. Similarly, principal coordinate analysis and structure analysis clearly separated the AA genome type from the rest of the genome types. Through the analysis of molecular variance, more variance (51%) was observed among the individual, whereas less (14%) was observed among the population. Thus, our findings may offer a valuable resource for studying the genetic diversity and relationship to facilitate the understanding of the complex mechanism of the origin and evolutionary processes of different Oryza species and wild relative rice.

Development and Characterization of EST-SSR Markers From RNA-Seq Data in Phyllostachys violascens

Bamboo are woody grass species containing important economic and ecological values. Lei bamboo (Phyllostachys violascens) is a kind of shoot-producing bamboo species with the highest economic yield per unit area. However, identifying different varieties of Lei bamboo based on morphological characteristics is difficult. Microsatellites play an important role in plant identification and genetic diversity analysis and are superior to other molecular markers. In this study, we identified 18,356 expressed sequence tag-simple sequence repeat (EST-SSR) loci in Lei bamboo transcriptome data. A total of 11,264 primer pairs were successfully designed from unigenes of all EST-SSR loci, and 96 primer pairs were randomly selected and synthesized. A total of 54 primer pairs were used for classifying 16 Lei bamboo varieties and 10 different Phyllostachys species. The number of polymorphism alleles among the 54 primer pairs ranged from 3 to 12 for P. violascens varieties and 3 to 20 for Phyllostachys. The phylogenetic tree based on polymorphism alleles successfully distinguished 16 P. violascens varieties and 10 Phyllostachys species. Our study provides abundant EST-SSR resources that are useful for genetic diversity analysis and molecular verification of bamboo and suggests that SSR markers developed from Lei bamboo are more efficient and reliable than ISSR, SRAP or AFLP markers.

Fine mapping of a major quantitative trait locus, qgnp7(t), controlling grain number per panicle in African rice (Oryza glaberrima S.)

Grain number per panicle is a major component of rice yield that is typically controlled by many quantitative trait loci (QTLs). The identification of genes controlling grain number per panicle in rice would be valuable for the breeding of high-yielding rice. The Oryza glaberrima chromosome segment substitution line 9IL188 had significantly smaller panicles compared with the recurrent parent 9311. QTL analysis in an F₂ population derived from a cross between 9IL188 and 9311 revealed that qgnp7(t), a major QTL located on the short arm of chromosome 7, was responsible for this phenotypic variation. Fine mapping was conducted using a large F₃ population containing 2250 individuals that were derived from the F₂ heterozygous plants. Additionally, plant height, panicle length, and grain number per panicle of the key F₄ recombinant families were examined. Through two-step substitution mapping, qgnp7(t) was finally localized to a 41 kb interval in which eight annotated genes were identified according to available sequence annotation databases. Phenotypic evaluation of near isogenic lines (NIL-qgnp7 and NIL-qGNP7) indicated that qgnp7(t) has pleiotropic effects on rice plant architecture and panicle structure. In addition, yield estimation of NILs indicated that qGNP7(t) derived from 9311 is the favorable allele. Our results provide a foundation for isolating qgnp7(t). Markers flanking this QTL will be a useful tool for the marker-assisted selection of favorable alleles in O. glaberrima improvement programs.

Tandem Repeats Contribute to Coding Sequence Variation in Bumblebees (Hymenoptera: Apidae)

Tandem repeats (TRs) are highly dynamic regions of the genome. Mutations at these loci represent a significant source of genetic variation and can facilitate rapid adaptation. Bumblebees are important pollinating insects occupying a wide range of habitats. However, to date, molecular mechanisms underlying the potential adaptation of bumblebees to diverse habitats are largely unknown. In the present study, we investigate how TRs contribute to genetic variation in bumblebees, thus potentially facilitating adaptation. We identified 26,595 TRs from the assembled 18 chromosome sequences of the buff-tailed bumblebee (Bombus terrestris), 66.7% of which reside in genic regions. We also compared TRs found in B. terrestris with those present in the assembled genome sequence of a congener, B. impatiens. We found that a total of 1,137 TRs were variable in length between the two sequenced bumblebee species, and further analysis reveals that 101 of them are located within coding regions. These 101 TRs are responsible for coding sequence variation and correspond to protein sequence length variation between the two bumblebee species. The variability of identified TRs in coding regions between bumblebees was confirmed by PCR amplification of a subset of loci. Functional classification of bumblebee genes where coding sequences include variable-length TRs suggests that a majority of genes (87%) that could be assigned to a protein class are related to transcriptional regulation. Our results show that TRs contribute to coding sequence variation in bumblebees, and thus may facilitate the adaptation of bumblebees through diversifying proteins involved in controlling gene expression.

Genetic analysis of roots and shoots in rice seedling by association mapping

The vigorous shoots and roots help to improve drought resistance and post-transplanting recovery in rice seedlings (Oryza sativa L.). Hundreds of loci related to root system have been identified recently, but little research has been done on shoot traits, and the relationship between roots and shoots development is also still unclear. The objective of this study was to identify associated loci for roots and shoots in rice seedlings as well as to screen pleiotropic QTLs involved in coordinated development of roots and shoots. Using mini core collection of 273 cultivated rice accessions and 280 simple-sequence repeat markers, we investigated six traits [root length (RL), root thickness (RT), root weight (RW), shoot length (SL), shoot weight (SW) and ratio of root-to-shoot mass] in seedlings. Study was performed in hydroponic medium and genetic analysis was performed by association mapping using general linear model (GLM) with population structure (Q) and mixed linear model (MLM) involving Q and familial relatedness (K). Two subgroups indica and japonica showed significant differences in RT, RW and SW. Maximum correlation was observed between RW and SW. Using GLM 65 QTLs for root and 43 QTLs associated with shoot traits were detected. Among them, seven QTLs were present between RL and RW and five common QTLs were detected between SL and SW with high phenotypic variation effects (PVEs). Two key pleiotropic QTLs were also identified involved in collaborative development of roots and shoots in rice seedlings. Importantly, 17 and 10 QTLs were identified for root and shoot traits respectively in both studies of GLM and MLM. More common QTLs with high PVEs between root and shoot traits suggested that longitudinal growth (RL and SL) played an important role in accumulation of biomass (RW and SW). Considering the obvious phenotypic differences and fewer common QTLs between indica and japonica, we suggested that there could be different mechanisms of seedling development between both subpopulations. Key pleiotropic QTLs and QTLs identified for root and shoot traits in both studies of GLM and MLM could be preferentially used in marker-assisted breeding for strong rice seedling.

Combining QTL-seq and linkage mapping to fine map a wild soybean allele characteristic of greater plant height

BACKGROUND

Plant height (PH) is an important agronomic trait and is closely related to yield in soybean [Glycine max (L.) Merr.]. Previous studies have identified many QTLs for PH. Due to the complex genetic background of PH in soybean, there are few reports on its fine mapping.

RESULTS

In this study, we used a mapping population derived from a cross between a chromosome segment substitution line CSSL3228 (donor N24852 (G. Soja), a receptor NN1138-2 (G. max)) and NN1138-2 to fine map a wild soybean allele of greater PH by QTL-seq and linkage mapping. We identified a QTL for PH in a 1.73 Mb region on soybean chromosome 13 through QTL-seq, which was confirmed by SSR marker-based classical QTL mapping in the mapping population. The linkage analysis showed that the QTLs of PH were located between the SSR markers BARCSOYSSR_13_1417 and BARCSOYSSR_13_1421 on chromosome 13, and the physical distance was 69.3 kb. RT-PCR and sequence analysis of possible candidate genes showed that Glyma.13 g249400 revealed significantly higher expression in higher PH genotypes, and the gene existed 6 differences in the amino acids encoding between the two parents.

CONCLUSIONS

Data presented here provide support for Glyma.13 g249400 as a possible candidate genes for higher PH in wild soybean line N24852.

Genetic dissection of QTLs and differentiation analysis of alleles for heading date genes in rice

Heading date is an important agronomic trait in rice (Oryza sativa L.); it determines the geographical and seasonal adaptability of the crop. Single segment substitution lines (SSSLs) have become the preferred experimental materials in mapping functional genetic variations as the particular chromosome segments from donor genotypes can be evaluated for their impact on the phenotype in a recurrent recipient background. The phenotypic differences can be attributed to the control of quantitative trait loci (QTLs). Here, we evaluated a library consisting of 1,123 SSSLs in the same genetic background of an elite rice variety, Huajingxian74 (HJX74), and revealed four SSSLs, W05-1-11-2-7-6 (W05), W08-16-3-2 (W08), W12-28-58-03-19-1 (W12), and W22-9-5-2-4-9-3 (W22), which had a significantly different heading date compared to HJX74. To further genetically dissect the QTLs controlling heading date on chromosomes 3, 6, and 10, four SSSLs were used to develop 15 secondary SSSLs with the smaller substituted segments. The qHD-3 heading date QTL detected in W05 and W08 was delimited to an interval of 4.15 cM, whereas qHD-6-1 and qHD-6-2 heading date QTLs dissected from the substituted segments in W12 were mapped to the intervals of 2.25-cM and 2.55-cM, respectively. The qHD-10 QTL detected on the substituted segment in W22 was mapped to an interval of 6.85-cM. The nucleotide and amino acid sequence changes for those genes in the secondary SSSLs were also revealed. The allele variations of those genes might contribute to the heading date QTLs on chromosome 3 (DTH3, OsDof12, and EHD4), chromosome 6 (Hd3a, Hd17, and RFT1), and chromosome 10 (Ehd1 and Ehd2). These sequence variations in heading date genes would be useful resources for further studying the function of genes, and would be important for rice breeding. Overall, our results indicate that secondary SSSLs were powerful tools for genetic dissection of QTLs and identification of differentiation in the genes.

Mapping and verification of grain shape QTLs based on an advanced backcross population in rice

Grain shape is a key factor that influences both the appearance quality and grain yield of rice. To clarify the genetic basis of grain shape, an advanced backcross population was developed from the cross of a slender indica variety Jin23B and a round japonica variety QingGuAi, and a total of 10 quantitative trait loci (QTLs) for grain shape were detected over 2 years. Three QTLs, qGW1, qGS3 and qGS7 have large effects on grain shape and were detected in both years. To further validate their effects, the BC4F2 populations of the three QTLs were constructed. The alleles from QingGuAi of qGW1 and qGS7 both greatly increase grain width and the alleles from Jin23B of qGS3 greatly increase grain length. The allele variations of the three QTLs lead to the totally different grain shape of the two rice varieties. Our study provides insights into the genetic bases of grain shape and will promote the improvements of grain quality and yield in rice.

Detection of quantitative trait loci controlling grain zinc concentration using Australian wild rice, Oryza meridionalis, a potential genetic resource for biofortification of rice

Zinc (Zn) is one of the essential mineral elements for both plants and humans. Zn deficiency in human is one of the major causes of hidden hunger, a serious health problem observed in many developing countries. Therefore, increasing Zn concentration in edible part is an important issue for improving human Zn nutrition. Here, we found that an Australian wild rice O. meridionalis showed higher grain Zn concentrations compared with cultivated and other wild rice species. The quantitative trait loci (QTL) analysis was then performed to identify the genomic regions controlling grain Zn levels using backcross recombinant inbred lines derived from O. sativa 'Nipponbare' and O. meridionalis W1627. Four QTLs responsible for high grain Zn were detected on chromosomes 2, 9, and 10. The QTL on the chromosome 9 (named qGZn9), which showed the largest effect on grain Zn concentration was confirmed with the introgression line, which had a W1627 chromosomal segment covering the qGZn9 region in the genetic background of O. sativa 'Nipponbare'. Fine mapping of this QTL resulted in identification of two tightly linked loci, qGZn9a and qGZn9b. The candidate regions of qGZn9a and qGZn9b were estimated to be 190 and 950 kb, respectively. Furthermore, we also found that plants having a wild chromosomal segment covering qGZn9a, but not qGZn9b, is associated with fertility reduction. qGZn9b, therefore, provides a valuable allele for breeding rice with high Zn in the grains.

Genetic Loci Governing Grain Yield and Root Development under Variable Rice Cultivation Conditions

Drought is the major abiotic stress to rice grain yield under unpredictable changing climatic scenarios. The widely grown, high yielding but drought susceptible rice varieties need to be improved by unraveling the genomic regions controlling traits enhancing drought tolerance. The present study was conducted with the aim to identify quantitative trait loci (QTLs) for grain yield and root development traits under irrigated non-stress and reproductive-stage drought stress in both lowland and upland situations. A mapping population consisting of 480 lines derived from a cross between Dular (drought-tolerant) and IR64-21 (drought susceptible) was used. QTL analysis revealed three major consistent-effect QTLs for grain yield (qDTY1.1, qDTY1.3 , and qDTY8.1 ) under non-stress and reproductive-stage drought stress conditions, and 2 QTLs for root traits (qRT9.1 for root-growth angle and qRT5.1 for multiple root traits, i.e., seedling-stage root length, root dry weight and crown root number). The genetic locus qDTY1.1 was identified as hotspot for grain yield and yield-related agronomic and root traits. The study identified significant positive correlations among numbers of crown roots and mesocotyl length at the seedling stage and root length and root dry weight at depth at later stages with grain yield and yield-related traits. Under reproductive stage drought stress, the grain yield advantage of the lines with QTLs ranged from 24.1 to 108.9% under upland and 3.0-22.7% under lowland conditions over the lines without QTLs. The lines with QTL combinations qDTY1.3 +qDTY8.1 showed the highest mean grain yield advantage followed by lines having qDTY1.1 +qDTY8.1 and qDTY1.1 +qDTY8.1 +qDTY1.3 , across upland/lowland reproductive-stage drought stress. The identified QTLs for root traits, mesocotyl length, grain yield and yield-related traits can be immediately deployed in marker-assisted breeding to develop drought tolerant high yielding rice varieties.

Estimation of loci involved in non-shattering of seeds in early rice domestication

Rice (Oryza sativa L.) is widely cultivated around the world and is known to be domesticated from its wild form, O. rufipogon. A loss of seed shattering is one of the most obvious phenotypic changes selected for during rice domestication. Previously, three seed-shattering loci, qSH1, sh4, and qSH3 were reported to be involved in non-shattering of seeds of Japonica-type cultivated rice, O. sativa cv. Nipponbare. In this study, we focused on non-shattering characteristics of O. sativa Indica cv. IR36 having functional allele at qSH1. We produced backcross recombinant inbred lines having chromosomal segments from IR36 in the genetic background of wild rice, O. rufipogon W630. Histological and quantitative trait loci analyses of abscission layer formation were conducted. In the analysis of quantitative trait loci, a strong peak was observed close to sh4. We, nevertheless, found that some lines showed complete abscission layer formation despite carrying the IR36 allele at sh4, implying that non-shattering of seeds of IR36 could be regulated by the combination of mutations at sh4 and other seed-shattering loci. We also genotyped qSH3, a recently identified seed-shattering locus. Lines that have the IR36 alleles at sh4 and qSH3 showed inhibition of abscission layer formation but the degree of seed shattering was different from that of IR36. On the basis of these results, we estimated that non-shattering of seeds in early rice domestication involved mutations in at least three loci, and these genetic materials produced in this study may help to identify novel seed-shattering loci.

Genome-Wide Association Mapping Reveals Multiple QTLs Governing Tolerance Response for Seedling Stage Chilling Stress in Indica Rice

Rice crop is sensitive to cold stress at seedling stage. A panel of population representing 304 shortlisted germplasm lines was studied for seedling stage chilling tolerance in indica rice. Six phenotypic classes were exposed to six low temperature stress regimes under control phenotyping facility to investigate response pattern. A panel of 66 genotypes representing all phenotypic classes was used for ensuring genetic diversity, population structure and association mapping for the trait using 58 simple sequence repeat (SSR) and 2 direct trait linked markers. A moderate level of genetic diversity was detected in the panel population for the trait. Deviation of Hardy-Weinberg's expectation was detected in the studied population using Wright's F statistic. The panel showed 30% variation among population and 70% among individuals. The entire population was categorized into three sub-populations through STRUCTURE analysis. This revealed tolerance for the trait had a common primary ancestor for each sub-population with few admix individuals. The panel population showed the presence of many QTLs for cold stress tolerance in the individuals representing like genome-wide expression of the trait. Nineteen SSR markers were significantly associated at chilling stress of 8°C to 4°C for 7-21 days duration. Thus, the primers linked to the seedling stage cold tolerance QTLs namely qCTS9, qCTS-2, qCTS6.1, qSCT2, qSCT11, qSCT1a, qCTS-3.1, qCTS11.1, qCTS12.1, qCTS-1b, and CTB2 need to be pyramided for development of strongly chilling tolerant variety.

Integrated RNA Sequencing and QTL Mapping to Identify Candidate Genes from Oryza rufipogon Associated with Salt Tolerance at the Seedling Stage

Soil salinity is a common abiotic stress affecting crop productivity. To identify favorable alleles from wild rice (Oryza rufipogon Griff.) that enhance salinity tolerance of rice (O. sativa L.), a set of introgression lines (ILs) were developed. The ILs were derived from an O. rufipogon accession collected from Chaling (Hunan Province, China) as the donor, and a widely grown O. sativa indica cultivar 93-11 as the recipient. Through evaluating the salt tolerance of 285 ILs at the seedling stage, a total of 10 quantitative trait loci (QTLs) related to salt tolerance were identified on chromosomes 1, 5, 7 and 9-12, with individual QTLs explaining 2-8% of phenotypic variance. The O. rufipogon-derived alleles at four QTLs improved salt tolerance in the 93-11 background. At the same time, a salt-tolerant IL, 9L136, was identified and characterized. Compared with the recipient parent 93-11, a total of 1,391 differentially expressed genes (DEGs) were detected specifically in 9L136 between salt stress and normal condition through genome-wide expression analysis. Of these, four DEGs located in the QTL regions carried by 9L136, suggesting that the four genes might be candidates associated with salt tolerance. Both the highly salt-tolerant ILs and the favorable O. rufipogon-derived QTLs identified in the present study will provide new genetic resources for improving the resistance of cultivated rice against salinity stress using molecular breeding strategies in the future.

Development of genome-wide PCR-based markers from insertion, deletion and single nucleotide polymorphisms for closely related Japanese rice cultivars and identification of QTLs for the appearance of cooked rice and polished rice

Appearance of rice grain is an important property, affecting its acceptance by consumers. Moreover, appearance is a complex characteristic involving many components, including glossiness and whiteness. The genetic bases for the glossiness of cooked rice and the whiteness of polished rice (WPR) were determined using 133 recombinant inbred lines (RILs) derived from a cross between two closely related cultivars from Hokkaido, Joiku462, with high glossiness and whiteness, and Yukihikari, an ancestor of Joiku462 with low glossiness and whiteness. Analyses identified 167 genome-wide InDel markers, five cleaved amplified polymorphic sequences (CAPS) and eight derived CAPS markers differentiating the parental lines. The glossiness area (GLA) and glossiness strength (GLS) of cooked rice and WPR were determined for RILs in two locations, Pippu and Sapporo, Hokkaido. Four QTLs were detected. qGLA10 and qGLS9 were detected on chromosomes 10 and 9, respectively, with both being significant at both geographic locations. qWPR1 on chromosome 1 was significant at Pippu, and qWPR4 on chromosome 4 was significant at Sapporo. The Joiku462 alleles at all QTLs increased each trait. The PCR-based markers flanking these four QTLs may be useful for improvement of GLA, GLS and WPR.

Phylogenetic and Diversity Analysis of Dactylis glomerata Subspecies Using SSR and IT-ISJ Markers

The genus Dactylis, an important forage crop, has a wide geographical distribution in temperate regions. While this genus is thought to include a single species, Dactylis glomerata, this species encompasses many subspecies whose relationships have not been fully characterized. In this study, the genetic diversity and phylogenetic relationships of nine representative Dactylis subspecies were examined using SSR and IT-ISJ markers. In total, 21 pairs of SSR primers and 15 pairs of IT-ISJ primers were used to amplify 295 polymorphic bands with polymorphic rates of 100%. The average polymorphic information contents (PICs) of SSR and IT-ISJ markers were 0.909 and 0.780, respectively. The combined data of the two markers indicated a high level of genetic diversity among the nine D. glomerata subspecies, with a Nei’s gene diversity index value of 0.283 and Shannon’s diversity of 0.448. Preliminarily phylogenetic analysis results revealed that the 20 accessions could be divided into three groups (A, B, C). Furthermore, they could be divided into five clusters, which is similar to the structure analysis with K = 5. Phylogenetic placement in these three groups may be related to the distribution ranges and the climate types of the subspecies in each group. Group A contained eight accessions of four subspecies, originating from the west Mediterranean, while Group B contained seven accessions of three subspecies, originating from the east Mediterranean.

Influence of ethnic traditional cultures on genetic diversity of rice landraces under on-farm conservation in southwest China

BACKGROUND

Crop genetic resources are important components of biodiversity. However, with the large-scale promotion of mono-cropping, genetic diversity has largely been lost. Ex-situ conservation approaches were widely used to protect traditional crop varieties worldwide. However, this method fails to maintain the dynamic evolutionary processes of crop genetic resources in their original habitats, leading to genetic diversity reduction and even loss of the capacity of resistance to new diseases and pests. Therefore, on-farm conservation has been considered a crucial complement to ex-situ conservation. This study aimed at clarifying the genetic diversity differences between ex-situ conservation and on-farm conservation and to exploring the influence of traditional cultures on genetic diversity of rice landraces under on-farm conservation.

METHODS

The conservation status of rice landrace varieties, including Indica and Japonica, non-glutinous rice (Oryza sativa) and glutinous rice (Oryza sativa var. glutinosa Matsum), was obtained through ethno-biology investigation method in 12 villages of ethnic groups from Guizhou, Yunnan and Guangxi provinces of China. The genetic diversity between 24 pairs of the same rice landraces from different times were compared using simple sequence repeat (SSR) molecular markers technology. The landrace paris studied were collected in 1980 and maintained ex-situ, while 2014 samples were collected on-farm in southwest of China.

RESULTS

The results showed that many varieties of rice landraces have been preserved on-farm by local farmers for hundreds or thousands of years. The number of alleles (Na), effective number of alleles (Ne), Nei genetic diversity index (He) and Shannon information index (I) of rice landraces were significantly higher by 12.3-30.4 % under on-farm conservation than under ex-situ conservation. Compared with the ex-situ conservation approach, rice landraces under on-farm conservation programs had more alleles and higher genetic diversity. In every site we investigated, ethnic traditional cultures play a positive influence on rice landrace variety diversity and genetic diversity.

CONCLUSION

Most China's rice landraces were conserved in the ethnic areas of southwest China. On-farm conservation can effectively promote the allelic variation and increase the genetic diversity of rice landraces over the past 35 years. Moreover, ethnic traditional culture practices are a crucial foundation to increase genetic diversity of rice landraces and implement on-farm conservation.

Association Mapping for Aluminum Tolerance in a Core Collection of Rice Landraces

Trivalent aluminum (Al³⁺) has drastic effect on the rice production in acidic soils. Elite genes for aluminum (Al) tolerance might exist in rice landraces. Therefore, the purpose of this research is to mine the elite genes within rice landraces. Association mapping for Al tolerance traits [i.e., relative root elongation (RRE)] was performed by using a core collection of 150 accessions of rice landraces (i.e., Ting's rice core collection). Our results showed that the Ting's rice core collection possessed a wide-range of phenotypic variation for Al tolerance, and the index of Al tolerance (RRE) was ranged from 0.22 to 0.89. Moreover, the groups with different origins and compositions of indica and japonica rice showed different degrees of tolerance to varying levels of Al. These rice landraces were further screened with 274 simple sequence repeat markers, and association mapping was performed using a mixed linear model approach. The mapping results showed that a total of 23 significant (P < 0.05) trait-marker associations were detected for Al tolerance. Of these, three associations (13%) were identical to the quantitative trait loci reported previously, and other 20 associations were reported for the first time in this study. The proportion of phenotypic variance (R²) explained by 23 significant associations ranged from 5.03 to 20.03% for Al tolerance. We detected several elite alleles for Al tolerance based on multiple comparisons of allelic effects, which could be used to develop Al tolerant rice cultivars through marker-assisted breeding.

Identification and linkage analysis of a new rice bacterial blight resistance gene from XM14, a mutant line from IR24

Bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo) is a chief factor limiting rice productivity worldwide. XM14, a rice mutant line resistant to Xoo, has been obtained by treating IR24, which is susceptible to six Philippine Xoo races and six Japanese Xoo races, with N-methyl-N-nitrosourea. XM14 showed resistance to six Japanese Xoo races. The F₂ population from XM14 × IR24 clearly showed 1 resistant : 3 susceptible segregation, suggesting control of resistance by a recessive gene. The approximate chromosomal location of the resistance gene was determined using 10 plants with shortest lesion length in the F₂ population from XM14 × Koshihikari, which is susceptible to Japanese Xoo races. DNA marker-assisted analysis revealed that the gene was located on chromosome 3. IAS16 line carries IR24 genetic background with a Japonica cultivar Asominori segment of chromosome 3, on which the resistance gene locus was thought to be located. The F₂ population from IAS16 × XM14 showed a discrete distribution. Linkage analysis indicated that the gene is located around the centromeric region. The resistance gene in XM14 was a new gene, named XA42. This gene is expected to be useful for resistance breeding programs and for genetic analysis of Xoo resistance.

Development and characterization of chromosome segment substitution lines derived from Oryza rufipogon in the genetic background of O. sativa spp. indica cultivar 9311

Background

Wild rice (Oryza rufipogon) constitutes a primary gene source for rice breed improvement. Chromosome segment substitution line (CSSL) for O. rufipogon is a powerful tool for fine mapping of quantitative traits, new gene discovery, and marker-assisted breeding. Thus, they provide a basis for a wide range of genomic and genetic studies.

Results

In this study, a set of 198 CSSLs were developed from a cross between recurrent parent indica var. 9311 and an O. rufipogon donor parent; these were then genotyped using 313 polymorphic SSR markers evenly distributed across the 12 rice chromosomes. On average, each CSSL carried 2.16 introgressed segments, and the genetic distance of each segment was about 6 cM. The segments collectively covered 84.9 % of the wild rice genome. Based on these CSSLs, 25 QTLs involved in 10 agronomic traits were identified. Seven CSSLs were subjected to a whole-genome single nucleotide polymorphism chip assay and two QTLs, qSH4-1 and qDTH10-1, detected. In addition, a new QTL associated with the heading date was detected in a 78-Kb region on chromosome 10, thus proving the ability of these CSSLs to identify new QTLs and genes.

Conclusions

The newly developed CSSL population proved a useful tool for both gene identification and whole-genome research of wild rice. These CSSL materials will provide a foundation for rice variety improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2987-5) contains supplementary material, which is available to authorized users.

Population Structure, Genetic Diversity and Molecular Marker-Trait Association Analysis for High Temperature Stress Tolerance in Rice

Rice exhibits enormous genetic diversity, population structure and molecular marker-traits associated with abiotic stress tolerance to high temperature stress. A set of breeding lines and landraces representing 240 germplasm lines were studied. Based on spikelet fertility percent under high temperature, tolerant genotypes were broadly classified into four classes. Genetic diversity indicated a moderate level of genetic base of the population for the trait studied. Wright’s F statistic estimates showed a deviation of Hardy-Weinberg expectation in the population. The analysis of molecular variance revealed 25 percent variation between population, 61 percent among individuals and 14 percent within individuals in the set. The STRUCTURE analysis categorized the entire population into three sub-populations and suggested that most of the landraces in each sub-population had a common primary ancestor with few admix individuals. The composition of materials in the panel showed the presence of many QTLs representing the entire genome for the expression of tolerance. The strongly associated marker RM547 tagged with spikelet fertility under stress and the markers like RM228, RM205, RM247, RM242, INDEL3 and RM314 indirectly controlling the high temperature stress tolerance were detected through both mixed linear model and general linear model TASSEL analysis. These markers can be deployed as a resource for marker-assisted breeding program of high temperature stress tolerance.

Development and validation of cross-transferable and polymorphic DNA markers for detecting alien genome introgression in Oryza sativa from Oryza brachyantha

African wild rice Oryza brachyantha (FF), a distant relative of cultivated rice Oryza sativa (AA), carries genes for pests and disease resistance. Molecular marker assisted alien gene introgression from this wild species to its domesticated counterpart is largely impeded due to the scarce availability of cross-transferable and polymorphic molecular markers that can clearly distinguish these two species. Availability of the whole genome sequence (WGS) of both the species provides a unique opportunity to develop markers, which are cross-transferable. We observed poor cross-transferability (~0.75 %) of O. sativa specific sequence tagged microsatellite (STMS) markers to O. brachyantha. By utilizing the genome sequence information, we developed a set of 45 low cost PCR based co-dominant polymorphic markers (STS and CAPS). These markers were found cross-transferrable (84.78 %) between the two species and could distinguish them from each other and thus allowed tracing alien genome introgression. Finally, we validated a Monosomic Alien Addition Line (MAAL) carrying chromosome 1 of O. brachyantha in O. sativa background using these markers, as a proof of concept. Hence, in this study, we have identified a set molecular marker (comprising of STMS, STS and CAPS) that are capable of detecting alien genome introgression from O. brachyantha to O. sativa.

Overcoming inter-subspecific hybrid sterility in rice by developing indica-compatible japonica lines

Rice (Oryza sativa L.) is an important staple crop. The exploitation of the great heterosis that exists in the inter-subspecific crosses between the indica and japonica rice has long been considered as a promising way to increase the yield potential. However, the male and female sterility frequently occurred in the inter-subspecific hybrids hampered the utilization of the heterosis. Here we report that the inter-subspecific hybrid sterility in rice is mainly affected by the genes at Sb, Sc, Sd and Se loci for F1 male sterility and the gene at S5 locus for F1 female sterility. The indica-compatible japonica lines (ICJLs) developed by pyramiding the indica allele (S-i) at Sb, Sc, Sd and Se loci and the neutral allele (S-n) at S5 locus in japonica genetic background through marker-assisted selection are compatible with indica rice in pollen fertility and in spikelet fertility. These results showed a great promise of overcoming the inter-subspecific hybrid sterility and exploiting the heterosis by developing ICJLs.

An efficient protocol for isolation of inhibitor-free nucleic acids even from recalcitrant plants

For fast and easy isolation of inhibitor-free genomic DNA even from the toughest plant leaf samples, including those high in polyphenols and polysaccharides, a protocol has been developed. To prevent the solubility of polysaccharides in the DNA extract, high salt concentration (1.4 M) was used in the extraction buffer. Polyvinylpyrrolidone (PVP) was used for the removal of polyphenols as polymerase chain reaction (PCR) inhibitors. Proteins like various enzymes were degraded by proteinase K and removed by centrifugation from plant extracts during the isolation process resulting in pure DNA and RNA ready to use in downstream applications including PCR, quantitative polymerase chain reaction (qPCR), ligation, restriction and sequencing. This protocol yielded a high molecular weight DNA and RNA isolated from leaves and roots of recalcitrant plants which was free from contamination and color. The average yields of total RNA from roots and shoot of Betula and Grape ranged from 285 to 364 ng/µl with A260/A280 between 1.9 and 2.08. The RNA isolated with this protocol was verified to be suitable for PCR, quantitative real-time PCR, semi-quantitative reverse transcription polymerase chain reaction, cDNA synthesis and expression analysis. This protocol shown here is reproducible and can be used for a broad spectrum of plant species which have polyphenols and polysaccharide compounds.

Lowered Diversity and Increased Inbreeding Depression within Peripheral Populations of Wild Rice Oryza rufipogon

BACKGROUND

The distribution of genetic variability from the interior towards the periphery of a species' range is of great interest to evolutionary biologists. Although it has been long presumed that population genetic variation should decrease as a species' range is approached, results of empirical investigations still remain ambiguous. Knowledge regarding patterns of genetic variability as well as affected factors is particularly not conclusive in plants.

METHODOLOGY/PRINCIPAL FINDINGS

To determine genetic divergence in peripheral populations of the wild rice Oryza rufipogon Griff. from China, genetic diversity and population structure were studied in five northern & northeastern peripheral and 16 central populations using six microsatellite loci. We found that populations resided at peripheries of the species possessed markedly decreased microsatellite diversity than those located in its center. Population size was observed to be positively correlated with microsatellite diversity. Moreover, there are significantly positive correlations between levels of microsatellite diversity and distances from the northern and northeastern periphery of this species. To investigate genetic structure and heterozygosity variation between generations of O. rufipogon, a total of 2382 progeny seeds from 186 maternal families were further assayed from three peripheral and central populations, respectively. Peripheral populations exhibited significantly lower levels of heterozygosities than central populations for both seed and maternal generations. In comparisons with maternal samples, significantly low observed heterozygosity (HO) and high heterozygote deficit within populations (FIS) values were detected in seed samples from both peripheral and central populations. Significantly lower observed heterozygosity (HO) and higher FIS values were further observed in peripheral populations than those in central populations for seed samples. The results indicate an excess of homozygotes and thus high inbreeding depression in peripheral populations.

CONCLUSIONS/SIGNIFICANCE

Our results together suggest that historical contraction of geographical range, demographic changes, and environmental conditions near the northern and northeastern margins of O. rufipogon favor inbreeding and possibly selfing, leading to the rapidly decreased effective population size. Genetic drift, reduced gene flow, and possible local selection, consequently lead to lowered gene diversity, accelerated genetic divergence and increased inbreeding depression found in peripheral populations of O. rufipogon. Given these characteristics observed, northern and northeastern peripheral populations deserve relatively different conservation strategies for either germplasm sampling of ex situ conservation or setting in situ reserves for the adaptation to possible environmental changes and the future germplasm utilization of wild rice.

QTL Location and Epistatic Effect Analysis of 100-Seed Weight Using Wild Soybean (Glycine soja Sieb. & Zucc.) Chromosome Segment Substitution Lines

Increasing the yield of soybean (Glycine max L. Merrill) is a main aim of soybean breeding. The 100-seed weight is a critical factor for soybean yield. To facilitate genetic analysis of quantitative traits and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population consisting of 194 chromosome segment substitution lines (CSSLs) was developed. In these lines, different chromosomal segments of the Chinese cultivar Suinong 14 were substituted into the genetic background of wild soybean (Glycine soja Sieb. & Zucc.) ZYD00006. Based on these CSSLs, a genetic map covering the full genome was generated using 121 simple sequence repeat (SSR) markers. In the quantitative trait loci (QTL) analysis, twelve main effect QTLs (qSW-B1-1/2/3, qSW-D1b-1/2, qSW-D2-1/2, qSW-G-1/2/3, qSW-M-2 and qSW-N-2) underlying 100-seed weight were identified in 2011 and 2012. The epistatic effects of pairwise interactions between markers were analyzed in 2011 and 2012. The results clearly demonstrated that these CSSLs could be used to identify QTLs, and that an epistatic analysis was able to detect several sites with important epistatic effects on 100-seed weight. Thus, we identified loci that will be valuable for improving soybean 100-seed weight. These results provide a valuable foundation for identifying the precise location of genes of interest, and for designing cloning and marker-assisted selection breeding strategies targeting the 100-seed weight of soybean.

Marker assisted pyramiding of drought yield QTLs into a popular Malaysian rice cultivar, MR219

Background

Three drought yield QTLs, qDTY_2.2, qDTY_3.1, and qDTY_12.1 with consistent effect on grain yield under reproductive stage drought stress were pyramided through marker assisted breeding with the objective of improving the grain yield of the elite Malaysian rice cultivar MR219 under reproductive stage drought stress. Foreground selection using QTL specific markers, recombinant selection using flanking markers, and background selection were performed. BC₁F₃-derived lines with different combinations of qDTY_2.2, qDTY_3.1, and qDTY_12.1 were evaluated under both reproductive stage drought stress and non-stress during the dry seasons of 2013 and 2014 at IRRI.

Results

The grain yield reductions in the stress trials compared to non-stress trials ranged from 79 to 93 %. In the stress trials, delay in days to flowering and reduction in plant height were observed. In both seasons, MR219 did not produce any yield under stress, however it produced a yield of 5917 kg ha⁻¹ during the 2013 dry season and 8319 kg ha⁻¹ during the 2014 dry season under non-stress. Selected introgressed lines gave a yield advantage of 903 to 2500 kg ha⁻¹ over MR219 under reproductive stage drought stress and a yield of more than 6900 kg ha⁻¹ under non-stress during the 2014 dry season. Among lines with single qDTY, lines carrying qDTY_2.2 provided a higher yield advantage under reproductive stage drought stress in the MR219 background. Two-qDTY combinations (qDTY_3.1+qDTY_2.2 and qDTY_3.1+qDTY_12.1) performed better than lines with three qDTY combinations, indicating the absence of positive interactions between the three qDTYs.

Conclusion

We successfully developed drought-tolerant MR219 pyramided lines with a yield advantage of more than 1500 kg ha⁻¹. Differential yield advantages of different combinations of the qDTYs indicate a differential synergistic relationship among qDTYs. This is the first report on the successful effect of qDTYs in increasing the yield under drought in genetic backgrounds other than those in which the qDTYs were earlier identified.

Inconsistent diversities between nuclear and plastid genomes of AA genome species in the genus Oryza

AA genome species in the genus Oryza are valuable resources for improvement of cultivated rice. Oryza rufipogon and O. barthii were progenitors of two domesticated rice species, O. sativa and O. glaberrima, respectively. We used chloroplast single-nucleotide repeats (RCt1-10) to evaluate genetic diversity among AA genome species. Higher diversity was detected in the American species O. glumaepatula and the Asian species O. rufipogon. Other chloroplast sequences indicated that O. glumaepatula shares high similarity with O. longistaminata. Insertions of retrotransposable elements, however, showed a close relation between O. barthii and O. glumaepatula. To clarify phylogenetic relationships among AA genomes, whole-genome sequences obtained from different species were used to develop chloroplast INDEL markers. The INDEL patterns clearly showed multiple maternal origins of O. glumaepatula. The complicated origins have resulted in high genetic diversity in this species. In contrast, the Australian endemic species O. meridionalis tended to show narrower diversity than the other species. High variation in O. rufipogon, reconfirmed using the chloroplast INDELs, covered the variation in O. meridionalis and part of the variation in O. glumaepatula. Maternal lineages including O. barthii, O. longistaminata and the remainder of O. glumaepatula were phylogenetically close to each other and carried low genetic diversity. They were separated from independent lineages, suggesting that they had diverged from a single ancestral maternal lineage, but diverged later to keep gene flow within respective species, as SSR compositions suggested. Genetic relationships among AA genome species indicate how these species have evolved and become distributed across four continents.

Evaluation of genetic variation among wild rice populations in Cambodia

A total of 448 samples in five natural populations of wild rice (Oryza rufipogon) were collected in Cambodia. They were examined using 12 SSR and two chloroplast markers to evaluate the degree of variation among populations and the genetic structure within populations. In the two annual populations, the number of plants with homozygous alleles at all 12 SSR loci were high (66.3% and 79.5%), suggesting that these plants propagate mainly through self-pollination. In the three perennial populations, no individuals had all homozygous genotypes, but redundant genotypes resulted from clonal propagation were observed. Percentages of the redundant genotypes were highly varied (3.6%, 29.2% and 86.0%). This may be due to the different stable levels of environmental conditions. As for chloroplast genome, most of the wild plants showed the same chloroplast types as most Indica-type cultivars have. However, plants with different chloroplast types were maintained, even in the same population. In tropical Asian countries, many wild rice populations were observed under similar ecological conditions examined in this study. Therefore, the present results concerning population structure will be important to further elucidate genetic features of wild rice, and will also give strong clues to utilize and conserve wild natural genetic resources.

Action of multiple intra-QTL genes concerted around a co-localized transcription factor underpins a large effect QTL

Sub-QTLs and multiple intra-QTL genes are hypothesized to underpin large-effect QTLs. Known QTLs over gene families, biosynthetic pathways or certain traits represent functional gene-clusters of genes of the same gene ontology (GO). Gene-clusters containing genes of different GO have not been elaborated, except in silico as coexpressed genes within QTLs. Here we demonstrate the requirement of multiple intra-QTL genes for the full impact of QTL qDTY12.1 on rice yield under drought. Multiple evidences are presented for the need of the transcription factor 'no apical meristem' (OsNAM12.1) and its co-localized target genes of separate GO categories for qDTY12.1 function, raising a regulon-like model of genetic architecture. The molecular underpinnings of qDTY12.1 support its effectiveness in further improving a drought tolerant genotype and for its validity in multiple genotypes/ecosystems/environments. Resolving the combinatorial value of OsNAM12.1 with individual intra-QTL genes notwithstanding, identification and analyses of qDTY12.1has fast-tracked rice improvement towards food security.

Population Structure of and Conservation Strategies for Wild Pyrus ussuriensis Maxim. in China

Pyrus ussriensis Maxim. is native to the northern part of China, but whose habitats are currently being destroyed by environmental changes and human deforestation. An investigation of population structure and genetic diversity of wild Ussurian pear is a priority in order to acquire fundamental knowledge for conservation. A total of 153 individuals of wild Ussurian pear from the main habitats, Heilongjiang, Jilin, and Inner Mongolia in China, possessed low genetic diversity as a result of habitat fragmentation. The genetic diversity of the populations in Inner Mongolia and north east of Heilongjiang was especially low and there was the possibility of inbreeding. Wild Ussurian pears were divided into 5 groups based on the Bayesian clustering method using 20 nuclear SSRs (nSSRs) and 5 groups by haplotype distributions using 16 chloroplast SSRs (cpSSRs), and the populations in Inner Mongolia and north east of Heilongjiang represented unique genotypes. AMOVA indicated there was a 20.05% variation in nSSRs and a 44.40% variation in cpSSRs among populations. These values are relatively high when compared to those of other tree species. Haplotype E, positioned in the center of the cpSSR analysis network and showed the largest number of connections with other haplotypes, represented the most important haplotype. Inner Mongolia and the north east of Heilongjiang are two areas that need urgent conservation because of their genetic vulnerability and peculiarity. We determined 4 conservation units based on the clustering by nSSRs and cpSSRs, and geographic factor. This information is helpful in deciding the conservation strategies for wild Ussurian pear in China.