Sequencing and alignment of mitochondrial genomes of Tibetan chicken and two lowland chicken breeds

Tibetan chicken lives in high-altitude area and has adapted well to hypoxia genetically. Shouguang chicken and Silky chicken are both lowland chicken breeds. In the present study, the complete mitochondrial genome sequences of the three chicken breeds were all sequenced. The results showed that the mitochondrial DNAs (mtDNAs) of Shouguang chicken and Silky chicken consist of 16784 bp and 16785 bp respectively, and Tibetan chicken mitochondrial genome varies from 16784 bp to 16786 bp. After sequence analysis, 120 mutations, including 4 single nucleotide polymorphisms (SNPs) in tRNA genes, 9 SNPs and 1 insertion in rRNA genes, 38 SNPs and 1 deletion in D-LOOP, 66 SNPs in protein-coding genes, were found. This work will provide clues for the future study on the association between mitochondrial genes and the adaptation to hypoxia.     

Sequencing and alignment of mitochondrial genomes of Tibetan chicken and two lowland chicken breeds

Tibetan chicken lives in high-altitude area and has adapted well to hypoxia genetically. Shouguang chicken and Silky chicken are both lowland chicken breeds. In the present study, the complete mito-chondrial genome sequences of the three chicken breeds were all sequenced. The results showed that the mitochondrial DNAs (mtDNAs) of Shouguang chicken and Silky chicken consist of 16784 bp and 16785 bp respectively, and Tibetan chicken mitochondrial genome varies from 16784 bp to 16786 bp. After sequence analysis, 120 mutations, including 4 single nucleotide polymorphisms (SNPs) in tRNA genes, 9 SNPs and 1 insertion in rRNA genes, 38 SNPs and 1 deletion in D-LOOP, 66 SNPs in pro-tein-coding genes, were found. This work will provide clues for the future study on the association between mitochondrial genes and the adaptation to hypoxia.Tibetan chicken, lowland chicken, mitochondrial genome, hypoxia.     

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F _st values against genome location. The highest F _st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.     

A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland

Background

Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals’ phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Pu?awska, Z?otnicka White and Z?otnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations.

Results

A two largely complementary statistical methods were used for signatures detection, including: pairwise F_ST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo.

Conclusions

The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.

     

Genetic relationships between Japanese native and commercial breeds using 70 chicken autosomal SNP genotypes by the DigiTag2 assay

Recently, single nucleotide polymorphisms (SNPs) have been used to identify genes or genomic regions responsible for economic traits, including genetic diseases in domestic animals, and to examine genetic diversity of populations. In this study, we genotyped 70 chicken autosomal SNPs using DigiTag2 assay to understand the genetic structure of the Japanese native chicken breeds Satsumadori and Ingie, and the relationship of these breeds with other established breeds, Rhode Island Red (RIR), commercial broiler and layer. Five breeds, each consisting of approximately 20 chickens, were subjected to the assay, revealing the following: Average expected heterozygosities of broiler, Satsumadori, RIR, layer and Ingie were 0.265, 0.254, 0.244, 0.179 and 0.176, respectively. Phylogenetic analysis using the concatenated 70 autosomal SNP genotypes distinguished all chickens and formed clusters of chickens belonging to the respective breeds. In addition, the 2-D scatter plot of the first two principal components was consistent with the phylogenic tree. Taken together with the pairwise F(st) distances, broiler and RIR were closely positioned near each other, while Ingie was positioned far from the other breeds. Structure analysis revealed that the probable number of genetic clusters (K) was six and four with maximum likelihood and ΔK values, respectively. The clustering with maximum likelihood revealed that, in addition to the clustering of the other five breeds, the Satsumadori was subdivided into two genetic clusters. The clustering with ΔK value indicated that the broiler and Rhode Island Red were assigned to the same genetic cluster.     

Evolutionary relationships of Red Jungle Fowl and chicken breeds

Published results were reassessed and original data are provided regarding the origin and relatedness of four postulated chicken breed lineages, egg-type, game, meat-type and Bantam, to each other and to the basic ancestral species of jungle fowls, Gallus gallus. A system approach was employed concerning the planning of the experiments. One element of the system approach is the choice of the breeds to be compared with G. gallus. These breeds were supposed to represent major evolutionary branches of chickens. Four experiments on genetic relationships were conducted using different estimation criteria including morphological discrete characters, body measurements, biochemical markers, and the activity of serum esterase-1. The greatest similarity was found between G. gallus and the egg-type breeds of Mediterranean roots and/or true Bantams. This fact might testify that the indicated chicken groups occupied earlier stages in the evolution from the wild progenitor to the present biodiversity of chickens in the world.     

Population genetic analyses of seven Chinese indigenous chicken breeds in a context of global breeds

Jiangxi province in China is rich in indigenous chicken breeds, which have diverse phenotypes and represent a valuable genetic resource for further genetic improvement of modern breeds. Here, we conducted a series of analyses to reveal genetic diversity, phylogenetic relationships and population structure of seven chicken breeds in Jiangxi province in the context of nine non‐local chicken breeds, using 600K SNP data. We show that Jiangxi indigenous breeds have more abundant nucleotide diversity than do European local and commercial breeds. Among Jiangxi breeds, Dongxiang Blue‐eggshell (DX) and Chongren Partride (CR) display remarkably reduced genetic diversity, as the two breeds exhibit increased inbreeding coefficients, runs of homozygosity, extent of linkage disequilibrium and reduced expected heterozygosity. DX, CR and Taihe Silkie (TH) represent three ancestral lineages of the Jiangxi chicken and display genetic differentiation from the other four Jiangxi breeds, which show a signature of admixture with European commercial breeds. These findings provide insight for the establishment of an efficient conservation program for Jiangxi chicken breeds. Considering the current status of genetic diversity and ancestral representativeness, particular attention should be paid to DX, CR and TH chickens.     

Maternal genealogical patterns of chicken breeds sampled in Europe

The aim of this study was to investigate the maternal genealogical pattern of chicken breeds sampled in Europe. Sequence polymorphisms of 1256 chickens of the hypervariable region (D‐loop) of mitochondrial DNA (mtDNA) were used. Median‐joining networks were constructed to establish evolutionary relationships among mtDNA haplotypes of chickens, which included a wide range of breeds with different origin and history. Chicken breeds which have had their roots in Europe for more than 3000 years were categorized by their founding regions, encompassing Mediterranean type, East European type and Northwest European type. Breeds which were introduced to Europe from Asia since the mid‐19th century were classified as Asian type, and breeds based on crossbreeding between Asian breeds and European breeds were classified as Intermediate type. The last group, Game birds, included fighting birds from Asia. The classification of mtDNA haplotypes was based on Liu et al.'s (2006) nomenclature. Haplogroup E was the predominant clade among the European chicken breeds. The results showed, on average, the highest number of haplotypes, highest haplotype diversity, and highest nucleotide diversity for Asian type breeds, followed by Intermediate type chickens. East European and Northwest European breeds had lower haplotype and nucleotide diversity compared to Mediterranean, Intermediate, Game and Asian type breeds. Results of our study support earlier findings that chicken breeds sampled in Europe have their roots in the Indian subcontinent and East Asia. This is consistent with historical and archaeological evidence of chicken migration routes to Europe.     

Signatures of contemporary selection in the Israeli Holstein dairy cattle

Strong selection in the Israeli Holstein dairy cattle population over the last three decades should have left clear signatures of selection. Two experimental approaches were applied to detect evidence of contemporary selection based on the 54K BeadChip genotypes of ~1000 Israeli Holstein bulls: (i) the long-range haplotype test, which searches for structural evidence resulting from selective sweep, and (ii) direct analysis of the changes in haplotypes frequencies over time combined with linkage disequilibrium blocks haplotype-based association analysis. Ten traits were analyzed: the PD07 Israeli selection index, milk, milk fat, % fat, milk protein, % protein, somatic cell score, female fertility, milk production persistency and herd life. The long-range haplotype test detected ~15% of the 3288 haplotypes that showed significant positive frequency trends (P < 0.05) and was significantly correlated with the substitution effects of the haplotypes and the selection intensities for the different traits. Thirty signatures of recent selection, which correspond to both approaches and affect the Israeli PD07 selection index, were identified on 17 of the 29 autosomes. The second experimental approach also was used to estimate the selection intensity of the different traits. The correlation between the selection intensities for the traits analyzed, derived from changes in haplotype frequencies in the population of bulls, and those derived from trait-based analysis of the cow population was 0.93 over all traits. Thus, the changes in haplotypes frequencies in the bulls' population accurately estimate genetic trends in the general cow population and can be used to detect signatures of recent selection.     

Using 2-DE for the differentiation of local chicken breeds

The aim of this study was to apply a proteomic approach for the characterisation of local chicken breeds. The experiment involved a total of 29 males of three local Italian chicken breeds: Pépoi, Padovana and Ermellinata di Rovigo. Sarcoplasmic protein fractions of breast muscle were analysed by 2-DE. Image analysis followed by statistical analysis enabled to differentiate groups of individuals based on the similarities of protein expression. Individuals were distinguished into clusters and groups, corresponding to the breed of origin. Distances among individuals, calculated using data on spot volumes, were used to draw a neighbour-joining tree, showing clear individual and breed grouping. The most relevant spots regarding breed differentiation were detected; 11 were identified by MS revealing preliminary evidences on the mechanisms of the breed differentiation process. The results evidenced the ability of proteomic analyses for the characterisation of chicken breeds.     

Annotation of selection strengths in viral genomes

MOTIVATION: Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses. RESULTS: We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.     

Signatures of natural selection in the mitochondrial genomes of Tachycineta swallows and their implications for latitudinal patterns of the ‘pace of life’

Latitudinal variation in avian life histories can be summarized as a slow–fast continuum, termed the ‘pace of life’, that encompasses patterns in life span, reproduction, and rates of development among tropical and temperate species. Much of the variation in avian pace of life is tied to differences in rates of long-term metabolic energy expenditure. Given the vital role of the mitochondrion in metabolic processes, studies of variation in the mitochondrial genome may offer opportunities to establish mechanistic links between genetic variation and latitudinal ‘pace of life’ patterns. Using comparative genomic analyses, we examined complete mitochondrial genome sequences obtained from nine, broadly distributed Tachycineta swallow species to test for signatures of natural selection across the mitogenome within a phylogenetic framework. Our results show that although purifying selection is the dominant selective force acting on the mitochondrial genome in Tachycineta, three mitochondrial genes (ND2, ND5, and CYTB) contain regions that exhibit signatures of diversifying selection. Two of these genes (ND2 and ND5) encode interacting subunits of NADH dehydrogenase, and amino residues that were inferred to be targets of positive selection were disproportionately concentrated in these genes. Moreover, the positively selected sites exhibited a phylogenetic pattern that could be indicative of adaptive divergence between “fast” and “slow” lineages. These results suggest that functional variation in cytochrome b and NADH dehydrogenase could mechanistically contribute to latitudinal ‘pace of life’ patterns in Tachycineta.     

The origin and genetic diversity of Chinese native chicken breeds

The first 539 bases of mitochondrial DNA D-loop region of six Chinese native chicken breeds (Gallus gallus domesticus) were sequenced and compared to those of the red junglefowl (Gallus gallus), the gray junglefowl (Gallus sonneratii), the green junglefowl (Gallus varius) and Lafayette's junglefowl (Gallus lafayettei) reported in GenBank, and the phylogenetic trees for the chickens were constructed based on the D-loop sequences. The results showed that the four species of the genus Gallus had great differences among each other, the G. g. domesticus was closest to the red junglefowl in Thailand and its adjacent regions, suggesting the Chinese domestic fowl probably originated from the red junglefowl in these regions. The two subs pecies of Thailand, G. g. gallus and G. g. spadiceus, should belong to one subspecies because of their resemblance. In the case of native breeds, there existed a great difference between the egg breeds and general purpose breeds, which suggested different maternal origins of the two types.     

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-934) contains supplementary material, which is available to authorized users.     

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype–phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP‐EHH (cross‐population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori‐Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.     

Untargeted metabolomics reveals differences between commercial and non-commercial Camellia sinensis cultivars used in black tea production

Journal of Plant Biochemistry and Biotechnology - Tea (Camellia sinensis) has enthralled both consumers and researchers, due to its taste, aroma and its medicinal attributes. Tea consumers concern...     

Signatures of natural selection among lineages and habitats in Oncorhynchus mykiss

Recent advances in molecular interrogation techniques now allow unprecedented genomic inference about the role of adaptive genetic divergence in wild populations. We used high-throughput genotyping to screen a genome-wide panel of 276 single nucleotide polymorphisms (SNPs) for the economically and culturally important salmonid Oncorhynchus mykiss. Samples included 805 individuals from 11 anadromous and resident populations from the northwestern United States and British Columbia, and represented two major lineages including paired populations of each life history within single drainages of each lineage. Overall patterns of variation affirmed clear distinctions between lineages and in most instances, isolation by distance within them. Evidence for divergent selection at eight candidate loci included significant landscape correlations, particularly with temperature. High diversity of two nonsynonymous mutations within the peptide-binding region of the major histocompatibility complex (MHC) class II (DAB) gene provided signatures of balancing selection. Weak signals for potential selection between sympatric resident and anadromous populations were revealed from genome scans and allele frequency comparisons. Our results suggest an important adaptive role for immune-related functions and present a large genomic resource for future studies.     

Signatures of positive selection in LY96 gene in vertebrates

As a secreted glycoprotein that binds to the extracellular domain of Toll-like receptor 4 (TLR4), Lymphocyte Antigen 96 (LY96), also called myeloid differentiation 2 (MD2), is required for the activation of TLR4 by lipopolysaccharide (LPS) and plays an important role in innate immunity, which is the first line of defence against microbial infections. Previous studies have proposed that mammalian toll-like receptors (TLRs) have evolved under diversifying selection due to their role in pathogen detection. Given the fact that LY96 is highly functionally linked to TLR4, it would be interesting to test whether LY96 is under the intense pressure of natural selection. To investigate the natural selection hypothesis, we compared the coding sequences from 13 vertebrates and evaluated the molecular evolution of LY96 gene in these species. Result shows that natural selection at exon 4 has indeed played a role in shaping the function of LY96 in the course of evolution. In addition to the study of Nakajima, we found the two branch nodes with Ka/Ks ratios greater than 1: the one leading to cow and pig and the other to rabbit and the primates.