First genome-wide CNV mapping in FELIS CATUS using next generation sequencing data.

BACKGROUND: Copy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium. RESULTS: Reads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated. CONCLUSION: This study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage.

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA-originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last. A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein-coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.

Looking at genetic structure and selection signatures of the Mexican chicken population using single nucleotide polymorphism markers.

Genetic variation enables both adaptive evolutionary changes and artificial selection. Genetic makeup of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. The aim of this study was to characterize the genetic variability of México's chicken population to reveal any underlying population structure. A total of 213 chickens were sampled in different rural production units located in 25 states of México. Genotypes were obtained using the Affymetrix Axiom® 600 K Chicken Genotyping Array. The Identity by Descent (IBD) and the principal components analysis (PCA) were performed by SVS software on pruned single nucleotide polymorphisms (SNPs).ADMIXTURE analyses identified 3 ancestors and the proportion of the genetic contribution of each of them has been determined in each individual. The results of the Neighbor-Joining (NJ) analysis resulted consistent with those obtained by the PCA. All methods utilized in this study did not allow a classification of Mexican chicken in distinct clusters or groups. A total of 3,059 run of homozygosity (ROH) were identified and, being mainly short in length (<4 Mb), these regions are indicative of a low inbreeding level in the population. Finally, findings from the ROH analysis indicated the presence of natural selective pressure in the population of Mexican chicken.The study indicates that the Mexican chicken clearly appear to be a unique creole chicken population that was not subjected to a specific artificial selection. Results provide a genetic knowledge that can be used as a basis for the genetic management of a unique and very large creole population, especially in the view of using it in production of hybrids to increase the productivity and economic revenue of family farming agriculture, which is widely present in México.

Preselection statistics and Random Forest classification identify population informative single nucleotide polymorphisms in cosmopolitan and autochthonous cattle breeds.

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.

Genomic variability in Mexican chicken population using copy number variants.

BACKGROUND: Copy number variations are genome polymorphism that influence phenotypic variation and are an important source of genetic variation in populations. The aim of this study was to investigate genetic variability in the Mexican Creole chicken population using CNVs. RESULTS: The Hidden Markov Model of the PennCNV software detected a total of 1924 CNVs in the genome of the 256 samples processed with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix). The mapped CNVs comprised 1538 gains and 386 losses, resulting at population level in 1216 CNV regions (CNVRs), of which 959 gains, 226 losses and 31 complex (i.e. containing both losses and gains). The CNVRs covered a total of 47 Mb of the whole genome sequence length, corresponding to 5.12% of the chicken galGal4 autosome assembly. CONCLUSIONS: This study allowed a deep insight into the structural variation in the genome of unselected Mexican chicken population, which up to now has not been genetically characterized. The genomic study disclosed that the population, even if presenting extreme morphological variation, cannot be organized in differentiated genetic subpopulations. Finally this study provides a chicken CNV map based on the 600 K SNP chip array jointly with a genome-wide gene copy number estimates in a native unselected for more than 500 years chicken population.

Genome-wide association study for milk somatic cell score in holstein cattle using copy number variation as markers.

Mastitis, the most common and expensive disease in dairy cows, implies significant losses in the dairy industry worldwide. Many efforts have been made to improve genetic mastitis resistance in dairy populations, but low heritability of this trait made this process not as effective as desired. The purpose of this study was to identify genomic regions explaining genetic variation of somatic cell count using copy number variations (CNVs) as markers in the Holstein population, genotyped with the Illumina BovineHD BeadChip. We found 24 and 47 copy number variation regions significantly associated with estimated breeding values for somatic cell score (SCS_EBVs) using SVS 8.3.1 and PennCNV-CNVRuler software, respectively. The association analysis performed with these two software allowed the identification of 18 candidate genes (TERT, NOTCH1, SLC6A3, CLPTM1L, PPARα, BCL-2, ABO, VAV2, CACNA1S, TRAF2, RELA, ELF3, DBH, CDK5, NF2, FASN, EWSR1 and MAP3K11) that result classified in the same functional cluster. These genes are also part of two gene networks, whose genes share the 'stress', 'cell death', 'inflammation' and 'immune response' GO terms. Combining CNV detection/association analysis based on two different algorithms helps towards a more complete identification of genes linked to phenotypic variation of the somatic cell count.

Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers.

Genomic and genetic variation among six Italian chicken native breeds (Livornese, Mericanel della Brianza, Milanino, Bionda Piemontese, Bianca di Saluzzo and Siciliana) were studied using single nucleotide polymorphism (SNP) and copy number variants (CNV) as markers. A total of 94 DNA samples genotyped with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix) were used in the analyses. The results showed the genetic and genomic variability occurring among the six Italian chicken breeds. The genetic relationship among animals was established with a principal component analysis. The genetic diversity within breeds was calculated using heterozygosity values (expected and observed) and with Wright's F-statistics. The individual-based CNV calling, based on log R ratio and B-allele frequency values, was done by the Hidden-Markov Model (HMM) of PennCNV software on autosomes. A hierarchical agglomerative clustering was applied in each population according to the absence or presence of definite CNV regions (CNV were grouped by overlapping of at least 1 bp). The CNV map was built on a total of 1003 CNV found in individual samples, after grouping by overlaps, resulting in 564 unique CNV regions (344 gains, 213 losses and 7 complex), for a total of 9.43 Mb of sequence and 1.03% of the chicken assembly autosome. All the approaches using SNP data showed that the Siciliana breed clearly differentiate from other populations, the Livornese breed separates into two distinct groups according to the feather colour (i.e. white and black) and the Bionda Piemontese and Bianca di Saluzzo breeds are closely related. The genetic variability found using SNP is comparable with that found by other authors in the same breeds using microsatellite markers. The CNV markers analysis clearly confirmed the SNP results.

Quantitative trait loci mapping for conjugated linoleic acid, vaccenic acid and ∆(9) -desaturase in Italian Brown Swiss dairy cattle using selective DNA pooling.

A selective DNA pooling approach was applied to identify QTL for conjugated linoleic acid (CLA), vaccenic acid (VA) and Δ(9) -desaturase (D9D) milk content in Italian Brown Swiss dairy cattle. Milk samples from 60 animals with higher values (after correction for environmental factors) and 60 animals with lower values for each of these traits from each of five half-sib families were pooled separately. The pools were genotyped using the Illumina BovineSNP50 BeadChip. Sire allele frequencies were compared between high and low tails at the sire and marker level for SNPs for which the sires were heterozygous. An r procedure was implemented to perform data analysis in a selective DNA pooling design. A correction for multiple tests was applied using the proportion of false positives among all test results. BTA 19 showed the largest number of markers in association with CLA. Associations between SNPs and the VA and Δ(9) -desaturase traits were found on several chromosomes. A bioinformatics survey identified genes with an important role in pathways for milk fat and fatty acids metabolism within 1 Mb of SNP markers associated with fatty acids contents.

An international parentage and identification panel for the domestic cat (Felis catus).

Seventeen commercial and research laboratories participated in two comparison tests under the auspices of the International Society for Animal Genetics to develop an internationally tested, microsatellite-based parentage and identification panel for the domestic cat (Felis catus). Genetic marker selection was based on the polymorphism information content and allele ranges from seven random-bred populations (n = 261) from the USA, Europe and Brazil and eight breeds (n = 200) from the USA. Nineteen microsatellite markers were included in the comparison test and genotyped across the samples. Based on robustness and efficiency, nine autosomal microsatellite markers were ultimately selected as a single multiplex 'core' panel for cat identification and parentage testing. Most markers contained dinucleotide repeats. In addition to the autosomal markers, the panel included two gender-specific markers, amelogenin and zinc-finger XY, which produced genotypes for both the X and Y chromosomes. This international cat parentage and identification panel has a power of exclusion comparable to panels used in other species, ranging from 90.08% to 99.79% across breeds and 99.47% to 99.87% in random-bred cat populations.

Short communication: quantitative trait loci affecting the somatic cell score on chromosomes 4 and 26 in Italian Holstein cattle.

This work aimed to confirm previously reported quantitative trait loci (QTL) affecting the somatic cell score (SCS) in dairy cattle on Bos taurus autosomes (BTA) 4 and 26. A granddaughter design with selective genotyping was implemented that included half-sib families from 12 male lines of Italian Holstein cattle. The animals were genotyped for 5 microsatellite markers each on regions of BTA 4 (average marker spacing 9.42 cM) and BTA 26 (average marker spacing 5.26 cM), previously reported by other authors as carrying QTL for somatic cell count. Quantitative trait loci analyses were performed using interval mapping by regressing sire breeding values for SCS onto genotype probabilities at 1-cM intervals along the 2 chromosome regions. Breeding values for SCS were estimated for the whole population using a test-day repeatability animal model. Results were not significant on a chromosome basis, but a possible QTL was found at BM4505 on BTA 26, confirming this region for further studies of QTL affecting SCS in the Italian Holstein population.