Niche specificity, polygeny, and pleiotropy in herbivorous insects

Why do herbivorous insects tend to be host specialists? Population genetic models predict specialization when there is antagonistic pleiotropy at a gene for host-use performance. But empirically, host-use performance is governed by many genetic regions, and antagonistic pleiotropy is rare. Here, we use individual-based quantitative genetic simulations to investigate the role of pleiotropy in the evolution of host-use specialization when host-use performance is polygenic. We find that if host-preference is allowed to evolve without cost, parasite populations tend to evolve host-use specialization even without pleiotropy; thus, it would seem that for a polygenic trait, the benefit of maintaining adaptive combinations of conditionally-neutral alleles suffices to drive specialization. But if there is a fecundity cost for host-preference, or if host patches are demographically volatile, host-use generalists evolve, even with high levels of pleiotropy. In sum, if pleiotropy is much more pervasive than has been observed in nature, our simulations show that it could play a role in driving the evolution of polygenic specialization. But pleiotropy is not necessary, and even when it is extensive, selection can favor generalist genotypes.

Investigating mechanisms underlying genetic resistance to Salmon Rickettsial Syndrome in Atlantic salmon using RNA sequencing

Background Salmon Rickettsial Syndrome (SRS), caused by Piscirickettsia salmonis, is one of the primary causes of morbidity and mortality in Atlantic salmon aquaculture, particularly in Chile. Host resistance is a heritable trait, and functional genomic studies have highlighted genes and pathways important in the response of salmon to the bacteria. However, the functional mechanisms underpinning genetic resistance are not yet well understood. In the current study, a large population of salmon pre-smolts were challenged with P. salmonis, with mortality levels recorded and samples taken for genotyping. In parallel, head kidney and liver samples were taken from animals of the same population with high and low genomic breeding values for resistance, and used for RNA-Sequencing to compare their transcriptome profile both pre and post infection. Results A significant and moderate heritability (h2 = 0.43) was shown for the trait of binary survival. Genome-wide association analyses using 38 K imputed SNP genotypes across 2265 animals highlighted that resistance is a polygenic trait. Several thousand genes were identified as differentially expressed between controls and infected samples, and enriched pathways related to the host immune response were highlighted. In addition, several networks with significant correlation with SRS resistance breeding values were identified, suggesting their involvement in mediating genetic resistance. These included apoptosis, cytoskeletal organisation, and the inflammasome. Conclusions While resistance to SRS is a polygenic trait, this study has highlighted several relevant networks and genes that are likely to play a role in mediating genetic resistance. These genes may be future targets for functional studies, including genome editing, to further elucidate their role underpinning genetic variation in host resistance.

Characterization, Genetic Analyses, and Identification of QTLs Conferring Metabolic Resistance to a 4-Hydroxyphenylpyruvate Dioxygenase Inhibitor in Sorghum (Sorghum bicolor)

Postemergence grass weed control continues to be a major challenge in grain sorghum [Sorghum bicolor (L.) Moench], primarily due to lack of herbicide options registered for use in this crop. The development of herbicide-resistant sorghum technology to facilitate broad-spectrum postemergence weed control can be an economical and viable solution. The 4-hydroxyphenylpyruvate dioxygenase-inhibitor herbicides (e.g., mesotrione or tembotrione) can control a broad spectrum of weeds including grasses, which, however, are not registered for postemergence application in sorghum due to crop injury. In this study, we identified two tembotrione-resistant sorghum genotypes (G-200, G-350) and one susceptible genotype (S-1) by screening 317 sorghum lines from a sorghum association panel (SAP). These tembotrione-resistant and tembotrione-susceptible genotypes were evaluated in a tembotrione dose–response [0, 5.75, 11.5, 23, 46, 92 (label recommended dose), 184, 368, and 736 g ai ha–1] assay. Compared with S-1, the genotypes G-200 and G-350 exhibited 10- and seven fold more resistance to tembotrione, respectively. To understand the inheritance of tembotrione-resistant trait, crosses were performed using S-1 and G-200 or G-350 to generate F1 and F2 progeny. The F1 and F2 progeny were assessed for their response to tembotrione treatment. Genetic analyses of the F1 and F2 progeny demonstrated that the tembotrione resistance in G-200 and G-350 is a partially dominant polygenic trait. Furthermore, cytochrome P450 (CYP)-inhibitor assay using malathion and piperonyl butoxide suggested possible CYP-mediated metabolism of tembotrione in G-200 and G-350. Genotype-by-sequencing based quantitative trait loci (QTL) mapping revealed QTLs associated with tembotrione resistance in G-200 and G-350 genotypes. Overall, the genotypes G-200 and G-350 confer a high level of metabolic resistance to tembotrione and controlled by a polygenic trait. There is an enormous potential to introgress the tembotrione resistance into breeding lines to develop agronomically desirable sorghum hybrids.

Investigating mechanisms underlying genetic resistance to Salmon Rickettsial Syndrome in Atlantic salmon using RNA sequencing

ABSTRACTBackgroundSalmon Rickettsial Syndrome (SRS), caused by Piscirickettsia salmonis, is one of the primary causes of morbidity and mortality in Atlantic salmon aquaculture, particularly in Chile. Host resistance is a heritable trait, and functional genomic studies have highlighted genes and pathways important in the response of salmon to the bacteria. However, the functional mechanisms underpinning genetic resistance are not yet well understood. In the current study, a large population of salmon pre-smolts were challenged with P. salmonis, with mortality levels recorded and samples taken for genotyping. In parallel, head kidney and liver samples were taken from animals of the same population with high and low genomic breeding values for resistance, and used for RNA-Sequencing to compare their transcriptome profile both pre and post infection.ResultsA significant and moderate heritability (h2 = 0.43) was shown for the trait of binary survival. Genome-wide association analyses using 38K imputed SNP genotypes across 2,251 animals highlighted that resistance is a polygenic trait. Several thousand genes were identified as differentially expressed between controls and infected samples, and enriched pathways related to the host immune response were highlighted. In addition, several networks with significant correlation with SRS resistance breeding values were identified, suggesting their involvement in mediating genetic resistance. These included apoptosis, cytoskeletal organisation, and the inflammasome.ConclusionsWhile resistance to SRS is a polygenic trait, this study has highlighted several relevant networks and genes that are likely to play a role in mediating genetic resistance. These genes may be future targets for functional studies, including genome editing, to further elucidate their role underpinning genetic variation in host resistance.

Resistance to tembotrione, a 4- Hydroxyphenylpyruvate Dioxygenase-Inhibitor in Sorghum bicolor

Postemergence grass weed control continues to be a big challenge in grain sorghum (Sorghum bicolor L. Moench), primarily due to a lack of herbicide options registered for use in this crop. The development of herbicide-resistant sorghum technology to facilitate broad-spectrum postemergence weed control is an economical and viable solution. The 4-hydroxyphenylpyruvate dioxygenase-inhibitor herbicides (e.g. mesotrione or tembotrione) can control broad-spectrum of weeds including grasses, which, however, is not registered for postemergence application in sorghum due to crop injury. In this study we identified two tembotrione-resistant sorghum genotypes (G-200, G-350) and one highly susceptible genotype (S-1) through screening 317 sorghum lines from the sorghum association panel (SAP). Compared to S-1, G-200 and G-350 exhibited 10- and 7-fold more resistance to tembotrione, respectively. Genetic analyses of the F1 and F2 progeny generated from a cross between tembotrione-resistant and -susceptible genotypes demonstrated that the resistance is a semi-dominant polygenic trait. Furthermore, cytochrome P450 (CYP)-inhibitor assay using malathion and piperonyl butoxide suggested possible CYP-mediated metabolism of tembotrione in G-200 and G-350. Genotype-by-sequencing based quantitative trait loci (QTL) mapping revealed eight QTLs associated with tembotrione resistance in grain sorghum. Sorghum genotypes G-200 and G-350 confer a high level of metabolic resistance to tembotrione and controlled by a polygenic trait. There is an enormous potential to introgress the tembotrione resistance into breeding lines to develop agronomically desirable sorghum hybrids.One-sentence summaryThis research focuses on characterization, genetic analyses, identification of QTLs-linked to metabolic resistance to tembotrione in Sorghum bicolor, for improved weed control and increased yield

Disentangling selection on genetically correlated polygenic traits using whole-genome genealogies

We present a full-likelihood method to estimate and quantify polygenic adaptation from contemporary DNA sequence data. The method combines population genetic DNA sequence data and GWAS summary statistics from up to thousands of nucleotide sites in a joint likelihood function to estimate the strength of transient directional selection acting on a polygenic trait. Through population genetic simulations of polygenic trait architectures and GWAS, we show that the method substantially improves power over current methods. We examine the robustness of the method under uncorrected GWAS stratification, uncertainty and ascertainment bias in the GWAS estimates of SNP effects, uncertainty in the identification of causal SNPs, allelic heterogeneity, negative selection, and low GWAS sample size. The method can quantify selection acting on correlated traits, fully controlling for pleiotropy even among traits with strong genetic correlation (|rg| = 80%; c.f. schizophrenia and bipolar disorder) while retaining high power to attribute selection to the causal trait. We apply the method to study 56 human polygenic traits for signs of recent adaptation. We find signals of directional selection on pigmentation (tanning, sunburn, hair, P=5.5e-15, 1.1e-11, 2.2e-6, respectively), life history traits (age at first birth, EduYears, P=2.5e-4, 2.6e-4, respectively), glycated hemoglobin (HbA1c, P=1.2e-3), bone mineral density (P=1.1e-3), and neuroticism (P=5.5e-3). We also conduct joint testing of 137 pairs of genetically correlated traits. We find evidence of widespread correlated response acting on these traits (2.6-fold enrichment over the null expectation, P=1.5e-7). We find that for several traits previously reported as adaptive, such as educational attainment and hair color, a significant proportion of the signal of selection on these traits can be attributed to correlated response, vs direct selection (P=2.9e-6, 1.7e-4, respectively). Lastly, our joint test uncovers antagonistic selection that has acted to increase type 2 diabetes (T2D) risk and decrease HbA1c (P=1.5e-5).

Polygenic adaptation: From sweeps to subtle frequency shifts

1AbstractEvolutionary theory has produced two conflicting paradigms for the adaptation of a polygenic trait. While population genetics views adaptation as a sequence of selective sweeps at single loci underlying the trait, quantitative genetics posits a collective response, where phenotypic adaptation results from subtle allele frequency shifts at many loci. Yet, a synthesis of these views is largely missing and the population genetic factors that favor each scenario are not well understood. Here, we study the architecture of adaptation of a binary polygenic trait (such as resistance) with negative epistasis among the loci of its basis. The genetic structure of this trait allows for a full range of potential architectures of adaptation, ranging from sweeps to small frequency shifts. By combining computer simulations and a newly devised analytical framework based on Yule branching processes, we gain a detailed understanding of the adaptation dynamics for this trait. Our key analytical result is an expression for the joint distribution of mutant alleles at the end of the adaptive phase. This distribution characterizes the polygenic pattern of adaptation at the underlying genotype when phenotypic adaptation has been accomplished. We find that a single compound parameter, the population-scaled background mutation rate Θbg, explains the main differences among these patterns. For a focal locus, Θbg measures the mutation rate at all redundant loci in its genetic background that offer alternative ways for adaptation. For adaptation starting from mutation-selection-drift balance, we observe different patterns in three parameter regions. Adaptation proceeds by sweeps for small Θbg ≾ 0.1, while small polygenic allele frequency shifts require large Θbg ≿ 100. In the large intermediate regime, we observe a heterogeneous pattern of partial sweeps at several interacting loci.2Author summaryIt is still an open question how complex traits adapt to new selection pressures. While population genetics champions the search for selective sweeps, quantitative genetics proclaims adaptation via small concerted frequency shifts. To date the empirical evidence of clear sweep signals is more scarce than expected, while subtle shifts remain notoriously hard to detect. In the current study we develop a theoretical framework to predict the expected adaptive architecture of a simple polygenic trait, depending on parameters such as mutation rate, effective population size, size of the trait basis, and the available genetic variability at the onset of selection. For a population in mutation-selection-drift balance we find that adaptation proceeds via complete or partial sweeps for a large set of parameter values. We predict adaptation by small frequency shifts for two main cases. First, for traits with a large mutational target size and high levels of genetic redundancy among loci, and second if the starting frequencies of mutant alleles are more homogeneous than expected in mutation-selection-drift equilibrium, e.g. due to population structure or balancing selection.