A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable. RESULTS: Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq). CONCLUSIONS: This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

A terpene synthase supergene locus determines chemotype in Melaleuca alternifolia (tea tree).

Leaf oil terpenes vary categorically in many plant populations, leading to discrete phenotypes of adaptive and economic significance, but for most species, a genetic explanation for the concerted fluctuation in terpene chemistry remains unresolved. To uncover the genetic architecture underlying multi-component terpene chemotypes in Melaleuca alternifolia (tea tree), a genome-wide association study was undertaken for 148 individuals representing all six recognised chemotypes. A number of single nucleotide polymorphisms in a genomic region of c. 400 kb explained large proportions of the variation in key monoterpenes of tea tree oil. The region contained a cluster of 10 monoterpene synthase genes, including four genes predicted to encode synthases for 1,8-cineole, terpinolene, and the terpinen-4-ol precursor, sabinene hydrate. Chemotype-dependent null alleles at some sites suggested structural variants within this gene cluster, providing a possible basis for linkage disequilibrium in this region. Genotyping in a separate domesticated population revealed that all alleles surrounding this gene cluster were fixed after artificial selection for a single chemotype. These observations indicate that a supergene accounts for chemotypes in M. alternifolia. A genetic model with three haplotypes, encompassing the four characterised monoterpene synthase genes, explained the six terpene chemotypes, and was consistent with available biparental cross-segregation data.

Molecular dissection of connected rice populations revealed important genomic regions for agronomic and biofortification traits.

Breeding staple crops with increased micronutrient concentration is a sustainable approach to address micronutrient malnutrition. We carried out Multi-Cross QTL analysis and Inclusive Composite Interval Mapping for 11 agronomic, yield and biofortification traits using four connected RILs populations of rice. Overall, MC-156 QTLs were detected for agronomic (115) and biofortification (41) traits, which were higher in number but smaller in effects compared to single population analysis. The MC-QTL analysis was able to detect important QTLs viz: qZn5.2, qFe7.1, qGY10.1, qDF7.1, qPH1.1, qNT4.1, qPT4.1, qPL1.2, qTGW5.1, qGL3.1 , and qGW6.1 , which can be used in rice genomics assisted breeding. A major QTL (qZn5.2 ) for grain Zn concentration has been detected on chromosome 5 that accounted for 13% of R2. In all, 26 QTL clusters were identified on different chromosomes. qPH6.1 epistatically interacted with qZn5.1 and qGY6.2 . Most of QTLs were co-located with functionally related candidate genes indicating the accuracy of QTL mapping. The genomic region of qZn5.2 was co-located with putative genes such as OsZIP5, OsZIP9, and LOC_OS05G40490 that are involved in Zn uptake. These genes included polymorphic functional SNPs, and their promoter regions were enriched with cis-regulatory elements involved in plant growth and development, and biotic and abiotic stress tolerance. Major effect QTL identified for biofortification and agronomic traits can be utilized in breeding for Zn biofortified rice varieties.

Critical disease burdens of Australian adults with cystic fibrosis: Results from an online survey.

BACKGROUND: The objective of this study was to conduct a web-based questionnaire to investigate self-reported phenotypes and disease burdens of individuals living in Australia and diagnosed with cystic fibrosis (CF) using a case-control study design. METHODS: An online questionnaire was distributed to individuals with CF and healthy control subjects. Overall health rating, medications, family history, education, clinical indicators of disease, and symptoms, including their severity and frequency, were evaluated. RESULTS: There was a total of 119 respondents consisting of 59 people living with CF and 60 controls. The CF cohort had significantly lower tertiary educational levels compared to controls. The analysis specific to the CF cohort depicted a significant correlation between the frequency of hospitalizations and the level of education in the CF cohort. Of the 26 self-reported symptoms of CF that were analyzed, 14 were significantly higher in the people living with CF. The CF cohort reporting symptoms of chronic pain (25%) described an increase in the burden of disease, depicting a 30% longer mean hospitalization, increased consumption of medications and significant relationships with four other symptoms, including muscle aches, digestive issues, pancreatic insufficiency, and abdominal swelling. CONCLUSIONS: The nationwide survey identified a diverse range of clinical manifestations experienced by the Australian CF population. Chronic pain, linked to aging and the changing landscape of disease, was a significant indicator of the burden of disease. A comprehensive understanding of the phenotypic profiles and symptom variability will contribute to future research and provide insights into the impacts of disease and the burden of therapy, particularly in children, at the start of their health journey.

Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379.

The cotton chromosome substitution line, CS-B15sh, exhibits 41% lower injury from 2,4-D when applied at the field recommended rate of 1.12 kg ae ha-1 (1×) than does Texas Marker-1 (TM-1). CS-B15sh was developed in the genetic background of Gossypium hirsutum L. cv TM-1 and has chromosome introgression on the short arm of chromosome 15 from Gossypium barbadense L. cv. Pima 379. In a previous experiment, we observed reduced translocation of [14C]2,4-D outside the treated leaf tissue in CS-B15sh, which contrasted with an increased translocation of the herbicide in the tissues above and below the treated leaf in TM-1. Our results indicate a potential 2,4-D tolerance mechanism in CS-B15sh involving altered movement of 2,4-D. Here, we used RNA sequencing (RNA-seq) to determine the differential expression of genes between 2,4-D-challenged and control plants of the tolerant (CS-B15sh) and susceptible lines (TM-1 and Pima 379). Several components of the 2,4-D/auxin-response pathway-including ubiquitin E3 ligase, PB1|AUX/IAA, ARF transcription factors, and F-box proteins of the SCFTIR1/AFB complex-were upregulated with at least threefold higher expression in TM-1 compared with CS-B15sh, while both Pima 379 and TM-1 showed the same fold change expression for PB1|AUX/IAA mRNA. Some genes associated with herbicide metabolism, including flavin monooxygenase (Gohir.A01G174100) and FAD-linked oxidase (Gohir.D06G002600), exhibited at least a twofold increase in CS-B15sh than in TM-1 (the gene was not expressed in Pima 379), suggesting a potential relationship between the gene's expression and 2,4-D tolerance. It is interesting to note that glutathione S-transferase was differentially expressed in both CS-B15sh and Pima 379 but not in TM-1, while cytochrome P450 and other genes involved in the oxidation-reduction process were significantly expressed only in CS-B15sh in response to 2,4-D. Gene set enrichment analysis on the union DEGs of the three cotton genotypes revealed the depletion of transcripts involved in photosynthesis and enrichment of transcripts involved in ABA response and signaling.

Establishing a Common Nutritional Vocabulary - From Food Production to Diet.

Informed policy and decision-making for food systems, nutritional security, and global health would benefit from standardization and comparison of food composition data, spanning production to consumption. To address this challenge, we present a formal controlled vocabulary of terms, definitions, and relationships within the Compositional Dietary Nutrition Ontology (CDNO, www.cdno.info) that enables description of nutritional attributes for material entities contributing to the human diet. We demonstrate how ongoing community development of CDNO classes can harmonize trans-disciplinary approaches for describing nutritional components from food production to diet.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.).

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Shotgun proteomics of Brassica rapa seed proteins identifies vicilin as a major seed storage protein in the mature seed.

Proteins make up a large percentage of the Brassica seed and are second only to the oil in economic importance with uses for both animal and human nutrition. The most abundant proteins reported in the seeds of Brassica are the seed storage proteins cruciferin and napin, belonging to the 12S globulin and 2S albumin families of proteins, respectively. To gain insight into the Brassica rapa seed proteome and to confirm the presence and relative quantity of proteins encoded by candidate seed storage genes in the mature seed, shotgun proteomics was carried out on protein extracts from seeds of B. rapa inbred line R-o-18. Following liquid chromatography tandem mass spectrometry, a total of 34016 spectra were mapped to 323 proteins, where 233 proteins were identified in 3 out of 4 biological replicates by at least 2 unique peptides. 2S albumin like napin seed storage proteins (SSPs), 11/12S globulin like cruciferin SSPs and 7S globulin like vicilin SSPs were identified in the samples, along with other notable proteins including oil body proteins, namely ten oleosins and two oil body-associated proteins. The identification of vicilin like proteins in the mature B. rapa seed represents the first account of these proteins in the Brassicaceae and analysis indicates high conservation of sequence motifs to other 7S vicilin-like allergenic proteins as well as conservation of major allergenic epitopes in the proteins. This study enriches our existing knowledge on rapeseed seed proteins and provides a robust foundation and rational basis for plant bioengineering of seed storage proteins.

Challenges for FAIR-compliant description and comparison of crop phenotype data with standardized controlled vocabularies.

Crop phenotypic data underpin many pre-breeding efforts to characterize variation within germplasm collections. Although there has been an increase in the global capacity for accumulating and comparing such data, a lack of consistency in the systematic description of metadata often limits integration and sharing. We therefore aimed to understand some of the challenges facing findable, accesible, interoperable and reusable (FAIR) curation and annotation of phenotypic data from minor and underutilized crops. We used bambara groundnut (Vigna subterranea) as an exemplar underutilized crop to assess the ability of the Crop Ontology system to facilitate curation of trait datasets, so that they are accessible for comparative analysis. This involved generating a controlled vocabulary Trait Dictionary of 134 terms. Systematic quantification of syntactic and semantic cohesiveness of the full set of 28 crop-specific COs identified inconsistencies between trait descriptor names, a relative lack of cross-referencing to other ontologies and a flat ontological structure for classifying traits. We also evaluated the Minimal Information About a Phenotyping Experiment and FAIR compliance of bambara trait datasets curated within the CropStoreDB schema. We discuss specifications for a more systematic and generic approach to trait controlled vocabularies, which would benefit from representation of terms that adhere to Open Biological and Biomedical Ontologies principles. In particular, we focus on the benefits of reuse of existing definitions within pre- and post-composed axioms from other domains in order to facilitate the curation and comparison of datasets from a wider range of crops. Database URL: https://www.cropstoredb.org/cs_bambara.html.

Characterization of the Cannabis sativa glandular trichome proteome.

Cannabis sativa has been cultivated since antiquity as a source of fibre, food and medicine. The recent resurgence of C. sativa as a cash crop is mainly driven by the medicinal and therapeutic properties of its resin, which contains compounds that interact with the human endocannabinoid system. Compared to other medicinal crops of similar value, however, little is known about the biology of C. sativa. Glandular trichomes are small hair-like projections made up of stalk and head tissue and are responsible for the production of the resin in C. sativa. Trichome productivity, as determined by C. sativa resin yield and composition, is only beginning to be understood at the molecular level. In this study the proteomes of glandular trichome stalks and heads, were investigated and compared to the proteome of the whole flower tissue, to help further elucidate C. sativa glandular trichome biochemistry. The data suggested that the floral tissue acts as a major source of carbon and energy to the glandular trichome head sink tissue, supplying sugars which drive secondary metabolite biosynthesis. The trichome stalk seems to play only a limited role in secondary metabolism and acts as both source and sink.

A high-quality draft genome for Melaleuca alternifolia (tea tree): a new platform for evolutionary genomics of myrtaceous terpene-rich species.

The economically important Melaleuca alternifolia (tea tree) is the source of a terpene-rich essential oil with therapeutic and cosmetic uses around the world. Tea tree has been cultivated and bred in Australia since the 1990s. It has been extensively studied for the genetics and biochemistry of terpene biosynthesis. Here, we report a high quality de novo genome assembly using Pacific Biosciences and Illumina sequencing. The genome was assembled into 3128 scaffolds with a total length of 362 Mb (N50  = 1.9 Mb), with significantly higher contiguity than a previous assembly (N50  = 8.7 Kb). Using a homology-based, RNA-seq evidence-based and ab initio prediction approach, 37,226 protein-coding genes were predicted. Genome assembly and annotation exhibited high completeness scores of 98.1% and 89.4%, respectively. Sequence contiguity was sufficient to reveal extensive gene order conservation and chromosomal rearrangements in alignments with Eucalyptus grandis and Corymbia citriodora genomes. This new genome advances currently available resources to investigate the genome structure and gene family evolution of M. alternifolia. It will enable further comparative genomic studies in Myrtaceae to elucidate the genetic foundations of economically valuable traits in this crop.

An extreme-phenotype genome-wide association study identifies candidate cannabinoid pathway genes in Cannabis.

Cannabis produces a class of isoprenylated resorcinyl polyketides known as cannabinoids, a subset of which are medically important and exclusive to this plant. The cannabinoid alkyl group is a critical structural feature that governs therapeutic activity. Genetic enhancement of the alkyl side-chain could lead to the development of novel chemical phenotypes (chemotypes) for pharmaceutical end-use. However, the genetic determinants underlying in planta variation of cannabinoid alkyl side-chain length remain uncharacterised. Using a diversity panel derived from the Ecofibre Cannabis germplasm collection, an extreme-phenotype genome-wide association study (XP-GWAS) was used to enrich for alkyl cannabinoid polymorphic regions. Resequencing of chemotypically extreme pools revealed a known cannabinoid synthesis pathway locus as well as a series of chemotype-associated genomic regions. One of these regions contained a candidate gene encoding a ß-keto acyl carrier protein (ACP) reductase (BKR) putatively associated with polyketide fatty acid starter unit synthesis and alkyl side-chain length. Association analysis revealed twenty-two polymorphic variants spanning the length of this gene, including two nonsynonymous substitutions. The success of this first reported application of XP-GWAS for an obligate outcrossing and highly heterozygote plant genus suggests that this approach may have generic application for other plant species.

Characterization of the Leaf Microbiome from Whole-Genome Sequencing Data of the 3000 Rice Genomes Project.

BACKGROUND: The crop microbial communities are shaped by interactions between the host, microbes and the environment, however, their relative contribution is beginning to be understood. Here, we explore these interactions in the leaf bacterial community across 3024 rice accessions. FINDINGS: By using unmapped DNA sequencing reads as microbial reads, we characterized the structure of the rice bacterial microbiome. We identified central bacteria taxa that emerge as microbial "hubs" and may have an influence on the network of host-microbe interactions. We found regions in the rice genome that might control the assembly of these microbial hubs. To our knowledge this is one of the first studies that uses raw data from plant genome sequencing projects to characterize the leaf bacterial communities. CONCLUSION: We showed, that the structure of the rice leaf microbiome is modulated by multiple interactions among host, microbes, and environment. Our data provide insight into the factors influencing microbial assemblage in the rice leaf and also opens the door for future initiatives to modulate rice consortia for crop improvement efforts.

Comparative Transcriptomics and Co-Expression Networks Reveal Tissue- and Genotype-Specific Responses of qDTYs to Reproductive-Stage Drought Stress in Rice (Oryza sativa L.).

Rice (Oryza sativa L.) is more sensitive to drought stress than other cereals. To dissect molecular mechanisms underlying drought-tolerant yield in rice, we applied differential expression and co-expression network approaches to transcriptomes from flag-leaf and emerging panicle tissues of a drought-tolerant yield introgression line, DTY-IL, and the recurrent parent Swarna, under moderate reproductive-stage drought stress. Protein turnover and efficient reactive oxygen species scavenging were found to be the driving factors in both tissues. In the flag-leaf, the responses further included maintenance of photosynthesis and cell wall reorganization, while in the panicle biosynthesis of secondary metabolites was found to play additional roles. Hub genes of importance in differential drought responses included an expansin in the flag-leaf and two peroxidases in the panicle. Overlaying differential expression data with allelic variation in DTY-IL quantitative trait loci allowed for the prioritization of candidate genes. They included a differentially regulated auxin-responsive protein, with DTY-IL-specific amino acid changes in conserved domains, as well as a protein kinase with a DTY-IL-specific frameshift in the C-terminal region. The approach highlights how the integration of differential expression and allelic variation can aid in the discovery of mechanism and putative causal contribution underlying quantitative trait loci for drought-tolerant yield.

Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741).

Macadamia integrifolia is a representative of the large basal eudicot family Proteaceae and the main progenitor species of the Australian native nut crop macadamia. Since its commercialisation in Hawaii fewer than 100 years ago, global production has expanded rapidly. However, genomic resources are limited in comparison to other horticultural crops. The first draft assembly of M. integrifolia had good coverage of the functional gene space but its high fragmentation has restricted its use in comparative genomics and association studies. Here we have generated an improved assembly of cultivar HAES 741 (4,094 scaffolds, 745 Mb, N50 413 kb) using a combination of Illumina paired and PacBio long read sequences. Scaffolds were anchored to 14 pseudo-chromosomes using seven genetic linkage maps. This assembly has improved contiguity and coverage, with >120 Gb of additional sequence. Following annotation, 34,274 protein-coding genes were predicted, representing 90% of the expected gene content. Our results indicate that the macadamia genome is repetitive and heterozygous. The total repeat content was 55% and genome-wide heterozygosity, estimated by read mapping, was 0.98% or an average of one SNP per 102 bp. This is the first chromosome-scale genome assembly for macadamia and the Proteaceae. It is expected to be a valuable resource for breeding, gene discovery, conservation and evolutionary genomics.

Identification, characterization and epitope mapping of proteins encoded by putative allergenic napin genes from Brassica rapa.

BACKGROUND: Brassica rapeseed crops contain high concentrations of oil in the seed. The remaining meal, following oil extraction, has a high protein content, but is of low value due to the presence of high amounts of napin seed storage proteins. These 2S albumin-like proteins are difficult to digest and have been identified as major allergens in humans. OBJECTIVE: To comprehensively characterize the napin gene (NG) family in Brassica rapa and to gain an understanding of the structural basis of allergenicity of the expressed proteins. METHODS: To identify candidate napin genes in B rapa, 2S albumin-like napin genes of Arabidopsis thaliana were used as query sequences to search for similarity against the B rapa var. pekinensis Chiifu-401 v2 and the var. trilocularis R-o-18 v1.5 genomes. Multiple sequence alignment (MSA) and epitope modelling was carried out to determine structural and evolutionary relationships of NGs and their potential allergenicity. RESULTS: Four candidate napin genes in R-o-18 and ten in Chiifu-401 were identified with high sequence similarity to A thaliana napin genes. Multiple sequence alignment revealed strong conservation among the candidate genes. An epitope survey indicated high conservation of allergenic epitope motifs with known 2S albumin-like allergens. CONCLUSION: Napin is thought to be responsible for a  high prevalence of food allergies. Characterization of the napin gene family in B rapa will give important insight into the protein structure, and epitope modelling will help to advance studies into allergenicity including the development of precise diagnostic screenings and therapies for this potential food allergy as well as the possible manipulation of napin levels in the seed by gene editing technology.

A platinum standard pan-genome resource that represents the population structure of Asian rice.

As the human population grows from 7.8 billion to 10 billion over the next 30 years, breeders must do everything possible to create crops that are highly productive and nutritious, while simultaneously having less of an environmental footprint. Rice will play a critical role in meeting this demand and thus, knowledge of the full repertoire of genetic diversity that exists in germplasm banks across the globe is required. To meet this demand, we describe the generation, validation and preliminary analyses of transposable element and long-range structural variation content of 12 near-gap-free reference genome sequences (RefSeqs) from representatives of 12 of 15 subpopulations of cultivated Asian rice. When combined with 4 existing RefSeqs, that represent the 3 remaining rice subpopulations and the largest admixed population, this collection of 16 Platinum Standard RefSeqs (PSRefSeq) can be used as a template to map resequencing data to detect virtually all standing natural variation that exists in the pan-genome of cultivated Asian rice.

Maximising recombination across macadamia populations to generate linkage maps for genome anchoring.

The Proteaceae genus Macadamia has a recent history of domestication as a commercial nut crop. We aimed to establish the first sequence-based haploid-correlated reference genetic linkage maps for this primarily outcrossing perennial tree crop, with marker density suitable for genome anchoring. Four first generation populations were used to maximise the segregation patterns available within full-sib, biparental and self-pollinated progeny. This allowed us to combine segregation data from overlapping subsets of >4,000 informative sequence-tagged markers to increase the effective coverage of the karyotype represented by the recombinant crossover events detected. All maps had 14 linkage groups, corresponding to the Macadamia haploid chromosome number, and enabled the anchoring and orientation of sequence scaffolds to construct a pseudo-chromosomal genome assembly for macadamia. Comparison of individual maps indicated a high level of congruence, with minor discrepancies satisfactorily resolved within the integrated maps. The combined set of maps significantly improved marker density and the proportion (70%) of the genome sequence assembly anchored. Overall, increasing our understanding of the genetic landscape and genome for this nut crop represents a substantial advance in macadamia genetics and genomics. The set of maps, large number of sequence-based markers and the reconstructed genome provide a toolkit to underpin future breeding that should help to extend the macadamia industry as well as provide resources for the long term conservation of natural populations in eastern Australia of this unique genus.

Exploring genetic architecture of grain yield and quality traits in a 16-way indica by japonica rice MAGIC global population.

Identification of Quantitative Trait Loci (QTL) has been a challenge for complex traits due to the use of populations with narrow genetic base. Most of QTL mapping studies were carried out from crosses made within the subspecies, either indica × indica or japonica × japonica. In this study we report advantages of using Multi-parent Advanced Generation Inter-Crosses global population, derived from a combination of eight indica and eight japonica elite parents, in QTL discovery for yield and grain quality traits. Genome-wide association study and interval mapping identified 38 and 34 QTLs whereas Bayesian networking detected 60 QTLs with 22 marker-marker associations, 32 trait-trait associations and 65 marker-trait associations. Notably, nine known QTLs/genes qPH1/OsGA20ox2, qDF3/OsMADS50, PL, QDg1, qGW-5b, grb7-2, qGL3/GS3, Amy6/Wx gene and OsNAS3 were consistently identified by all approaches for nine traits whereas qDF3/OsMADS50 was co-located for both yield and days-to-flowering traits on chromosome 3. Moreover, we identified a number of candidate QTLs in either one or two analyses but further validations will be needed. The results indicate that this new population has enabled identifications of significant QTLs and interactions for 16 traits through multiple approaches. Pyramided recombinant inbred lines provide a valuable source for integration into future breeding programs.