A genome-wide association study of Cloninger's temperament scales: implications for the evolutionary genetics of personality

Variation in personality traits is 30-60% attributed to genetic influences. Attempts to unravel these genetic influences at the molecular level have, so far, been inconclusive. We performed the first genome-wide association study of Cloninger's temperament scales in a sample of 5117 individuals, in order to identify common genetic variants underlying variation in personality. Participants' scores on Harm Avoidance, Novelty Seeking, Reward Dependence, and Persistence were tested for association with 1,252,387 genetic markers. We also performed gene-based association tests and biological pathway analyses. No genetic variants that significantly contribute to personality variation were identified, while our sample provides over 90% power to detect variants that explain only 1% of the trait variance. This indicates that individual common genetic variants of this size or greater do not contribute to personality trait variation, which has important implications regarding the genetic architecture of personality and the evolutionary mechanisms by which heritable variation is maintained.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Common genetic influences underlie comorbidity of migraine and endometriosis

We examined the co-occurrence of migraine and endometriosis within the largest known collection of families containing multiple women with surgically confirmed endometriosis and in an independent sample of 815 monozygotic and 457 dizygotic female twin pairs. Within the endometriosis families, a significantly increased risk of migrainous headache was observed in women with endometriosis compared to women without endometriosis (odds ratio [OR] 1.57, 95% confidence interval [CI]: 1.12-2.21, P=0.009). Bivariate heritability analyses indicated no evidence for common environmental factors influencing either migraine or endometriosis but significant genetic components for both traits, with heritability estimates of 69 and 49%, respectively. Importantly, a significant additive genetic correlation (r(G) = 0.27, 95% CI: 0.06-0.47) and bivariate heritability (h(2)=0.17, 95% CI: 0.08-0.27) was observed between migraine and endometriosis. Controlling for the personality trait neuroticism made little impact on this association. These results confirm the previously reported comorbidity between migraine and endometriosis and indicate common genetic influences completely explain their co-occurrence within individuals. Given pharmacological treatments for endometriosis typically target hormonal pathways and a number of findings provide support for a relationship between hormonal variations and migraine, hormone-related genes and pathways are highly plausible candidates for both migraine and endometriosis. Therefore, taking into account the status of both migraine and endometriosis may provide a novel opportunity to identify the genes underlying them. Finally, we propose that the analysis of such genetically correlated comorbid traits can increase power to detect genetic risk loci through the use of more specific, homogenous and heritable phenotypes.

Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10−9 to P = 1.8 × 10−40) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10−3 to P = 1.2 × 10−13). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

A genome-wide linkage scan provides evidence for both new and previously reported loci influencing common migraine

Latent class analysis was performed on migraine symptom data collected in a Dutch population sample (N = 12,210, 59% female) in order to obtain empirical groupings of individuals suffering from symptoms of migraine headache. Based on these heritable groupings (h(2) = 0.49, 95% CI: 0.41-0.57) individuals were classified as affected (migrainous headache) or unaffected. Genome-wide linkage analysis was performed using genotype data from 105 families with at least 2 affected siblings. In addition to this primary phenotype, linkage analyses were performed for the individual migraine symptoms. Significance levels, corrected for the analysis of multiple traits, were determined empirically via a novel simulation approach. Suggestive linkage for migrainous headache was found on chromosomes 1 (LOD = 1.63; pointwise P = 0.0031), 13 (LOD = 1.63; P = 0.0031), and 20 (LOD = 1.85; P = 0.0018). Interestingly, the chromosome 1 peak was located close to the ATP1A2 gene, associated with familial hemiplegic migraine type 2 (FHM2). Individual symptom analysis produced a LOD score of 1.97 (P = 0.0013) on chromosome 5 (photo/phonophobia), a LOD score of 2.13 (P = 0.0009) on chromosome 10 (moderate/severe pain intensity) and a near significant LOD score of 3.31 (P = 0.00005) on chromosome 13 (pulsating headache). These peaks were all located near regions previously reported in migraine linkage studies. Our results provide important replication and support for the presence of migraine susceptibility genes within these regions, and further support the utility of an LCA-based phenotyping approach and analysis of individual symptoms in migraine genetic research. Additionally, our novel "2-step" analysis and simulation approach provides a powerful means to investigate linkage to individual trait components.

Trait components provide tools to dissect the genetic susceptibility of migraine

The commonly used "end diagnosis" phenotype that is adopted in linkage and association studies of complex traits is likely to represent an oversimplified model of the genetic background of a disease. This is also likely to be the case for common types of migraine, for which no convincingly associated genetic variants have been reported. In headache disorders, most genetic studies have used end diagnoses of the International Headache Society (IHS) classification as phenotypes. Here, we introduce an alternative strategy; we use trait components--individual clinical symptoms of migraine--to determine affection status in genomewide linkage analyses of migraine-affected families. We identified linkage between several traits and markers on chromosome 4q24 (highest LOD score under locus heterogeneity [HLOD] 4.52), a locus we previously reported to be linked to the end diagnosis migraine with aura. The pulsation trait identified a novel locus on 17p13 (HLOD 4.65). Additionally, a trait combination phenotype (IHS full criteria) revealed a locus on 18q12 (HLOD 3.29), and the age at onset trait revealed a locus on 4q28 (HLOD 2.99). Furthermore, suggestive or nearly suggestive evidence of linkage to four additional loci was observed with the traits phonophobia (10q22) and aggravation by physical exercise (12q21, 15q14, and Xp21), and, interestingly, these loci have been linked to migraine in previous studies. Our findings suggest that the use of symptom components of migraine instead of the end diagnosis provides a useful tool in stratifying the sample for genetic studies.

Phosphine Resistance in the Rust Red Flour Beetle, Tribolium castaneum (Coleoptera: Tenebrionidae): Inheritance, Gene Interactions and Fitness Costs

The recent emergence of heritable high level resistance to phosphine in stored grain pests is a serious concern among major grain growing countries around the world. Here we describe the genetics of phosphine resistance in the rust red flour beetle Tribolium castaneum (Herbst), a pest of stored grain as well as a genetic model organism. We investigated three field collected strains of T. castaneum viz., susceptible (QTC4), weakly resistant (QTC1012) and strongly resistant (QTC931) to phosphine. The dose-mortality responses of their test- and inter-cross progeny revealed that most resistance was conferred by a single major resistance gene in the weakly (3.2x) resistant strain. This gene was also found in the strongly resistant (431x) strain, together with a second major resistance gene and additional minor factors. The second major gene by itself confers only 12-206x resistance, suggesting that a strong synergistic epistatic interaction between the genes is responsible for the high level of resistance (431x) observed in the strongly resistant strain. Phosphine resistance is not sex linked and is inherited as an incompletely recessive, autosomal trait. The analysis of the phenotypic fitness response of a population derived from a single pair inter-strain cross between the susceptible and strongly resistant strains indicated the changes in the level of response in the strong resistance phenotype; however this effect was not consistent and apparently masked by the genetic background of the weakly resistant strain. The results from this work will inform phosphine resistance management strategies and provide a basis for the identification of the resistance genes.

The Relationship Between the Stay-Green Trait and Grain Yield in Elite Sorghum Hybrids Grown in a Range of Environments

The stay-green drought adaptation mechanism has been widely promoted as a way of improving grain yield and lodging resistance in sorghum [Sorghum bicolor (L.) Moench] and as a result has been the subject of many physiological and genetic studies. The relevance of these studies to elite sorghum hybrids is not clear given that they sample a limited number of environments and were conducted using inbred lines or relatively small numbers of experimental F-1 hybrids. In this study we investigated the relationship between stay-green and yield using data from breeding trials that sampled 1668 unique hybrid combinations and 23 environments whose mean yields varied from 2.3 to 10.5 t ha(-1). The strength and direction of the association between stay-green and grain yield varied with both environment and genetic background (male tester). The majority of associations were positive, particularly in environments with yields below 6 t ha(-1). As trial mean yield increased above 6 t ha(-1) there was a trend toward an increased number of negative associations; however, the number and magnitude of the positive associations were larger. Given that post-flowering drought is very commonly experienced by sorghum crops world wide and average yields are 1.2 and 2.5 t ha(-1) for the world and Australia, respectively, our results indicate that selection for stay-green in elite sorghum hybrids may be broadly beneficial for increasing yield in a wide range of environments.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60-103 days or 94-136 days, and those for AGECL were 95-125 or 141-176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

Quantitative effects of in-line operations on Campylobacter and Escherichia coli through two Australian broiler processing plants

Campylobacter is an important food borne pathogen, mainly associated with poultry. A lack of through-chain quantitative Campylobacter data has been highlighted within quantitative risk assessments. The aim of this study was to quantitatively and qualitatively measure Campylobacter and Escherichia coli concentration on chicken carcasses through poultry slaughter. Chickens (n = 240) were sampled from each of four flocks along the processing chain, before scald, after scald, before chill, after chill, after packaging and from individual caeca. The overall prevalence of Campylobacter after packaging was 83% with a median concentration of 0.8 log10 CFU/mL. The processing points of scalding and chilling had significant mean reductions of both Campylobacter (1.8 and 2.9 log10 CFU/carcase) and E. coli (1.3 and 2.5 log10 CFU/carcase). The concentration of E. coli and Campylobacter was significantly correlated throughout processing indicating that E. coli may be a useful indicator organism for reductions in Campylobacter concentration. The carriage of species varied between flocks, with two flocks dominated by Campylobacter coli and two flocks dominated by Campylobacter jejuni. Current processing practices can lead to significant reductions in the concentration of Campylobacter on carcasses. Further understanding of the variable effect of processing on Campylobacter and the survival of specific genotypes may enable more targeted interventions to reduce the concentration of this poultry associated pathogen.

Quantitative methods to measure pigmentation variation in farmed Giant Tiger Prawns, Penaeus monodon, and the effects of different harvest methods on cooked colour

Cooked prawn colour is known to be a driver of market price and a visual indicator of product quality for the consumer. Although there is a general understanding that colour variation exists in farmed prawns, there has been no attempt to quantify this variation or identify where this variation is most prevalent. The objectives of this study were threefold: firstly to compare three different quantitative methods to measure prawn colour or pigmentation, two different colorimeters and colour quantification from digital images. Secondly, to quantify the amount of pigmentation variation that exists in farmed prawns within ponds, across ponds and across farms. Lastly, to assess the effects of ice storage or freeze-thawing of raw product prior to cooking. Each method was able to detect quantitative differences in prawn colour, although conversion of image based quantification of prawn colour from RGB to Lab was unreliable. Considerable colour variation was observed between prawns from different ponds and different farms, and this variation potentially affects product value. Different post-harvest methods prior to cooking were also shown to have a profound detrimental effect on prawn colour. Both long periods of ice storage and freeze thawing of raw product were detrimental to prawn colour. However, ice storage immediately after cooking was shown to be beneficial to prawn colour. Results demonstrated that darker prawn colour was preserved by holding harvested prawns alive in chilled seawater, limiting the time between harvesting and cooking, and avoiding long periods of ice storage or freeze thawing of uncooked product.

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.

Highly heritable resistance to root-lesion nematode (Pratylenchus thornei) in Australian chickpea germplasm observed using an optimised glasshouse method and multi-environment trial analysis

Pratylenchus thornei is a root-lesion nematode (RLN) of economic significance in the grain growing regions of Australia. Chickpea (Cicer arietinum) is a significant legume crop grown throughout these regions, but previous testing found most cultivars were susceptible to P. thornei. Therefore, improved resistance to P. thornei is an important objective of the Australian chickpea breeding program. A glasshouse method was developed to assess resistance of chickpea lines to P. thornei, which requires relatively low labour and resource input, and hence is suited to routine adoption within a breeding program. Using this method, good differentiation of chickpea cultivars for P. thornei resistance was measured after 12 weeks. Nematode multiplication was higher for all genotypes than the unplanted control, but of the 47 cultivars and breeding lines tested, 17 exhibited partial resistance, allowing less than two fold multiplication. The relative differences in resistance identified using this method were highly heritable (0.69) and were validated against P. thornei data from seven field trials using a multi-environment trial analysis. Genetic correlations for cultivar resistance between the glasshouse and six of the field trials were high (>0.73). These results demonstrate that resistance to P. thornei in chickpea is highly heritable and can be effectively selected in a limited set of environments. The improved resistance found in a number of the newer chickpea cultivars tested shows that some advances have been made in the P. thornei resistance of Australian chickpea cultivars, and that further targeted breeding and selection should provide incremental improvements.

Living on voles - plastic life of the Ural owl

Individuals face variable environmental conditions during their life. This may be due to migration, dispersion, environmental changes or, for example, annual variation in weather conditions. Genetic adaptation to a novel environment happens through natural selection. Phenotypic plasticity allows, however, a quick individual response to a new environment. Phenotypic plasticity may also be beneficial for individual if the environment is highly variable. For example, eggs are costly to produce. If the food conditions vary significantly between breeding seasons it is useful to be able to adjust the clutch and egg size according to the food abundance. In this thesis I use Ural owl vole system to study phenotypic plasticity and natural selection using a number of reproduction related traits. The Ural owl (Strix uralensis) is a long-lived and sedentary species. The reproduction and survival of the Ural owl, in fact their whole life, is tied to the dramatically fluctuating vole densities. Ural owls do not cause vole cycles but they have to adjust their behaviour to the rather predictable population fluctuations of these small mammals. Earlier work with this system has shown that Ural owl laying date and clutch size are plastic in relation to vole abundance. Further, individual laying date clutch size reaction norms have been shown to vary in the amount of plasticity. My work extends the knowledge of natural selection and phenotypic plasticity in traits related to reproduction. I show that egg size, timing of the onset of incubation and nest defense aggressiveness are plastic traits with fitness consequences for the Ural owl. Although egg size is in general thought to be a fixed characteristic of an individual, this highly heritable trait in the Ural owl is also remarkably plastic in relation to the changes in vole numbers, Ural owls are laying the largest eggs when their prey is most abundant. Timing of the onset of incubation is an individual-specific property and plastic in relation to clutch size. Timing of incubation is an important underlying cause for asynchronous hatching in birds. Asynchronous hatching is beneficial to offspring survival in Ural owl. Hence, timing of the onset of incubation may also be under natural selection. Ural owl females also adjust their nest defense aggressiveness according to the vole dynamics, being most aggressive in years when they produce the largest broods. Individual females show different levels of nest defense aggressiveness. Aggressiveness is positively correlated with the phenotypic plasticity of aggressiveness. As elevated nest defense aggressiveness is selected for, it may promote the plasticity of aggressive nest defense behaviour. All the studied traits are repeatable or heritable on individual level, and their expression is either directly or indirectly sensitive to changes in vole numbers. My work considers a number of important fitness-related traits showing phenotypic plasticity in all of them. Further, in two chapters I show that there is individual variation in the amount of plasticity exhibited. These findings on plasticity in reproduction related traits suggest that variable environments indeed promote plasticity.