Detection and use of QTL for complex traits in multiple environments.

QTL mapping methods for complex traits are challenged by new developments in marker technology, phenotyping platforms, and breeding methods. In meeting these challenges, QTL mapping approaches will need to also acknowledge the central roles of QTL by environment interactions (QEI) and QTL by trait interactions in the expression of complex traits like yield. This paper presents an overview of mixed model QTL methodology that is suitable for many types of populations and that allows predictive modeling of QEI, both for environmental and developmental gradients. Attention is also given to multi-trait QTL models which are essential to interpret the genetic basis of trait correlations. Biophysical (crop growth) model simulations are proposed as a complement to statistical QTL mapping for the interpretation of the nature of QEI and to investigate better methods for the dissection of complex traits into component traits and their genetic controls.

A new semi-dwarfing gene identified by molecular mapping of quantitative trait loci in barley.

Semi-dwarfing genes have been widely used in spring barley (Hordeum vulgare L.) breeding programs in many parts of the world, but the success in developing barley cultivars with semi-dwarfing genes has been limited in North America. Exploiting new semi-dwarfing genes may help in solving this dilemma. A recombinant inbred line population was developed by crossing ZAU 7, a semi-dwarf cultivar from China, to ND16092, a tall breeding line from North Dakota. To identify quantitative trait loci (QTL) controlling plant height, a linkage map comprised of 111 molecular markers was constructed. Simple interval mapping was performed for each of the eight environments. A consistent QTL for plant height was found on chromosome 7HL. This QTL is not associated with maturity and rachis internode length. We suggest the provisional name Qph-7H for this QTL. Qph-7H from ZAU 7 reduced plant height to about 3/4 of normal; thus, Qph-7H is considered a semi-dwarfing gene. Other QTLs for plant height were found, but their expression was variable across the eight environments tested.

Leaf trait co-ordination in relation to construction cost, carbon gain and resource-use efficiency in exotic invasive and native woody vine species

Background and Aims: Success of invasive plant species is thought to be linked with their higher leaf carbon fixation strategy, enabling them to capture and utilize resources better than native species, and thus pre-empt and maintain space. However, these traits are not well-defined for invasive woody vines. Methods: In a glass house setting, experiments were conducted to examine how leaf carbon gain strategies differ between non-indigenous invasive and native woody vines of south-eastern Australia, by investigating their biomass gain, leaf structural, nutrient and physiological traits under changing light and moisture regimes. Key Results: Leaf construction cost (CC), calorific value and carbon : nitrogen (C : N) ratio were lower in the invasive group, while ash content, N, maximum photosynthesis, light-use efficiency, photosynthetic energyuse efficiency (PEUE) and specific leaf area (SLA) were higher in this group relative to the native group. Trait plasticity, relative growth rate (RGR), photosynthetic nitrogen-use efficiency and water-use efficiency did not differ significantly between the groups. However, across light resource, regression analyses indicated that at a common (same) leaf CC and PEUE, a higher biomass RGR resulted for the invasive group; also at a common SLA, a lower CC but higher N resulted for the invasive group. Overall, trait co-ordination (using pair-wise correlation analyses) was better in the invasive group. Ordination using 16 leaf traits indicated that the major axis of invasive-native dichotomy is primarily driven by SLA and CC (including its components and/or derivative of PEUE) and was significantly linked with RGR. Conclusions: These results demonstrated that while not all measures of leaf resource traits may differ between the two groups, the higher level of trait correlation and higher revenue returned (RGR) per unit of major resource need (CC) and use (PEUE) in the invasive group is in line with their rapid spread where introduced.

Understanding invasion history: genetic structure and diversity of two globally invasive plants and implications for their management

Aim: Resolving the origin of invasive plant species is important for understanding the introduction histories of successful invaders and aiding strategies aimed at their management. This study aimed to infer the number and origin(s) of introduction for the globally invasive species, Macfadyena unguis-cati and Jatropha gossypiifolia using molecular data. Location: Native range: Neotropics; Invaded range: North America, Africa, Europe, Asia, Pacific Islands and Australia. Methods: We used chloroplast microsatellites (cpSSRs) to elucidate the origin(s) of introduced populations and calculated the genetic diversity in native and introduced regions. Results: Strong genetic structure was found within the native range of M. unguis-cati, but no genetic structuring was evident in the native range of J. gossypiifolia. Overall, 27 haplotypes were found in the native range of M. unguis-cati. Only four haplotypes were found in the introduced range, with more than 96% of introduced specimens matching a haplotype from Paraguay. In contrast, 15 haplotypes were found in the introduced range of J. gossypiifolia, with all invasive populations, except New Caledonia, comprising multiple haplotypes. Main conclusions: These data show that two invasive plant species from the same native range have had vastly different introduction histories in their non-native ranges. Invasive populations of M. unguis-cati probably came from a single or few independent introductions, whereas most invasive J. gossypiifolia populations arose from multiple introductions or alternatively from a representative sample of genetic diversity from a panmictic native range. As introduced M. unguis-cati populations are dominated by a single haplotype, locally adapted natural enemies should make the best control agents. However, invasive populations of J. gossypiifolia are genetically diverse and the selection of bio-control agents will be considerably more complex.

A genome-wide association study of tick burden and milk composition in cattle.

To study the genetic basis of tick burden and milk production and their interrelationship, we collected a sample of 1961 cattle with multiple tick counts from northern Australia of which 973 had dairy production data in the Australian Dairy Herd Information Service database. We calculated heritabilities, genetic and phenotypic correlations for these traits and showed a negative relationship between tick counts and milk and milk component yield. Tests of polymorphisms of four genes associated with milk yield, ABCG2, DGAT1, GHR and PRLR, showed no statistically significant effect on tick burden but highly significant associations to milk component yield in these data and we confirmed separate effects for GHR and PRLR on bovine chromosome 20. To begin to identify some of the molecular genetic bases for these traits, we genotyped a sample of 189 of these cattle for 7397 single nucleotide polymorphisms in a genome-wide association study. Although the allele effects for adjusted milk fat and protein yield were highly correlated (r = 0.66), the correlations of allele effects of these milk component yields and tick burden were small (|r| <= 0.10). These results agree in general with the phenotypic correlations between tick counts and milk component yield and suggest that selection on markers for tick burden or milk component yield may have no undesirable effect on the other trait.

QTL mapping of multiple foliar disease and root-lesion nematode resistances in wheat.

A genetic linkage map, based on a cross between the synthetic hexaploid CPI133872 and the bread wheat cultivar Janz, was established using 111 F1-derived doubled haploid lines. The population was phenotyped in multiple years and/or locations for seven disease resistance traits, namely, Septoria tritici blotch (Mycosphaeralla graminicola), yellow leaf spot also known as tan spot (Pyrenophora tritici-repentis), stripe rust (Puccinia striiformis f. sp. tritici), leaf rust (Puccinia triticina), stem rust (Puccinia graminis f. sp. tritici) and two species of root-lesion nematode (Pratylenchyus thornei and P. neglectus). The DH population was also scored for coleoptile colour and the presence of the seedling leaf rust resistance gene Lr24. Implementation of a multiple-QTL model identified a tightly linked cluster of foliar disease resistance QTL in chromosome 3DL. Major QTL each for resistance to Septoria tritici blotch and yellow leaf spot were contributed by the synthetic hexaploid parent CPI133872 and linked in repulsion with the coincident Lr24Sr24/ locus carried by parent Janz. This is the first report of linked QTL for Septoria tritici blotch and yellow leaf spot contributed by the same parent. Additional QTL for yellow leaf spot were detected in 5AS and 5BL. Consistent QTL for stripe rust resistance were identified in chromosomes 1BL, 4BL and 7DS, with the QTL in 7DS corresponding to the Yr18Lr34/ region. Three major QTL for P. thornei resistance (2BS, 6DS, 6DL) and two for P. neglectus resistance (2BS, 6DS) were detected. The recombinants combining resistance to Septoria tritici blotch, yellow leaf spot, rust diseases and root-lesion nematodes from parents CPI133872 and Janz constitute valuable germplasm for the transfer of multiple disease resistance into new wheat cultivars.

Quantitative trait loci and epistatic interactions in barley conferring resistance to net type net blotch (Pyrenophora teres f. teres) isolates.

Net type net blotch (NTNB) is an important barley disease in Australia and elsewhere, with significant yield reduction. This trait is important in selection along with other traits of quality and agronomic value. Two-hundred doubled-haploid lines were generated through anther culture from a cross between 'Pompadour' and 'Stirling'. Quantitative trait loci (QTL) were identified against five isolates of Pyrenophora teres f. teres, which represent virulences across Australia. QTL were mapped on chromosomes 3H and 6H using simple sequence repeat (SSR) markers. The resistance locus on 6H was detected with all isolates while the 3H locus was detected with two isolates. The 6H QTL from 'Pompadour' contributed resistance to isolates 97NB1, 95NB100 and NB81, whereas 6H QTL from 'Stirling' contributed resistance to isolates NB50 and NB52B. The 3H QTL from 'Pompadour' contributed resistance to NB50 and NB52B. Significant epistatic interactions were detected between QTL on 3H and 6H. These resistance QTL are a useful resource and identifying closely linked SSR markers with allelic combinations will facilitate in marker-assisted selection to develop NTNB resistant breeding lines.

Quantitative and molecular genetics of juvenile wood traits in radiata and slash/Caribbean pines.

The Juvenile Wood Initiative (JWI) project has been running successfully since July 2003 under a Research Agreement with FWPA and Letters of Association with the consortium partners STBA (Southern Tree Breeding Association), ArborGen and FPQ (Forestry Plantations Queensland). Over the last five and half years, JWI scientists in CSIRO, FPQ, and STBA have completed all 12 major milestones and 28 component milestones according to the project schedule. We have made benchmark progress in understanding the genetic control of wood formation and interrelationships among wood traits. The project has made 15 primary scientific findings and several results have been adopted by industry as summarized below. This progress was detailed in 10 technical reports to funding organizations and industry clients. Team scientists produced 16 scientific manuscripts (8 published, 1 in press, 2 submitted, and several others in the process of submission) and 15 conference papers or presentations. Primary Scientific Findings. The 15 major scientific findings related to wood science, inheritance and the genetic basis of juvenile wood traits are: 1. An optimal method to predict stiffness of standing trees in slash/Caribbean pine is to combine gravimetric basic density from 12 mm increment cores with a standing tree prediction of MoE using a time of flight acoustic tool. This was the most accurate and cheapest way to rank trees for breeding selection for slash/Caribbean hybrid pine. This method was also recommended for radiata pine. 2. Wood density breeding values were predicted for the first time in the STBA breeding population using a large sample of 7,078 trees (increment cores) and it was estimated that selection of the best 250 trees for deployment will produce wood density gains of 12.4%. 3. Large genetic variation for a suite of wood quality traits including density, MFA, spiral grain, shrinkage, acoustic and non-acoustic stiffness (MoE) for clear wood and standing trees were observed. Genetic gains of between 8 and 49% were predicted for these wood quality traits with selection intensity between 1 to 10% for radiata pine. 4. Site had a major effect on juvenile-mature wood transition age and the effect of selective breeding for a shorter juvenile wood formation phase was only moderate (about 10% genetic gain with 10% selection intensity, equivalent to about 2 years reduction of juvenile wood). 5. The study found no usable site by genotype interactions for the wood quality traits of density, MFA and MoE for both radiata and slash/Caribbean pines, suggesting that assessment of wood properties on one or two sites will provide reliable estimates of the genetic worth of individuals for use in future breeding. 6. There were significant and sizable genotype by environment interactions between the mainland and Tasmanian regions and within Tasmania for DBH and branch size. 7. Strong genetic correlations between rings for density, MFA and MoE for both radiata and slash/Caribbean pines were observed. This suggests that selection for improved wood properties in the innermost rings would also result in improvement of wood properties in the subsequent rings, as well as improved average performance of the entire core. 8. Strong genetic correlations between pure species and hybrid performance for each of the wood quality traits were observed in the hybrid pines. Parental performance can be used to identify the hybrid families which are most likely to have superior juvenile wood properties of the slash/Caribbean F1 hybrid in southeast Queensland. 9. Large unfavourable genetic correlations between growth and wood quality traits were a prominent feature in radiata pine, indicating that overcoming this unfavourable genetic correlation will be a major technical issue in progressing radiata pine breeding. 10. The project created the first radiata pine 18 k cDNA microarray and generated 5,952 radiata pine xylogenesis expressed sequence tags (ESTs) which assembled into 3,304 unigenes. 11. A total of 348 genes were identified as preferentially expressed genes in earlywood or latewood while a total of 168 genes were identified as preferentially expressed genes in either juvenile or mature wood. 12. Juvenile earlywood has a distinct transcriptome relative to other stages of wood development. 13. Discovered rapid decay of linkage disequilibrium (LD) in radiata pine with LD decaying to approximately 50% within 1,700 base pairs (within a typical gene). A total of 913 SNPS from sequencing 177,380 base pairs were identified for association genetic studies. 14. 149 SNPs from 44 genes and 255 SNPs from a further 51 genes (total 95 genes) were selected for association analysis with 62 wood traits, and 30 SNPs were shortlisted for their significant association with variation of wood quality traits (density, MFA and MoE) with individual significant SNPs accounting for between 1.9 and 9.7% of the total genetic variation in traits. 15. Index selection using breeding objectives was the most profitable selection method for radiata pine, but in the long term it may not be the most effective in dealing with negative genetic correlations between wood volume and quality traits. A combination of economic and biological approaches may be needed to deal with the strong adverse correlation.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.

Location of major effect genes in sorghum (Sorghum bicolor (L.) Moench).

Major effect genes are often used for germplasm identification, for diversity analyses and as selection targets in breeding. To date, only a few morphological characters have been mapped as major effect genes across a range of genetic linkage maps based on different types of molecular markers in sorghum (Sorghum bicolor (L.) Moench). This study aims to integrate all available previously mapped major effect genes onto a complete genome map, linked to the whole genome sequence, allowing sorghum breeders and researchers to link this information to QTL studies and to be aware of the consequences of selection for major genes. This provides new opportunities for breeders to take advantage of readily scorable morphological traits and to develop more effective breeding strategies. We also provide examples of the impact of selection for major effect genes on quantitative traits in sorghum. The concepts described in this paper have particular application to breeding programmes in developing countries where molecular markers are expensive or impossible to access.

Negligible evidence for regional genetic population structure for two shark species Rhizoprionodon acutus (Rüppell, 1837) and Sphyrna lewini (Griffith & Smith, 1834) with contrasting biology

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n = 237) and the milk shark (Rhizoprionodon acutus, n = 207) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751–0.903, respectively; microsatellite loci, 0.038–0.047 respectively). Our results support the spatially homogeneous monitoring and management plan for shark species in Queensland, Australia.

Migration of green turtles (Chelonia mydas) from Australasian feeding grounds inferred from genetic analyses

Coastal seagrass habitats in tropical and subtropical regions support aggregations of resident green turtles (Chelonia mydas) from several genetically distinct breeding populations. Migration of individuals to their respective dispersed breeding sites provides a complex pattern of migratory connectivity among nesting and feeding habitats of this species. An understanding of this pattern is important in regions where the persistence of populations is under threat from anthropogenic impacts. The present study uses mitochondrial DNA and mixed-stock analyses to assess the connectivity among seven feeding grounds across the north Australian coast and adjacent areas and 17 genetically distinct breeding populations from the Indo-Pacific region. It was hypothesised that large and geographically proximate breeding populations would dominate at nearby feeding grounds. As expected, each sampled feeding area appears to support multiple breeding populations, with two aggregations dominated by a local breeding population. Geographic distance between breeding and feeding habitat strongly influenced whether a breeding population contributed to a feeding ground (wi = 0.654); however, neither distance nor size of a breeding population was a good predictor of the extent of their contribution. The differential proportional contributions suggest the impact of anthropogenic mortality at feeding grounds should be assessed on a case-by-case basis.

Towards better management of Australia’s shark fishery: genetic analyses reveal unexpected ratios of cryptic blacktip species Carcharhinus tilstoni and C. limbatus

The common blacktip shark (Carcharhinus limbatus) and the Australian blacktip shark (C. tilstoni) are morphologically similar species that co-occur in subtropical and tropical Australia. In striking contrast to what has been previously reported, we demonstrate that the common blacktip shark is not rare in northern Australia but occurs in approximately equal frequencies with the Australian blacktip shark. Management of shark resources in northern Australia needs to take account of this new information. Species identification was performed using nucleotide sequences of the control, NADH dehydrogenase subunit 4 (ND4) and cytochrome oxidase I (COI) regions in the mitochondrial genome. The proportion of overall genetic variation (FST) between the two species was small (0.042, P < 0.01) based on allele frequencies at five microsatellite loci. We confirm that a third blacktip species (C. amblyrhynchoides, graceful shark) is closely related to C. tilstoni and C. limbatus and can be distinguished from them on the basis of mtDNA sequences from two gene regions. The Australian blacktip shark (C. tilstoni) was not encountered among 20 samples from central Indonesia that were later confirmed to be common blacktip and graceful sharks. Fisheries regulators urgently need new information on life history, population structure and morphological characters for species identification of blacktip shark species in Australia.

Population genetic evidence for the east–west division of the narrow-barred Spanish mackerel (Scomberomorus commerson, Perciformes: Teleostei) along Wallace’s Line

Mitochondrial DNA D-loop (control) region (426-bp) was used to infer the genetic structure of Spanish mackerel (Scomberomorus commerson) from populations in Southeast Asia (Brunei, East and West Malaysia, Philippines, Thailand, Singapore, and China) and northern Australia (including western Timor). An east–west division along Wallace’s Line was strongly supported by a significant AMOVA, with 43% of the total sequence variation partitioned among groups of populations. Phylogenetic and network analyses supported two clades: clade A and clade B. Members of clade A were found in Southeast Asia and northern Australia, but not in locations to the west (Gulf of Thailand) or north (China). Clade B was found exclusively in Southeast Asia. Genetic division along Wallace’s Line suggests that co-management of S. commerson populations for future sustainability may not be necessary between Southeast Asian nations and Australia, however all countries should share the task of management of the species in Southeast Asia equally. More detailed genetic studies of S. commerson populations in the region are warranted.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: ‘Least Concern’ and ‘Vulnerable’. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric FST analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.