Identification of genomic regions of Sorghum bicolor (L.) Moench linked to biofuel-related traits in grain x sweet sorghum recombinant inbred lines

Abstract

The increasing demand of food and the rising concerns about climate change and energy security pushes Sorghum bicolor to the top of global agendas. The main attraction is its ability to provide food and animal feed as well as renewable energy products, and industrial commodities. The species comprises grain and sweet sorghum types. Sweet sorghum has major advantages compared to other sugar and biomass crops. Firstly, sorghum is a species showing extensive genetic variation, including drought and heat tolerant genotypes, which enable the usage of marginal land that is not suitable for cultivating other crops. Secondly, competition between the use of land for food or for energy is less because the grains can be used for food or feed while stems can be used for biofuel production. In the current study, field trials and related genetic analysis were carried out to identify major chromosomal regions that are linked to biofuel-related traits. An experimental population of 213 RILs from a cross between grain sorghum (M71) and sweet sorghum (SS79) was planted in 5 environments in Germany in the vegetation periods 2007 and 2008. This population segregates for alleles controlling biomass yield, juice content, sugar content, fibre content, and many other traits. The approach to identify these chromosomal regions was (1) to phenotype the RIL population for the mentioned traits, (2) to construct a genetic map from the RILs, and (3) to correlate the phenotypic data and genetic map informations. The phenotypic data was analysed for individual location and environments (years*location) using mixed model from SAS® 9.1 version. Genotypes were treated as fixed effects while all other components were treated as random effects. Regarding, sugar-related traits, we analysed brix, glucose content, sucrose content, and sugar content, measured stem diameter, stem juice weight, plant height and fresh panicle weight, and counted the days to anthesis and the number of tillers per plant. We partitioned the variance components as genetic, location, year, genetic x location, genetic x year, location x year, genetic x location x year, and harvest date. The results showed that although location contributed highly to the phenotype in terms of percentage, it was not a significant source of variance. The interaction between genotypes and environment was observed in most traits except for stem juice weight, fresh panicle weight and glucose content. Genetic variance contributed small but significant to the phenotypic variance. And thus the heritability estimates were moderate to high in all traits except stem juice weight (H2=0.18). Pearson correlation estimates showed that these traits were significantly correlated with each other. Flowering dates correlated with all the traits while number of tillers only positively correlated with stem juice weight, and negatively correlated with sugar content and flowering dates. In fibre-related traits, we quantified fibre content traits (acidic detergent fibre, acidic detergent lignin, neutral detergent fibre, cellulose, hemicelluloses), and measured related agronomic traits (fresh leaf mass, stripped stalk mass, dry stalk mass, fresh biomass and dry biomass). We partitioned our variance components as genetic, environment, harvest date and genotype x environment. Genetic variance was significant in all traits except for dry biomass while environment was not significant for all traits. There was a significant genotype x environment interaction in all traits except hemicelluloses. Dry stalk mass, dry biomass and hemicelluloses had low heritability estimates of <=0.24. Pearson correlation estimates showed that fibre content traits significantly and positively correlated with each other as well as agronomic traits, while there were negative correlations among agronomic and fibre content traits.A total of 213 RILs were used to construct the genetic linkage map using 157 markers, i.e. simple sequence repeats (SSR), amplified fragment length polymorphism (AFLP) and expressed sequence tag-simple sequence repeats (EST-SSR). The linkage map was compared to other sorghum published maps and linkage groups were assigned according to the recent sorghum nomenclature. The phenotypic data of each trait was averaged to obtain the combined environmental means to detect QTLs. Composite interval mapping (CIM) of the PLABQTL software was used for analysis. The model AA was used to estimate additive x additive QTL interaction, and the command environ was used to estimate significant QTL in each environment and to detect QTL x environmental interaction. CIM detected a few additive QTLs per trait across the 10 sorghum chromosomes. The additive QTLs co-localized and showed clusters on some chromosome, where the most prominent clusters were observed on chromosomes 1, 2, 6, and 7. QTLs showing multiple effects were observed on all chromosomes except chromosome 5 and 10 and most of these pleiotropic QTLs showed the contribution of the positive allele in each trait. For example, a QTL on SBI-06, position 14, showed the allele contribution of M71 (high cellulose and compact panicle parent) for cellulose and fresh panicle weight, along with the allele contribution of SS79 (high sucrose and biomass parent) for sucrose and fresh biomass. However, QTLs that showed pleiotropic effect also exhibited QTL x environment interaction. Most of the QTLs detected for agronomic biomass (i.e. stripped stalk mass, fresh leaf mass and fresh biomass) interacted significantly with the environment. Additive x additive interaction was detected in 14 out of 20 traits studied totalling 41 digenic pairs. This indicates the abundance of digenic epistasis across the genome. However, failure to detect QTL x QTL interaction can never be interpreted as being absent but can be attributed to QTL x environment interaction, sampling error or accuracy of statistic methods. Among other characteristics, stem sugar components and lignocellulosic biomass are the important characteristics of Sorghum bicolor to determine it as one of the major bionergy crops. There is limited literature on sugar and fibre related traits in sorghum at the moment, and this study among other few, serves as one of the base investigations on the topic. The phenotypic variation of the lines, the moderate to high heritability of the desired traits, the significant correlation and pleiotropic QTLs to the positive direction identified, the additive stable QTLs and their interactions detected suggest that it will be possible to develop improved lines and hybrids by breeding and to ultimately enhance the sorghum plant as a bioenergy crop. Die steigende Nachfrage nach Nahrungsmitteln und regenerativen Energieträgern sowie die Sorge um den fortschreitenden Klimawandel machen Sorghum zu einer hochinteressanten landwirtschaftlichen Pflanze. Seine hauptsächlichen Vorzüge liegen in seiner gleichzeitigen Verwendbarkeit als Nahrungs- und Futtermittel, als Energieträger und als Rohstofflieferant für industrielle Erzeugnisse. Die Art Sorghum umfaßt sowohl Korn- als auch Zuckertypen, wobei sich Zuckersorghum sich gegenüber anderen Zucker- und Biomassekulturen durch einige wesentliche Eigenschaften auszeichnet: Erstens zeigt Sorghum eine enorme genetische Variabilität, die sogar dürre- und hitzetolerate Genotypen umfaßt. Zweitens verringert die gleichzeitige Nutzung des Sorghums für Nahrungs- und Futtermittelgewinnung aus dem Korn und Biokraftstoff aus den Stängeln den Wettbewerb um die Landnutzung durch verschiedene Pflanzen für unterschiedliche Zwecke. In der vorliegenden Arbeit wurden Feldversuche und genetische Analysen mit dem Ziel durchgeführt, hautverantwortliche chromosomale Regionen zu identifizieren, die mit Eigenschaften gekoppelt sind, welche die Eignung zur Gewinnung von Biokraftstoff beeinflussen. Eine Versuchspopulation von 213 RILs (rekombinante Inzuchtlinien) aus einer Kreuzung zwischen Kornsorghum (M71) und Zuckersorghum (SS79) wurde in fünf Umwelten in Deutschland während der Vegetationsperioden 2007 und 2008 gepflanzt. Diese Population spaltet in Allele, die den Biomasseertrag, den Saftgehalt, den Zuckergehalt, den Fasergehalt, sowie viele weitere Eigenschaften kontrollieren. Der Ansatz, diese chromosomalen Regionen zu identifizieren, umfaßte (1) die Phänotypisierung der RIL-Population für die genannten Merkmale, (2) die Erstellung einer genetische Karte der RIL-Linien und (3) die Korrelation der phänotyosichen Daten mit den genotypischen Daten.Die phänotypischen Daten wurden für einzelne Standorte und Umwelten mit SAS® Version 9.1 im gemischten Model (mixed model) analysiert. Die Genotypen wurden dabei als fixe Effekte betrachtet, während die übrigen Komponenten als zufällige Effekte behandelt wurden.In Hinblick auf zuckerrelevante Eigenschaften analysierten wir den Brix, den Sucrose- und Glucosegehalt, den Zuckergehalt, den Stängeldurchmesser, die Stängelsaftmenge, die Pflanzenhöhe, das Frischgewicht der Rispen, die Zeitspanne bis zur Blüte und die Anzahl der Triebe pro Pflanze. Die Varianz-Komponenten unterteilten wir in Genotyp, Standort, Jahr, Genotyp x Standort, Standort x Jahr, Genotyp x Standort x Jahr und Erntedaten. Die Ergebnisse zeigen, daß der Standort trotz seines großen Einflusses auf den Phänotyp keine signifikante Quelle für die Varianz ist. Mit Ausnahme von Stängelsaftmenge, Frischgewicht der Rispe und Glucosegehalt wurden Wechselwirkungen zwischen Genotyp und Umwelt für alle Merkmale beobachtet. Die genetische Varianz beinflußt die phänotypische Varianz zwar wenig, jedoch signifikant. Die geschätzten Heritabilitäten waren daher mittel bis hoch in allen Merkmalen mit Ausnahme der Stängelsaftmenge (H2=0.18). Pearsons Korrelationsschätzungen zeigen, daß diese Merkmale auch untereinander signifikant korreliert sind. So korrelieren die Blütedaten mit allen Merkmalen. Die Anzahl der Stängel pro Pflanze ist dagegen zwar positiv mit der Stängelsaftmenge korreliert, aber negativ mit dem Zuckergehalt und Blütedaten.Bezüglich der Fasermerkmale wurden Fasergehalte, wie Säure-detergenfaser, Neutral-detergenzfaser, Säure-detergenzlignin, Zellulose und Hemizellulose und pflanzenbauliche Merkmale, wie Blattfrischmasse, Geschälte-Stängel-Masse, Stängeltrockenmasse, Gesamtfrischbio¬masse, Gesamttrockenbiomasse bestimmt. Die Varianz-Komponenten unterteilten wir in Genotyp, Umwelt, Erntezeitpunkt, Genotyp x Umwelt. Die genetische Varianz ist signifikant in allen Merkmalen mit Ausnahme der Gesamttrockenbiomasse. Die Umwelt ist hingegen in keinem der Merkmale signifikant. Die Genotyp-Umwelt-Wechselwirkung ist für alle Merkmale außer Hemizellulose signifikant. Stängeltrockenmasse, Gesamttrockenbiomasse und Hemizellulose weisen niedrige Heritabilitätsschätzungen von <=0.24 auf. Pearsons Korrelationsschätzungen zeigen, daß sowohl jeweils Fasergehaltmerkmale als auch pflanzenbauliche Merkmale signifikant und positiv untereinander korreliert sind. Zwischen diesen Merkmalsgruppen werden allerdings negative Korrelationen beobachtet.Eine genetische Kopplungskarte wurde aus den genotypischen Daten von 213 RILs und 157 Markern, wie simple sequence repeats (SSR), amplified fragment length polymorphism (AFLP) and expressed sequence tag-simple sequence repeats (EST-SSR) erstellt. Diese Karte wurde mit anderen publizierten Sorghum-Karten verglichen und die Kopplungsgruppen entsprechend der aktuellen Sorghum-Nomenklatur bezeichnet. Die phänotypischen Daten wurden gemittelt, um mittels der kombinierten Umweltdaten QTLs zu detektieren. Die QTL-Analyse wurde mit der Software PLABQTL unter Nutzung der Composite interval mapping - Methode (CIM) durchgeführt. Additive x additive QTL-Wechselwirkungen wurden mit dem Model AA geschätzt. Der Befehl environ wurde angewandt, um signifikante QTLs in jeder Umwelt zu ermitteln und QTL-Umwelt-Wechselwirkungen zu ermitteln. Es konnten einige additive QTLs pro Merkmal in den 10 Sorghumchromosomen detektiert werden. Die additiven QTL ko-lokalisieren und bilden auf einigen Chromosomen Cluster. Die bedeutensten Cluster finden sich auf den Chromosomen 1,2,6 und 7. Pleitrope QTLs konnten auf allen Chromosomen außer Chromosom 5 und 10 beobachtet werden. Die meisten von ihnen beeinflussen das jeweils positive Allel der betroffenen Merkmale. Zum Bespiel zeigt ein QTL auf SBI-06, Position 14, den allelischen Einfluß des Elters M71 (Hoch-Zellulose und kompakte Rispe) auf Zellulose und Rispenfrischmasse, sowie den allelischen Einfluß des Elters SS79 (Hoch-Sucrose, Hoch-Biomasse) auf Sucrose und Gesamtfrischbiomasse. QTLs mit pleitropen Effekten zeigen darüberhinaus QTL-Umwelt-Wechselwirkungen. Die meisten ermittelten QTL für pflanzenbauliche Merkmale (z.B. Geschälte-Stängel-Masse, Blattfrischmasse, Gesamtfrischbiomasse) zeigen signifikante Wechselwirkungen mit der Umwelt. In 14 von 20 untersuchten Merkmalen wurden Additive x Additive Wechselwirkungen beobachtet, insgesamt sind es 41 digenische Paare. Dies verdeutlicht die Verbreitung digenischer Epistasie über das Genom. Dessenungeachtet darf eine nicht nachgewiesene QTL-QTL-Wechselwirkung nicht ignoriert werden. Sie kann auch der Genauigkeit oder Fehlerhaftigkeit der statistischen Methoden und Probennahme, oder aber einer möglichen QTL-Umwelt-Wechselwirkung zugeschreiben werden.Neben anderen Charakteristika sind Stängelzuckerkomponenten und Lignozellulose-Biomasse diejenigen Merkmale von Sorghum bicolor, die zur seiner Bedeutung als Bioenergieträger beitragen. Gegenwärtig ist nur wenig zu Zucker- und Fasermerkmalen in Sorghum publiziert. Die vorliegende Arbeit dient neben anderen wenigen der Grundlagenforschung auf diesem Gebiet. Die phänotypische Varianz der Linien, die mittlere bis hohe Heritabilität der gewünschten Merkmale, die signifikante Korrelation, sowie pleitrope QTL bezüglich der positiven Allele zusammen mit stabilen QTLs und ihren Wechselwirkungen lassen erwarten, daß die Entwicklung verbesserter Linien und Hybriden durch Züchtung möglich ist, wodurch Sorhum als Bioenergiepflanze gestärkt wird.