Population differentiation and behavioural association of the two ‘personality’ genes DRD4 and SERT in dunnocks (Prunella modularis)

Quantifying the variation in behaviour‐related genes within and between populations provides insight into how evolutionary processes shape consistent behavioural traits (i.e. personality). Deliberate introductions of non‐native species offer opportunities to investigate how such genes differ between native and introduced populations and how polymorphisms in the genes are related to variation in behaviour. Here, we compared the genetic variation of the two ‘personality’ genes, DRD4 and SERT, between a native (United Kingdom, UK) and an introduced (New Zealand, NZ) population of dunnocks, Prunella modularis. The NZ population showed a significantly lower number of single nucleotide polymorphisms (SNPs) compared to the UK population. Standardized F’_st estimates of the personality genes and neutral microsatellites indicate that selection (anthropogenic and natural) probably occurred during and post the introduction event. Notably, the largest genetic differentiation was found in the intronic regions of the genes. In the NZ population, we also examined the association between polymorphisms in DRD4 and SERT and two highly repeatable behavioural traits: flight‐initiation distance and mating status (promiscuous females and cobreeding males). We found 38 significant associations (for different allele effect models) between the two behavioural traits and the studied genes. Further, 22 of the tested associations showed antagonistic allele effects for males and females. Our findings illustrate how introduction events and accompanying ecological changes could influence the genetic diversity of behaviour‐related genes.     

Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild‐derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]‐element insertion mutants in Sema‐5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild‐type control. We observed significant variation in olfactory responses to benzaldehyde among F₁ genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome‐wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry‐go‐round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema‐5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome‐wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.     

Genic rather than genome‐wide differences between sexually deceptive Ophrys orchids with different pollinators

High pollinator specificity and the potential for simple genetic changes to affect pollinator attraction make sexually deceptive orchids an ideal system for the study of ecological speciation, in which change of flower odour is likely important. This study surveys reproductive barriers and differences in floral phenotypes in a group of four closely related, coflowering sympatric Ophrys species and uses a genotyping‐by‐sequencing (GBS) approach to obtain information on the proportion of the genome that is differentiated between species. Ophrys species were found to effectively lack postpollination barriers, but are strongly isolated by their different pollinators (floral isolation) and, to a smaller extent, by shifts in flowering time (temporal isolation). Although flower morphology and perhaps labellum coloration may contribute to floral isolation, reproductive barriers may largely be due to differences in flower odour chemistry. GBS revealed shared polymorphism throughout the Ophrys genome, with very little population structure between species. Genome scans for F_ST outliers identified few markers that are highly differentiated between species and repeatable in several populations. These genome scans also revealed highly differentiated polymorphisms in genes with putative involvement in floral odour production, including a previously identified candidate gene thought to be involved in the biosynthesis of pseudo‐pheromones by the orchid flowers. Taken together, these data suggest that ecological speciation associated with different pollinators in sexually deceptive orchids has a genic rather than a genomic basis, placing these species at an early phase of genomic divergence within the ‘speciation continuum’.     

Dissection of the genetic architecture of three seed‐quality traits and consequences for breeding in Brassica napus

Genome‐wide association studies (GWASs) combining high‐throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed‐quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well‐characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed‐quality traits, which would be useful for genetic improvement of B. napus.     

Genome sequence of M6, a diploid inbred clone of the high‐glycoalkaloid‐producing tuber‐bearing potato species Solanum chacoense,reveals residual heterozygosity

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self‐compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high‐confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome‐enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.     

Odour‐mediated orientation behaviour of the lily beetle Lilioceris lilii

The lily beetle Lilioceris lilii (Scopoli) (Coleoptera: Chrysomelidae) feeds on Lilium, Fritillaria and Cardiocrinum plants and is a serious pest in gardens and amenity plantings in parts of Northern Europe and North America. In the present study, the odour‐mediated behaviour of L. lilii is investigated by behavioural bioassays using a linear‐track olfactometer. The behavioural responses of L. lilii to hosts and conspecifics are, at least in part, odour‐mediated and the responses differ with respect to the physiological (reproductive) state of the adult beetle (i.e. pre‐ or post‐diapause). Significantly more diapaused female L. lilii move into air streams containing the odour of intact host plants than into clean air, and move into air streams containing odour of host plants and beetles combined in preference to odour from manually‐damaged host plants. Diapaused females also move into air streams containing odours from intact plants over those from larval‐infested plants. Pre‐diapause males move into the air streams of intact plants rather than L. lilii‐infested plants. Pre‐diapause females show no significant response in any experiment. The data indicate that the odour‐mediated responses of L. lilii are consistent with those known for other chrysomelids that produce a male aggregation pheromone to which reproductive individuals of both sexes respond.     

Identification of quantitative trait loci and candidate genes for an anxiolytic‐like response to ethanol in BXD recombinant inbred strains

Genetic differences in acute behavioral responses to ethanol contribute to the susceptibility to alcohol use disorder and the reduction of anxiety is a commonly reported motive underlying ethanol consumption among alcoholics. Therefore, we studied the genetic variance in anxiolytic‐like responses to ethanol across the BXD recombinant inbred (RI) mouse panel using the light–dark transition model of anxiety. Strain‐mean genetic mapping and a mixed‐model quantitative trait loci (QTL) analysis replicated several previously published QTL for locomotor activity and identified several novel anxiety‐related loci. Significant loci included a chromosome 11 saline anxiety‐like QTL (Salanq1) and a chromosome 12 locus (Etanq1) influencing the anxiolytic‐like response to ethanol. Etanq1 was successfully validated by studies with BXD advanced intercross strains and fine‐mapped to a region comprising less than 3.5 Mb. Through integration of genome‐wide mRNA expression profiles of the mesocorticolimbic reward circuit (prefrontal cortex, nucleus accumbens and ventral midbrain) across the BXD RI panel, we identified high priority candidate genes within Etanq1, the strongest of which was Ninein (Nin), a Gsk3β‐interacting protein that is highly expressed in the brain.     

Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome‐wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As‐rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non‐redundant quantitative trait loci (QTL s) associated with As concentrations were identified (P  ≤ 2.04 × 10^?6) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%–24.65% of the phenotypic variation for each QTL , on average. Additionally, four QTL s [involving 15 single nucleotide polymorphisms (SNP s)] were detected in the single and the combined locations, indicating that these loci/SNP s might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non‐redundant QTL s (6 SNP s), including a QTL that was detected in multiple locations according to the genome‐wide association study, were found to co‐localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.     

Aim‐then‐shoot anemotaxis involved in the hopping approach of potato tuberworm moth Phthorimaea operculella toward a sex pheromone source

Male moths locate conspecific females by pheromone‐induced upwind flight maintained by detecting a visual flow, termed optomotor anemotaxis. Their behavioural pattern is characterized by an upwind surge in response to a pheromone stimulus and crosswind casting after odour loss, which is considered to be reset and restarted on receipt of another pheromone pulse. However, pheromone‐stimulated males of the potato tuberworm moth Phthorimaea operculella exhibit a series of short and straight intermittent flights, or hops, when moving upwind. It is unclear whether they navigate by employing the same behavioural pattern and wind detection mechanism as that used by flying moths. To analyze odour‐modulated anemotaxis in male potato tuberworm moths, a flat wind tunnel is constructed to give regular odour stimuli to an insect regardless of its location. Moths are subjected to pheromone pulses of different frequencies to test whether they show a behavioural pattern that is reset and restarted by a pheromone pulse. Moths on the ground are also subjected to crosswind shear to examine their detection of wind direction. Path analyses reveal that males surge upwind when they receive a pheromone pulse and exhibit casting by successive hops when they lose odour. This behavioural pattern appears to be similar to that of flying moths. When the direction of the airflow is switched orthogonally, males adjust their course angle accordingly when they are on the ground. It is suggested that, instead of optomotor anemotaxis, this ‘aim‐then‐shoot’ system aids the detection of wind direction, possibly by mechanosensory means.     

Sex differences in olfaction‐mediated visual acuity in blowflies and its consequences for gender‐specific trapping

The behaviour of 650 female Calliphora vicina Robineau‐Desvoidy (Diptera: Calliphoridae) was examined in a wind tunnel using odour, in combination with six artificial visual stimuli, ranging from a simple black square to a three‐dimensional model of a dead mouse. The carcasses of laboratory mice were used to provide a natural odour and visual source, and a blend consisting of dimethyl trisulphide, mercaptoethanol, and o‐cresol was used to provide a synthetic lure. Significant differences were found in attraction to these odour sources: 90% of the flies oriented upwind to the natural source and 62% to the synthetic lure. No significant differences were found in upwind orientation towards different visual stimuli, but flies showed significantly more landings if the visual cues provided a vertical contrast against the background. A horizontal contrast gave no difference in landing rate compared to treatments without visual cues. In a field study, the blowfly genera Pollenia, Calliphora, and Lucilia were caught. The overall blowfly catch was significantly higher when an odour lure was present (Pollenia: 3×, Calliphora: 15×, Lucilia: >79×). A significant three‐way interaction between visual cue, genus, and gender was found. The saprophagous Lucilia and Calliphora showed a gender‐specific response to visual stimuli, whereas the parasitic Pollenia did not. A 2:1 female:male sex ratio was found for Calliphora spp. and a 12:1 sex ratio for Lucilia spp. The data suggest that landing responses of male and female saprophagous blowflies, and consequently trap catches, result from olfaction, but also from gender‐specific visual responses when under the influence of odour.     

Genome‐wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome‐wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome‐wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome‐wide associations identified known (Twin Sister of FT, FRIGIDA‐like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage‐Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.     

Whole‐genome resequencing‐based QTL‐seq identified candidate genes and molecular markers for fresh seed dormancy in groundnut

The subspecies fastigiata of cultivated groundnut lost fresh seed dormancy (FSD) during domestication and human‐made selection. Groundnut varieties lacking FSD experience precocious seed germination during harvest imposing severe losses. Development of easy‐to‐use genetic markers enables early‐generation selection in different molecular breeding approaches. In this context, one recombinant inbred lines (RIL) population (ICGV 00350 × ICGV 97045) segregating for FSD was used for deploying QTL‐seq approach for identification of key genomic regions and candidate genes. Whole‐genome sequencing (WGS) data (87.93 Gbp) were generated and analysed for the dormant parent (ICGV 97045) and two DNA pools (dormant and nondormant). After analysis of resequenced data from the pooled samples with dormant parent (reference genome), we calculated delta‐SNP index and identified a total of 10,759 genomewide high‐confidence SNPs. Two candidate genomic regions spanning 2.4 Mb and 0.74 Mb on the B05 and A09 pseudomolecules, respectively, were identified controlling FSD. Two candidate genes—RING‐H2 finger protein and zeaxanthin epoxidase—were identified in these two regions, which significantly express during seed development and control abscisic acid (ABA) accumulation. QTL‐seq study presented here laid out development of a marker, GMFSD1, which was validated on a diverse panel and could be used in molecular breeding to improve dormancy in groundnut.     

Natural variation involving deletion alleles of FRIGIDA modulate temperature‐sensitive flowering responses in Arabidopsis thaliana

Ambient temperature is one of the major environmental factors that modulate plant growth and development. There is extensive natural genetic variation in thermal responses of plants exemplified by the variation exhibited by the accessions of Arabidopsis thaliana. In this work we have studied the enhanced temperature response in hypocotyl elongation and flowering shown by the Tsu‐0 accession in long days. Genetic mapping in the Col‐0 × Tsu‐0 recombinant inbred line (RIL) population identified several QTLs for thermal response including three major effect loci encompassing candidate genes FRIGIDA (FRI), FLOWERING LOCUS C (FLC) and FLOWERING LOCUS T (FT). We confirm and validate these QTLs. We show that the Tsu‐0 FRI allele, which is the same as FRI‐Ler is associated with late flowering but only at lower temperatures in long days. Using transgenic lines and accessions, we show that the FRI‐Ler allele confers temperature‐sensitive late flowering confirming a role for FRI in photoperiod‐dependent thermal response. Through quantitative complementation with heterogeneous inbred families, we further show that cis‐regulatory variation at FT contributes to the observed hypersensitivity of Tsu‐0 to ambient temperature. Overall our results suggest that multiple loci that interact epistatically govern photoperiod‐dependent thermal responses of A. thaliana.     

Quantitative Genetics of Food Intake in Drosophila melanogaster

Food intake is an essential animal activity, regulated by neural circuits that motivate food localization, evaluate nutritional content and acceptance or rejection responses through the gustatory system, and regulate neuroendocrine feedback loops that maintain energy homeostasis. Excess food consumption in people is associated with obesity and metabolic and cardiovascular disorders. However, little is known about the genetic basis of natural variation in food consumption. To gain insights in evolutionarily conserved genetic principles that regulate food intake, we took advantage of a model system, Drosophila melanogaster, in which food intake, environmental conditions and genetic background can be controlled precisely. We quantified variation in food intake among 182 inbred, sequenced lines of the Drosophila melanogaster Genetic Reference Panel (DGRP). We found significant genetic variation in the mean and within-line environmental variance of food consumption and observed sexual dimorphism and genetic variation in sexual dimorphism for both food intake traits (mean and variance). We performed genome wide association (GWA) analyses for mean food intake and environmental variance of food intake (using the coefficient of environmental variation, CV _E, as the metric for environmental variance) and identified molecular polymorphisms associated with both traits. Validation experiments using RNAi-knockdown confirmed 24 of 31 (77%) candidate genes affecting food intake and/or variance of food intake, and a test cross between selected DGRP lines confirmed a SNP affecting mean food intake identified in the GWA analysis. The majority of the validated candidate genes were novel with respect to feeding behavior, and many had mammalian orthologs implicated in metabolic diseases.     

Next‐generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan)

To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing‐based bulked segregant analysis (Seq‐BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R‐ and S‐bulks with the help of draft genome sequence and reference‐guided assembly of ICPL 20096 (resistant parent). Seq‐BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re‐sequenced and their combined analysis with R‐ and S‐bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2‐Mb flanking regions of seven candidate SNPs identified through Seq‐BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re‐sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics‐assisted breeding in pigeonpea.