The Evolution of Innate Immune Genes: Purifying and Balancing Selection on ß-Defensins in Waterfowl.

In disease dynamics, high immune gene diversity can confer a selective advantage to hosts in the face of a rapidly evolving and diverse pathogen fauna. This is supported empirically for genes involved in pathogen recognition and signalling. In contrast, effector genes involved in pathogen clearance may be more constrained. ß-Defensins are innate immune effector genes; their main mode of action is via disruption of microbial membranes. Here, five ß-defensin genes were characterized in mallards (Anas platyrhynchos) and other waterfowl; key reservoir species for many zoonotic diseases. All five genes showed remarkably low diversity at the individual-, population-, and species-level. Furthermore, there was widespread sharing of identical alleles across species divides. Thus, specific ß-defensin alleles were maintained not only spatially but also over long temporal scales, with many amino acid residues being fixed across all species investigated. Purifying selection to maintain individual, highly efficacious alleles was the primary evolutionary driver of these genes in waterfowl. However, we also found evidence for balancing selection acting on the most recently duplicated ß-defensin gene (AvBD3b). For this gene, we found that amino acid replacements were more likely to be radical changes, suggesting that duplication of ß-defensin genes allows exploration of wider functional space. Structural conservation to maintain function appears to be crucial for avian ß-defensin effector molecules, resulting in low tolerance for new allelic variants. This contrasts with other types of innate immune genes, such as receptor and signalling molecules, where balancing selection to maintain allelic diversity has been shown to be a strong evolutionary force.

DNA double-strand breaks in incubating female common eiders (Somateria mollissima): Comparison between a low and a high polluted area.

Alterations in the genetic material may have severe consequences for individuals and populations. Hence, genotoxic effects of environmental exposure to pollutants are of great concern. We assessed the impact of blood concentrations of persistent organic pollutants (POPs) and mercury (Hg) on DNA double-strand break (DSB) frequency, in blood cells of a high-exposed Baltic, and lower exposed Arctic population of common eiders (Somateria mollissima). Furthermore, we examined whether the genotoxic response was influenced by antioxidant concentration (plasma total glutathione (tGSH) and total antioxidant capacity) and female body mass. The DNA DSB frequency did not differ between the two populations. We found significant positive relationships between Hg and DNA DSB frequency in Baltic, but not in Arctic eiders. Although both p,p'-DDE and PCB 118 had a lesser effect than Hg, they exhibited a positive association with DNA DSB frequency in Baltic eiders. Antioxidant levels were not important for the genotoxic effect, suggesting alternative mechanisms other than GSH depletion for the relationship between Hg and DNA DSBs. Hence, the Baltic population, which is considered to be endangered and is under the influence of several environmental stressors, may be more susceptible to genotoxic effects of environmental exposure to Hg than the Arctic population.

Variation at the DRD4 locus is associated with wariness and local site selection in urban black swans.

BACKGROUND: Interactions between wildlife and humans are increasing. Urban animals are often less wary of humans than their non-urban counterparts, which could be explained by habituation, adaptation or local site selection. Under local site selection, individuals that are less tolerant of humans are less likely to settle in urban areas. However, there is little evidence for such temperament-based site selection, and even less is known about its underlying genetic basis. We tested whether site selection in urban and non-urban habitats by black swans (Cygnus atratus) was associated with polymorphisms in two genes linked to fear in animals, the dopamine receptor D4 (DRD4) and serotonin transporter (SERT) genes. RESULTS: Wariness in swans was highly repeatable between disturbance events (repeatability = 0.61) and non-urban swans initiated escape from humans earlier than urban swans. We found no inter-individual variation in the SERT gene, but identified five DRD4 genotypes and an association between DRD4 genotype and wariness. Individuals possessing the most common DRD4 genotype were less wary than individuals possessing rarer genotypes. As predicted by the local site selection hypothesis, genotypes associated with wary behaviour were over three times more frequent at the non-urban site. This resulted in moderate population differentiation at DRD4 (FST = 0.080), despite the sites being separated by only 30 km, a short distance for this highly-mobile species. Low population differentiation at neutrally-selected microsatellite loci and the likely occasional migration of swans between the populations reduces the likelihood of local site adaptations. CONCLUSION: Our results suggest that wariness in swans is partly genetically-determined and that wary swans settle in less-disturbed areas. More generally, our findings suggest that site-specific management strategies may be necessary that consider the temperament of local animals.

Avian sarcoma and leukosis virus gag gene in the Anser anser domesticus genome.

Endogenous retroviruses are regarded as ideal genetic markers for evolutionary analyses. Birds were some of the initial vertebrates found to contain endogenous retroviruses. However, few studies have investigated the presence and distribution of endogenous retroviruses in goose. In this study, we detected the avian sarcoma and leukosis virus gag gene in the genomic DNA of 8 Chinese native breeds using polymerase chain reaction method. The results indicated that a 1.2-kb avian sarcoma and leukosis virus gag sequence was integrated into all 8 goose breeds. The mean genetic pairwise distance was 0.918% among the investigated geese. To the best of our knowledge, this is the first report demonstrating the presence of the endogenous retroviruses in the domestic goose genome. The genetic structure should be further examined in the domestic goose.

Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds.

BACKGROUND: Toll-like receptors (TLR) are essential activators of the innate part of the vertebrate immune system. In this study, we analysed the interspecific variability of three TLR (bacterial-sensing TLR4 and TLR5 and viral-sensing TLR7) within the Galloanserae bird clade, investigated their phylogeny, assessed their structural conservation and estimated site-specific selection pressures. RESULTS: Physiochemical properties varied according to the TLR analysed, mainly with regards to the surface electrostatic potential distribution. The predicted ligand-binding features (mainly in TLR4 and TLR5) differed between the avian proteins and their fish and mammalian counterparts, but also varied within the Galloanserae birds. We identified 20 positively selected sites in the three TLR, among which several are topologically close to ligand-binding sites reported for mammalian and fish TLR. We described 26, 28 and 25 evolutionarily non-conservative sites in TLR4, TLR5 and TLR7, respectively. Thirteen of these sites in TLR4, and ten in TLR5 were located in functionally relevant regions. The variability appears to be functionally more conserved for viral-sensing TLR7 than for the bacterial-sensing TLR. Amino-acid positions 268, 270, 343, 383, 444 and 471 in TLR4 and 180, 183, 209, 216, 264, 342 and 379 in TLR5 are key candidates for further functional research. CONCLUSIONS: Host-pathogen co-evolution has a major effect on the features of host immune receptors. Our results suggest that avian and mammalian TLR may be differentially adapted to pathogen-derived ligand recognition. We have detected signatures of positive selection even within the Galloanserae lineage. To our knowledge, this is the first study to depict evolutionary pressures on Galloanserae TLR and to estimate the validity of current knowledge on TLR function (based on mammalian and chicken models) for non-model species of this clade.

DNA double-strand breaks in relation to persistent organic pollutants in a fasting seabird.

Lipophilic persistent organic pollutants (POPs) are released from fat reserves during fasting, causing increased blood concentrations. Thus, POPs represent a potential anthropogenic stressor during fasting periods. We analysed the blood of female common eiders (Somateria mollissima) by using agarose gel electrophoresis and image data analysis to quantify the DNA-fraction, of total DNA, that migrated into the gel (DNA-FTM) as a relative measure of DNA double strand-breaks (DSBs) during the fasting incubation period in the high arctic. In 2008 and in 2009 blood samples were obtained for analysis of 9 POPs and DNA-FTM at day 5 of the incubation period, and then in the same individuals at day 20. This unique study design gave us the opportunity to analyse the same individuals throughout two points in time, with low and high stress burdens. During the incubation period the body mass (BM) decreased by 21-24%, whereas the POP levels increased by 148-639%. The DNA-FTM increased by 61-67% (being proportional to the increase in DSBs). At day 5, but not day 20, DNA-FTM was positively correlated with most analysed POPs. The increase in DNA-FTM was positively correlated with the decrease in BM (g) during incubation. Thus, we suggest that fasting stress (BM loss) decreases DNA integrity and that stress caused by fasting on BM loss appeared to override the additional stress caused by concurrent increase in levels of the analysed POPs in the eiders. Blood levels of POPs in the eiders in Svalbard were relatively low, and additive and/or synergistic genotoxic effects of fasting stress and POP exposure may occur in populations with higher POP levels.

Karyotype and C-banding pattern of the domestic geese Anser anser populations (Aves: Anatidae) in Egypt.

The karyotype and C-banding pattern of domestic Greylag geese Anser anser anser populations collected from five localities in El Minia, Egypt, that have either whitish grey or white feather color patterns were described. All populations have a diploid number of 2n= 80 chromosomes. Of the 80 chromosomes, 10 pairs, including ZW chromosomes, were macrochromosomes and the remaining 30 pairs were microchromosomes. Slight variation in the size of macrochromosomes was observed amongst populations. However, obvious variation of C-banding distribution was found and attributed to variation of euchromatin content and its correlation with chromosome size and arrangement of constitutive heterochromatin. Nevertheless, significant variation in the mean number of C-heterochromatin blocks in microchromosomes was attributed to either transformation of heterochromatin into euchromatin and vice versa or to involvement of structural chromosomal aberrations during karyotype evolution. The present results show that A. anser populations common in Egypt could be distinguished from those of A. anser and A. cygnoides occurring elsewhere in Europe and Asiaviavariability in chromosome morphology of pairs nos. 2, 3 and 4.

Genetic diversity of the endangered scaly-sided merganser (Mergus squamatus) in the wintering habitat of central-southern China.

The scaly-sided merganser (Mergus squamatus), found in temperate East Asia, has been reduced to a very small population. Central and southern China are its main wintering habitat. However, populations have declined greatly since the 1980s due to habitat loss and degradation, and poaching. To meet the urgent need for up-to-date conservation information, we examined RAPD DNA markers from 156 specimens in 6 populations in Jiangxi Province. We found that genetic diversity (based on individual similarities) is in fact low; molecular variance between populations ranged from 0.137 to 0.347. Genetic similarity ranged from 0.683 to 0.866. In conclusion, the geographical pattern of genetic diversity supports the long-term refugial status of the scaly-sided merganser in central-southern China; strong conservation measures should be taken to maintain the merganser in this region.

A high-quality chromosomal-level genome assembly of Greater Scaup (Aythya marila).

Aythya marila is one of the few species of Anatidae, and the only Aythya to live in the circumpolar. However, there is a relative lack of research on genetics of this species. In this study, we reported and assembled the first high-quality chromosome-level genome assembly of A. marila. This genome was assembled using Nanopore long reads, and errors corrected using Illumina short reads, with a final genome size of 1.14 Gb, scaffold N50 of 85.44 Mb, and contig N50 of 32.46 Mb. 106 contigs were clustered and ordered onto 35 chromosomes based on Hi-C data, covering approximately 98.28% of the genome. BUSCO assessment showed that 97.0% of the highly conserved genes in aves_odb10 were present intact in the genome assembly. In addition, a total of 154.94 Mb of repetitive sequences were identified. 15,953 protein-coding genes were predicted in the genome, and 98.96% of genes were functionally annotated. This genome will be a valuable resource for future genetic diversity and genomics studies of A. marila.

Bacterial composition along the digestive tract of the Horned Screamer (Anhima cornuta), a tropical herbivorous bird.

Background: The Horned Screamer (Anhima cornuta) is an herbivorous bird that inhabits wetlands of the South American tropical region. We hypothesize that due to its herbivorous niche, its digestive tract compartments may have bacteria specialized in fermenting complex plant carbohydrates. To test this hypothesis, we compared the bacterial communities along the gastrointestinal tract (GIT) of a Horned Screamer captured in Venezuela. Methods: Samples were taken from tissues and content of the proventriculus and the small intestine (considered for this study as upper GIT), and the large intestine and cecum (lower GIT). The bacterial community was characterized by sequencing the V4 region of the 16S rRNA gene. Bioinformatic analysis was performed using QIIME, QIITA and Microbiome Analyst. The association between microbial taxonomy and function was analyzed using their Greengenes OTU IDs and a custom KEGG BRITE hierarchical tree and visualized with BURRITO. Results: The Screamer's gastrointestinal microbiota was composed by seven phyla being Firmicutes and Bacteroidetes the most predominant. The dominant taxa in the upper GIT were Helicobacter, Vibrio, Enterobacter, Acinetobacter and Staphylococcus. The dominant taxa in the lower GIT were Oribacterium, Blautia, Roseburia, Ruminococcus, Desulfovibrio, Intestinimonas, Marvinbryantia and Parabacteroides. Complete degradation of cellulose to the end-products acetate, propanoate, butanoate and acetoacetate was found in the upper and lower GIT without significant differences. Conclusion: Our study confirmed changes in bacterial community composition throughout the GIT of the Horned Screamer primarily associated with the production of metabolic end-products of carbohydrate digestion essential for the fermentation of the herbivorous diet.

Kin association during brood care in a facultatively social bird: active discrimination or by-product of partner choice and demography?

Intra-group relatedness does not necessarily imply kin selection, a leading explanation for social evolution. An overlooked mechanism for generating population genetic structure is variation in longevity and fecundity, referred to as individual quality, affecting kin structure and the potential for cooperation. Individual quality also affects choosiness in partner choice, a key process explaining cooperation through direct fitness benefits. Reproductive skew theory predicts that relatedness decreases with increasing group size, but this relationship could also arise because of quality-dependent demography and partner choice, without active kin association. We addressed whether brood-rearing eider (Somateria mollissima) females preferentially associated with kin using a 6-year data set with individuals genotyped at 19 microsatellite loci and tested whether relatedness decreased with increasing female group size. We also determined the relationship between local relatedness and indices of female age and body condition. We further examined whether the level of female intracoalition relatedness differed from background relatedness in any year. As predicted, median female intra-group relatedness decreased with increasing female group size. However, the proportion of related individuals increased with advancing female age, and older females prefer smaller brood-rearing coalitions, potentially producing a negative relationship between group size and relatedness. There were considerable annual fluctuations in the level of relatedness between coalition-forming females, and in 1year this level exceeded that expected by random association. Thus, both passive and active mechanisms govern kin associations in brood-rearing eiders. Eiders apparently can discriminate between kin, but the benefits of doing so may vary over time.

Breeding site fidelity and winter admixture in a long-distance migrant, the tufted duck (Aythya fuligula).

Long-distance migrants are, by definition, highly mobile but it is poorly understood if this leads to high rates of gene flow and an essentially panmictic global population structure. Genetic divergence in migratory species could be promoted, for example, by fidelity to distinct migratory pathways. In this study, we investigate the population genetic structure of tufted duck (Aythya fuligula), a long-distance migrant with a largely continuous breeding distribution across Eurasia. Distinct, longitudinally oriented flyways have been postulated based on geographically disjunct wintering areas and are supported by evidence from ringing data. We generated sequences of the mitochondrial control region and multi-locus microsatellite genotypes for several hundreds of samples from the European and Asian breeding and wintering grounds including some individuals infected with highly pathogenic avian influenza virus H5N1. Significant differentiation between breeding sites was observed for both marker types, but F(ST) values were approximately 10 times higher for maternally inherited mitochondrial DNA than for biparentally transmitted nuclear markers. The genetic differentiation between the postulated European and Asian flyways was similar to that observed within continents and, in general, genetic divergence was not associated with geographic distance. Neither marker type showed evidence of genetic substructure among aggregations on the European wintering grounds. Our results suggest some breeding site fidelity, especially in females, but extensive population admixture on the wintering grounds. Several scenarios may explain the observed lack of genetic divergence between Europe and Asia including non-equilibrium conditions following a recent range expansion or contemporary gene flow across the postulated migratory divides.

Molecular tools for species and sex identification in the mixed-species flocks of bean geese and white-fronted geese.

Genetic studies on protected species can be difficult, particularly when they form a mixed-species flock with other species. The bean goose (Anser fabalis), which is internationally recognized as a threatened species, was observed to form overwintering foraging flocks with white-fronted goose (Anser albifrons) at agricultural lands to feed on the grains in Korea. Non-invasive samples such as feces and feather that are readily available in their foraging ground can be useful for understanding the structure and composition of populations, but they often require specific experimental conditions due to small amount or low quality of DNA. In this study, we designed sets of primers that would allow efficient molecular identification of species and sex of individuals of bean geese and white-fronted geese. Species-specific primers (WFG-F/BG-F and G-R), developed from ND2 region of mitochondrial DNA, produced PCR products with different sizes which allow easy species identification without further sequencing. Based on published CHD 1 sequences, we designed internal primers (Gsex-F and Gsex-R) for sex determination that can be used in nested PCR after applying P2/P8 primers, and our methods clearly showed high success rate of molecular sexing from non-invasive samples. These molecular tools open the possibilities for genetic studies using non-invasive samples collected from a mixed-species aggregation containing bean geese and white-fronted geese.

[Phylogeography of the mallard Anas platyrhynchos from Eurasia inferred from sequencing of the mtDNA control region].

Phylogenetic relationships, demographic history, and geographic distribution of the mtDNA haplotypes of the mallard Anas platyrhynchos were examined in three populations, Indian, Northern European, and Far Eastern. Two divergent haplotype groups, A and B, were found in the Far Eastern population, while haplotypes identified in Northern European and Far Eastern populations were exclusively of the A group. The presence of B group haplotypes in the Far Eastern population can be explained either in terms of hybridization of the mallard with spot-billed duck Anas zonorhyncha at the south of the Russian Far East, or by the mtDNA paraphyly in mallards. In general, mallards from Eurasia were characterized by low genetic population differentiation along with slightly expressed phylogeographic structure. The most differentiated was the population from India (Phi(st) = 0.076-0.077), while the difference between Northern European and Far Eastern populations was extremely low (Phi(st) = 0.0029). Differentiation of Anas platyrhynchos Indian population was determined by the fact that a part of the population, inhabiting southern and eastern coasts of the Hindustan Peninsula, was resident.

Synteny conservation of chicken macrochromosomes 1-10 in different avian lineages revealed by cross-species chromosome painting.

Cross-species chromosome painting can directly visualize syntenies between diverged karyotypes and, thus, increase our knowledge on avian genome evolution. DNA libraries of chicken (Gallus gallus, GGA) macrochromosomes 1 to 10 were hybridized to metaphase spreads of 9 different species from 3 different orders (Anseriformes, Gruiformes and Passeriformes). Depending on the analyzed species, GGA1-10 delineated 11 to 13 syntenic chromosome regions, indicating a high degree of synteny conservation. No exchange between the GGA macrochromosome complement and microchromosomes of the analyzed species was observed. GGA1 and GGA4 were distributed on 2 or 3 chromosomes each in some of the analyzed species, indicating rare evolutionary rearrangements between macrochromosomes. In all 6 analyzed species of Passeriformes, GGA1 was diverged on 2 macrochromosomes, representing a synapomorphic marker for this order. GGA4 was split on 2 chromosomes in most karyotypes, but syntenic to a single chromosome in blackcap (Passeriformes). GGA5/10 and also GGA8/9 associations on chromosomes were found to be important cytogenetic features of the Eurasian nuthatch (Passeriformes) karyotype. Fusion of GGA4 and GGA5 segments and of entire GGA6 and GGA7, respectively, was seen in the 2 analyzed species of Gruiformes. Consistent with the literature, our inter-species chromosome painting demonstrates remarkable conservation of macrochromosomal synteny over 100 million years of avian evolution. The low rate of rearrangements between macrochromosomes and the absence of detectable macrochromosome-microchromosome exchanges suggests a predominant role for rearrangements within the gene-dense microchromosome complement in karyotypic diversification.

Discord reigns among nuclear, mitochondrial and phenotypic estimates of divergence in nine lineages of trans-Beringian birds.

Proposals for genetic thresholds for species delimitation assume that simple genetic data sets (e.g. mitochondrial sequence data) are correlated with speciation; i.e. such data sets accurately reflect organismal lineage divergence. We used taxonomically stratified phenotypic levels of differentiation (populations, subspecies and species) among nine avian lineages using paired, trans-Beringian samples from three lineages each in three orders (Anseriformes, Charadriiformes, and Passeriformes) to test this assumption. Using mitochondrial DNA sequence data and nuclear genomic data (amplified fragment length polymorphisms), we found a lack of concordance between these two genomes in their respective estimates of divergence and little or no relationship between phenotype (taxonomic relatedness) and genetic differentiation between taxon pairs. There are several possible reasons for the discord observed (e.g. selection on one of the genomes or perhaps lineage sorting), but the implications are that genetic estimates of lineage divergence may not be correlated with estimates from other parts of the genome, are not well correlated with the speciation process and are thus not reliable indicators of species limits.

Health monitoring in birds using bio-loggers and whole blood transcriptomics.

Monitoring and early detection of emerging infectious diseases in wild animals is of crucial global importance, yet reliable ways to measure immune status and responses are lacking for animals in the wild. Here we assess the usefulness of bio-loggers for detecting disease outbreaks in free-living birds and confirm detailed responses using leukocyte composition and large-scale transcriptomics. We simulated natural infections by viral and bacterial pathogens in captive mallards (Anas platyrhynchos), an important natural vector for avian influenza virus. We show that body temperature, heart rate and leukocyte composition change reliably during an acute phase immune response. Using genome-wide gene expression profiling of whole blood across time points we confirm that immunostimulants activate pathogen-specific gene regulatory networks. By reporting immune response related changes in physiological and behavioural traits that can be studied in free-ranging populations, we provide baseline information with importance to the global monitoring of zoonotic diseases.

The complete mitochondrial genomes of the whistling duck (Dendrocygna javanica) and black swan (Cygnus atratus): dating evolutionary divergence in Galloanserae.

Galloanserae is an ancient and diverse avian group, for which comprehensive molecular evidence relevant to phylogenetic analysis in the context of molecular chronology is lacking. In this study, we present two additional mitochondrial genome sequences of Galloanserae (the whistling duck, Dendrocygna javanica, and the black swan, Cygnus atratus) to broaden the scope of molecular phylogenetic reconstruction. The lengths of the whistling duck's and black swan's mitochondrial genomes are 16,753 and 16,748 bases, respectively. Phylogenetic analyses suggest that Dendrocygna is more likely to be in a basal position of the branch consisting of Anatinae and Anserinae, an affiliation that does not conform to its traditional classification. Bayesian approaches were employed to provide a rough timescale for Galloanserae evolution. In general, a narrow range of 95% confidence intervals gave younger estimates than those based on limited genes and estimated that at least two lineages originated before the Coniacian epoch around 90 MYA, well before the Cretaceous-Tertiary boundary. In addition, these results, which were compatible with estimates from fossil evidence, also imply that the origin of numerous genera in Anseriformes took place in the late Oligocene to early Miocene. Taken together, the results presented here provide a working framework for future research on Galloanserae evolution, and they underline the utility of whole mitochondrial genome sequences for the resolution of deep divergence.

Effects of dietary cadmium contamination on bird Anas platyrhynchos--comparison with species Cairina moschata.

This study aimed to assess the effect of two dietary cadmium (Cd) levels (C1: 1 mgkg(-1); C10: 10 mgkg(-1)) on bird Anas platyrhynchos exposed for 10, 20 and 40 days (5 animals per experimental condition). Ducks were able to accumulate high amounts of Cd, especially in kidneys (after 40 days: C1 8.1 ± 1 mgkg(-1), C10 37.7 ± 4.3 mgkg(-1)). After 40 days, the lowest Cd level triggered oxidative stress and stimulated mitochondrial metabolism. At the same time, highest amounts of Cd (C10 group) only triggered repression of genes encoding for catalase and acetyl-CoA carboxylase, with repression factors of 1/50 and 1/5, respectively. High dose exposures were then associated with the repression of genes encoding for antioxidant, whereas low dose exposure triggered their induction. In contrast, the onset of MT gene expression appeared quickly for the C10 group even if a time delay was observed between gene expression and protein accumulation. Through the comparison of A. platyrhynchos and Cairina moschata, the response to Cd toxicity appeared species-dependent. Discrepancies between species could be explained by differential utilization of MT. This pathway of detoxification seemed sufficient to counter Cd toxicity.

[Genetic structure of the Far Eastern population of Eurasian wigeon Anas penelope inferred from sequencing of the mitochondrial DNA control region].

Sequence variation of the 5' end of the mitochondrial DNA control region (600 bp) was examined in the population samples of Eurasian wigeon Anas penelope from Anadyr' and Primorye. A total of 11 different mtDNA haplotypes were identified, with one of these belonging to American wigeon Anas americana. The presence of the mtDNA haplotype from the species closely relative to A. penelope in the Anadyr' sample can be considered as the genetic evidence in favor of interspecific hybridization. This suggestion is in the good agreement with ornithological data. Genetic differentiation of the Primorye and Anadyr' populations was low (phi(ST) = 0.096). The phylogeographic structure was not pronounced.