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.

Ecology of uncultivated oscillospira species in the rumen of cattle, sheep, and reindeer as assessed by microscopy and molecular approaches

The ecology of the uncultured, but large and morphologically conspicuous, rumen bacterium Oscillospira spp. was studied. Oscillospira-specific 16S rRNA gene sequences were detected in North American domestic cattle, sheep from Australia and Japan, and Norwegian reindeer. Phylogenetic analysis of the sequences obtained allowed definition of three operational taxonomic units within the Oscillospira clade. Consistent with this genetic diversity, we observed atypical smaller morphotypes by using an Oscillospira-specific fluorescence in situ hybridization probe. Despite the visual disappearance of typical large Oscillospira morphotypes, the presence of Oscillospira spp. was still detected by Oscillospira-specific PCR in the rumen of cattle and sheep. These observations suggest the broad presence of Oscillospira species in various rumen ecosystems with the level, and most likely the morphological form, dependent on diet. An ecological analysis based on enumeration of the morphologically conspicuous, large-septate form confirms that the highest counts are associated with the feeding of fresh forage diets to cattle and sheep and in two different subspecies of reindeer investigated.

Molecular tools for studying puberty in the grey mullet, Mugil cephalus: cloning of cDNAs that regulate reproductive function

The cDNAs coding for the brain GnRHs (AY373449-51), pituitary GH, SL and PRL, and liver IGFs (AY427954-5) were isolated. Partial cDNA sequences of the brain (Cyp19b) and gonadal (Cyp19a) aromatases have also been obtained. These tools would be utilized to study the endocrine regulation of puberty in the grey mullet.

Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Molecular characterization and pathogenicity of Pythium species associated with damping-off in greenhouse cucumber (Cucumis sativus) in Oman

A study was undertaken in 2004 and 2005 to characterize pathogens associated with damping-off of greenhouse-grown cucumber seedlings in 13 districts in Oman. Identification of Pythium to the species level was based on sequences of the internal transcribed spacer (ITS) of the ribosomal DNA. Of the 98 Pythium isolates collected during the survey, Pythium aphanidermatum, P. spinosum, P. splendens and P. oligandrum accounted for 76%, 22%, 1% and 1%, respectively. Pythium aphanidermatum was isolated from all of the districts, while P. spinosum was isolated from seven districts. Pathogenicity tests showed inter- and intraspecific variation in aggressiveness between Pythium species. Pythium aphanidermatum, P. spinosum and P. splendens were found to be highly aggressive at 25°C. However, the aggressiveness of P. spinosum decreased when the temperature was raised to 30°C, which was found to correspond to the lower frequency of isolation of P. spinosum in the warmer seasons, compared to the cooler time of the year. Pythium aphanidermatum exhibited limited intraspecific variation in the sequences of the ITS region of the rDNA and showed 100% similarity to the corresponding P. aphanidermatum sequences from GenBank. The ITS sequence data, as well as morphological characteristics of P. spinosum isolates, showed a high level of similarity within and between P. spinosum and P. kunmingense, and suggested that the two species were synonymous. This study represents the first report of P. spinosum, P. splendens and P. oligandrum in Oman.

Molecular discrimination of Perna (Mollusca: Bivalvia) species using the polymerase chain reaction and species-specific mitochondrial primers

This work was prompted by the need to be able to identify the invasive mussel species, Perna viridis, in tropical Australian seas using techniques that do not rely solely on morphology. DNA-based molecular methods utilizing a polymerase chain reaction (PCR) approach were developed to distinguish unambiguously between the three species in the genus Perna. Target regions were portions of two mitochondrial genes, cox1 and nad4, and the intergenic spacer between these that occurs in at least two Perna species. Based on interspecific sequence comparisons of the nad4 gene, a conserved primer has been designed that can act as a forward primer in PCRs for any Perna species. Four reverse primers have also been designed, based on nad4 and intergenic spacer sequences, which yield species-specific products of different lengths when paired with the conserved forward primer. A further pair of primers has been designed that will amplify part of the cox1 gene of any Perna species, and possibly other molluscs, as a positive control to demonstrate that the PCR is working.

Recent advances in the molecular biology of Leifsonia xyli subsp. xyli, causal organism of ratoon stunting disease

Twelve years ago our understanding of ratoon stunting disease (RSD) was confined almost exclusively to diagnosis of the disease and control via farm hygiene, with little understanding of the biology of the interaction between the causal agent (Leifsonia xyli subsp. xyli) and the host plant sugarcane (Saccharum spp. hybrids). Since then, research has focused on developing the molecular tools to dissect L. xyli subsp. xyli, so that better control strategies can be developed to prevent losses from RSD. Within this review, we give a brief overview of the progression in research on L. xyli subsp. xyli and highlight future challenges. After a brief historical background on RSD, we discuss the development of molecular tools such as transformation and transposon mutagenesis and discuss the apparent lack of genetic diversity within the L. xyli subsp. xyli world population. We go on to discuss the sequencing of the genome of L. xyli subsp. xyli, describe the key findings and suggest some future research based on known deficiencies that will capitalise on this tremendous knowledge base to which we now have access.

Assessing the relationship between patch type and soil mites: A molecular approach

An urgent need exists for indicators of soil health and patch functionality in extensive rangelands that can be measured efficiently and at low cost. Soil mites are candidate indicators, but their identification and handling is so specialised and time-consuming that their inclusion in routine monitoring is unlikely. The aim of this study was to measure the relationship between patch type and mite assemblages using a conventional approach. An additional aim was to determine if a molecular approach traditionally used for soil microbes could be adapted for soil mites to overcome some of the bottlenecks associated with soil fauna diversity assessment. Soil mite species abundance and diversity were measured using conventional ecological methods in soil from patches with perennial grass and litter cover (PGL), and compared to soil from bare patches with annual grasses and/or litter cover (BAL). Soil mite assemblages were also assessed using a molecular method called terminal-restriction fragment length polymorphism (T-RFLP) analysis. The conventional data showed a relationship between patch type and mite assemblage. The Prostigmata and Oribatida were well represented in the PGL sites, particularly the Aphelacaridae (Oribatida). For T-RFLP analysis, the mite community was represented by a series of DNA fragment lengths that reflected mite sequence diversity. The T-RFLP data showed a distinct difference in the mite assemblage between the patch types. Where possible, T-RFLP peaks were matched to mite families using a reference 18S rDNA database, and the Aphelacaridae prevalent in the conventional samples at PGL sites were identified, as were prostigmatids and oribatids. We identified limits to the T-RFLP approach and this included an inability to distinguish some species whose DNA sequences were similar. Despite these limitations, the data still showed a clear difference between sites, and the molecular taxonomic inferences also compared well with the conventional ecological data. The results from this study indicated that the T-RFLP approach was effective in measuring mite assemblages in this system. The power of this technique lies in the fact that species diversity and abundance data can be obtained quickly because of the time taken to process hundreds of samples, from soil DNA extraction to data output on the gene analyser, can be as little as 4 days.

Rice molecular breeding laboratories in the genomics era: Current status and future considerations

Using DNA markers in plant breeding with marker-assisted selection (MAS) could greatly improve the precision and efficiency of selection, leading to the accelerated development of new crop varieties. The numerous examples of MAS in rice have prompted many breeding institutes to establish molecular breeding labs. The last decade has produced an enormous amount of genomics research in rice, including the identification of thousands of QTLs for agronomically important traits, the generation of large amounts of gene expression data, and cloning and characterization of new genes, including the detection of single nucleotide polymorphisms. The pinnacle of genomics research has been the completion and annotation of genome sequences for indica and japonica rice. This information-coupled with the development of new genotyping methodologies and platforms, and the development of bioinformatics databases and software tools-provides even more exciting opportunities for rice molecular breeding in the 21st century. However, the great challenge for molecular breeders is to apply genomics data in actual breeding programs. Here, we review the current status of MAS in rice, current genomics projects and promising new genotyping methodologies, and evaluate the probable impact of genomics research. We also identify critical research areas to "bridge the application gap" between QTL identification and applied breeding that need to be addressed to realize the full potential of MAS, and propose ideas and guidelines for establishing rice molecular breeding labs in the postgenome sequence era to integrate molecular breeding within the context of overall rice breeding and research programs.

Concordant molecular and phenotypic data delineate new taxonomy and conservation priorities for the endangered mountain yellow-legged frog

The mountain yellow-legged frog Rana muscosa sensu lato, once abundant in the Sierra Nevada of California and Nevada, and the disjunct Transverse Ranges of southern California, has declined precipitously throughout its range, even though most of its habitat is protected. The species is now extinct in Nevada and reduced to tiny remnants in southern California, where as a distinct population segment, it is classified as Endangered. Introduced predators (trout), air pollution and an infectious disease (chytridiomycosis) threaten remaining populations. A Bayesian analysis of 1901 base pairs of mitochondrial DNA confirms the presence of two deeply divergent clades that come into near contact in the Sierra Nevada. Morphological studies of museum specimens and analysis of acoustic data show that the two major mtDNA clades are readily differentiated phenotypically. Accordingly, we recognize two species, Rana sierrae, in the northern and central Sierra Nevada, and R. muscosa, in the southern Sierra Nevada and southern California. Existing data indicate no range overlap. These results have important implications for the conservation of these two species as they illuminate a profound mismatch between the current delineation of the distinct population segments (southern California vs. Sierra Nevada) and actual species boundaries. For example, our study finds that remnant populations of R. muscosa exist in both the southern Sierra Nevada and the mountains of southern California, which may broaden options for management. In addition, despite the fact that only the southern California populations are listed as Endangered, surveys conducted since 1995 at 225 historic (1899-1994) localities from museum collections show that 93.3% (n=146) of R. sierrae populations and 95.2% (n=79) of R. muscosa populations are extinct. Evidence presented here underscores the need for revision of protected population status to include both species throughout their ranges.

A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and agricultural traits

Background: Molecular marker technologies are undergoing a transition from largely serial assays measuring DNA fragment sizes to hybridization-based technologies with high multiplexing levels. Diversity Arrays Technology (DArT) is a hybridization-based technology that is increasingly being adopted by barley researchers. There is a need to integrate the information generated by DArT with previous data produced with gel-based marker technologies. The goal of this study was to build a high-density consensus linkage map from the combined datasets of ten populations, most of which were simultaneously typed with DArT and Simple Sequence Repeat (SSR), Restriction Enzyme Fragment Polymorphism (RFLP) and/or Sequence Tagged Site (STS) markers. Results: The consensus map, built using a combination of JoinMap 3.0 software and several purpose-built perl scripts, comprised 2,935 loci (2,085 DArT, 850 other loci) and spanned 1,161 cM. It contained a total of 1,629 'bins' (unique loci), with an average inter-bin distance of 0.7 ± 1.0 cM (median = 0.3 cM). More than 98% of the map could be covered with a single DArT assay. The arrangement of loci was very similar to, and almost as optimal as, the arrangement of loci in component maps built for individual populations. The locus order of a synthetic map derived from merging the component maps without considering the segregation data was only slightly inferior. The distribution of loci along chromosomes indicated centromeric suppression of recombination in all chromosomes except 5H. DArT markers appeared to have a moderate tendency toward hypomethylated, gene-rich regions in distal chromosome areas. On the average, 14 ± 9 DArT loci were identified within 5 cM on either side of SSR, RFLP or STS loci previously identified as linked to agricultural traits. Conclusion: Our barley consensus map provides a framework for transferring genetic information between different marker systems and for deploying DArT markers in molecular breeding schemes. The study also highlights the need for improved software for building consensus maps from high-density segregation data of multiple populations.

Molecular and serological detection of A. centrale- and A. marginale-infected cattle grazing within an endemic area

A reverse line blot hybridization (RLB) one-stage nested PCR (nPCR) for Anaplasma centrale and a nested PCR for Anaplasma marginale were used to detect infected cattle grazing within an endemic region in Israel. A novel set of PCR primers and oligonucleotide probes based on a 16S ribosomal RNA gene was designed for RLB detection of both Anaplasma species, and the performance of the molecular assays compared. The immunofluorescent antibody test (IFA) was used to detect antibodies to both Anaplasma species, whereas, a highly sensitive and specific competitive enzyme-linked immunosorbent assay (cELISA) was used to detect antibodies in A. centrale-vaccinated cattle. The RLB and the nested PCR procedures showed bacteremia with sensitivity of 50 infected erythrocytes per milliliter. Up to 93% of the A. centrale vaccinates carried specific antibodies that were detected by cELISA, and up to 71% of the vaccinated cattle were found to be naturally infected with A. marginale according to the PCR and the RLB assays. Nevertheless, no severe outbreaks of A. marginale infection occurred among vaccinated herds in this endemic region. It appears that both, molecular tools and serology are useful for evaluation of the vaccine efficacy. In the light of wide natural field infection with A. marginale, strong recommendations to continue the A. centrale vaccination program regime will continue until a new generation of non-blood-based vaccine will be developed.

Molecular characterisation, pathogenesis and fungicide sensitivity of Pythium spp. from table beet (Beta vulgaris var. vulgaris) grown in the Lockyer Valley, Queensland

Table beet production in the Lockyer Valley of south-eastern Queensland is known to be adversely affected by soilborne root disease from infection by Pythium spp. However, little is known regarding the species or genotypes that are the causal agents of both pre- and post-emergence damping off. Based on RFLP analysis with HhaI, HinfI and MboI of the PCR amplified ITS region DNA from soil and diseased plant samples, the majority of 130 Pythium isolates could be grouped into three genotypes, designated LVP A, LVP B and LVP C. These groups comprised 43, 41 and 7% of all isolates, respectively. Deoxyribonucleic acid sequence analysis of the ITS region indicated that LVP A was a strain of Pythium aphanidermatum, with greater than 99% similarity to the corresponding P. aphanidermatum sequences from the publicly accessible databases. The DNA sequences from LVP B and LVP C were most closely related to P. ultimum and P. dissotocum, respectively. Lower frequencies of other distinct isolates with unique RFLP patterns were also obtained with high levels of similarity (>97%) to P. heterothallicum, P. periplocum and genotypes of P. ultimum other than LVP B. Inoculation trials of 1- and 4-week-old beet seedlings indicated that compared with isolates of the LVP B genotype, a higher frequency of LVP A isolates caused disease. Isolates with the LVP A, LVP B and LVP C genotypes were highly sensitive to the fungicide Ridomil MZ, which suppressed radial growth on V8 agar between approximately four and thirty fold at 5 μg/mL metalaxyl and 40 μg/mL mancozeb, a concentration far lower than the recommended field application rate.

Identification of quantitative trait loci for resistance to two species of root-lesion nematode (Pratylenchus thornei and P. neglectus) in wheat

Pratylenchus thornei and P. neglectus are two species of root-lesion nematode that cause substantial yield losses in wheat. No commercially available wheat variety has resistance to both species. A doubled-haploid population developed from a cross between the synthetic hexaploid wheat line CPI133872 and the bread wheat Janz was used to locate and tag quantitative trait loci (QTLs) associated with resistance to both P. thornei and P. neglectus. Wheat plants were inoculated with both species of nematode in independent replicated glasshouse trials repeated over 2 years. Known locations of wheat microsatellite markers were used to construct a framework map. After an initial single-marker analysis to detect marker-trait linkages, chromosome regions associated with putative QTLs were targetted with microsatellite markers to increase map density in the chromosome regions of interest. In total, 148 wheat microsatellite markers and 21 amplified fragment length polymorphism markers were mapped. The codominant microsatellite marker Xbarc183 on the distal end of chromosome 6DS was allelic for resistance to both P. thornei and P. neglectus. The QTL were designated QRlnt.lrc-6D.1 and QRlnn.lrc-6D.1, for the 2 traits, respectively. The allele inherited from CPI133872 explained 22.0-24.2% of the phenotypic variation for P. thornei resistance, and the allele inherited from Janz accounted for 11.3-14.0% of the phenotypic variation for P. neglectus resistance. Composite interval mapping identified markers that flank a second major QTL on chromosome 6DL (QRlnt.lrc-6D.2) that explained 8.3-13.4% of the phenotypic variation for P. thornei resistance. An additional major QTL associated with P. neglectus resistance was detected on chromosome 4DS (QRlnn.lrc-4D.1) and explained a further 10.3-15.4% of the phenotypic variation. The identification and tagging of nematode resistance genes with molecular markers will allow appropriate allele combinations to be selected, which will aid the successful breeding of wheat with dual nematode resistance.

Application of phytotoxicity data to a new Australian soil quality guideline framework for biosolids

To protect terrestrial ecosystems and humans from contaminants many countries and jurisdictions have developed soil quality guidelines (SQGs). This study proposes a new framework to derive SQGs and guidelines for amended soils and uses a case study based on phytotoxicity data of copper (Cu) and zinc (Zn) from field studies to illustrate how the framework could be applied. The proposed framework uses normalisation relationships to account for the effects of soil properties on toxicity data followed by a species sensitivity distribution (SSD) method to calculate a soil added contaminant limit (soil ACL) for a standard soil. The normalisation equations are then used to calculate soil ACLs for other soils. A soil amendment availability factor (SAAF) is then calculated as the toxicity and bioavailability of pure contaminants and contaminants in amendments can be different. The SAAF is used to modify soil ACLs to ACLs for amended soils. The framework was then used to calculate soil ACLs for copper (Cu) and zinc (Zn). For soils with pH of 4-8 and OC content of 1-6%, the ACLs range from 8 mg/kg to 970 mg/kg added Cu. The SAAF for Cu was pH dependant and varied from 1.44 at pH 4 to 2.15 at pH 8. For soils with pH of 4-8 and OC content of 1-6%, the ACLs for amended soils range from 11 mg/kg to 2080 mg/kg added Cu. For soils with pH of 4-8 and a CEC from 5-60, the ACLs for Zn ranged from 21 to 1470 mg/kg added Zn. A SAAF of one was used for Zn as it concentrations in plant tissue and soil to water partitioning showed no difference between biosolids and soluble Zn salt treatments, indicating that Zn from biosolids and Zn salts are equally bioavailable to plants.