Using RNA-seq for genomic scaffold placement, correcting assemblies, and genetic map creation in a commonBrassica rapamapping population

Abstract

Brassica rapa is a model species for agronomic, ecological, evolutionary, and translational studies. Here, we describe high-density SNP discovery and genetic map construction for a B. rapa recombinant inbred line (RIL) population derived from field collected RNA sequencing (RNA-Seq) data. This high-density genotype data enables the detection and correction of putative genome misassemblies and accurate assignment of scaffold sequences to their likely genomic locations. These assembly improvements represent 7.1-8.0% of the annotated B. rapa genome. We demonstrate how using this new resource leads to a significant improvement for QTL analysis over the current low-density genetic map. Improvements are achieved by the increased mapping resolution and by having known genomic coordinates to anchor the markers for candidate gene discovery. These new molecular resources and improvements in the genome annotation will benefit the Brassicaceae genomics community and may help guide other communities in fine-tuning genome annotations. KEYWORDS

Multi-Omics Driven Assembly and Annotation of the Sandalwood (Santalum album) Genome

“…In the same plant, RNA-seq analysis allowed for the identification of several novel protein-coding genes together with novel transcribed regions, noncoding RNAs, and small RNA loci that were not annotated earlier (Cheng et al, 2017). RNA-seq-assisted genome annotation of two other members of the Brassicaceae (Arabidopsis lyrata and Brassica rapa) also has been reported (Rawat et al, 2015;Markelz et al, 2017). Among monocots, proteogenomics studies have been carried out in maize (Zea mays), wherein a total of 165 novel protein-coding genes were reported and gene models for 741 additional genes were corrected (Castellana et al, 2014).…”

Allocation to male vs female floral function varies by currency and responds differentially to density and moisture stress

“…QTL mapping and candidate gene analyses. We used a dense linkage map consisting of 1273 SNPs that were identified from RNAseq analyses and were mapped across the 10 B. rapa chromosomes in the BraIRRI population (Iniguez-Luy et al, 2009;Brock et al, 2016;Markelz et al, 2017. To map QTL, we performed composite interval mapping as implemented in QTL Cartographer (ver. 2.5; Wang et al, 2007).…”

“…We mapped components of female (average seed per fruit and total seed set) and male fitness (average viable pollen per flower and total viable pollen production), organ mass, and CNP content, concentration and proportion in each organ. Given the refined mapping resolution of a recent high-density linkage map (Markelz et al, 2017) and large number of traits, we mapped a total of 200 QTL. Individual QTL explained 9.3% (range: 4-23%) of the genetic variance on average.…”

A framework for genomics-informed ecophysiological modeling in plants

“…Genomic prediction was performed on the BraIRRi RIL population (Iniguez-Luy et al , 2009) for the median, shape, and rate estimates of each leaf growth parameter, K , N 0 , and r . A published genome-wide 1482 single nucleotide polymorphism (SNP) data set (Markelz et al , 2017) was utilized for 123 individuals that also had phenotypic data. After quality control, 1481 SNPs with minor allele frequency >0.01 were retained for genomic prediction.…”