A mixed model QTL analysis for sugarcane multiple-harvest-location trial data

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.

Identificazione di un QTL principale per resistenza a ruggine bruna sul cromosoma 7B di frumento duro

Leaf rust caused by Puccinia triticina is a serious disease of durum wheat (Triticum durum) worldwide. However, genetic and molecular mapping studies aimed at characterizing leaf rust resistance genes in durum wheat have been only recently undertaken. The Italian durum wheat cv. Creso shows a high level of resistance to P. triticina that has been considered durable and that appears to be due to a combination of a single dominant gene and one or more additional factors conferring partial resistance. In this study, the genetic basis of leaf rust resistance carried by Creso was investigated using 176 recombinant inbred lines (RILs) from the cross between the cv. Colosseo (C, leaf rust resistance donor) and Lloyd (L, susceptible parent). Colosseo is a cv. directly related to Creso with the leaf rust resistance phenotype inherited from Creso, and was considered as resistance donor because of its better adaptation to local (Emilia Romagna, Italy) cultivation environment. RILs have been artificially inoculated with a mixture of 16 Italian P. triticina isolates that were characterized for virulence to seedlings of 22 common wheat cv. Thatcher isolines each carrying a different leaf rust resistance gene, and for molecular genotypes at 15 simple sequence repeat (SSR) loci, in order to determine their specialization with regard to the host species. The characterization of the leaf rust isolates was conducted at the Cereal Disease Laboratory of the University of Minnesota (St. Paul, USA) (Chapter 2). A genetic linkage map was constructed using segregation data from the population of 176 RILs from the cross CL. A total of 662 loci, including 162 simple sequence repeats (SSRs) and 500 Diversity Arrays Technology markers (DArTs), were analyzed by means of the package EasyMap 0.1. The integrated SSR-DArT linkage map consisted of 554 loci (162 SSR and 392 DArT markers) grouped into 19 linkage blocks with an average marker density of 5.7 cM/marker. The final map spanned a total of 2022 cM, which correspond to a tetraploid genome (AABB) coverage of ca. 77% (Chapter 3). The RIL population was phenotyped for their resistance to leaf rust under artificial inoculation in 2006; the percentage of infected leaf area (LRS, leaf rust susceptibility) was evaluated at three stages through the disease developmental cycle and the area under disease progress curve (AUDPC) was then calculated. The response at the seedling stage (infection type, IT) was also investigated. QTL analysis was carried out by means of the Composite Interval Mapping method based on a selection of markers from the CL map. A major QTL (QLr.ubo-7B.2) for leaf rust resistance controlling both the seedling and the adult plant response, was mapped on the distal region of chromosome arm 7BL (deletion bin 7BL10-0.78-1.00), in a gene-dense region known to carry several genes/QTLs for resistance to rusts and other major cereal fungal diseases in wheat and barley. QLr.ubo-7B.2 was identified within a supporting interval of ca. 5 cM tightly associated with three SSR markers (Xbarc340.2, Xgwm146 e Xgwm344.2), and showed an R2 and an LOD peak value for the AUDPC equal to 72.9% an 44.5, respectively. Three additional minor QTLs were also detected (QLr.ubo-7B.1 on chr. 7BS; QLr.ubo-2A on chr. 2AL and QLr.ubo-3A on chr. 3AS) (Chapter 4). The presence of the major QTL (QLr.ubo-7B.2) was validated by a linkage disequilibrium (LD)-based test using field data from two different plant materials: i) a set of 62 advanced lines from multiple crosses involving Creso and his directly related resistance derivates Colosseo and Plinio, and ii) a panel of 164 elite durum wheat accessions representative of the major durum breeding program of the Mediterranean basin. Lines and accessions were phenotyped for leaf rust resistance under artificial inoculation in two different field trials carried out at Argelato (BO, Italy) in 2006 and 2007; the durum elite accessions were also evaluated in two additional field experiments in Obregon (Messico; 2007 and 2008) and in a green-house experiment (seedling resistance) at the Cereal Disease Laboratory (St. Paul, USA, 2008). The molecular characterization involved 14 SSR markers mapping on the 7BL chromosome region found to harbour the major QTL. Association analysis was then performed with a mixed-linear-model approach. Results confirmed the presence of a major QTL for leaf rust resistance, both at adult plant and at seedling stage, located between markers Xbarc340.2, Xgwm146 and Xgwm344.2, in an interval that coincides with the supporting interval (LOD-2) of QLr.ubo-7B.2 as resulted from the RIL QTL analysis. (Chapter 5). The identification and mapping of the major QTL associated to the durable leaf rust resistance carried by Creso, together with the identification of the associated SSR markers, will enhance the selection efficiency in durum wheat breeding programs (MAS, Marker Assisted Selection) and will accelerate the release of cvs. with durable resistance through marker-assisted pyramiding of the tagged resistance genes/QTLs most effective against wheat fungal pathogens.

Genetische Kartierung und QTL-Analyse agronomischer Eigenschaften der Weinrebe unter besonderer Berücksichtigung der Pilzresistenz

153 Nachkommen einer Kreuzung aus der pilzresistenten Rebsorte ‘Regent‘ und ‘Lemberger‘ als klassischer pilzsensitiver Sorte zeigen quantitative Merkmalsvariation bezüglich der Resistenz gegen Plasmopara viticola und Uncinula necator sowie für weitere Eigenschaften, die z.B. das Eintreten der Beerenreife betreffen. Auf dem Weg über die genetische Kartierung mit molekularen Markern und der Lokalisierung von QTL-Effekten konnten Hinweise auf weinbaulich relevante Genomregionen gewonnen werden; dies liefert z.B. die Basis für markergestützte Selektion bei Zuchtvorhaben mit dem Resistenzträger ‘Regent’ (vgl. auch FISCHER et al., 2004). Ein Major-QTL für die Resistenz gegen den Echten Mehltau Uncinula necator sowie zwei Major QTL für die Resistenz gegen den Erreger des Falschen Mehltau, Plasmopara viticola, traten mit hoher Signifikanz auf drei verschiedenen Kopplungsgruppen von ‘Regent‘ auf. Auch Regionen mit Relevanz für das Eintreten der Beerenreife wurden beschrieben. Über die Isolierung, Sequenzierung und anschließende Analyse einzelner Markerfragmente mit Methoden der Bioinformatik ist es gelungen, ein putatives T10P12.4-Ortholog der Weinrebe (ein thioredoxinähnliches Protein) in enger Kopplung zu einem Major-QTL-Maximum für Plasmopara viticola-Resistenz zu identifizieren, das als Kandidat für die Beteiligung an der Pathogenantwort in Frage kommt. Es konnte exemplarisch gezeigt werden, dass die eingesetzten Methoden der Kartierung und QTL-Analyse unter Verwendung PCR-basierter Markertypen wie SSR und AFLP und einer beschleunigten Analyse über computergestützte Kapillargelelektrophorese in vertretbarem Zeitrahmen bis zur Isolation potentieller Schlüsselgene führen können. Die grundsätzliche Eignung der QTL-Analyse als effizientes Werkzeug gezielter Züchtungsplanung für den Weinbau bestätigte sich. Ihre Anwendung im Rahmen der vorliegenden Dissertation hat die Basis für die Nutzung von QTL-Information bei dem Vergleich etablierter und der Entwicklung neuer Sorten gelegt und zum Verständnis von Prozessen beigetragen, die den betrachteten Eigenschaften wie der Pilzresistenz möglicherweise zu Grunde liegen. Ein großer Teil der gewonnenen Daten bringt auch die Untersuchungen anderer Kultivare voran und ist intervarietal übertragbar. Darüber hinaus haben sich Chancen für vergleichende Studien zwischen der Weinrebe einerseits und der Modellpflanze Arabidopsis thaliana sowie weiteren Kulturpflanzen andererseits abgezeichnet. Die Hinweise auf die zentrale Rolle und universelle Natur des Redox-Signalling haben interessante Perspektiven zum Verständnis organismenübergreifender physiologischer Zusammenhänge eröffnet. Dies betrifft z.B. auch die Reaktion auf Verwundung oder die Pathogenantwort.

A multiparental cross population for mapping QTL for relevant agronomic traits in durum wheat (Triticum durum Desf.)

Multiparental cross designs for mapping quantitative trait loci (QTL) in crops are efficient alternatives to conventional biparental experimental populations because they exploit a broader genetic basis and higher mapping resolution. We describe the development and deployment of a multiparental recombinant inbred line (RIL) population in durum wheat (Triticum durum Desf.) obtained by crossing four elite cultivars characterized by different traits of agronomic value. A linkage map spanning 2,663 cM and including 7,594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs with a wheat-dedicated 90k SNP chip. A cluster file was developed for correct allele calling in the framework of the tetraploid durum wheat genome. Based on phenotypic data collected over four field experiments, a multi-trait quantitative trait loci (QTL) analysis was carried out for 18 traits of agronomic relevance (including yield, yield-components, morpho-physiological and seed quality traits). Across environments, a total of 63 QTL were identified and characterized in terms of the four founder haplotypes. We mapped two QTL for grain yield across environments and 23 QTL for grain yield components. A novel major QTL for number of grain per spikelet/ear was mapped on chr 2A and shown to control up to 39% of phenotypic variance in this cross. Functionally different QTL alleles, in terms of direction and size of genetic effect, were distributed among the four parents. Based on the occurrence of QTL-clusters, we characterized the breeding values (in terms of effects on yield) of most of QTL for heading and maturity as well as yield component and quality QTL. This multiparental RIL population provides the wheat community with a highly informative QTL mapping resource enabling the dissection of the genetic architecture of multiple agronomic relevant traits in durum wheat.

Estudio genético y molecular de la respuesta a la vernalización en zanahoria (Daucus carota L.)

La zanahoria es una planta bienal de estación fría que requiere de un período de vernalización para florecer. Sin embargo, algunos cultivares adaptados a zonas más cálidas requieren de menor vernalización y son clasificados como anuales o de floración temprana. El objetivo de esta tesis fue determinar la base genética e identificar los genes y/o regiones cromosómicas involucradas en los requerimientos de vernalización en la zanahoria. Para ello fueron evaluadas a campo familias segregantes F1, F2, F3, RC1 y RC2 obtenidas a partir de un cruzamiento entre una planta anual y una bienal. En base a los patrones de segregación observados se concluyó que la anualidad, o bajos requerimientos de vernalización, estaría determinada por un gen simple dominante. Al evaluar introducciones de zanahoria anuales y bienales de diversos orígenes geográficos y sus cruzamientos, se volvió a observar la total dominancia de la anualidad y se encontró variabilidad en el ciclo entre materiales anuales y entre materiales bienales. Utilizando un método molecular que se basa en similitud completa (BLAST), no se encontró en el genoma de la zanahoria secuencias homólogas al gen FLC, el cual juega un rol central en la respuesta a la vernalización en Arabidopsis y otras especies como las Brassicas. Mediante la técnica de mapeo se encontró una región cromosómica ligada a la respuesta a la vernalización en zanahoria. La misma se localizó en un grupo de ligamiento con 78 marcadores moleculares a una distancia de 0,69 cM y de 0,79 cM de los marcadores más cercanos. Este mapa servirá como base para el desarrollo de marcadores moleculares ligados al carácter y en un futuro para el mapeo físico y secuenciación de la región de interés utilizando una librería génica de BACs de zanahoria.

Detection of growth-related QTLs in turbot (Scophtalmus maximux)

Background The turbot (Scophthalmus maximus) is a highly appreciated European aquaculture species. Growth related traits constitute the main goal of the ongoing genetic breeding programs of this species. The recent construction of a consensus linkage map in this species has allowed the selection of a panel of 100 homogeneously distributed markers covering the 26 linkage groups (LG) suitable for QTL search. In this study we addressed the detection of QTL with effect on body weight, length and Fulton's condition factor. Results Eight families from two genetic breeding programs comprising 814 individuals were used to search for growth related QTL using the panel of microsatellites available for QTL screening. Two different approaches, maximum likelihood and regression interval mapping, were used in order to search for QTL. Up to eleven significant QTL were detected with both methods in at least one family: four for weight on LGs 5, 14, 15 and 16; five for length on LGs 5, 6, 12, 14 and 15; and two for Fulton's condition factor on LGs 3 and 16. In these LGs an association analysis was performed to ascertain the microsatellite marker with the highest apparent effect on the trait, in order to test the possibility of using them for marker assisted selection. Conclusions The use of regression interval mapping and maximum likelihood methods for QTL detection provided consistent results in many cases, although the high variation observed for traits mean among families made it difficult to evaluate QTL effects. Finer mapping of detected QTL, looking for tightly linked markers to the causative mutation, and comparative genomics are suggested to deepen in the analysis of QTL in turbot so they can be applied in marker assisted selection programs.

Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid

The genetic basis of heterosis was investigated in an elite rice hybrid by using a molecular linkage map with 150 segregating loci covering the entire rice genome. Data for yield and three traits that were components of yield were collected over 2 years from replicated field trials of 250 F2:3 families. Genotypic variations explained from about 50% to more than 80% of the total variation. Interactions between genotypes and years were small compared with the main effects. A total of 32 quantitative trait loci (QTLs) were detected for the four traits; 12 were observed in both years and the remaining 20 were detected in only one year. Overdominance was observed for most of the QTLs for yield and also for a few QTLs for the component traits. Correlations between marker heterozygosity and trait expression were low, indicating that the overall heterozygosity made little contribution to heterosis. Digenic interactions, including additive by additive, additive by dominance, and dominance by dominance, were frequent and widespread in this population. The interactions involved large numbers of marker loci, most of which individually were not detectable on single-locus basis; many interactions among loci were detected in both years. The results provide strong evidence that epistasis plays a major role as the genetic basis of heterosis.

Comparative organization of cattle chromosome 5 revealed by comparative mapping by annotation and sequence similarity and radiation hybrid mapping

A whole genome cattle-hamster radiation hybrid cell panel was used to construct a map of 54 markers located on bovine chromosome 5 (BTA5). Of the 54 markers, 34 are microsatellites selected from the cattle linkage map and 20 are genes. Among the 20 mapped genes, 10 are new assignments that were made by using the comparative mapping by annotation and sequence similarity strategy. A LOD-3 radiation hybrid framework map consisting of 21 markers was constructed. The relatively low retention frequency of markers on this chromosome (19%) prevented unambiguous ordering of the other 33 markers. The length of the map is 398.7 cR, corresponding to a ratio of ≈2.8 cR5,000/cM. Type I genes were binned for comparison of gene order among cattle, humans, and mice. Multiple internal rearrangements within conserved syntenic groups were apparent upon comparison of gene order on BTA5 and HSA12 and HSA22. A similarly high number of rearrangements were observed between BTA5 and MMU6, MMU10, and MMU15. The detailed comparative map of BTA5 should facilitate identification of genes affecting economically important traits that have been mapped to this chromosome and should contribute to our understanding of mammalian chromosome evolution.

A locus on chromosome 7 determines myocardial cell necrosis and calcification (dystrophic cardiac calcinosis) in mice.

Dystrophic cardiac calcinosis, an age-related cardiomyopathy that occurs among certain inbred strains of mice, involves myocardial injury, necrosis, and calcification. Using a complete linkage map approach and quantitative trait locus analysis, we sought to identify genetic loci determining dystrophic cardiac calcinosis in an F2 intercross of resistant C57BL/6J and susceptible C3H/HeJ inbred strains. We identified a single major locus, designated Dyscalc, located on proximal chromosome 7 in a region syntenic with human chromosomes 19q13 and 11p15. The statistical significance of Dyscalc (logarithm of odds score 14.6) was tested by analysis of permuted trait data. Analysis of BxH recombinant inbred strains confirmed the mapping position. The inheritance pattern indicated that this locus influences susceptibility of cells both to enter necrosis and to subsequently undergo calcification.

Identification of a major locus conferring resistance to powdery mildew (Erysiphe polygoni DC) in mungbean (Vigna radiata L. Wilczek) by QTL analysis

A major locus conferring resistance to the causal organism of powdery mildew, Erysiphe polygoni DC,, in mungbean (Vigna radiata L. Wilczek) was identified using QTL analysis with a population of 147 recombinant inbred individuals. The population was derived from a cross between 'Berken', a highly susceptible variety, and ATF 3640, a highly resistant line. To test for response to powdery mildew, F-7 and F-8 lines were inoculated by dispersing decaying mungbean leaves with residual conidia of E. polygoni amongst the young plants to create an artificial epidemic and assayed in a glasshouse facility. To generate a linkage map, 322 RFLP clones were tested against the two parents and 51 of these were selected to screen the mapping population. The 51 probes generated 52 mapped loci, which were used to construct a linkage map spanning 350 cM of the mungbean genome over 10 linkage groups. Using these markers, a single locus was identified that explained up to a maximum of 86% of the total variation in the resistance response to the pathogen.

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class (Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue (Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.

Resistance gene analogues in sugarcane and sorghum and their association with quantitative trait loci for rust resistance

Fifty-four different sugarcane resistance gene analogue (RGA) sequences were isolated, characterized, and used to identify molecular markers linked to major disease-resistance loci in sugarcane. Ten RGAs were identified from a sugarcane stem expressed sequence tag (EST) library; the remaining 44 were isolated from sugarcane stem, leaf, and root tissue using primers designed to conserved RGA motifs. The map location of 31 of the RGAs was determined in sugarcane and compared with the location of quantitative trait loci (QTL) for brown rust resistance. After 2 years of phenotyping, 3 RGAs were shown to generate markers that were significantly associated with resistance to this disease. To assist in the understanding of the complex genetic structure of sugarcane, 17 of the 31 RGAs were also mapped in sorghum. Comparative mapping between sugarcane and sorghum revealed syntenic localization of several RGA clusters. The 3 brown rust associated RGAs were shown to map to the same linkage group (LG) in sorghum with 2 mapping to one region and the third to a region previously shown to contain a major rust-resistance QTL in sorghum. These results illustrate the value of using RGAs for the identification of markers linked to disease resistance loci and the value of simultaneous mapping in sugarcane and sorghum.

Analysing diversity in sugarcane resistance gene analogues

As resistance genes have been shown to contain conserved motifs and cluster in many plant genomes, the identification of resistance gene analogues can be used as a strategy for both the discovery of DNA markers linked to disease resistance loci and the map-based cloning of disease resistance genes. Sugarcane suffers from many important diseases and an analysis of resistance gene analogues offers a means to identify DNA markers linked to resistance loci. However, sugarcane has the most complex genome of any crop plant and initially it is important to understand the extent of resistance gene analogue diversity in the sugarcane genome before genetic analysis. We review herein how more than 100 expressed sequence tags with homology to different resistance genes have been identified in sugarcane with many mapped as single-dose restriction fragment length polymorphism markers. Importantly, some of these resistance gene analogues have been shown to be linked to disease resistance genes or disease quantitative trait loci. In an attempt to more efficiently analyse additional resistance gene analogues in sugarcane, we report on experiments aimed at investigating the molecular diversity of several resistance gene analogue families using a modified form of a technique termed Ecotilling. Using Ecotilling, we were able to rapidly detect single nucleotide polymorphisms in fragments amplified by PCR from four different resistance gene analogue families, SoRP1D, SoPTO, SoXa21 and SoHs1pro-1. An analysis of a diverse set of sugarcane varieties, including modern sugarcane cultivars and several S. officinarum and S. spontaneum clones, indicated that all amplicons, apart from SoHs1pro-1, contained significant polymorphism within the gene region studied. However, a comparison among these sugarcane clones, including between the parents of two sugarcane mapping populations, indicated that most polymorphisms were multi-dose, not single-dose, preventing their genetic map location or association with disease susceptibility or resistance from being determined.

Assessing the taxonomic status of dingoes Canis familiaris dingo for conservation

1. The conservation status of the dingo Canis familiaris dingo is threatened by hybridization with the domestic dog C. familiaris familiaris. A practical method that can estimate the different levels of hybridization in the field is urgently required so that animals below a specific threshold of dingo ancestry (e.g. 1/4 or 1/2 dingoes) can reliably be identified and removed from dingo populations. 2. Skull morphology has been traditionally used to assess dingo purity, but this method does not discriminate between the different levels of dingo ancestry in hybrids. Furthermore, measurements can only be reliably taken from the skulls of dead animals. 3. Methods based on the analysis of variation in DNA are able to discriminate between the different levels of hybridization, but the validity of this method has been questioned because the materials currently used as a reference for dingoes are from captive animals of unproven genetic purity. The use of pre-European materials would improve the accuracy of this method, but suitable material has not been found in sufficient quantity to develop a reliable reference population. Furthermore, current methods based on DNA are impractical for the field-based discrimination of hybrids because samples require laboratory analysis. 4. Coat colour has also been used to estimate the extent of hybridization and is possibly the most practical method to apply in the field. However, this method may not be as powerful as genetic or morphological analyses because some hybrids (e.g. Australian cattle dog x dingo) are similar to dingoes in coat colour and body form. This problem may be alleviated by using additional visual characteristics such as the presence/absence of ticking and white markings.

DNA markers linked to yield, yield components, and morphological traits in autotetraploid lucerne (Medicago sativa L.)

We have mapped and identifed DNA markers linked to morphology, yield, and yield components of lucerne, using a backcross population derived from winter-active parents. The high-yielding and recurrent parent, D, produced individual markers that accounted for up to 18% of total yield over 6 harvests, at Gatton, south-eastern Queensland. The same marker, AC/TT8, was consistently identified at each individual harvest, and in individual harvests accounted for up to 26% of the phenotypic variation for yield. This marker was located in linkage group 2 of the D map, and several other markers positively associated with yield were consistently identified in this linkage group. Similarly, markers negatively associated with yield were consistently identified in the W116 map, W116 being the low-yielding parent. Highly significant positive correlations were observed between total yield and yield for harvests 1-6, and between total yield and stem length, tiller number, leaf yield/plant, leaf yield/5 stems, stem yield/plant, and stem yield/5 stems. Highly significant QTL were located for all these characters as well as for leaf shape and pubescence.