A Consensus Map in Cultivated Hexaploid Oat Reveals Conserved Grass Synteny with Substantial Subgenome Rearrangement

Citation: Chaffin, A. S., Huang, Y. F., Smith, S., Bekele, W. A., Babiker, E., Gnanesh, B. N., . . . Tinker, N. A. (2016). A Consensus Map in Cultivated Hexaploid Oat Reveals Conserved Grass Synteny with Substantial Subgenome Rearrangement. Plant Genome, 9(2), 21. doi:10.3835/plantgenome2015.10.0102Hexaploid oat (Avena sativa L., 2n = 6x = 42) is a member of the Poaceae family and has a large genome (similar to 12.5 Gb) containing 21 chromosome pairs from three ancestral genomes. Physical rearrangements among parental genomes have hindered the development of linkage maps in this species. The objective of this work was to develop a single high-density consensus linkage map that is representative of the majority of commonly grown oat varieties. Data from a cDNA-derived single-nucleotide polymorphism (SNP) array and genotyping-by-sequencing (GBS) were collected from the progeny of 12 biparental recombinant inbred line populations derived from 19 parents representing oat germplasm cultivated primarily in North America. Linkage groups from all mapping populations were compared to identify 21 clusters of conserved collinearity. Linkage groups within each cluster were then merged into 21 consensus chromosomes, generating a framework consensus map of 7202 markers spanning 2843 cM. An additional 9678 markers were placed on this map with a lower degree of certainty. Assignment to physical chromosomes with high confidence was made for nine chromosomes. Comparison of homeologous regions among oat chromosomes and matches to orthologous regions of rice (Oryza sativa L.) reveal that the hexaploid oat genome has been highly rearranged relative to its ancestral diploid genomes as a result of frequent translocations among chromosomes. Heterogeneous chromosome rearrangements among populations were also evident, probably accounting for the failure of some linkage groups to match the consensus. This work contributes to a further understanding of the organization and evolution of hexaploid grass genomes

Population Genomics Related to Adaptation in Elite Oat Germplasm

Six hundred thirty five oat ( L.) lines and 4561 single-nucleotide polymorphism (SNP) loci were used to evaluate population structure, linkage disequilibrium (LD), and genotype–phenotype association with heading date. The first five principal components (PCs) accounted for 25.3% of genetic variation. Neither the eigenvalues of the first 25 PCs nor the cross-validation errors from = 1 to 20 model-based analyses suggested a structured population. However, the PC and = 2 model-based analyses supported clustering of lines on spring oat vs. southern United States origin, accounting for 16% of genetic variation ( < 0.0001). Single-locus -statistic () in the highest 1% of the distribution suggested linkage groups that may be differentiated between the two population subgroups. Population structure and kinship-corrected LD of = 0.10 was observed at an average pairwise distance of 0.44 cM (0.71 and 2.64 cM within spring and southern oat, respectively). On most linkage groups LD decay was slower within southern lines than within the spring lines. A notable exception was found on linkage group Mrg28, where LD decay was substantially slower in the spring subpopulation. It is speculated that this may be caused by a heterogeneous translocation event on this chromosome. Association with heading date was most consistent across location-years on linkage groups Mrg02, Mrg12, Mrg13, and Mrg24

SNP Discovery and Chromosome Anchoring Provide the First Physically-Anchored Hexaploid Oat Map and Reveal Synteny with Model Species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources

Barley grain (1,3;1,4)-β-glucan content:effects of transcript and sequence variation in genes encoding the corresponding synthase and endohydrolase enzymes

The composition of plant cell walls is important in determining cereal end uses. Unlike other widely consumed cereal grains barley is comparatively rich in (1,3;1,4)-β-glucan, a source of dietary fibre. Previous work showed Cellulose synthase-like genes synthesise (1,3;1,4)-β-glucan in several tissues. HvCslF6 encodes a grain (1,3;1,4)-β-glucan synthase, whereas the function of HvCslF9 is unknown. Here, the relationship between mRNA levels of HvCslF6, HvCslF9, HvGlbI (1,3;1,4)-β-glucan endohydrolase, and (1,3;1,4)-β-glucan content was studied in developing grains of four barley cultivars. HvCslF6 was differentially expressed during mid (8-15 DPA) and late (38 DPA) grain development stages while HvCslF9 transcript was only clearly detected at 8-10 DPA. A peak of HvGlbI expression was detected at 15 DPA. Differences in transcript abundance across the three genes could partially explain variation in grain (1,3;1,4)-β-glucan content in these genotypes. Remarkably narrow sequence variation was found within the HvCslF6 promoter and coding sequence and does not explain variation in (1,3;1,4)-β-glucan content. Our data emphasise the genotype-dependent accumulation of (1,3;1,4)-β-glucan during barley grain development and a role for the balance between hydrolysis and synthesis in determining (1,3;1,4)-β-glucan content, and suggests that other regulatory sequences or proteins are likely to be involved in this trait in developing grain.Guillermo Garcia-Gimenez, Joanne Russell, Matthew K. Aubert, Geoffrey B. Fincher, Rachel A. Burton, Robbie Waugh, Matthew R. Tucker, Kelly Housto

A new genetic linkage map of barley (\u3ci\u3eHordeum vulgare\u3c/i\u3e L.) facilitates genetic dissection of height and spike length and angle

Plant height and spike length and angle are important agronomic traits in the production of barley (Hordeum vulgare L.) due to strong correlations with lodging and disease. The objective of this study was to use QTL analysis to identify genetic regions associated with each trait in a recombinant inbred line (RIL) mapping population derived from a cross of Falcon by Azhul. Falcon is a spring six-row hulless feed barley with long spikes displaying obtuse angles, while Azhul is a spring dwarf, six-row hulless food barley with short spikes displaying acute angles. The population was genotyped using SNP, DArT and SSR markers and quantitative trait loci (QTL) were detected on chromosomes 2H (102.8 cM, spikelength), 3H (89.2 cM, plant height and 38.2, spike angle and length), 4H (19.0 cM, spike length), and 5H(106.7 cM, spike angle). In conclusion, we developed a barley genetic map, which incorporated SNP, DArT and SSR markers, for detection of height and spike length and angle QTL. Three spike angle, one spike length and one plant height QTL were novel and by using comparative genomics we identified possible candidate genes involved in gibberellic acid signaling and auxin- and ethylene-responsive pathways. This knowledge can be used to generate suitable markers for barley breeding improvement

Impact of Food Matrices on Digestibility of Allergens and Poorly Allergenic Homologs

BACKGROUND: Protease resistance is considered a risk factor for allergenicity of proteins, although the correlation is low. It is nonetheless a part of the weight-of-evidence approach, proposed by Codex, for assessing the allergenicity risk of novel food proteins. Susceptibility of proteins to pepsin is commonly tested with purified protein in solution. OBJECTIVE: Food proteins are rarely consumed in purified form. Our aim was to evaluate the impact of experimental and endogenous food matrices on protease susceptibility of homologous protein pairs with different degrees of allergenicity. METHODS: Porcine and shrimp tropomyosin (ST) were subjected to sequential exposure to amylase, pepsin, and pancreatin in their respective endogenous matrix (pork tenderloin/boiled shrimp) and in three different experimental matrices (dessert mousse [DM], soy milk [SM], and chocolate bar [CB]). Digestion was monitored by immunoblotting using tropomyosin-specific antibodies. Recombinant peach and strawberry lipid transfer protein were biotinylated, spiked into both peach and strawberry fruit pulp, and subjected to the same sequential digestion protocol. Digestion was monitored by immunoblotting using streptavidin for detection. RESULTS: Chocolate bar, and to a lesser extent SM, had a clear protective effect against pepsin digestion of porcine tropomyosin (PT) and to a lesser extent of ST. Increased resistance was associated with increased protein content. Spiking experiments with bovine serum albumin (BSA) confirmed the protective effect of a protein-rich matrix. The two tropomyosins were both highly resistant to pepsin in their protein-rich and lean native food matrix. Pancreatin digestion remained rapid and complete, independent of the matrix. The fat-rich environment did not transfer protection against pepsin digestion. Spiking of recombinant peach and strawberry lipid transfer proteins into peach and strawberry pulp did not reveal any differential protective effect that could explain differences in allergenicity of both fruits. CONCLUSIONS: Protein-rich food matrices delay pepsin digestion by saturating the protease. This effect is most apparent for proteins that are highly pepsin susceptible in solution. The inclusion of food matrices does not help in understanding why some proteins are strong primary sensitizers while homologs are very poor allergens. Although for induction of symptoms in food allergic patients (elicitation), a protein-rich food matrix that may contribute to increased risk, our results indicate that the inclusion of food matrices in the weight-of-evidence approach for estimating the potential risks of novel proteins to become allergens (sensitization), is most likely of very limited value

Comparative Systems Biology Reveals Allelic Variation Modulating Tocochromanol Profiles in Barley (Hordeum vulgare L.)

Tocochromanols are recognized for nutritional content, plant stress response, and seed longevity. Here we present a systems biological approach to characterize and develop predictive assays for genes affecting tocochromanol variation in barley. Major QTL, detected in three regions of a SNP linkage map, affected multiple tocochromanol forms. Candidate genes were identified through barley/rice orthology and sequenced in genotypes with disparate tocochromanol profiles. Gene-specific markers, designed based on observed polymorphism, mapped to the originating QTL, increasing R2 values at the respective loci. Polymorphism within promoter regions corresponded to motifs known to influence gene expression. Quantitative PCR analysis revealed a trend of increased expression in tissues grown at cold temperatures. These results demonstrate utility of a novel method for rapid gene identification and characterization, and provide a resource for efficient development of barley lines with improved tocochromanol profiles

Tocochromanol means of Falcon and Azhul parents and two barley checks grown over four location years.

<p>Means and heritability estimates were calculated for Falcon x Azhul (FA) RILs using combined years at each location.</p>a<p>Concentrations given in µg/g. Values followed by the same letter or number of asterisks are not significantly different within a column (p<0.05). For means of genotypes within an form, a single asterisk is equivalent to A, two asterisks to B, and three asterisks to C.</p>+<p>indicates an intermediate value.</p>b<p>Irrigated field trial in Aberdeen.</p>c<p>Non-irrigated field trial at Tetonia.</p>d<p>Broad sense heritability calculated as genotype variance divided by cumulative variance including error.</p