On systems thinking, systems biology and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

The potential of microarrays to assist shrimp breeding and production: a review

The shrimp aquaculture industry is a relatively new livestock industry, having developed over the past 30 years. Thus, it is poised to take advantage of new technologies from the outset of selective breeding programs. This contrasts with long established livestock industries, where there are already highly specialised breeds. This review focuses specifically on the potential application of microarrays to shrimp breeding. Potential applications of microarrays in selective breeding programs are summarised. Microarrays can be used as a rapid means to generate molecular markers for genetic linkage mapping, and genetic maps have been constructed for yeast, Arabidopsis and barley using microarray technology. Microarrays can also be used in the hunt for candidate genes affecting particular traits, leading to development of perfect markers for these traits (i.e. causative mutations). However, this requires that microarray analysis be combined with genetic linkage mapping, and that substantial genomic information is available for the species in question. A novel application of microarrays is to treat gene expression as a quantitative trait in itself and to combine this with linkage mapping to identify quantitative trait loci controlling the levels of gene expression; this approach may identify higher level regulatory genes in specific pathways. Finally, patterns of gene expression observed using microarrays may themselves be treated as phenotypic traits in selection programs (e.g. a particular pattern of gene expression might be indicative of a disease tolerant individual). Microarrays are now being developed for a number of shrimp species in laboratories around the world, primarily with a focus on identifying genes involved in the immune response. However, at present, there is no central repository of shrimp genomic information, which limits the rate at which shrimp genomic research can be progressed. The application of microarrays to shrimp breeding will be extremely limited until there is a shared repository of genomic information for shrimp, and the collective will and resources to develop comprehensive genomic tools for shrimp.

Relationships between hard-seededness and seed weight in mungbean (Vigna radiata) assessed by QTL analysis

Weather damage reduces the value of commercial mungbean, but hard-seededness can reduce the level of damage. However, attempts to breed large- and hard-seeded mungbean varieties have been unsuccessful. To understand the relationship between seed weight and hard-seededness, these traits were investigated using a quantitative trait loci (QTL) mapping approach with a recombinant inbred population derived from a cross between a completely soft-seeded variety and a completely hard-seeded genotype. The two parental genotypes also had a sixfold difference in seed weight. QTL analyses revealed four loci for hard-seededness and I I loci for seed weight. Two of the hardseededness loci co-localized with seed weight QTL. When seed weight was used as a covariate in the analysis of hard-seededness from the field data, two of the four hard-seeded QTL remained significant with the effect at one of these remaining unchanged. These results explain why retaining hard-seededness in large seeded mungbean lines has been unsuccessful. The existence of a persistent locus, however, indicated that breeding large and persistently hard-seeded varieties of mungbean may be possible.

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within- family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

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.

Congruence in QTL for adventitious rooting in Pinus elliottii X Pinus caribaea hybrids resolves between and within-species effects

Targeting between-species effects for improvement in synthetic hybrid populations derived from outcrossing parental tree species may be one way to increase the efficacy and predictability of hybrid breeding. We present a comparative analysis of the quantitative trait loci (QTL) which resolved between from within-species effects for adventitious rooting in two populations of hybrids between Pinus elliottii and P. caribaea, an outbred F-1 (n=287) and an inbred-like F-2 family (n=357). Most small to moderate effect QTL (each explaining 2-5% of phenotypic variation, PV) were congruent (3 out of 4 QTL in each family) and therefore considered within-species effects as they segregated in both families. A single large effect QTL (40% PV) was detected uniquely in the F-2 family and assumed to be due to a between-species effect, resulting from a genetic locus with contrasting alleles in each parental species. Oligogenic as opposed to polygenic architecture was supported in both families (60% and 20% PV explained by 4 QTL in the F-2 and F-1 respectively). The importance of adventitious rooting for adaptation to survive water-logged environments was thought in part to explain oligogenic architecture of what is believed to be a complex trait controlled by many hundreds of genes.

Dissecting the complex genetic basis of mate choice

The genetic analysis of mate choice is fraught with difficulties. Males produce complex signals and displays that can consist of a combination of acoustic, visual, chemical and behavioural phenotypes. Furthermore, female preferences for these male traits are notoriously difficult to quantify. During mate choice, genes not only affect the phenotypes of the individual they are in, but can influence the expression of traits in other individuals. How can genetic analyses be conducted to encompass this complexity? Tighter integration of classical quantitative genetic approaches with modern genomic technologies promises to advance our understanding of the complex genetic basis of mate choice.

Genetic improvement of lucerne for anthracnose (Colletotrichum trifolii) resistance

Anthracnose, caused by Colletotrichum trifolii, is one of the most serious diseases influencing lucerne persistence and productivity in eastern Australia. The disease is largely controlled by plant resistance; however, new pathotypes of C. trifolii have developed in Australia, seriously limiting the productive life of susceptible cultivars. This paper describes an incompletely recessive and quantitatively inherited resistance to C. trifolii identified in a clone (W116) from cv. Sequel. S-1, F-1, F-2 and backcross populations of W116 and D (highly susceptible clone) were studied for their reaction to C. trifolii race 1. Resistance was found to be quantitatively inherited, and quantitative trait loci associated with resistance and susceptibility were identified in a backcross population (D x W116) x D using random amplified polymorphic DNA and amplified fragment length polymorphic markers. A multi-locus region on linkage group 4 was found to contribute significantly to the resistance phenotype. The application of DNA markers to allow exploitation of this quantitatively inherited resistance in lucerne breeding is discussed.

Trans-ethnic fine mapping highlights kidney-function genes linked to salt sensitivity

We analyzed genome-wide association studies (GWASs), including data from 71,638 individuals from four ancestries, for estimated glomerular filtration rate (eGFR), a measure of kidney function used to define chronic kidney disease (CKD). We identified 20 loci attaining genome-wide-significant evidence of association (p < 5 × 10(-8)) with kidney function and highlighted that allelic effects on eGFR at lead SNPs are homogeneous across ancestries. We leveraged differences in the pattern of linkage disequilibrium between diverse populations to fine-map the 20 loci through construction of "credible sets" of variants driving eGFR association signals. Credible variants at the 20 eGFR loci were enriched for DNase I hypersensitivity sites (DHSs) in human kidney cells. DHS credible variants were expression quantitative trait loci for NFATC1 and RGS14 (at the SLC34A1 locus) in multiple tissues. Loss-of-function mutations in ancestral orthologs of both genes in Drosophila melanogaster were associated with altered sensitivity to salt stress. Renal mRNA expression of Nfatc1 and Rgs14 in a salt-sensitive mouse model was also reduced after exposure to a high-salt diet or induced CKD. Our study (1) demonstrates the utility of trans-ethnic fine mapping through integration of GWASs involving diverse populations with genomic annotation from relevant tissues to define molecular mechanisms by which association signals exert their effect and (2) suggests that salt sensitivity might be an important marker for biological processes that affect kidney function and CKD in humans.

Evidence for different QTL underlying the immune and hypersensitive responses of Eucalyptus globulus to the rust pathogen Puccinia psidii

The rust Puccinia psidii infects many species in the family Myrtaceae. Native to South America, the pathogen has recently entered Australia which has a rich Myrtaceous flora, including trees of the ecologically and economically important genus Eucalyptus. We studied the genetic basis of variation in rust resistance in Eucalyptus globulus, the main plantation eucalypt in Australia. Quantitative trait loci (QTL) analysis was undertaken using 218 genotypes of an outcross F2 mapping family, phenotyped by controlled inoculation of their open pollinated progeny with the strain of P. psidii found in Australia. QTL analyses were conducted using a binary classification of individuals with no symptoms (immune) versus those with disease symptoms, and in a separate analysis dividing plants with disease symptoms into those exhibiting the hypersensitive response versus those with more severe symptoms. Four QTL were identified, two influencing whether a plant exhibited symptoms (Ppr2 and Ppr3), and two influencing the presence or absence of a hypersensitive reaction (Ppr4 and Ppr5). These QTL mapped to four different linkage groups, none of which overlap with Ppr1, the major QTL previously identified for rust resistance in Eucalyptus grandis. Candidate genes within the QTL regions are presented and possible mechanisms discussed. Together with past findings, our results suggest that P. psidii resistance in eucalypts is quantitative in nature and influenced by the complex interaction of multiple loci of variable effect.

Markers for seedling and adult plant crown rot resistance in four partially resistant bread wheat sources

QTL identified for seedling and adult plant crown rot resistance in four partially resistant hexaploid wheat sources. PCR-based markers identified for use in marker-assisted selection. Crown rot, caused by Fusarium pseudograminearum, is an important disease of wheat in many wheat-growing regions globally. Complete resistance to infection by F. pseudograminearum has not been observed in a wheat host, but germplasm with partial resistance to this pathogen has been identified. The partially resistant wheat hexaploid germplasm sources 2-49, Sunco, IRN497 and CPI133817 were investigated in both seedling and adult plant field trials to identify markers associated with the resistance which could be used in marker-assisted selection programs. Thirteen different quantitative trait loci (QTL) conditioning crown rot resistance were identified in the four different sources. Some QTL were only observed in seedling trials whereas others appeared to be adult plant specific. For example while the QTL on chromosomes 1AS, 1BS, and 4BS contributed by 2-49 and on 2BS contributed by Sunco were detected in both seedling and field trials, the QTL on 1DL present in 2-49 and the QTL on 3BL in IRN497 were only detected in seedling trials. Genetic correlations between field trials of the same population were strong, as were correlations between seedling trials of the same population. Low to moderate correlations were observed between seedling and field trials. Flanking markers, most of which are less than 10 cM apart, have now been identified for each of the regions associated with crown rot resistance.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Genetic diversity of stay-green sorghums and their derivatives revealed by microsatellites

The genetic variability of 28 sorghum genotypes of known senescence phenotype was investigated using 66 SSR markers well-distributed across the sorghum genome. The genotypes of a number of lines from breeding programmes for stay-green were also determined. This included lines selected phenotypically for stay-green and also RSG 03123, a marker-assisted backcross progeny of R16 (recurrent parent) and B35 (stay-green donor). A total of 419 alleles were detected with a mean of 6.2 per locus. The number of alleles ranged from one for Xtxp94 to 14 for Xtxp88. Chromosome SBI-10 had the highest mean number of alleles (8.33), while SBI-05 had the lowest (4.17). The PIC values obtained ranged from zero to 0.89 in Xtxp94 and Xtxp88, respectively, with a mean of 0.68. On a chromosome basis, mean PIC values were highest in SBI-10 (0.81) and lowest in SBI-05 (0.53). Most of the alleles from B35 in RSG 03123 were found on chromosomes SBI-01, SBI-02 and SBI-03, confirming the successful introgression of quantitative trait loci associated with stay-green from B35 into the senescent background R16. However, the alternative stay-green genetic sources were found to be distinct based on either all the SSRs employed or using only those associated with the stay-green trait in B35. Therefore, the physiological and biochemical basis of each stay-green source should be evaluated in order to enhance the understanding of the functioning of the trait in the various backgrounds. These genetic sources of stay-green could provide a valuable resource for improving this trait in sorghum breeding programmes.

Identification of regulatory polymorphisms associated with breast cancer risk

Dissertação de Mestrado, Ciências Biomédicas, Departamento de Ciências Biomédicas e Medicina, Universidade do Algarve, 2014

Genome selection in fruit breeding: application to table grapes.

Genomic selection (GS) has recently been proposed as a new selection strategy which represents an innovative paradigm in crop improvement, now widely adopted in animal breeding. Genomic selection relies on phenotyping and high-density genotyping of a sufficiently large and representative sample of the target breeding population, so that the majority of loci that regulate a quantitative trait are in linkage disequilibrium with one or more molecular markers and can thus be captured by selection. In this study we address genomic selection in a practical fruit breeding context applying it to a breeding population of table grape obtained from a cross between the hybrid genotype D8909-15 (Vitis rupestris × Vitis arizonica/girdiana), which is resistant to dagger nematode and Pierce?s disease (PD), and ?B90-116?, a susceptible Vitis vinifera cultivar with desirable fruit characteristics. Our aim was to enhance the knowledge on the genomic variation of agronomical traits in table grape populations for future use in marker-assisted selection (MAS) and GS, by discovering a set of molecular markers associated with genomic regions involved in this variation. A number of Quantitative Trait Loci (QTL) were discovered but this method is inaccurate and the genetic architecture of the studied population was better captured by the BLasso method of genomic selection, which allowed for efficient inference about the genetic contribution of the various marker loci. The technology of genomic selection afforded greater efficiency than QTL analysis and can be very useful in speeding up the selection procedures for agronomic traits in table grapes.