Rapid sperm length divergence in a polygynandrous passerine: a mechanism of cryptic speciation?

When populations become geographically isolated, they begin to diverge in various traits and at variable rates. The dynamics of such trait divergences are relevant for understanding evolutionary processes such as local adaptation and speciation. Here we examine divergences in sperm and body structures in a polygynandrous songbird, the alpine accentor (Prunella collaris) between two allopatric high-altitude populations, in Morocco and Spain. The populations diverged around 82,000 years ago, as estimated with a coalescence-based phylogenetic analysis of genome-wide single-nucleotide polymorphisms. We found that birds in the two areas had nonoverlapping sperm lengths, which suggests adaptation to divergent female reproductive tract environments. Sperm length also showed an exceptionally low coefficient of among-male variation, a signal of strong stabilizing selection imposed by sperm competition. The evolutionary rate of sperm length was almost twice the rates for the most divergent morphological traits and more than three times higher than expected from literature data over a similar generational timescale. This rapid evolution of a key reproductive trait has implications for reproductive isolation and ultimately for speciation. Strong selection for different sperm length optima in allopatry predicts conspecific sperm precedence and disruptive selection in sympatry, hence a possible postcopulatory prezygotic barrier to gene flow.

Guanylate Cyclase C Activation Shapes the Intestinal Microbiota in Patients with Familial Diarrhea and Increased Susceptibility for Crohn's Disease.

BACKGROUND: With 25% prevalence of Crohn's disease, Familial GUCY2C diarrhea syndrome (FGDS) is a monogenic disorder potentially suited to study initiating factors in inflammatory bowel disease (IBD). We aimed to characterize the impact of an activating GUCY2C mutation on the gut microbiota in patients with FGDS controlling for Crohn's disease status and to determine whether changes share features with those observed in unrelated patients with IBD. METHODS: Bacterial DNA from fecal samples collected from patients with FGDS (N = 20), healthy relatives (N = 11), unrelated healthy individuals (N = 263), and IBD controls (N = 46) was subjected to sequencing of the V3-V4 region of the 16S rRNA gene to determine gut microbiota composition. Food frequency questionnaires were obtained from patients with FGDS and their relatives. RESULTS: Compared with healthy controls, FGDS displayed prominent changes in many microbial lineages including increase in Enterobacteriaceae, loss of Bifidobacterium and Faecalibacterium prausnitzii but an unchanged intraindividual (alpha) diversity. The depletion of F. prausnitzii is in line with what is typically observed in Crohn's disease. There was no significant difference in the dietary profile between the patients and related controls. The gut microbiota in related and unrelated healthy controls was also similar, suggesting that diet and familial factors do not explain the gut microbiota alterations in FGDS. CONCLUSIONS: The findings support that the activating mutation in GUCY2C creates an intestinal environment with a major influence on the microbiota, which could contribute to the increased susceptibility to IBD in patients with FGDS.

Complete mitochondrial genomes of eleven extinct or possibly extinct bird species.

Natural history museum collections represent a vast source of ancient and historical DNA samples from extinct taxa that can be utilized by high-throughput sequencing tools to reveal novel genetic and phylogenetic information about them. Here, we report on the successful sequencing of complete mitochondrial genome sequences (mitogenomes) from eleven extinct bird species, using de novo assembly of short sequences derived from toepad samples of degraded DNA from museum specimens. For two species (the Passenger Pigeon Ectopistes migratorius and the South Island Piopio Turnagra capensis), whole mitogenomes were already available from recent studies, whereas for five others (the Great Auk Pinguinis impennis, the Imperial Woodpecker Campehilus imperialis, the Huia Heteralocha acutirostris, the Kauai Oo Moho braccathus and the South Island Kokako Callaeas cinereus), there were partial mitochondrial sequences available for comparison. For all seven species, we found sequence similarities of >98%. For the remaining four species (the Kamao Myadestes myadestinus, the Paradise Parrot Psephotellus pulcherrimus, the Ou Psittirostra psittacea and the Lesser Akialoa Akialoa obscura), there was no sequence information available for comparison, so we conducted blast searches and phylogenetic analyses to determine their phylogenetic positions and identify their closest extant relatives. These mitogenomes will be valuable for future analyses of avian phylogenetics and illustrate the importance of museum collections as repositories for genomics resources.

Development of a high-resolution NGS-based HLA-typing and analysis pipeline.

The human leukocyte antigen (HLA) complex contains the most polymorphic genes in the human genome. The classical HLA class I and II genes define the specificity of adaptive immune responses. Genetic variation at the HLA genes is associated with susceptibility to autoimmune and infectious diseases and plays a major role in transplantation medicine and immunology. Currently, the HLA genes are characterized using Sanger- or next-generation sequencing (NGS) of a limited amplicon repertoire or labeled oligonucleotides for allele-specific sequences. High-quality NGS-based methods are in proprietary use and not publicly available. Here, we introduce the first highly automated open-kit/open-source HLA-typing method for NGS. The method employs in-solution targeted capturing of the classical class I (HLA-A, HLA-B, HLA-C) and class II HLA genes (HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). The calling algorithm allows for highly confident allele-calling to three-field resolution (cDNA nucleotide variants). The method was validated on 357 commercially available DNA samples with known HLA alleles obtained by classical typing. Our results showed on average an accurate allele call rate of 0.99 in a fully automated manner, identifying also errors in the reference data. Finally, our method provides the flexibility to add further enrichment target regions.

Small duct primary sclerosing cholangitis without inflammatory bowel disease is genetically different from large duct disease.

BACKGROUND & AIMS: Small duct primary sclerosing cholangitis (PSC) is phenotypically a mild version of large duct PSC, but it is unknown whether these phenotypes share aetiology. We aimed to characterize their relationship by investigating genetic associations in the human leucocyte antigen (HLA) complex, which represent the strongest genetic risk factors in large duct PSC. METHODS: Four classical HLA loci (HLA-A, HLA-B, HLA-C and HLA-DRB1) were genotyped in 87 small duct PSC patients, 485 large duct PSC patients and 1117 controls across three geographical regions. RESULTS: HLA-DRB1*13:01 (OR = 2.0, 95% CI 1.2-3.4, P = 0.01) and HLA-B*08 (OR = 1.6, 95% CI 1.1-2.4, P = 0.02) were significantly associated with small duct PSC compared with healthy controls. Based on the observed frequency of HLA-B*08 in small duct PSC, the strongest risk factor in large duct PSC, an estimated 32% (95% CI 4-65%) of this population can be hypothesized to represent early stages or mild variants of large duct PSC. This subgroup may be constituted by small duct PSC patients with inflammatory bowel disease (IBD), which greatly resembled large duct PSC in its HLA association. In contrast, small duct PSC without IBD was only associated with HLA-DRB1*13:01(P = 0.03) and was otherwise distinctly dissimilar from large duct PSC. CONCLUSIONS: Small duct PSC with IBD resembles large duct PSC in its HLA association and may represent early stages or mild variants of large duct disease. Different HLA associations in small duct PSC without IBD could indicate that this subgroup is a different entity. HLA-DRB1*13:01 may represent a specific risk factor for inflammatory bile duct disease.

Female promiscuity is positively associated with neutral and selected genetic diversity in passerine birds.

Passerine birds show large interspecific variation in extrapair paternity rates. There is accumulating evidence that such promiscuous behavior is driven by indirect, genetic benefits to females. Sexual selection theory distinguishes between two types of genetic benefits, additive and nonadditive effects, mediated by preferences for good and compatible genes, respectively. Good genes preferences should imply directional selection and mating skew among males, and thus reduced genetic diversity in the population. In contrast, compatible genes preferences should give balancing selection that retains genetic diversity. Here, we test how well these predictions fit with patterns of variation in genetic diversity and promiscuity levels among passerine birds. We found that more promiscuous species had higher nucleotide diversity at autosomal introns, but not at Z-chromosome introns. We also found that major histocompatibility complex (MHC) class IIB alleles had higher sequence diversity, and therefore should recognize a broader spectrum of pathogens, in more promiscuous species. Our results suggest that female promiscuity targets a multitude of autosomal genes for their nonadditive, compatibility benefits. Also, as immunity genes seem to be of particular importance, we hypothesize that interspecific variation in female promiscuity among passerine birds has arisen in response to the strength of pathogen-mediated selection.

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae).

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca.

BACKGROUND: Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts. FINDINGS: Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (ß chain). Taking the pied flycatcher Ficedula hypoleuca as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species. CONCLUSIONS: Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.

Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes.

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide-binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.

Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10.

BACKGROUND: Microsatellites are frequently used genetic markers in a wide range of applications, primarily due to their high length polymorphism levels that can easily be genotyped by fragment length analysis. However, the mode of microsatellite evolution is yet not fully understood, and the role of interrupting motifs for the stability of microsatellites remains to be explored in more detail. Here we present a sequence analysis of mutation events and a description of the structure of repeated regions in the hypervariable, pentanucleotide microsatellite locus HrU10 in barn swallows (Hirundo rustica) and tree swallows (Tachycineta bicolor). RESULTS: In a large-scale parentage analysis in barn swallows and tree swallows, broods were screened for mutations at the HrU10 locus. In 41 cases in the barn swallows and 15 cases in the tree swallows, mutations corresponding to the loss or gain of one or two repeat units were detected. The parent and mutant offspring alleles were sequenced for 33 of these instances (26 in barn swallows and 7 in tree swallows). Replication slippage was considered the most likely mutational process. We tested the hypothesis that HrU10, a microsatellite with a wide allele size range, has an increased probability of introductions of interruptive motifs (IMs) with increasing length of the repeated region. Indeed, the number and length of the IMs was strongly positively correlated with the total length of the microsatellite. However, there was no significant correlation with the length of the longest stretch of perfectly repeated units, indicating a threshold level for the maximum length of perfectly repeated pentanucleotide motifs in stable HrU10 alleles. The combination of sequence and pedigree data revealed that 15 barn swallow mutations (58%) produced alleles that were size homoplasic to other alleles in the data set. CONCLUSION: Our results give further insights into the mode of microsatellite evolution, and support the assumption of increased slippage rate with increased microsatellite length and a stabilizing effect of interrupting motifs for microsatellite regions consisting of perfect repeats. In addition, the observed extent of size homoplasy may impose a general caution against using hypervariable microsatellites in genetic diversity measures when alleles are identified by fragment length analysis only.