Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus

Time trees and clock genes: a systematic review and comparative analysis of contemporary avian migration genetics

Timing is a crucial aspect for survival and reproduction in seasonal environments leading to carefully scheduled annual programs of migration in many species. But what are the exact mechanisms through which birds (class: Aves) can keep track of time, anticipate seasonal changes, and adapt their behaviour? One proposed mechanism regulating annual behaviour is the circadian clock, controlled by a highly conserved set of genes, collectively called 'clock genes' which are well established in controlling the daily rhythmicity of physiology and behaviour. Due to diverse migration patterns observed within and among species, in a seemingly endogenously programmed manner, the field of migration genetics has sought and tested several candidate genes within the clock circuitry that may underlie the observed differences in breeding and migration behaviour. Among others, length polymorphisms within genes such as Clock and Adcyap1 have been hypothesised to play a putative role, although association and fitness studies in various species have yielded mixed results. To contextualise the existing body of data, here we conducted a systematic review of all published studies relating polymorphisms in clock genes to seasonality in a phylogenetically and taxonomically informed manner. This was complemented by a standardised comparative re-analysis of candidate gene polymorphisms of 76 bird species, of which 58 are migrants and 18 are residents, along with population genetics analyses for 40 species with available allele data. We tested genetic diversity estimates, used Mantel tests for spatial genetic analyses, and evaluated relationships between candidate gene allele length and population averages for geographic range (breeding- and non-breeding latitude), migration distance, timing of migration, taxonomic relationships, and divergence times. Our combined analysis provided evidence (i) of a putative association between Clock gene variation and autumn migration as well as a putative association between Adcyap1 gene variation and spring migration in migratory species; (ii) that these candidate genes are not diagnostic markers to distinguish migratory from sedentary birds; and (iii) of correlated variability in both genes with divergence time, potentially reflecting ancestrally inherited genotypes rather than contemporary changes driven by selection. These findings highlight a tentative association between these candidate genes and migration attributes as well as genetic constraints on evolutionary adaptation.

The structure of songbird MHC class I reveals antigen binding that is flexible at the N-terminus and static at the C-terminus

Long-distance migratory animals such as birds and bats have evolved to withstand selection imposed by pathogens across the globe, and pathogen richness is known to be particularly high in tropical regions. Immune genes, so-called Major Histocompatibility Complex (MHC) genes, are highly duplicated in songbirds compared to other vertebrates, and this high MHC diversity has been hypothesised to result in a unique adaptive immunity. To understand the rationale behind the evolution of the high MHC genetic diversity in songbirds, we determined the structural properties of an MHC class I protein, Acar3, from a long-distance migratory songbird, the great reed warbler Acrocephalus arundinaceus (in short: Acar). The structure of Acar3 was studied in complex with pathogen-derived antigens and shows an overall antigen presentation similar to human MHC class I. However, the peptides bound to Acar3 display an unusual conformation: Whereas the N-terminal ends of the peptides display enhanced flexibility, the conformation of their C-terminal halves is rather static. This uncommon peptide-binding mode in Acar3 is facilitated by a central Arg residue within the peptide-binding groove that fixes the backbone of the peptide at its central position, and potentially permits successful interactions between MHC class I and innate immune receptors. Our study highlights the importance of investigating the immune system of wild animals, such as birds and bats, to uncover unique immune mechanisms which may neither exist in humans nor in model organisms.

Genetic Control of Avian Migration: Insights from Studies in Latitudinal Passerine Migrants

Twice-a-year, large-scale movement of billions of birds across latitudinal gradients is one of the most fascinating behavioral phenomena seen among animals. These seasonal voyages in autumn southwards and in spring northwards occur within a discrete time window and, as part of an overall annual itinerary, involve close interaction of the endogenous rhythm at several levels with prevailing photoperiod and temperature. The overall success of seasonal migrations thus depends on their close coupling with the other annual sub-cycles, namely those of the breeding, post-breeding recovery, molt and non-migratory periods. There are striking alterations in the daily behavior and physiology with the onset and end of the migratory period, as shown by the phase inversions in behavioral (a diurnal passerine bird becomes nocturnal and flies at night) and neural activities. Interestingly, there are also differences in the behavior, physiology and regulatory strategies between autumn and spring (vernal) migrations. Concurrent molecular changes occur in regulatory (brain) and metabolic (liver, flight muscle) tissues, as shown in the expression of genes particularly associated with 24 h timekeeping, fat accumulation and the overall metabolism. Here, we present insights into the genetic basis of migratory behavior based on studies using both candidate and global gene expression approaches in passerine migrants, with special reference to Palearctic-Indian migratory blackheaded and redheaded buntings.

Inversions maintain differences between migratory phenotypes of a songbird

Structural rearrangements have been shown to be important in local adaptation and speciation, but have been difficult to reliably identify and characterize in non-model species. Here we combine long reads, linked reads and optical mapping to characterize three divergent chromosome regions in the willow warbler Phylloscopus trochilus, of which two are associated with differences in migration and one with an environmental gradient. We show that there are inversions (0.4-13 Mb) in each of the regions and that the divergence times between inverted and non-inverted haplotypes are similar across the regions (~1.2 Myrs), which is compatible with a scenario where inversions arose in either of two allopatric populations that subsequently hybridized. The improved genomes allow us to detect additional functional differences in the divergent regions, providing candidate genes for migration and adaptations to environmental gradients.

Clock-linked genes underlie seasonal migratory timing in a diurnal raptor

Seasonal migration is a dynamic natural phenomenon that allows organisms to exploit favourable habitats across the annual cycle. While the morphological, physiological and behavioural changes associated with migratory behaviour are well characterized, the genetic basis of migration and its link to endogenous biological time-keeping pathways are poorly understood. Historically, genome-wide research has focused on genes of large effect, whereas many genes of small effect may work together to regulate complex traits like migratory behaviour. Here, we explicitly relax stringent outlier detection thresholds and, as a result, discover how multiple biological time-keeping genes are important to migratory timing in an iconic raptor species, the American kestrel (Falco sparverius). To validate the role of candidate loci in migratory timing, we genotyped kestrels captured across autumn migration and found significant associations between migratory timing and genetic variation in metabolic and light-input pathway genes that modulate biological clocks (top1, phlpp1, cpne4 and peak1). Further, we demonstrate that migrating individuals originated from a single panmictic source population, suggesting the existence of distinct early and late migratory genotypes (i.e. chronotypes). Overall, our results provide empirical support for the existence of a within-population-level polymorphism in genes underlying migratory timing in a diurnally migrating raptor.

Sexual-biased gene expression of olfactory-related genes in the antennae of Conogethes pinicolalis (Lepidoptera: Crambidae)

BACKGROUND

Conogethes pinicolalis (Lepidoptera: Crambidae), is similar to Conogethes punctiferalis (yellow peach moth) and its host plant is gymnosperms, especially for masson pine. So far, less literature was reported on this pest. In the present study, we sequenced and characterized the antennal transcriptomes of male and female C. pinicolalis for the first time.

RESULTS

Totally, 26 odorant-binding protein (OBP) genes, 19 chemosensory protein (CSP) genes, 55 odorant receptor (OR) genes and 20 ionotropic receptor (IR) genes were identified from the C. pinicolalis antennae transcriptome and amino sequences were annotated against homologs of C. punctiferalis. The neighbor-joining tree indicated that the amino acid sequence of olfactory related genes is highly homologous with C. punctiferalis. Furthermore, the reference genes were selected, and we recommended the phosphate dehydrogenase gene (GAPDH) or ribosomal protein 49 gene (RP49) to verify the target gene expression during larval development stages and RP49 or ribosomal protein L13 gene (RPL13) for adult tissues.

CONCLUSIONS

Our study provides a starting point on the molecular level characterization between C. pinicolalis and C. punctiferalis, which might be supportive for pest management studies in future.

Population Genetics and Speciation of Yellow-Bellied, Red-Naped, and Red-Breasted Sapsuckers (Sphyrapicus varius, S. nuchalis, and S. ruber)

The root of understanding speciation lies in determining the forces which drive it. In many closely-related species, including Sphyrapicus varius, S. nuchalis, and S. ruber, it is assumed that speciation occurred due to isolation in multiple Pleistocene refugia. We used genetic data from 457 samples at the control region (CR), cytochrome oxidase I (COI), and chromo-helicase DNA binding protein (CHD1Z) to examine range-wide population genetic structure and differentiation amongst these 3 species across each species' breeding range. In addition, we modeled these species' ecological niches for the Holocene (~6000 years ago), Last Glacial Maximum (~22000 years ago), and Last Interglacial (~120000-140000 years ago) to determine if Pleistocene glaciations could have contributed to allopatric distributions, therefore allowing these groups to differentiate. Population genetic data show a potential Pleistocene refugium in Haida Gwaii, an east-west split among S. varius, and low genetic differentiation within each species. Our CR data show some polyphyly, while COI and CHD1Z data show differentiation among species using composite genotypes. Ecological Niche Modeling shows a large amount of niche overlap at each time period suggesting that S. varius, S. nuchalis, and S. ruber may not have been completely allopatric, and these species likely had repeated intermittent contact. Our data support the growing body of research that suggests differentiation despite gene flow.

Length polymorphisms at two candidate genes explain variation of migratory behaviors in blackpoll warblers (Setophaga striata)

Migratory behaviors such as the timing and duration of migration are genetically inherited and can be under strong natural selection, yet we still know very little about the specific genes or molecular pathways that control these behaviors. Studies in candidate genes Clock and Adcyap1 have revealed that both of these loci can be significantly correlated with migratory behaviors in birds, though observed relationships appear to vary across species. We investigated geographic genetic structure of Clock and Adcyap1 in four populations of blackpoll warblers (Setophaga striata), a Neotropical-Nearctic migrant that exhibits geographic variation in migratory timing and duration across its boreal breeding distribution. Further, we used data on migratory timing and duration, obtained from light-level geolocator trackers to investigate candidate genotype-phenotype relationships at the individual level. While we found no geographic structure in either candidate gene, we did find evidence that candidate gene lengths are correlated with five of the six migratory traits. Maximum Clock allele length was significantly and negatively associated with spring arrival date. Minimum Adcyap1 allele length was significantly and negatively associated with spring departure date and positively associated with fall arrival date at the wintering grounds. Additionally, we found a significant interaction between Clock and Adcyap1 allele lengths on both spring and fall migratory duration. Adcyap1 heterozygotes also had significantly shorter migration duration in both spring and fall compared to homozygotes. Our results support the growing body of evidence that Clock and Adcyap1 allele lengths are correlated with migratory behaviors in birds.

Signatures of Selection on Standing Genetic Variation Underlie Athletic and Navigational Performance in Racing Pigeons

Racing pigeons have been selectively bred to find their way home quickly over what are often extremely long distances. This breed is of substantial commercial value and is also an excellent avian model to gain empirical insights into the evolution of traits associated with flying performance and spatial orientation. Here, we investigate the molecular basis of the superior athletic and navigational capabilities of racing pigeons using whole-genome and RNA sequencing data. We inferred multiple signatures of positive selection distributed across the genome of racing pigeons. The strongest signature overlapped the CASK gene, a gene implicated in the formation of neuromuscular junctions. However, no diagnostic alleles were found between racing pigeons and other breeds, and only a small proportion of highly differentiated variants were exclusively detected in racing pigeons. We can thus conclude that very few individual genetic changes, if any, are either strictly necessary or sufficient for superior athletics and navigation. Gene expression analysis between racing and nonracing breeds revealed modest differences in muscle (213) and brain (29). These transcripts, however, showed only slightly elevated levels of genetic differentiation between the two groups, suggesting that most differential expression is not causative but likely a consequence of alterations in regulatory networks. Our results show that the unique suite of traits that enable fast flight, long endurance, and accurate navigation in racing pigeons, do not result from few loci acting as master switches but likely from a polygenic architecture that leveraged standing genetic variation available at the onset of the breed formation.

Transcriptional Signatures in Liver Reveal Metabolic Adaptations to Seasons in Migratory Blackheaded Buntings

The molecular underpinnings of metabolic adaptation to seasons are poorly understood in long- distance migrants. We measured changes in physiology and performed de novo sequencing of RNA extracted from liver samples collected at 4-h intervals over a period of 24 h from a long-distance avian migrant, the blackheaded bunting (Emberiza melanocephala), during two states: photostimulated vernal migratory (M) state and photorefractory non-migratory (nM) state. The M state was differentiated from the nM state based on body fattening and weight gain, as well as on Zugunruhe, that is, nocturnal migratory restlessness in caged birds. We found that baseline blood glucose and triglyceride levels were significantly higher in the M state than the nM state; conversely, surface body temperature was higher in the nM state than the M state. In a total of 6 liver samples that were sequenced from each state, 11,246 genes were annotated, including 4448 genes that were cyclic over 24 h. We found 569 differentially expressed genes (DEGs) between the M and the nM state, and the M state showed 131 upregulated and 438 downregulated genes. These DEGs formed core gene hubs associated with specific biological processes in both the states. In addition, weighted gene coexpression network analysis revealed two discrete modules of coexpressed genes, with a significant difference in the expression pattern of metab olism-associated genes between M and nM states. These results demonstrate, for the first time, transcriptome-wide changes in the liver between two distinct physiological states and give molecular insights into seasonal metabolic adaptations in latitudinal migrants.

Hypothalamic and liver transcriptome from two crucial life-history stages in a migratory songbird

NEW FINDINGS

What is the central question of this study? What are the molecular underpinnings of the seasonal adaptation in a latitudinal migratory songbird? What is the main finding and its importance? We found changes in mRNA levels after a photoperiod-induced alteration of seasonal state in a captive long-distance latitudinal avian migrant. The hypothalamus and liver transcriptomes revealed genes involved in the regulatory and functional pathways between non-migratory and migratory states. Our results provide insights into mechanisms underlying homeostasis during seasonal changes that are conserved across most species, including humans.

ABSTRACT

Very little is understood about genetic mechanisms underlying the onset of spring migration in latitudinal avian migrants. To gain insight into the genetic architecture of the hypothalamus and liver tissues of a long-distance migrant, we examined and compared the transcriptome profile of captive night-migratory black-headed buntings (Emberiza melanocephala) between photoperiod-induced winter non-migratory (WnM) and spring migratory (SM) life-history states under short and long days, respectively. High-throughput 454 pyrosequenced transcripts were mapped initially with reference to the genome of two phylogenetically close species, Taeniopygia guttata and Ficedula albicollis. The F. albicollis genome gave higher annotation results and was used for further analysis. A total of 216 (78 in hypothalamus; 138 in liver) genes were found to be expressed differentially between the WnM and SM life-history states. These genes were enriched for physiological pathways that might be involved in the regulation of seasonal migrations in birds. For example, genes for the ATP binding pathway in the hypothalamus were expressed at a significantly higher level in SM than in the WnM life-history state. Likewise, upregulated genes associated with the myelin sheath and focal adhesion were enriched in the hypothalamus, and those with cell-to-cell junction, intracellular protein transport, calcium ion transport and small GTPase-mediated signal transduction were enriched in the liver. Many of these genes are a part of physiological pathways potentially involved in the regulation of seasonal migration in birds. These results show molecular changes at the regulatory and metabolic levels associated with seasonal transitions in a long-distance migrant and provide the basis for future studies aimed at unravelling the genetic control of migration in birds.

Genetic differences between willow warbler migratory phenotypes are few and cluster in large haplotype blocks

It is well established that differences in migratory behavior between populations of songbirds have a genetic basis but the actual genes underlying these traits remains largely unknown. In an attempt to identify such candidate genes we de novo assembled the genome of the willow warbler Phylloscopus trochilus, and used whole‐genome resequencing and a SNP array to associate genomic variation with migratory phenotypes across two migratory divides around the Baltic Sea that separate SW migrating P. t. trochilus wintering in western Africa and SSE migrating P. t. acredula wintering in eastern and southern Africa. We found that the genomes of the two migratory phenotypes lack clear differences except for three highly differentiated regions located on chromosomes 1, 3, and 5 (containing 146, 135, and 53 genes, respectively). Within each migratory phenotype we found virtually no differences in allele frequencies for thousands of SNPs, even when comparing geographically distant populations breeding in Scandinavia and Far East Russia (>6000 km). In each of the three differentiated regions, multidimensional scaling‐based clustering of SNP genotypes from more than 1100 individuals demonstrates the presence of distinct haplotype clusters that are associated with each migratory phenotype. In turn, this suggests that recombination is absent or rare between haplotypes, which could be explained by inversion polymorphisms. Whereas SNP alleles on chromosome 3 correlate with breeding altitude and latitude, the allele distribution within the regions on chromosomes 1 and 5 perfectly matches the geographical distribution of the migratory phenotypes. The most differentiated 10 kb windows and missense mutations within these differentiated regions are associated with genes involved in fatty acid synthesis, possibly representing physiological adaptations to the different migratory strategies. The ∼200 genes in these regions, of which several lack described function, will direct future experimental and comparative studies in the search for genes that underlie important migratory traits.

Animal tracking meets migration genomics: transcriptomic analysis of a partially migratory bird species

Seasonal migration is a widespread phenomenon, which is found in many different lineages of animals. This spectacular behaviour allows animals to avoid seasonally adverse environmental conditions to exploit more favourable habitats. Migration has been intensively studied in birds, which display astonishing variation in migration strategies, thus providing a powerful system for studying the ecological and evolutionary processes that shape migratory behaviour. Despite intensive research, the genetic basis of migration remains largely unknown. Here, we used state-of-the-art radio-tracking technology to characterize the migratory behaviour of a partially migratory population of European blackbirds (Turdus merula) in southern Germany. We compared gene expression of resident and migrant individuals using high-throughput transcriptomics in blood samples. Analyses of sequence variation revealed a nonsignificant genetic structure between blackbirds differing by their migratory phenotype. We detected only four differentially expressed genes between migrants and residents, which might be associated with hyperphagia, moulting and enhanced DNA replication and transcription. The most pronounced changes in gene expression occurred between migratory birds depending on when, in relation to their date of departure, blood was collected. Overall, the differentially expressed genes detected in this analysis may play crucial roles in determining the decision to migrate, or in controlling the physiological processes required for the onset of migration. These results provide new insights into, and testable hypotheses for, the molecular mechanisms controlling the migratory phenotype and its underlying physiological mechanisms in blackbirds and other migratory bird species.

Patterns of genetic, phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.

Antennal transcriptome and differential expression of olfactory genes in the yellow peach moth, Conogethes punctiferalis (Lepidoptera: Crambidae)

The yellow peach moth (YPM), Conogethes punctiferalis (Guenée), is a multivoltine insect pest of crops and fruits. Antennal-expressed receptors are important for insects to detect olfactory cues for host finding, mate attraction and oviposition site selection. However, few olfactory related genes were reported in YPM until now. In the present study, we sequenced and characterized the antennal transcriptomes of male and female YPM. In total, 15 putative odorant binding proteins (OBPs), 46 putative odorant receptors (ORs) and 7 putative ionotropic receptors (IRs) were annotated and identified as olfactory-related genes of C. punctiferalis. Further analysis of RT-qPCR revealed that all these olfactory genes are primarily or uniquely expressed in male and female antennae. Among which, 3 OBPs (OBP4, OBP8 and PBP2) and 4 ORs (OR22, OR26, OR44 and OR46) were specially expressed in male antennae, whereas 4 ORs (OR5, OR16, OR25 and OR42) were primarily expressed in female antennae. The predicted protein sequences were compared with homologs in other lepidopteran species and model insects, which showed high sequence homologies between C. punctiferalis and O. furnacalis. Our work allows for further functional studies of pheromone and general odorant detection genes, which might be meaningful targets for pest management.

Gene expression in the brain of a migratory songbird during breeding and migration

BACKGROUND

We still have limited knowledge about the underlying genetic mechanisms that enable migrating species of birds to navigate the globe. Here we make an attempt to get insight into the genetic architecture controlling this complex innate behaviour. We contrast the gene expression profiles of two closely related songbird subspecies with divergent migratory phenotypes. In addition to comparing differences in migratory strategy we include a temporal component and contrast patterns between breeding adults and autumn migrating juvenile birds of both subspecies. The two willow warbler subspecies, Phylloscopus trochilus trochilus and P. t. acredula, are remarkably similar both in phenotype and genotype and have a narrow contact zone in central Scandinavia. Here we used a microarray gene chip representing 23,136 expressed sequence tags (ESTs) from the zebra finch Taeniopygia guttata to identify mRNA level differences in willow warbler brain tissue in relation to subspecies and season.

RESULTS

Out of the 22,109 EST probe sets that remained after filtering poorly binding probes, we found 11,898 (51.8 %) probe sets that could be reliably and uniquely matched to a total of 6,758 orthologous zebra finch genes. The two subspecies showed very similar levels of gene expression with less than 0.1 % of the probe sets being significantly differentially expressed. In contrast, 3,045 (13.8 %) probe sets were found to be differently regulated between samples collected from breeding adults and autumn migrating juvenile birds. The genes found to be differentially expressed between seasons appeared to be enriched for functional roles in neuronal firing and neuronal synapse formation.

CONCLUSIONS

Our results show that only few genes are differentially expressed between the subspecies. This suggests that the different migration strategies of the subspecies might be governed by few genes, or that the expression patterns of those genes are time-structured or tissue-specific in ways, which our approach fails to uncover. Our findings will be useful in the planning of new experiments designed to unravel the genes involved in the migratory program of birds.

Characterization of early transcriptional responses to cadmium in the root and leaf of Cd-resistant Salix matsudana Koidz

BACKGROUND

Salix matsudana Koidz. is a fast growing tree species. It has a high cadmium (Cd) tolerance capacity, making it potentially suitable for phytoremediation. Presently, transcriptomic and physiological Cd response mechanisms are poorly understood. Transcriptomic analysis in early response to high (50 μM) Cd levels was investigated in leaf and root of Cd-resistant S. matsudana Koidz..

RESULTS

Analysis of the response profiles demonstrate the existence of a complex transcriptional network in the root and leaf when exposed to Cd. The main response in the root involved up-regulation of genes associated with defence response via callose deposition in the cell wall and cell wall thickening. In the leaf, transcripts related to biotic stress signalling and secondary metabolism were activated. Additionally, many lignin and brassinosteroids synthesis pathway genes were induced mainly in the leaf, indicating that gene response to Cd was tissue-specific. The Cd transcriptome results were consistent with observed physiological changes.

CONCLUSION

The sub-localization, transcriptional network, and physiological regulation demonstrate the tissue-specific manner of Cd response, and provide a novel insight into in early response of tree species to Cd exposure.

No evidence for assortative mating within a willow warbler migratory divide

Introduction

In contact zones, genetic mixing of two taxa can be restricted by prezygotic (e.g. assortative mating) or postzygotic (lower fitness of hybrid offspring) barriers, or a combination of the two. A hybrid zone between two willow warbler subspecies (Phylloscopus trochilus trochilus, P. t. acredula) with distinctive migratory strategies occurs in central Sweden. These subspecies exhibit differences in migratory direction and distance, resulting in geographically distinct wintering areas in Africa. The subspecies may have diverged from a common refuge after the last ice age, and neutral genetic markers are homogeneous across their range. By contrast, several phenotypic traits and genetic markers of two chromosomal regions previously identified show steep clines across the divide. The evolutionary forces that maintain this migratory divide remain unknown. Here we use plumage colour, morphology, genetic markers and feather stable nitrogen-isotopes (δ¹⁵N) to assess if assortative mating between migratory phenotypes could be acting as a possible mechanism for keeping the two forms genetically separate and maintaining the migratory divide. We colour-ringed a willow warbler breeding population in the central part of the hybrid zone and observed the breeding population to assess phenotypic and genotypic traits of social pairs.

Results

Our data suggest that wintering area and genetic ancestry had an effect on male arrival time to the breeding grounds which could contribute to assortment. However, evidence for assortative mating could not be detected based on a comparison of plumage colour, morphology and δ¹⁵N between social mates.

Conclusion

This finding was strengthened by analyses of subspecies-specific genetic markers, which allowed us to identify the presence of a large proportion of potential hybrids and backcrosses at the study site. Our results supported the hypothesis that pre-mating isolation in willow warblers is weak, resulting in extensive hybridisation across the migratory divide.

De novo analysis of the Nilaparvata lugens (Stål) antenna transcriptome and expression patterns of olfactory genes

We sequenced the antenna transcriptome of the brown planthopper (BPH), Nilaparvata lugens (Stål), a global rice pest, and performed transcriptome analysis on BPH antenna. We obtained about 40million 90bp reads that were assembled into 75,874 unigenes with a mean size of 456bp. Among the antenna transcripts, 32,856 (43%) showed significant similarity (E-value <1e(-5)) to known proteins in the NCBI database. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to classify functions of BPH antenna genes. We identified 10 odorant-binding proteins (OBPs), including 7 previously unidentified, and 11 chemosensory proteins (CSPs), including two new members. The expression profiles of 4 OBPs and 2 CSPs were determined by q-PCR for antenna, abdomen, leg and wing of insects of different age, gender, and mating status including two BPH adult wing-morphology types. NlugCSP10 and 4 OBPs appeared to be antenna-specific because they were highly and differentially expressed in male and female antennae. NlugCSP11 was expressed ubiquitously, with particularly high expression in wings. The transcript levels of several olfactory genes depended on adult wing form, age, gender, and mating status, although no clear expression patterns were determined.