Polygenic basis for adaptive morphological variation in a threatened Aotearoa | New Zealand bird, the hihi (Notiomystis cincta).

To predict if a threatened species can adapt to changing selective pressures, it is crucial to understand the genetic basis of adaptive traits, especially in species historically affected by severe bottlenecks. We estimated the heritability of three hihi (Notiomystis cincta) morphological traits known to be under selection (nestling tarsus length, body mass and head-bill length) using 523 individuals and 39 699 single nucleotide polymorphisms (SNPs) from a 50 K Affymetrix SNP chip. We then examined the genetic architecture of the traits via chromosome partitioning analyses and genome-wide association scans (GWAS). Heritabilities estimated using pedigree relatedness or genomic relatedness were low. For tarsus length, the proportion of genetic variance explained by each chromosome was positively correlated with its size, and more than one chromosome explained significant variation for body mass and head-bill length. Finally, GWAS analyses suggested many loci of small effect contributing to trait variation for all three traits, although one locus (an SNP within an intron of the transcription factor HEY2) was tentatively associated with tarsus length. Our findings suggest a polygenic nature for the morphological traits, with many small effect size loci contributing to the majority of the variation, similar to results from many other wild populations. However, the small effective population size, polygenic architecture and already low heritabilities suggest that both the total response and rate of response to selection are likely to be limited in hihi.

Overexpression of the Chromosome Partitioning Gene parA in Azorhizobium caulinodans ORS571 Alters the Bacteroid Morphotype in Sesbania rostrata Stem Nodules.

Azorhizobium caulinodans ORS571 is a diazotroph that forms N2-fixing nodules on the roots and stems of the tropical legume Sesbania rostrata. Deletion of the parA gene of this bacterium results in cell cycle defects, pleiomorphic cell shape, and formation of immature stem nodules on its host plant. In this study, we constructed a parA overexpression mutant (PnptII-parA) to complement a previous study and provide new insights into bacteroid formation. We found that overproduction of ParA did not affect growth, cell morphology, chromosome partitioning, or vegetative nitrogen fixation in the free-living state. Under symbiosis, however, distinctive features, such as a single swollen bacteroid in one symbiosome, relatively narrow symbiosome space, and polyploid cells were observed. The morphotype of the PnptII-parA bacteroid is reminiscent of terminal differentiation in some IRLC indeterminate nodules, but S. rostrata is not thought to produce the NCR peptides that induce terminal differentiation in rhizobia. In addition, the transcript patterns of many symbiosis-related genes elicited by PnptII-parA were different from those elicited by the wild type. Accordingly, we propose that the particular symbiosome formation in PnptII-parA stem-nodules is due to cell cycle disruption caused by excess ParA protein in the symbiotic cells during nodulation.

Accounting for heteroscedasticity and censoring in chromosome partitioning analyses.

A fundamental assumption in quantitative genetics is that traits are controlled by many loci of small effect. Using genomic data, this assumption can be tested using chromosome partitioning analyses, where the proportion of genetic variance for a trait explained by each chromosome (h2 c ), is regressed on its size. However, as h2 c -estimates are necessarily positive (censoring) and the variance increases with chromosome size (heteroscedasticity), two fundamental assumptions of ordinary least squares (OLS) regression are violated. Using simulated and empirical data we demonstrate that these violations lead to incorrect inference of genetic architecture. The degree of bias depends mainly on the number of chromosomes and their size distribution and is therefore specific to the species; using published data across many different species we estimate that not accounting for this effect overall resulted in 28% false positives. We introduce a new and computationally efficient resampling method that corrects for inflation caused by heteroscedasticity and censoring and that works under a large range of dataset sizes and genetic architectures in empirical datasets. Our new method substantially improves the robustness of inferences from chromosome partitioning analyses.

Inference of genetic architecture from chromosome partitioning analyses is sensitive to genome variation, sample size, heritability and effect size distribution.

Genomewide association studies have contributed immensely to our understanding of the genetic basis of complex traits. One major conclusion arising from these studies is that most traits are controlled by many loci of small effect, confirming the infinitesimal model of quantitative genetics. A popular approach to test for polygenic architecture involves so-called "chromosome partitioning" where phenotypic variance explained by each chromosome is regressed on the size of the chromosome. First developed for humans, this has now been repeatedly used in other species, but there has been no evaluation of the suitability of this method in species that can differ in their genome characteristics such as number and size of chromosomes. Nor has the influence of sample size, heritability of the trait, effect size distribution of loci controlling the trait or the physical distribution of the causal loci in the genome been examined. Using simulated data, we show that these characteristics have major influence on the inferences of the genetic architecture of traits we can infer using chromosome partitioning analyses. In particular, small variation in chromosome size, small sample size, low heritability, a skewed effect size distribution and clustering of loci can lead to a loss of power and consequently altered inference from chromosome partitioning analyses. Future studies employing this approach need to consider and derive an appropriate null model for their study system, taking these parameters into consideration. Our simulation results can provide some guidelines on these matters, but further studies examining a broader parameter space are needed.

Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations.

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long-term studies across Europe and consequently are considered an ecological 'model organism'. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long-term study populations of great tits--one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population-specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population-specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free-living populations.

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population.

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high-density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex-, age- and environment-specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.