Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.)

BACKGROUND

This study aims to decipher the genetic basis governing yield components and quality attributes of peanuts, a critical aspect for advancing molecular breeding techniques. Integrating genotype re-sequencing and phenotypic evaluations of seven yield components and two grain quality traits across four distinct environments allowed for the execution of a genome-wide association study (GWAS).

RESULTS

The nine phenotypic traits were all continuous and followed a normal distribution. The broad heritability ranged from 88.09 to 98.08%, and the genotype-environment interaction effects were all significant. There was a highly significant negative correlation between protein content (PC) and oil content (OC). The 10× genome re-sequencing of 199 peanut accessions yielded a total of 631,988 high-quality single nucleotide polymorphisms (SNPs), with 374 significant SNP loci identified in association with the nine traits of interest. Notably, 66 of these pertinent SNPs were detected in multiple environments, and 48 of them were linked to multiple traits of interest. Five loci situated on chromosome 16 (Chr16) exhibited pleiotropic effects on yield traits, accounting for 17.64-32.61% of the observed phenotypic variation. Two loci on Chr08 were found to be strongly associated with protein and oil contents, accounting for 12.86% and 14.06% of their respective phenotypic variations, respectively. Linkage disequilibrium (LD) block analysis of these seven loci unraveled five nonsynonymous variants, leading to the identification of one yield-related candidate gene and two quality-related candidate genes. The correlation between phenotypic variation and SNP loci in these candidate genes was validated by Kompetitive allele-specific PCR (KASP) marker analysis.

CONCLUSIONS

Overall, molecular markers were developed for genetic loci associated with yield and quality traits through a GWAS investigation of 199 peanut accessions across four distinct environments. These molecular tools can aid in the development of desirable peanut germplasm with an equilibrium of yield and quality through marker-assisted breeding.

The pear genomics database (PGDB): a comprehensive multi-omics research platform for Pyrus spp

BACKGROUND

Pears are among the most important temperate fruit trees in the world, with significant research efforts increasing over the last years. However, available omics data for pear cannot be easily and quickly retrieved to enable further studies using these biological data.

DESCRIPTION

Here, we present a publicly accessible multi-omics pear resource platform, the Pear Genomics Database (PGDB). We collected and collated data on genomic sequences, genome structure, functional annotation, transcription factor predictions, comparative genomics, and transcriptomics. We provide user-friendly functional modules to facilitate querying, browsing and usage of these data. The platform also includes basic and useful tools, including JBrowse, BLAST, phylogenetic tree building, and additional resources providing the possibility for bulk data download and quick usage guide services.

CONCLUSIONS

The Pear Genomics Database (PGDB, http://pyrusgdb.sdau.edu.cn ) is an online data analysis and query resource that integrates comprehensive multi-omics data for pear. This database is equipped with user-friendly interactive functional modules and data visualization tools, and constitutes a convenient platform for integrated research on pear.

Discovery of genetic loci and causal genes for seed germination via deep re-sequencing in soybean

High seed germination is crucial for mechanical sowing, seedling establishment, growth potential, multiple resistances, and the formation of yield and quality. However, few genetic loci and candidate genes conferring seed germination were explored in soybean at present. In view of this, a natural population containing 199 accessions was assessed for the germination potential (GP) and germination rate (GR) and also was re-sequenced at the average sequencing depth of 18.4 × for each accession. In total, 5,665,469 SNPs were obtained for association analysis, and 470 SNPs in 55 loci on 18 chromosomes were identified to associate with seed germination. Of them, 85 SNPs on chromosomes 1, 10, and 14 were associated with mean value and BLUP value for GP and GR, simultaneously. Moreover, 324 SNPs (68.9% of the total) in four loci were located on chromosome 14 for seed germination, of which 11 SNPs were located in the exons, 30 in introns, 17 in 5'UTR or 3'UTR, and 46 in upstream or downstream. Based on these, 131 candidate genes flanking the associated SNPs were analyzed for gene annotation, SNP mutation, and RNA expression, and three causal genes, Glyma.14G069800 (RNA-binding protein), Glyma.14G071400 (bZIP transcription factor), and Glyma.17G033200 (nucleic acid-binding protein), were screened out and might be responsible for the seed germination. The closely associated SNPs and causal genes provided an important resource and dissecting of genetic basis for seed germination improvement in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01316-6.

The mutated Bacillus amyloliquefaciens strain shows high resistance to Aeromonas hydrophila and Aeromonas veronii in grass carp

Bacillus amyloliquefaciens X030 (BaX030) has broad-spectrum antibacterial activity against the fish pathogens Aeromonas hydrophila and Aeromonas veronii. To improve its antibacterial effect, BaX030 was subjected to compound mutagenesis of atmospheric and room temperature plasma (ARTP) and nitrosoguanidine (NTG). The results showed that, compared with the original strain, the production of macrolactin A and oxydifficidin in mutated strain N-11 increased to 39 % and 268 %, respectively. The re-sequencing analysis suggested that there were SNPs and InDels in the gene clusters focused on the sucrose utilization pathway, glycolysis pathway and fatty acid synthesis pathway. Scanning electron microscopy revealed that strain N-11 became thin and long. The qRT-PCR results indicated that the expression of immune factors in the liver or kidney tissue of grass carp increased after feeding with N-11. H&E staining and protection experiments also showed that the mortality and surface symptoms of grass carp infected by the two pathogens were significantly reduced. The study identified a probiotic strain with potential application value in aquaculture production and provided a new strategy for the discovery of new strains with higher antibacterial biological activity.

Whole genome re-sequencing identifies unique adaption of single nucleotide polymorphism, insertion/deletion and structure variation related to hypoxia in Tibetan chickens

BACKGROUND

The adaptation to hypoxia in high altitude areas has great research value in the field of biological sciences. Tibetan chicken has unique adaptability to high-altitude, low pressure and anoxic conditions, and served as a biological model to search for genetic diversity of hypoxia adaption.

METHODS

The whole genome re-sequencing technology was conducted to investigate the genetic diversity.

RESULTS

In this study, we obtained quantity genetic resource, contained 5164926 single nucleotide polymorphisms (SNPs), 237504 Insertion/Deletion (InDel), 55606 structural variation types in all chromosomes of Tibetan chicken. Moreover, 17154 non-synonymous mutations, 45763 synonymous mutations, 258 InDel mutations and 9468 structural mutations were detected in coding sequencing (CDS) region. Furthermore, SNPs occur in 591 genes, including HIF1A, VEGF, MAPK 8/9/10/11, PPARA/D/G, NOTCH2, and ABCs, which were involved in 14 hypoxia-related pathways, such as VEGF signaling pathway, MAPK signaling pathway, PPAR signaling pathway and Notch signaling pathway. Among them, 19 genes with non-synonymous SNP variation in CDS were identified. Moreover, structure variation in CDS also occurred in the mentioned above genes with SNPs.

CONCLUSIONS

This study provides useful targets for clarifying the hypoxia adaptability of the domestication of chickens in Tibetan and may help breeding efforts to develop improved breeds for the highlands.

Genome-wide association studies detects candidate genes for wool traits by re-sequencing in Chinese fine-wool sheep

BACKGROUND

The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits.

RESULTS

A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263).

CONCLUSION

This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

Array-based dynamic allele specific hybridization (Array-DASH): Optimization-free microarray processing for multiple simultaneous genomic assays

We report proof-of-principle experiments regarding a dynamic microarray protocol enabling accurate and semi-quantitative DNA analysis for re-sequencing, fingerprinting and genotyping. Single-stranded target molecules hybridise to surface-bound probes during initial gradual cooling with high-fidelity. Real-time tracking of target denaturation (via fluorescence) during a 'dynamic' gradual heating phase permits 'melt-curve' analysis. The probe most closely matching the target sequence is identified based on the highest melting temperature. We demonstrated a >99% re-sequencing accuracy and a potential detection rate of 1% for SNPs. Experiments employing Hypericum ribosomal ITS regions and HIV genomes illustrated a reliable detection level of 5% plus simultaneous re-sequencing and genotyping. Such performance suggests a range of potential real-world applications involving rapid sequence interrogation, for example, in the Covid-19 pandemic. Guidance is offered towards the development of a commercial platform and dedicated software required to bring this technique into mainstream science.

Engineering of primary metabolic pathways for titer improvement of milbemycins in Streptomyces bingchenggensis

Milbemycins are used commercially as insect repellents and acaricides; however, their high cost remains a significant challenge to commercial production. Hence, improving the titer of milbemycins for commercial application is an urgent priority. The present study aimed to effectively increase the titer of milbemycins using a combination of genome re-sequencing and metabolic engineering. First, 133 mutation sites were identified by genome re-sequencing in the mutagenized high-yielding strain BC04. Among them, three modifiable candidate genes (sbi_04868 encoding citrate synthase, sbi_06921 and sbi_06922 encoding alpha and beta subunits of acetyl-CoA carboxylase, and sbi_04683 encoding carbon uptake system gluconate transporter) related to primary metabolism were screened and identified. Next, the DNase-deactivated Cpf1-based integrative CRISPRi system was used in S. bingchenggensis to downregulate the transcription level of gene sbi_04868. Then, overexpression of the potential targets sbi_06921-06922 and sbi_04683 further facilitated milbemycin biosynthesis. Finally, those candidate genes were engineered to produce strains with combinatorial downregulation and overexpression, which resulted in the titer of milbemycin A3/A4 increased by 27.6% to 3164.5 mg/L. Our research not only identified three genes in S. bingchenggensis that are closely related to the production of milbemycins, but also offered an efficient engineering strategy to improve the titer of milbemycins using genome re-sequencing. KEY POINTS: • We compared the genomes of two strains with different titers of milbemycins. • We found three genes belonging to primary metabolism influence milbemycin production. • We improved titer of milbemycins by a combinatorial engineering of three targets.

QTL mapping for pre-harvest sprouting resistance in japonica rice varieties utilizing genome re-sequencing

Pre-harvest sprouting (PHS) leads to serious economic losses because of reductions in yield and quality. To analyze the quantitative trait loci (QTLs) for PHS resistance in japonica rice, PHS rates on panicles were measured in 160 recombinant inbred lines (RILs) derived from a cross between the temperate japonica varieties Odae (PHS resistant) and Unbong40 (PHS susceptible) under two different environmental conditions—field (summer) and greenhouse (winter) environments. Genome re-sequencing of the parental varieties detected 266,773 DNA polymorphisms including 248,255 single nucleotide polymorphisms and 18,518 insertions/deletions. We constructed a genetic map comprising 239 kompetitive allele-specific PCR and 49 cleaved amplified polymorphic sequence markers. In the field environment, two major QTLs, qPHS-3^FD and qPHS-11^FD, were identified on chromosomes 3 and 11, respectively, whereas three major QTLs, qPHS-3^GH, qPHS-4^GH, and qPHS-11^GH, were identified on chromosomes 3, 4, and 11, respectively, in the greenhouse environment. qPHS-11^GH and qPHS-11^FD had similar locations on chromosome 11, suggesting the existence of a gene conferring stable PHS resistance effects under different environmental conditions. The QTLs identified in this study can be used to improve the PHS resistance of japonica varieties, and they may improve our understanding of the genetic basis of PHS resistance.

Whole genome sequencing analysis of horse populations inhabiting the Korean Peninsula and Przewalski's horse

BACKGROUND

The Jeju horse is an indigenous horse breed in Korea. However, there is a severe lack of genomic studies on Korean horse breeds.

OBJECTIVE

The objective of this study was to report genomic characteristics of domestic horse populations that inhabit South Korea (Jeju, Jeju crossbred, and Thoroughbred) and a wild horse breed (Przewalski's horse).

RESULTS

Using the equine reference genome assembly (EquCab 2.0), more than ~ 6.5 billion sequence reads were successfully mapped, which generated an average of 40.87-fold coverage throughout the genome. Using these data, we detected a total of 12.88 million SNPs, of which 73.7% were found to be novel. All the detected SNPs were deeply annotated to retrieve SNPs in gene regions using the RefSeq and Ensemble gene sets. Approximately 27% of the total SNPs were located within genes, whereas the remaining 73% were found in intergenic regions. Using 129,776 coding SNPs, we retrieved a total of 49,171 nonsynonymous SNPs in 12,351 genes. Furthermore, we identified a total of 10,770 deleterious nonsynonymous SNPs which are predicted to affect protein structure or function.

CONCLUSION

We showed numerous genomic variants from domestic and wild horse breeds. These results provide a valuable resource for further studies on functions of SNP-containing genes, and can aid in determining the molecular basis underlying variation in economically important traits of horses.

Re-sequencing Resources to Improve Starch and Grain Quality in Rice

Next-generation sequencing can identify differences in the rice genome that explain the genetic basis of grain quality variation. Differences in rice grain quality are mainly associated with differences in the major component of the grain, starch. Association of rice quality variation with rice genome variation can be conducted at the gene or whole-genome level. Re-sequencing of specific genes or whole genomes can be used depending on the extent to which candidate genes for the traits of interest are known. Amplicon sequencing of genes involved in starch metabolism can help in targeted discovery of the molecular genetic basis of differences in starch related quality attributes. Whole-genome re-sequencing can complement this, when the genetic basis of the trait is expected to be outside the coding region of starch metabolism genes. These approaches have been used successfully to understand the rice genome at specific loci and over the whole genome.

Genome-wide association study reveals candidate genes related to low temperature tolerance in rice (Oryza sativa) during germination

In this study, relative germination percentage (RGP) and delayed mean germination time (DMGT) were measured in various rice accessions at the germination stage and carried out association analysis to identify candidate genes related to low temperature germination (LTG) using a natural population comprising 137 rice cultivars and inbred lines selected from the Korean rice core set. Genome-wide association study using ~ 1.44 million high-quality SNPs, which were identified by re-sequencing all rice collections, revealed 48 candidate genes on chromosome 10 and 55 candidate genes on chromosome 11 in the high peak SNP sites of associated loci for RGP and DMGT, respectively. By detecting highly associated variations located inside genic regions and performing functional annotation of the genes, we detected 23 candidate genes for RGP and 18 genes for DMGT for LTG. In addition, the haplotype and sequence analysis of the candidate gene (Os10g0371100) with RGP trait and the candidate gene (Os11t0104240-00) with DMGT revealed correlation between sequences of functional variations and phenotypes. Several novel LTG-related candidate genes previously were known for the function during rice germination and uncovered their substantial natural variations. These candidate genes represent valuable resources for molecular breeding and genetic improvement of cold tolerance during rice germination.

Imputation from SNP chip to sequence: a case study in a Chinese indigenous chicken population

Background

Genome-wide association studies and genomic predictions are thought to be optimized by using whole-genome sequence (WGS) data. However, sequencing thousands of individuals of interest is expensive. Imputation from SNP panels to WGS data is an attractive and less expensive approach to obtain WGS data. The aims of this study were to investigate the accuracy of imputation and to provide insight into the design and execution of genotype imputation.

Results

We genotyped 450 chickens with a 600 K SNP array, and sequenced 24 key individuals by whole genome re-sequencing. Accuracy of imputation from putative 60 K and 600 K array data to WGS data was 0.620 and 0.812 for Beagle, and 0.810 and 0.914 for FImpute, respectively. By increasing the sequencing cost from 24X to 144X, the imputation accuracy increased from 0.525 to 0.698 for Beagle and from 0.654 to 0.823 for FImpute. With fixed sequence depth (12X), increasing the number of sequenced animals from 1 to 24, improved accuracy from 0.421 to 0.897 for FImpute and from 0.396 to 0.777 for Beagle. Using optimally selected key individuals resulted in a higher imputation accuracy compared with using randomly selected individuals as a reference population for re-sequencing. With fixed reference population size (24), imputation accuracy increased from 0.654 to 0.875 for FImpute and from 0.512 to 0.762 for Beagle as the sequencing depth increased from 1X to 12X. With a given total cost of genotyping, accuracy increased with the size of the reference population for FImpute, but the pattern was not valid for Beagle, which showed the highest accuracy at six fold coverage for the scenarios used in this study.

Conclusions

In conclusion, we comprehensively investigated the impacts of several key factors on genotype imputation. Generally, increasing sequencing cost gave a higher imputation accuracy. But with a fixed sequencing cost, the optimal imputation enhance the performance of WGP and GWAS. An optimal imputation strategy should take size of reference population, imputation algorithms, marker density, and population structure of the target population and methods to select key individuals into consideration comprehensively. This work sheds additional light on how to design and execute genotype imputation for livestock populations.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0241-5) contains supplementary material, which is available to authorized users.

Identification of deleterious rare variants in MTTP, PNPLA3, and TM6SF2 in Japanese males and association studies with NAFLD

Background

Non-alcoholic fatty liver disease (NAFLD) is a disorder characterized by excessive fat deposits in hepatocytes without excessive alcohol intake. NAFLD is influenced by genetic factors, and the heritability has been estimated at 0.35 to 0.6 by twin studies. We explored rare variants in known NAFLD-associated genes to investigate whether these rare variants are involved in the susceptibility to NAFLD.

Methods

The target genes for re-sequencing were PNPLA3, TM6SF2, and MTTP. All exons of these three genes were amplified from a discovery panel of 950 Japanese males, and the identified rare variants were further tested for genetic association in 3014 individuals from the Japanese general population and for in vitro functional evaluation.

Results

Target re-sequencing analysis using next-generation sequencing identified 29 rare variants in 65 Japanese males (6.84%), 12 of which were newly identified base substitutions. A splicing mutation in TM6SF2 that resulted in deletion of 31 amino acids was identified in an NAFLD case. Among eight genotyped rare single-nucleotide polymorphisms (SNPs; minor allele frequency < 0.02), rs143392071 (Tyr220Cys, PNPLA3) significantly increased (odds ratio = 3.52, P = 0.008) and rs756998920 (Val42Ile, MTTP) significantly decreased (odds ratio = 0.03, P = 0.019) the NAFLD risk. Functional assays showed that these two SNPs disrupted protein functions and supported the genetic association.

Conclusion

Collectively, 1.79% of individuals in our studied population were estimated carriers of rare variants that are potentially associated with NAFLD.

Electronic supplementary material

The online version of this article (10.1186/s12944-017-0570-y) contains supplementary material, which is available to authorized users.

Identification of a stable major-effect QTL (Parth 2.1) controlling parthenocarpy in cucumber and associated candidate gene analysis via whole genome re-sequencing

BACKGROUND

Parthenocarpy is an important trait for yield and quality in many plants. But due to its complex interactions with genetic and physiological factors, it has not been adequately understood and applied to breeding and production. Finding novel and effective quantitative trait loci (QTLs) is a critical step towards understanding its genetic mechanism. Cucumber (Cucumis sativus L.) is a typical parthenocarpic plant but the QTLs controlling parthenocarpy in cucumber were not mapped on chromosomes, and the linked markers were neither user-friendly nor confirmed by previous studies. Hence, we conducted a two-season QTL study of parthenocarpy based on the cucumber genome with 145 F2:3 families derived from a cross between EC1 (a parthenocarpic inbred line) and 8419 s-1 (a non-parthenocarpic inbred line) in order to map novel QTLs. Whole genome re-sequencing was also performed both to develop effective linked markers and to predict candidate genes.

RESULTS

A genetic linkage map, employing 133 Simple Sequence Repeats (SSR) markers and nine Insertion/Deletion (InDel) markers spanning 808.1 cM on seven chromosomes, was constructed from an F2 population. Seven novel QTLs were identified on chromosomes 1, 2, 3, 5 and 7. Parthenocarpy 2.1 (Parth2.1), a QTL on chromosome 2, was a major-effect QTL with a logarithm of odds (LOD) score of 9.0 and phenotypic variance explained (PVE) of 17.0 % in the spring season and with a LOD score of 6.2 and PVE of 10.2 % in the fall season. We confirmed this QTL using a residual heterozygous line97-5 (RHL97-5). Effectiveness of linked markers of the Parth2.1 was validated in F3:4 population and in 21 inbred lines. Within this region, there were 57 genes with nonsynonymous SNPs/InDels in the coding sequence. Based on further combined analysis with transcriptome data between two parents, CsARF19, CsWD40, CsEIN1, CsPPR, CsHEXO3, CsMDL, CsDJC77 and CsSMAX1 were predicted as potential candidate genes controlling parthenocarpy.

CONCLUSIONS

A major-effect QTL Parth2.1 and six minor-effect QTLs mainly contribute to the genetic architecture of parthenocarpy in cucumber. SSR16226 and Indel-T-39 can be used in marker-assisted selection (MAS) of cucumber breeding. Whole genome re-sequencing enhances the efficiency of polymorphic marker development and prediction of candidate genes.

Peeping into genomic architecture by re-sequencing of Ochrobactrum intermedium M86 strain during laboratory adapted conditions

Advances in de novo sequencing technologies allow us to track deeper insights into microbial genomes for restructuring events during the course of their evolution inside and outside the host. Bacterial species belonging to Ochrobactrum genus are being reported as emerging, and opportunistic pathogens in this technology driven era probably due to insertion and deletion of genes. The Ochrobactrumintermedium M86 was isolated in 2005 from a case of non-ulcer dyspeptic human stomach followed by its first draft genome sequence in 2009. Here we report re-sequencing of O. intermedium M86 laboratory adapted strain in terms of gain and loss of genes. We also attempted for finer scale genome sequence with 10 times more genome coverage than earlier one followed by comparative evaluation on Ion PGM and Illumina MiSeq. Despite their similarities at genomic level, lab-adapted strain mainly lacked genes encoding for transposase protein, insertion elements family, phage tail-proteins that were not detected in original strain on both chromosomes. Interestingly, a 5 kb indel was detected in chromosome 2 that was absent in original strain mapped with phage integrase gene of Rhizobium spp. and may be acquired and integrated through horizontal gene transfer indicating the gene loss and gene gain phenomenon in this genus. Majority of indel fragments did not match with known genes indicating more bioinformatic dissection of this fragment. Additionally we report genes related to antibiotic resistance, heavy metal tolerance in earlier and re-sequenced strain. Though SNPs detected, there did not span urease and flagellar genes. We also conclude that third generation sequencing technologies might be useful for understanding genomic architecture and re-arrangement of genes in the genome due to their ability of larger coverage that can be used to trace evolutionary aspects in microbial system.

Re-sequencing and genetic variation identification of a rice line with ideal plant architecture

BACKGROUND

The ideal plant architecture (IPA) includes several important characteristics such as low tiller numbers, few or no unproductive tillers, more grains per panicle, and thick and sturdy stems. We have developed an indica restorer line 7302R that displays the IPA phenotype in terms of tiller number, grain number, and stem strength. However, its mechanism had to be clarified.

FINDINGS

We performed re-sequencing and genome-wide variation analysis of 7302R using the Solexa sequencing technology. With the genomic sequence of the indica cultivar 9311 as reference, 307 627 SNPs, 57 372 InDels, and 3 096 SVs were identified in the 7302R genome. The 7302R-specific variations were investigated via the synteny analysis of all the SNPs of 7302R with those of the previous sequenced none-IPA-type lines IR24, MH63, and SH527. Moreover, we found 178 168 7302R-specific SNPs across the whole genome and 30 239 SNPs in the predicted mRNA regions, among which 8 517 were Non-syn CDS. In addition, 263 large-effect SNPs that were expected to affect the integrity of encoded proteins were identified from the 7302R-specific SNPs. SNPs of several important previously cloned rice genes were also identified by aligning the 7302R sequence with other sequence lines.

CONCLUSIONS

Our results provided several candidates account for the IPA phenotype of 7302R. These results therefore lay the groundwork for long-term efforts to uncover important genes and alleles for rice plant architecture construction, also offer useful data resources for future genetic and genomic studies in rice.