Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Concurrent invasions of European starlings in Australia and North America reveal population-specific differentiation in shared genomic regions.

A species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole-genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population-specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process.

Free-induction-decay magnetic field imaging with a microfabricated Cs vapor cell.

Magnetic field imaging is a valuable resource for signal source localization and characterization. This work reports an optically pumped magnetometer (OPM) based on the free-induction-decay (FID) protocol, that implements microfabricated cesium (Cs) vapor cell technology to visualize the magnetic field distributions resulting from various magnetic sources placed close to the cell. The slow diffusion of Cs atoms in the presence of a nitrogen (N2) buffer gas enables spatially independent measurements to be made within the same vapor cell by translating a 175 µm diameter probe beam over the sensing area. For example, the OPM was used to record temporal and spatial information to reconstruct magnetic field distributions in one and two dimensions. The optimal magnetometer sensitivity was estimated to be 0.43 pT/H z within a Nyquist limited bandwidth of 500 Hz. Furthermore, the sensor's dynamic range exceeds the Earth's field of approximately 50 µT, which provides a framework for magnetic field imaging in unshielded environments.

Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus Symbiodinium.

BACKGROUND: Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1-5 Gbp) and idiosyncratic genome features. RESULTS: Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species. CONCLUSIONS: Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Genomes of the dinoflagellate Polarella glacialis encode tandemly repeated single-exon genes with adaptive functions.

BACKGROUND: Dinoflagellates are taxonomically diverse and ecologically important phytoplankton that are ubiquitously present in marine and freshwater environments. Mostly photosynthetic, dinoflagellates provide the basis of aquatic primary production; most taxa are free-living, while some can form symbiotic and parasitic associations with other organisms. However, knowledge of the molecular mechanisms that underpin the adaptation of these organisms to diverse ecological niches is limited by the scarce availability of genomic data, partly due to their large genome sizes estimated up to 250 Gbp. Currently available dinoflagellate genome data are restricted to Symbiodiniaceae (particularly symbionts of reef-building corals) and parasitic lineages, from taxa that have smaller genome size ranges, while genomic information from more diverse free-living species is still lacking. RESULTS: Here, we present two draft diploid genome assemblies of the free-living dinoflagellate Polarella glacialis, isolated from the Arctic and Antarctica. We found that about 68% of the genomes are composed of repetitive sequence, with long terminal repeats likely contributing to intra-species structural divergence and distinct genome sizes (3.0 and 2.7 Gbp). For each genome, guided using full-length transcriptome data, we predicted > 50,000 high-quality protein-coding genes, of which ~40% are in unidirectional gene clusters and ~25% comprise single exons. Multi-genome comparison unveiled genes specific to P. glacialis and a common, putatively bacterial origin of ice-binding domains in cold-adapted dinoflagellates. CONCLUSIONS: Our results elucidate how selection acts within the context of a complex genome structure to facilitate local adaptation. Because most dinoflagellate genes are constitutively expressed, Polarella glacialis has enhanced transcriptional responses via unidirectional, tandem duplication of single-exon genes that encode functions critical to survival in cold, low-light polar environments. These genomes provide a foundational reference for future research on dinoflagellate evolution.

The quail genome: insights into social behaviour, seasonal biology and infectious disease response.

BACKGROUND: The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research. RESULTS: We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1. CONCLUSIONS: We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

Illuminating the dark side of the human transcriptome with long read transcript sequencing.

BACKGROUND: The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise. RESULTS: We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models. CONCLUSIONS: Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

Toll-Like Receptor Evolution in Birds: Gene Duplication, Pseudogenization, and Diversifying Selection.

Toll-like receptors (TLRs) are key sensor molecules in vertebrates triggering initial phases of immune responses to pathogens. The avian TLR family typically consists of ten receptors, each adapted to distinct ligands. To understand the complex evolutionary history of each avian TLR, we analyzed all members of the TLR family in the whole genome assemblies and target sequence data of 63 bird species covering all major avian clades. Our results indicate that gene duplication events most probably occurred in TLR1 before synapsids diversified from sauropsids. Unlike mammals, ssRNA-recognizing TLR7 has duplicated independently in several avian taxa, while flagellin-sensing TLR5 has pseudogenized multiple times in bird phylogeny. Our analysis revealed stronger positive, diversifying selection acting in TLR5 and the three-domain TLRs (TLR10 [TLR1A], TLR1 [TLR1B], TLR2A, TLR2B, TLR4) that face the extracellular space and bind complex ligands than in single-domain TLR15 and endosomal TLRs (TLR3, TLR7, TLR21). In total, 84 out of 306 positively selected sites were predicted to harbor substitutions dramatically changing the amino acid physicochemical properties. Furthermore, 105 positively selected sites were located in the known functionally relevant TLR regions. We found evidence for convergent evolution acting between birds and mammals at 54 of these sites. Our comparative study provides a comprehensive insight into the evolution of avian TLR genetic variability. Besides describing the history of avian TLR gene gain and gene loss, we also identified candidate positions in the receptors that have been likely shaped by direct molecular host-pathogen coevolutionary interactions and most probably play key functional roles in birds.

APOLLO-2: A Randomized, Placebo and Active-Controlled Phase III Study Investigating Oliceridine (TRV130), a G Protein-Biased Ligand at the µ-Opioid Receptor, for Management of Moderate to Severe Acute Pain Following Abdominoplasty.

OBJECTIVES: The clinical utility of conventional IV opioids is limited by the occurrence of opioid-related adverse events. Oliceridine is a novel G protein-biased µ-opioid receptor agonist designed to provide analgesia with an improved safety and tolerability profile. This phase III, double-blind, randomized trial (APOLLO-2 [NCT02820324]) evaluated the efficacy and safety of oliceridine for acute pain following abdominoplasty. METHODS: Patients received a loading dose of either placebo, oliceridine (1.5 mg), or morphine (4 mg), followed by demand doses via patient-controlled analgesia (0.1, 0.35, or 0.5 mg oliceridine; 1 mg morphine; or placebo) with a 6-minute lockout interval. The primary endpoint was the proportion of treatment responders over 24 hours for oliceridine regimens compared to placebo. Secondary outcomes included a predefined composite measure of respiratory safety burden (RSB, representing the cumulative duration of respiratory safety events) and the proportion of treatment responders vs. morphine. RESULTS: A total of 401 patients were treated with study medication. Effective analgesia was observed for all oliceridine regimens, with responder rates of 61.0%, 76.3%, and 70.0% for the 0.1-, 0.35-, and 0.5-mg regimens, respectively, compared with 45.7% for placebo (all P < 0.05) and 78.3% for morphine. Oliceridine 0.35- and 0.5-mg demand dose regimens were equi-analgesic to morphine using a noninferiority analysis. RSB showed a dose-dependent increase across oliceridine regimens (mean hours [standard deviation], 0.1 mg: 0.43 [1.56]; 0.35 mg: 1.48 [3.83]; 0.5 mg: 1.59 [4.26]; all comparisons not significant at P > 0.05 vs. placebo: 0.60 [2.82]). The RSB measure for morphine was 1.72 (3.86) (P < 0.05 vs. placebo). Gastrointestinal adverse events increased in a dose-dependent manner across oliceridine demand dose regimens (0.1 mg: 49.4%; 0.35 mg: 65.8%; 0.5 mg: 78.8%; vs. placebo: 47.0%; and morphine: 79.3%). In comparison to morphine, the proportion of patients experiencing nausea or vomiting was lower with the 2 equi-analgesic dose regimens of 0.35 and 0.5 mg oliceridine. CONCLUSIONS: Oliceridine is a safe and effective IV analgesic for the relief of moderate to severe acute postoperative pain in patients undergoing abdominoplasty. Since the low-dose regimen of 0.1 mg oliceridine was superior to placebo but not as effective as the morphine regimen, safety comparisons to morphine are relevant only to the 2 equi-analgesic dose groups of 0.35 and 0.5 mg, which showed a favorable safety and tolerability profile regarding respiratory and gastrointestinal adverse effects compared to morphine. These findings support that oliceridine may provide a new treatment option for patients with moderate to severe acute pain where an IV opioid is warranted.

Chicken anaemia virus evades host immune responses in transformed lymphocytes.

Chicken anaemia virus (CAV) is a lymphotropic virus that causes anaemia and immunosuppression in chickens. Previously, we proposed that CAV evades host antiviral responses in vivo by disrupting T-cell signalling, but the precise cellular targets and modes of action remain elusive. In this study, we examined gene expression in Marek's disease virus-transformed chicken T-cell line MSB-1 after infection with CAV using both a custom 5K immune-focused microarray and quantitative real-time PCR at 24, 48 and 72 h post-infection. The data demonstrate an intricate equilibrium between CAV and the host gene expression, displaying subtle but significant modulation of transcripts involved in the T-cell, inflammation and NF-κB signalling cascades. CAV efficiently blocked the induction of type-I interferons and interferon-stimulated genes at 72 h. The cell expression pattern implies that CAV subverts host antiviral responses and that the transformed environment of MSB-1 cells offers an opportunistic advantage for virus growth.

Normalized long read RNA sequencing in chicken reveals transcriptome complexity similar to human.

BACKGROUND: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences. RESULTS: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences. CONCLUSIONS: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

Visualisation of chicken macrophages using transgenic reporter genes: insights into the development of the avian macrophage lineage.

We have generated the first transgenic chickens in which reporter genes are expressed in a specific immune cell lineage, based upon control elements of the colony stimulating factor 1 receptor (CSF1R) locus. The Fms intronic regulatory element (FIRE) within CSF1R is shown to be highly conserved in amniotes and absolutely required for myeloid-restricted expression of fluorescent reporter genes. As in mammals, CSF1R-reporter genes were specifically expressed at high levels in cells of the macrophage lineage and at a much lower level in granulocytes. The cell lineage specificity of reporter gene expression was confirmed by demonstration of coincident expression with the endogenous CSF1R protein. In transgenic birds, expression of the reporter gene provided a defined marker for macrophage-lineage cells, identifying the earliest stages in the yolk sac, throughout embryonic development and in all adult tissues. The reporter genes permit detailed and dynamic visualisation of embryonic chicken macrophages. Chicken embryonic macrophages are not recruited to incisional wounds, but are able to recognise and phagocytose microbial antigens.

Mapping of leptin and its syntenic genes to chicken chromosome 1p.

BACKGROUND: Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly. RESULTS: A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = - 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes. CONCLUSION: Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Correction to: Mapping of leptin and its syntenic genes to chicken chromosome 1p.

CORRECTION: After the publication of this work [1] an error was noticed in one of the author surnames. The author name Leif Anderson should be spelt as Leif Andersson.

Prevalence of risk factors for noncommunicable diseases in an indigenous community in Santiago Atitlán, Guatemala.

OBJECTIVE: To describe the prevalence of noncommunicable disease (NCD) risk factors and assess knowledge of those risk factors in the indigenous community of Santiago Atitlán in Guatemala, a lower-middle income country. METHODS: A population-based, cross-sectional study was conducted using a modified version of the World Health Organization's STEPS protocol. Adults aged 20-65 years were surveyed regarding demographics and NCD risk factors, and the survey was followed by anthropometric and biochemical measurements. RESULTS: Out of 501 screened individuals, 350 respondents were enrolled. The mean age was 36.7 years, and 72.3% were women. Over 90% reported earning less than US$ 65 per month. Almost 80% were stunted. Among women, 37.3% were obese and over three-quarters had central obesity. Over three-quarters of the entire group had dyslipidemia and 18.3% had hypertension, but only 3.0% had diabetes. Overall, 36.0% of participants met criteria for metabolic syndrome. There was no significant association between participants' education and NCD risk factors except for an inverse association with obesity by percent body fat. CONCLUSIONS: Santiago Atitlán is a rural, indigenous Guatemalan community with high rates of poverty and stunting coexisting alongside high rates of obesity, particularly among women. Additionally, high rates of hypertension and dyslipidemia were found, but a low rate of diabetes mellitus. Knowledge of NCDs and their risk factors was low, suggesting that educational interventions may be a high-yield, low-cost approach to combating NCDs in this community.

A new look at the LTR retrotransposon content of the chicken genome.

BACKGROUND: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact. RESULTS: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor. CONCLUSIONS: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size.

Analysis of the early immune response to infection by infectious bursal disease virus in chickens differing in their resistance to the disease.

UNLABELLED: Chicken whole-genome gene expression arrays were used to analyze the host response to infection by infectious bursal disease virus (IBDV). Spleen and bursal tissue were examined from control and infected birds at 2, 3, and 4 days postinfection from two lines that differ in their resistance to IBDV infection. The host response was evaluated over this period, and differences between susceptible and resistant chicken lines were examined. Antiviral genes, including IFNA, IFNG, MX1, IFITM1, IFITM3, and IFITM5, were upregulated in response to infection. Evaluation of this gene expression data allowed us to predict several genes as candidates for involvement in resistance to IBDV. IMPORTANCE: Infectious bursal disease (IBD) is of economic importance to the poultry industry and thus is also important for food security. Vaccines are available, but field strains of the virus are of increasing virulence. There is thus an urgent need to explore new control solutions, one of which would be to breed birds with greater resistance to IBD. This goal is perhaps uniquely achievable with poultry, of all farm animal species, since the genetics of 85% of the 60 billion chickens produced worldwide each year is under the control of essentially two breeding companies. In a comprehensive study, we attempt here to identify global transcriptomic differences in the target organ of the virus between chicken lines that differ in resistance and to predict candidate resistance genes.

Evolution of the avian ß-defensin and cathelicidin genes.

BACKGROUND: ß-defensins and cathelicidins are two families of cationic antimicrobial peptides (AMPs) with a broad range of antimicrobial activities that are key components of the innate immune system. Due to their important roles in host defense against rapidly evolving pathogens, the two gene families provide an ideal system for studying adaptive gene evolution. In this study we performed phylogenetic and selection analyses on ß-defensins and cathelicidins from 53 avian species representing 32 orders to examine the evolutionary dynamics of these peptides in birds. RESULTS AND CONCLUSIONS: Avian ß-defensins are found in a gene cluster consisting of 13 subfamiles. Nine of these are conserved as one to one orthologs in all birds, while the others (AvBD1, AvBD3, AvBD7 and AvBD14) are more subject to gene duplication or pseudogenisation events in specific avian lineages. Avian cathelicidins are found in a gene cluster consisting of three subfamilies with species-specific duplications and gene loss. Evidence suggested that the propiece and mature peptide domains of avian cathelicidins are possibly co-evolving in such a way that the cationicity of the mature peptide is partially neutralised by the negative charge of the propiece prior to peptide secretion (further evidence obtained by repeating the analyses on primate cathelicidins). Negative selection (overall mean dN < dS) was detected in most of the gene domains examined, conserving certain amino acid residues that may be functionally crucial for the avian ß-defensins and cathelicidins, while episodic positive selection was also involved in driving the diversification of specific codon sites of certain AMPs in avian evolutionary history. These findings have greatly improved our understanding of the molecular evolution of avian AMPs and will be useful to understand their role in the avian innate immune response. Additionally, the large dataset of ß-defensin and cathelicidin peptides may also provide a valuable resource for translational research and development of novel antimicrobial agents in the future.

Detection and characterization of small insertion and deletion genetic variants in modern layer chicken genomes.

BACKGROUND: Small insertions and deletions (InDels) constitute the second most abundant class of genetic variants and have been found to be associated with many traits and diseases. The present study reports on the detection and characterisation of about 883 K high quality InDels from the whole-genome analysis of several modern layer chicken lines from diverse breeds. RESULTS: To reduce the error rates seen in InDel detection, this study used the consensus set from two InDel-calling packages: SAMtools and Dindel, as well as stringent post-filtering criteria. By analysing sequence data from 163 chickens from 11 commercial and 5 experimental layer lines, this study detected about 883 K high quality consensus InDels with 93% validation rate and an average density of 0.78 InDels/kb over the genome. Certain chromosomes, viz, GGAZ, 16, 22 and 25 showed very low densities of InDels whereas the highest rate was observed on GGA6. In spite of the higher recombination rates on microchromosomes, the InDel density on these chromosomes was generally lower relative to macrochromosomes possibly due to their higher gene density. About 43-87% of the InDels were found to be fixed within each line. The majority of detected InDels (86%) were 1-5 bases and about 63% were non-repetitive in nature while the rest were tandem repeats of various motif types. Functional annotation identified 613 frameshift, 465 non-frameshift and 10 stop-gain/loss InDels. Apart from the frameshift and stopgain/loss InDels that are expected to affect the translation of protein sequences and their biological activity, 33% of the non-frameshift were predicted as evolutionary intolerant with potential impact on protein functions. Moreover, about 2.5% of the InDels coincided with the most-conserved elements previously mapped on the chicken genome and are likely to define functional elements. InDels potentially affecting protein function were found to be enriched for certain gene-classes e.g. those associated with cell proliferation, chromosome and Golgi organization, spermatogenesis, and muscle contraction. CONCLUSIONS: The large catalogue of InDels presented in this study along with their associated information such as functional annotation, estimated allele frequency, etc. are expected to serve as a rich resource for application in future research and breeding in the chicken.