Targeted genotyping by sequencing: a new way to genome profile the cat

Targeted GBS is a recent approach for obtaining an effective characterization for hundreds to thousands of markers. The high throughput of next-generation sequencing technologies, moreover, allows sample multiplexing. The aims of this study were to (i) define a panel of single nucleotide polymorphisms (SNPs) in the cat, (ii) use GBS for profiling 16 cats, and (iii) evaluate the performance with respect to the inference using standard approaches at different coverage thresholds, thereby providing useful information for designing similar experiments. Probes for sequencing 230 variants were designed based on the Felis_catus_8.0. 8.0 genome. The regions comprised anonymous and non-anonymous SNPs. Sixteen cat samples were analysed, some of which had already been genotyped in a large group of loci and one having been whole-genome sequenced in the 99_Lives Cat Genome Sequencing Project. The accuracy of the method was assessed by comparing the GBS results with the genotypes already available. Overall, GBS achieved good performance, with 92-96% correct assignments, depending on the coverage threshold used to define the set of trustable genotypes. Analyses confirmed that (i) the reliability of the inference of each genotype depends on the coverage at that locus and (ii) the fraction of target loci whose genotype can be inferred correctly is a function of the total coverage. GBS proves to be a valid alternative to other methods. Data suggested a depth of less than 11× is required for greater than 95% accuracy. However, sequencing depth must be adapted to the total size of the targets to ensure proper genotype inference.

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.

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.

Gene expression of hepatic glucocorticoid receptor NR3C1 and correlation with plasmatic corticosterone in Italian chickens

This study examined breed-specific stress-related hormonal and gene expression profiles in three Italian chicken breeds (Valdarnese Bianca, Bionda Piemontese, Robusta Maculata) reared in controlled conditions. Glucocorticoids work through the glucocorticoid receptor (GR), which modulates target genes transcription. We investigated breed-specific changes in corticosterone (ELISA) and GR expression. GR mRNA levels were analyzed using one-tube, two-temperature real-time PCR for absolute quantification of the gene expression by the standard curve method. Our results show high expression of GR in hepatic tissue. Significant effect of the breed was recorded for plasma corticosterone concentration: Valdarnese Bianca 3.35 ng/mL, Bionda Piemontese 1.73 ng/mL, Robusta Maculata 2.02 ng/mL. Breed specific gene expression has been recorded with a GR ranging from 1.12E+04 (Robusta Maculata) to 1.00E+05 (Bionda Piemontese) mRNA copy number/100 ng total RNA. Negative correlation was found between gene expression and blood corticosterone level.