Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

CattleTickBase: An integrated Internet-based bioinformatics resource for Rhipicephalus (Boophilus) microplus

The Rhipicephalus micro plus genome is large and complex in structure, making it difficult to assemble a genome sequence and costly to resource the required bioinformatics. In light of this, a consortium of international collaborators was formed to pool resources to begin sequencing this genome. We have acquired and assembled genomic DNA into contigs that represent over 1.8 Gigabase pairs of DNA from gene-enriched regions of the R. micro plus genome. We also have several datasets containing transcript sequences from a number of gene expression experiments conducted by the consortium. A web-based resource was developed to enable the scientific community to access our datasets and conduct analysis through a web-based bioinformatics environment called YABI. The collective bioinformatics resource is termed CattleTickBase. Our consortium has acquired genomic and transcriptomic sequence data at approximately 0.9X coverage of the gene-coding regions of the R. microplus genome. The YABI tool will facilitate access and manipulation of cattle tick genome sequence data as the genome sequencing of R. microplus proceeds. During this process the CattleTickBase resource will continue to be updated. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology Inc.

The use of genetic correlations to evaluate associations between SNP markers and quantitative traits

Open-pollinated progeny of Corymbia citriodora established in replicated field trials were assessed for stem diameter, wood density, and pulp yield prior to genotyping single nucleotide polymorphisms (SNP) and testing the significance of associations between markers and assessment traits. Multiple individuals within each family were genotyped and phenotyped, which facilitated a comparison of standard association testing methods and an alternative method developed to relate markers to additive genetic effects. Narrow-sense heritability estimates indicated there was significant additive genetic variance within this population for assessment traits ( h ˆ 2 =0.28to0.44 ) and genetic correlations between the three traits were negligible to moderate (r G = 0.08 to 0.50). The significance of association tests (p values) were compared for four different analyses based on two different approaches: (1) two software packages were used to fit standard univariate mixed models that include SNP-fixed effects, (2) bivariate and multivariate mixed models including each SNP as an additional selection trait were used. Within either the univariate or multivariate approach, correlations between the tests of significance approached +1; however, correspondence between the two approaches was less strong, although between-approach correlations remained significantly positive. Similar SNP markers would be selected using multivariate analyses and standard marker-trait association methods, where the former facilitates integration into the existing genetic analysis systems of applied breeding programs and may be used with either single markers or indices of markers created with genomic selection processes.

Mapping Quantitative Trait Loci for Partial Resistance to Powdery Mildew in an Australian Barley Population

Genomic regions influencing resistance to powdery mildew [Blumeria graminis (DC.) E.O. Speer f. sp. hordei Em. Marchal] were detected in a doubled haploid (DH) barley (Hordeum vulgare L.) population derived from a cross between the breeding line ND24260 and cultivar Flagship when evaluated across four field environments in Australia and Uruguay. Significant quantitative trait loci (OIL) for resistance to B. graminis were detected on six of the seven chromosomes (1H, 2H, 3H, 4H, 5H, and 7H). A QTL with large effect donated by ND24260 mapped to the short arm of chromosome 1H (1 HS) conferring near immunity to B. graminis in Australia but was ineffective in Uruguay. Three OIL donated by Flagship contributed partial resistance to B. graminis and were detected in at least two environments. These OIL were mapped to chromosomes 3H, 4H, and 5H (5HS) accounting for up to 18.6, 3.4, and 8.8% phenotypic variation, respectively. The 5HS QTL contributed partial resistance to B. graminis in all field environments in both Australia and Uruguay and aligned with the genomic region of Rph20, a gene conferring adult plant resistance (APR) to leaf rust (Puccinia hordei Otth), which is found in some cultivars having Vada' or 'Emir' in their parentage. Selection for favorable marker haplotypes within the 3H, 4H, and 5H QTL regions can be performed even in the presence of single (major) gene resistance. Pyramiding such QTL may provide an effective and potentially durable form of resistance to B. graminis.

High resolution mapping of Dense spike-ar (dsp.ar) to the genetic centromere of barley chromosome 7H

Spike density in barley is under the control of several major genes, as documented previously by genetic analysis of a number of morphological mutants. One such class of mutants affects the rachis internode length leading to dense or compact spikes and the underlying genes were designated dense spike (dsp). We previously delimited two introgressed genomic segments on chromosome 3H (21 SNP loci, 35.5 cM) and 7H (17 SNP loci, 20.34 cM) in BW265, a BC7F3 nearly isogenic line (NIL) of cv. Bowman as potentially containing the dense spike mutant locus dsp.ar, by genotyping 1,536 single nucleotide polymorphism (SNP) markers in both BW265 and its recurrent parent. Here, the gene was allocated by high-resolution bi-parental mapping to a 0.37 cM interval between markers SC57808 (Hv_SPL14)-CAPSK06413 residing on the short and long arm at the genetic centromere of chromosome 7H, respectively. This region putatively contains more than 800 genes as deduced by comparison with the collinear regions of barley, rice, sorghum and Brachypodium, Classical map-based isolation of the gene dsp.ar thus will be complicated due to the infavorable relationship of genetic to physical distances at the target locus.

High global diversity of cycloviruses amongst dragonflies

Members of the family Circoviridae, specifically the genus Circovirus, were thought to infect only vertebrates; however, members of a sister group under the same family, the proposed genusCyclovirus, have been detected recently in insects. In an effort to explore the diversity of cycloviruses and better understand the evolution of these novel ssDNA viruses, here we present five cycloviruses isolated from three dragonfly species (Orthetrum sabina, Xanthocnemis zealandica and Rhionaeschna multicolor) collected in Australia, New Zealand and the USA, respectively. The genomes of these five viruses share similar genome structure to other cycloviruses, with a circular ~1.7 kb genome and two major bidirectionally transcribed ORFs. The genomic sequence data gathered during this study were combined with all cyclovirus genomes available in public databases to identify conserved motifs and regulatory elements in the intergenic regions, as well as determine diversity and recombinant regions within their genomes. The genomes reported here represent four different cyclovirus species, three of which are novel. Our results confirm that cycloviruses circulate widely in winged-insect populations; in eight different cyclovirus species identified in dragonflies to date, some of these exhibit a broad geographical distribution. Recombination analysis revealed both intra-and inter-species recombination events amongst cycloviruses, including genomes recovered from disparate sources (e.g. goat meat and human faeces). Similar to other well-characterized circular ssDNA viruses, recombination may play an important role in cyclovirus evolution. © 2013 SGM.

Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional-scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS-FLX Titanium platform. A total of 172837 raw reads (17.4Mbp) were retrieved, 95500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. FST=0.0025) can be resolved with a likely power 67% by using standard sample sizes of 50 and validated subsets of available markers.

Rhipicephalus (Boophilus) microplus tick in vitro feeding methods for functional (dsRNA) and vaccine candidate (antibody) screening

Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) ticks cause economic losses for cattle industries throughout tropical and subtropical regions of the world estimated at $US2.5 billion annually. Lack of access to efficacious long-lasting vaccination regimes and increases in tick acaricide resistance have led to the investigation of targets for the development of novel tick vaccines and treatments. In vitro tick feeding has been used for many tick species to study the effect of new acaricides on the transmission of tick-borne pathogens. Few studies have reported the use of in vitro feeding for functional genomic studies using RNA interference and/or the effect of specific anti-tick antibodies. In particular, in vitro feeding reports for the cattle tick are limited due to its relatively short hypostome. Previously published methods were further modified to broaden optimal tick sizes/weights, feeding sources including bovine and ovine serum, optimisation of commercially available blood anti-coagulant tubes, and IgG concentrations for effective antibody delivery. Ticks are fed overnight and monitored for ∼5–6 weeks to determine egg output and success of larval emergence using a humidified incubator. Lithium-heparin blood tubes provided the most reliable anti-coagulant for bovine blood feeding compared with commercial citrated (CPDA) and EDTA tubes. Although >30 mg semi-engorged ticks fed more reliably, ticks as small as 15 mg also fed to repletion to lay viable eggs. Ticks which gained less than ∼10 mg during in vitro feeding typically did not lay eggs. One mg/ml IgG from Bm86-vaccinated cattle produced a potent anti-tick effect in vitro (83% efficacy) similar to that observed in vivo. Alternatively, feeding of dsRNA targeting Bm86 did not demonstrate anti-tick effects (11% efficacy) compared with the potent effects of ubiquitin dsRNA. This study optimises R. microplus tick in vitro feeding methods which support the development of cattle tick vaccines and treatments.

Association mapping of resistance to Puccinia hordei in Australian barley breeding germplasm

Key message “To find stable resistance using association mapping tools, QTL with major and minor effects on leaf rust reactions were identified in barley breeding lines by assessing seedlings and adult plants.” Abstract Three hundred and sixty (360) elite barley (Hordeum vulgare L.) breeding lines from the Northern Region Barley Breeding Program in Australia were genotyped with 3,244 polymorphic diversity arrays technology markers and the results used to map quantitative trait loci (QTL) conferring a reaction to leaf rust (Puccinia hordei Otth). The F3:5 (Stage 2) lines were derived or sourced from different geographic origins or hubs of international barley breeding ventures representing two breeding cycles (2009 and 2011 trials) and were evaluated across eight environments for infection type at both seedling and adult plant stages. Association mapping was performed using mean scores for disease reaction, accounting for family effects using the eigenvalues from a matrix of genotype correlations. In this study, 15 QTL were detected; 5 QTL co-located with catalogued leaf rust resistance genes (Rph1, Rph3/19, Rph8/14/15, Rph20, Rph21), 6 QTL aligned with previously reported genomic regions and 4 QTL (3 on chromosome 1H and 1 on 7H) were novel. The adult plant resistance gene Rph20 was identified across the majority of environments and pathotypes. The QTL detected in this study offer opportunities for breeding for more durable resistance to leaf rust through pyramiding multiple genomic regions via marker-assisted selection.

Spot form of net blotch resistance in barley is under complex genetic control

Key message: Evaluation of resistance toPyrenophora teresf.maculatain barley breeding populations via association mapping revealed a complex genetic architecture comprising a mixture of major and minor effect genes. Abstract: In the search for stable resistance to spot form of net blotch (Pyrenophora teres f. maculata, SFNB), association mapping was conducted on four independent barley (Hordeum vulgare L.) breeding populations comprising a total of 898 unique elite breeding lines from the Northern Region Barley Breeding Program in Australia for discovery of quantitative trait loci (QTL) influencing resistance at seedling and adult plant growth stages. A total of 29 significant QTL were validated across multiple breeding populations, with 22 conferring resistance at both seedling and adult plant growth stages. The remaining 7 QTL conferred resistance at either seedling (2 QTL) or adult plant (5 QTL) growth stages only. These 29 QTL represented 24 unique genomic regions, of which five were found to co-locate with previously identified QTL for SFNB. The results indicated that SFNB resistance is controlled by a large number of QTL varying in effect size with large effects QTL on chromosome 7H. A large proportion of the QTL acted in the same direction for both seedling and adult responses, suggesting that phenotypic selection for SFNB resistance performed at either growth stage could achieve adequate levels of resistance. However, the accumulation of specific resistance alleles on several chromosomes must be considered in molecular breeding selection strategies.

Phylodynamic evidence of the migration of turnip mosaic potyvirus from Europe to Australia and New Zealand

Turnip mosaic virus (TuMV) is a potyvirus that is transmitted by aphids and infects a wide range of plant species. We investigated the evolution of this pathogen by collecting 32 isolates of TuMV, mostly from Brassicaceae plants, in Australia and New Zealand. We performed a variety of sequence-based phylogenetic and population genetic analyses of the complete genomic sequences and of three non-recombinogenic regions of those sequences. The substitution rates, divergence times and phylogeographical patterns of the virus populations were estimated. Six inter- and seven intralineage recombination-type patterns were found in the genomes of the Australian and New Zealand isolates, and all were novel. Only one recombination-type pattern has been found in both countries. The Australian and New Zealand populations were genetically different, and were different from the European and Asian populations. Our Bayesian coalescent analyses, based on a combination of novel and published sequence data from three nonrecombinogenic protein-encoding regions, showed that TuMV probably started to migrate from Europe to Australia and New Zealand more than 80 years ago, and that distinct populations arose as a result of evolutionary drivers such as recombination. The basal-B2 subpopulation in Australia and New Zealand seems to be older than those of the world-B2 and -B3 populations. To our knowledge, our study presents the first population genetic analysis of TuMV in Australia and New Zealand. We have shown that the time of migration of TuMV correlates well with the establishment of agriculture and migration of Europeans to these countries.

Two distinct classes of QTL determine rust resistance in sorghum

Background: Agriculture is facing enormous challenges to feed a growing population in the face of rapidly evolving pests and pathogens. The rusts, in particular, are a major pathogen of cereal crops with the potential to cause large reductions in yield. Improving stable disease resistance is an on-going major and challenging focus for many plant breeding programs, due to the rapidly evolving nature of the pathogen. Sorghum is a major summer cereal crop that is also a host for a rust pathogen which occurs in almost all sorghum growing areas of the world, causing direct and indirect yield losses in sorghum worldwide, however knowledge about its genetic control is still limited. In order to further investigate this issue, QTL and association mapping methods were implemented to study rust resistance in three bi-parental populations and an association mapping set of elite breeding lines in different environments. Results: In total, 64 significant or highly significant QTL and 21 suggestive rust resistance QTL were identified representing 55 unique genomic regions. Comparisons across populations within the current study and with rust QTL identified previously in both sorghum and maize revealed a high degree of correspondence in QTL location. Negative phenotypic correlations were observed between rust, maturity and height, indicating a trend for both early maturing and shorter genotypes to be more susceptible to rust. Conclusions: The significant amount of QTL co-location across traits, in addition to the consistency in the direction of QTL allele effects, has provided evidence to support pleiotropic QTL action across rust, height, maturity and stay-green, supporting the role of carbon stress in susceptibility to rust. Classical rust resistance QTL regions that did not co-locate with height, maturity or stay-green QTL were found to be significantly enriched for the defence-related NBS-encoding gene family, in contrast to the lack of defence-related gene enrichment in multi-trait effect rust resistance QTL. The distinction of disease resistance QTL hot-spots, enriched with defence-related gene families from QTL which impact on development and partitioning, provides plant breeders with knowledge which will allow for fast-tracking varieties with both durable pathogen resistance and appropriate adaptive traits.

SorGSD: a sorghum genome SNP database

Sorghum (Sorghum bicolor) is one of the most important cereal crops globally and a potential energy plant for biofuel production. In order to explore genetic gain for a range of important quantitative traits, such as drought and heat tolerance, grain yield, stem sugar accumulation, and biomass production, via the use of molecular breeding and genomic selection strategies, knowledge of the available genetic variation and the underlying sequence polymorphisms, is required.

Molecular Breeding for Complex Adaptive Traits: How Integrating Crop Ecophysiology and Modelling Can Enhance Efficiency

Progress in crop improvement is limited by the ability to identify favourable combinations of genotypes (G) and management practices (M) in relevant target environments (E) given the resources available to search among the myriad of possible combinations. To underpin yield advance we require prediction of phenotype based on genotype. In plant breeding, traditional phenotypic selection methods have involved measuring phenotypic performance of large segregating populations in multi-environment trials and applying rigorous statistical procedures based on quantitative genetic theory to identify superior individuals. Recent developments in the ability to inexpensively and densely map/sequence genomes have facilitated a shift from the level of the individual (genotype) to the level of the genomic region. Molecular breeding strategies using genome wide prediction and genomic selection approaches have developed rapidly. However, their applicability to complex traits remains constrained by gene-gene and gene-environment interactions, which restrict the predictive power of associations of genomic regions with phenotypic responses. Here it is argued that crop ecophysiology and functional whole plant modelling can provide an effective link between molecular and organism scales and enhance molecular breeding by adding value to genetic prediction approaches. A physiological framework that facilitates dissection and modelling of complex traits can inform phenotyping methods for marker/gene detection and underpin prediction of likely phenotypic consequences of trait and genetic variation in target environments. This approach holds considerable promise for more effectively linking genotype to phenotype for complex adaptive traits. Specific examples focused on drought adaptation are presented to highlight the concepts.