Effect of polymorphisms in the CAMKMT gene on growth traits in Ujumqin sheep

Our previous genome‐wide association study in sheep revealed that OAR3‐84073899.1 (SNP31) in intron 8 of the CAMKMT gene was significantly associated with post‐weaning gain at the genomic level. Herein, we performed a replication study to investigate single nucleotide polymorphisms (SNPs) within the CAMKMT gene exons, and 1000 bp of the 5′‐ and 3′‐intranslated regions (UTRs) and their associations with growth traits in Ujumqin sheep. Five SNPs were identified through DNA pool sequencing technology: SNP26 in the 5′‐UTR, SNP06 in exon 5, SNP07 in exon 8 and SNP27 and SNP28 in the 3′‐UTR. Six SNPs, including SNP31 in intron 8, were genotyped in the validation group of 343 Ujumqin sheep, and each SNP was classified into three genotypes. The chi‐square test suggested that all the variations were in Hardy–Weinberg equilibrium (P > 0.05) except for SNP28 and SNP31. Linkage disequilibrium analysis showed that SNP07 and SNP31 were strongly linked. An association analysis suggested that SNP06 was significantly associated with chest girth at 6 months of age (P < 0.05). SNP07 exhibited significant correlation with body weight and chest girth at 4 months of age and with body weight, chest girth and chest width at 6 months of age (P < 0.05). SNP27 was highly associated with body weight and chest girth at 4 months of age (P < 0.05), and SNP28 was extremely significantly associated with body weight and chest girth at 4 months of age and with chest girth at 6 months of age (P < 0.01). SNP31 was significantly associated with body weight and shin circumference at 4 months of age and with post‐weaning gain (P < 0.05). Association analysis of the combined effect of SNP07 and SNP31 showed significant correlation with body weight and chest girth at four of months of age (P < 0.05) and with body weight and chest girth at 6 months of age (P < 0.05). These results indicate that the SNPs could be used as meritorious and available genetic markers in growth traits breeding and that the CAMKMT gene may be one of the key candidate genes that affect Ujumqin economic traits.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Discovery of selection-driven genetic differences of Duroc,Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.     

Molecular characterization,expression patterns and subcellular localization of Myotrophin (MTPN) gene in porcine skeletal muscle

Myotrophin (MTPN) is an effective growth factor in promoting skeletal muscle growth in vitro and vivo and has been purified from porcine skeletal muscle. However, in pigs, the information on MTPN gene is very limited. In this study, we cloned cDNA sequences and analyzed the genomic structure of porcine MTPN gene. The deduced amino acid sequence of porcine MTPN contains two the ankyrin repeat domains. RT-PCR analysis revealed that porcine MTPN gene was widely expressed in many tissues, a high expression level was observed in the spleen, liver and uterus, and transient transfection indicated that porcine MTPN proteins was located in cytoplasms within Pig Kidney Epithelial cells (PK15). Quantitative real-time PCR (qRT-PCR) analyses showed that MTPN expression peaked at embryonic 65 day post conception (dpc). During postnatal muscle development, MTPN expression was down-regulated from the 3 day to the 180 day in Yorkshire pigs. This result suggests that the MTPN gene may be important gene for skeletal muscle growth and provides useful information for further studies on its roles in porcine skeletal muscle.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

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.     

Genome‐wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries)

Fat‐tailed sheep (Ovis aries) can survive in harsh environments and satisfy human's intake of dietary fat. However, the animals require more feed, which increases the cost of farming. Thus, most farmers currently prefer thin‐tailed, short‐tailed or docked sheep. To date, the molecular mechanism of the formation of fat tails in sheep has not been completely elucidated. Here, we conducted a genome‐wide association study using phenotypes and genotypes (the Ovine Infinium HD SNP BeadChip genotype data) of two breeds of contrasting tail types (78 Small‐tailed and 78 Large‐tailed Han sheep breeds) to identify functional genes and variants associated with fat deposition. We identified four significantly (rs416433540, rs409848439, rs408118325 and rs402128848) and three approximately associated autosomal SNPs (rs401248376, rs402445895 and rs416201901). Gene annotation indicated that the surrounding genes (CREB1, STEAP4, CTBP1 and RIP140, also known as NRIP1) function in lipid storage or fat cell regulation. Furthermore, through an X‐chromosome‐wide association analysis, we detected significantly associated SNPs in the OARX: 88–89 Mb region, which could be a strong candidate genomic region for fat deposition in tails of sheep. Our results represent a new genomic resource for sheep genetics and breeding. In addition, the findings provide novel insights into genetic mechanisms of fat deposition in the tail of sheep and other mammals.     

Bayesian genome‐wide association analysis for body weight in farmed Atlantic salmon (Salmo salar L.)

We performed a genome‐wide association study to detect markers associated with growth traits in Atlantic salmon. The analyzed traits included body weight at tagging (BWT) and body weight at 25 months (BW25M). Genotypes of 4662 animals were imputed from the 50K SNP chip to the 200K SNP chip using fimpute software. The markers were simultaneously modeled using Bayes C to identify genomic regions associated with the traits. We identified windows explaining a maximum of 3.71% and 3.61% of the genetic variance for BWT and BW25M respectively. We found potential candidate genes located within the top ten 1‐Mb windows for BWT and BW25M. For instance, the vitronectin (VTN) gene, which has been previously reported to be associated with cell growth, was found within one of the top ten 1‐Mb windows for BWT. In addition, the WNT1‐inducible‐signaling pathway protein 3, melanocortin 2 receptor accessory protein 2, myosin light chain kinase, transforming growth factor beta receptor type 3 and myosin light chain 1 genes, which have been reported to be associated with skeletal growth in humans, growth stimulation during the larval stage in zebrafish, body weight in pigs, feed conversion in chickens and growth rate of sheep skeletal muscle respectively, were found within some of the top ten 1‐Mb windows for BW25M. These results indicate that growth traits are most likely controlled by many variants with relatively small effects in Atlantic salmon. The genomic regions associated with the traits studied here may provide further insight into the functional regions underlying growth traits in this species.     

Genome‐wide association study for reproductive traits in a Large White pig population

Using the PorcineSNP80 BeadChip, we performed a genome‐wide association study for seven reproductive traits, including total number born, number born alive, litter birth weight, average birth weight, gestation length, age at first service and age at first farrowing, in a population of 1207 Large White pigs. In total, we detected 12 genome‐wide significant and 41 suggestive significant SNPs associated with six reproductive traits. The proportion of phenotypic variance explained by all significant SNPs for each trait ranged from 4.46% (number born alive) to 11.49% (gestation length). Among them, 29 significant SNPs were located within known QTL regions for swine reproductive traits, such as corpus luteum number, stillborn number and litter size, of which one QTL region associated with litter size contained the ALGA0098819 SNP for total number born. Subsequently, we found that 376 functional genes contained or were near these significant SNPs. Of these, 14 genes—BHLHA15, OCM2, IL1B2, GCK, SMAD2, HABP2, PAQR5, GRB10, PRELID2, DMKN, GPI, GPIHBP1, ADCY2 and ACVR2B—were considered important candidates for swine reproductive traits based on their critical roles in embryonic development, energy metabolism and growth development. Our findings contribute to the understanding of the genetic mechanisms for reproductive traits and could have a positive effect on pig breeding programs.     

Genome-wide association study of birth weight in sheep

Birth weight is the earliest available growth trait with considerable impacts on lamb survivability and growth performance traits. This study was conducted to perform a genome-wide association study of birth weight in a meat-type sheep. A total of 132 Lori-Bakhtiari sheep were selected based on estimated of breeding values (EBVs) for BW analyses. The selected animals were genotyped using Illumina Ovine SNP50 Bead Chip. After quality control, a total of 41 323 single-nucleotide polymorphisms (SNPs) and 130 sheep were used for subsequent analyses. Plink 1.90 beta software was used for the analyses. Seven SNPs on chromosomes 1, 16, 19 and 22 were detected based on genome-wide unadjusted P-values (P <10⁻⁶), which jointly accounted for 1.2% of total genetic variation. However, based on Bonferroni-adjusted P-values, only three SNPs on chromosome 1 had significant associations with EBVs for birth weight (P <0.05), which jointly explained 0.8% of total genetic variation. A total of seven genes were found in 50 kb intervals from the three significant SNPs on chromosome 1, but only three genes, including RAB6B (a member of RAS oncogene family), Tf serotransferrin and GIGYF2 (a GRB10 interacting GYF protein 2), could be considered as candidate genes for birth weight in future studies. The results of this study may facilitate potential use of the genes involving in growth and production traits for genetic improvement of productivity in sheep.     

Genome-Wide Association for Growth Traits in Canchim Beef Cattle

Studies are being conducted on the applicability of genomic data to improve the accuracy of the selection process in livestock, and genome-wide association studies (GWAS) provide valuable information to enhance the understanding on the genetics of complex traits. The aim of this study was to identify genomic regions and genes that play roles in birth weight (BW), weaning weight adjusted for 210 days of age (WW), and long-yearling weight adjusted for 420 days of age (LYW) in Canchim cattle. GWAS were performed by means of the Generalized Quasi-Likelihood Score (GQLS) method using genotypes from the BovineHD BeadChip and estimated breeding values for BW, WW, and LYW. Data consisted of 285 animals from the Canchim breed and 114 from the MA genetic group (derived from crossings between Charolais sires and ½ Canchim + ½ Zebu dams). After applying a false discovery rate correction at a 10% significance level, a total of 4, 12, and 10 SNPs were significantly associated with BW, WW, and LYW, respectively. These SNPs were surveyed to their corresponding genes or to surrounding genes within a distance of 250 kb. The genes DPP6 (dipeptidyl-peptidase 6) and CLEC3B (C-type lectin domain family 3 member B) were highlighted, considering its functions on the development of the brain and skeletal system, respectively. The GQLS method identified regions on chromosome associated with birth weight, weaning weight, and long-yearling weight in Canchim and MA animals. New candidate regions for body weight traits were detected and some of them have interesting biological functions, of which most have not been previously reported. The observation of QTL reports for body weight traits, covering areas surrounding the genes (SNPs) herein identified provides more evidence for these associations. Future studies targeting these areas could provide further knowledge to uncover the genetic architecture underlying growth traits in Canchim cattle.     

Genomic regions influencing intramuscular fat in divergently selected rabbit lines

Intramuscular fat (IMF) is one of the main meat quality traits for breeding programmes in livestock species. The main objective of this study was to identify genomic regions associated with IMF content comparing two rabbit populations divergently selected for this trait, and to generate a list of putative candidate genes. Animals were genotyped using the Affymetrix Axiom OrcunSNP Array (200k). After quality control, the data involved 477 animals and 93 540 SNPs. Two methods were used in this research: single marker regressions with the data adjusted by genomic relatedness, and a Bayesian multiple marker regression. Associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU8 and OCU13. The highest value for the percentage of the genomic variance explained by a genomic region was found in two consecutive genomic windows on OCU8 (7.34%). Genes in the associated regions of OCU1 and OCU8 presented biological functions related to the control of adipose cell function, lipid binding, transportation and localisation (APOLD1, PLBD1, PDE6H, GPRC5D and GPRC5A) and lipid metabolic processes (MTMR2). The EWSR1 gene, underlying the OCU13 region, is linked to the development of brown adipocytes. The findings suggest that there is a large component of polygenic effect behind the differences in IMF content in these two lines, as the variance explained by most of the windows was low. The genomic regions of OCU1, OCU8 and OCU13 revealed novel candidate genes. Further studies would be needed to validate the associations and explore their possible application in selection programmes.     

Genome-wide association study of milk production traits in a crossbred dairy sheep population using three statistical models

Milk production is one of the most important characteristics of dairy sheep, and the identification of genes affecting milk production traits is critical to understanding the genetics and improve milk production in future generations. Three statistical techniques, namely GWAS, ridge-regression BLUP and BayesC , were used to identify SNPs in significant association with three milk production traits (milk yield, fat yield and protein yield) in a crossbred dairy sheep population. The results suggested that chromosomes 1, 3, 4, 5, 7 and 11 were likely to harbor genes important to milk production because these chromosomes had the greatest top-100-SNP variance contributions on the three milk production traits. The GWAS analysis identified between 74 and 288 genome-wide significant SNP (P < 0.05) whereas the BayesCπ model revealed between six and 63 SNPs, each with >95% posterior probability of inclusion as having a non-zero association effect on at least one of the three milk production traits. Positional candidate genes for milk production in sheep were searched, based on the sheep genomic assembly OAR version 3.1, such as those which map position coincided with or was located within 0.1 Mbp of a genome-wide suggestive or significant SNP. These identified SNPs and candidate genes supported some previous findings and also added new information about genetic markers for genetic improvement of lactation in dairy sheep, but keeping in mind that the majority of these positional candidate genes are not necessarily true causative loci for these traits and future validations are thus necessary.     

The effect of divergent selection for intramuscular fat on the domestic rabbit genome

An experiment of divergent selection for intramuscular fat was carried out at Universitat Politècnica de València. The high response of selection in intramuscular fat content, after nine generations of selection, and a multidimensional scaling analysis showed a high degree of genomic differentiation between the two divergent populations. Therefore, local genomic differences could link genomic regions, encompassing selective sweeps, to the trait used as selection criterion. In this sense, the aim of this study was to identify genomic regions related to intramuscular fat through three methods for detection of selection signatures and to generate a list of candidate genes. The methods implemented in this study were Wright’s fixation index, cross population composite likelihood ratio and cross population – extended haplotype homozygosity. Genomic data came from the 9th generation of the two populations divergently selected, 237 from Low line and 240 from High line. A high single nucleotide polymorphism (SNP) density array, Affymetrix Axiom OrcunSNP Array (around 200k SNPs), was used for genotyping samples. Several genomic regions distributed along rabbit chromosomes (OCU) were identified as signatures of selection (SNPs having a value above cut-off of 1%) within each method. In contrast, 8 genomic regions, harbouring 80 SNPs (OCU1, OCU3, OCU6, OCU7, OCU16 and OCU17), were identified by at least 2 methods and none by the 3 methods. In general, our results suggest that intramuscular fat selection influenced multiple genomic regions which can be a consequence of either only selection effect or the combined effect of selection and genetic drift. In addition, 73 genes were retrieved from the 8 selection signatures. After functional and enrichment analyses, the main genes into the selection signatures linked to energy, fatty acids, carbohydrates and lipid metabolic processes were ACER2, PLIN2, DENND4C, RPS6, RRAGA (OCU1), ST8SIA6, VIM (OCU16), RORA, GANC and PLA2G4B (OCU17). This genomic scan is the first study using rabbits from a divergent selection experiment. Our results pointed out a large polygenic component of the intramuscular fat content. Besides, promising positional candidate genes would be analysed in further studies in order to bear out their contributions to this trait and their feasible implications for rabbit breeding programmes.     

Genome-Wide Association Studies for Growth and Meat Production Traits in Sheep

Background

Growth and meat production traits are significant economic traits in sheep. The aim of the study is to identify candidate genes affecting growth and meat production traits at genome level with high throughput single nucleotide polymorphisms (SNP) genotyping technologies.

Methodology and Results

Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep for 11 growth and meat production traits (birth weight, weaning weight, 6-month weight, eye muscle area, fat thickness, pre-weaning gain, post-weaning gain, daily weight gain, height at withers, chest girth, and shin circumference). After quality control, 319 sheep and 48,198 SNPs were analyzed by TASSEL program in a mixed linear model (MLM). 36 significant SNPs were identified for 7 traits, and 10 of them reached genome-wise significance level for post-weaning gain. Gene annotation was implemented with the latest sheep genome Ovis_aries_v3.1 (released October 2012). More than one-third SNPs (14 out of 36) were located within ovine genes, others were located close to ovine genes (878bp-398,165bp apart). The strongest new finding is 5 genes were thought to be the most crucial candidate genes associated with post-weaning gain: s58995.1 was located within the ovine genes MEF2B and RFXANK, OAR3_84073899.1, OAR3_115712045.1 and OAR9_91721507.1 were located within CAMKMT, TRHDE, and RIPK2 respectively. GRM1, POL, MBD5, UBR2, RPL7 and SMC2 were thought to be the important candidate genes affecting post-weaning gain too. Additionally, 25 genes at chromosome-wise significance level were also forecasted to be the promising genes that influencing sheep growth and meat production traits.

Conclusions

The results will contribute to the similar studies and facilitate the potential utilization of genes involved in growth and meat production traits in sheep in future.     

Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Hybrid zones allow the measurement of gene flow across the genome, producing insight into the genomic architecture of speciation. Such analysis is particularly powerful when applied to multiple pairs of hybridizing species, as patterns of genomic differentiation can then be related to age of the hybridizing species, providing a view into the build‐up of differentiation over time. We examined 33 809 single nucleotide polymorphisms (SNPs) in three hybridizing woodpecker species: Red‐breasted, Red‐naped and Yellow‐bellied sapsuckers (Sphyrapicus ruber, Sphyrapicus nuchalis and Sphyrapicus varius), two of which (ruber and nuchalis) are much more closely related than each is to the third (varius). To identify positions of SNPs on chromosomes, we developed a localization method based on comparative genomics. We found narrow clines, bimodal distributions of hybrid indices and genomic regions with decreased rates of introgression. These results suggest moderately strong reproductive isolation among species and selection against specific hybrid genotypes. We found 19 small regions of strong differentiation between species, partly shared among species pairs, but no large regions of differentiation. An association analysis revealed a single strong‐effect candidate locus associated with plumage, possibly explaining mismatch among the three species in genomic relatedness and plumage similarity. Our comparative analysis of species pairs of different age and their hybrid zones showed that moderately strong reproductive isolation can occur with little genomic differentiation, but that reproductive isolation is incomplete even with much greater genomic differentiation, implying there are long periods of time when hybridization is possible if diverging populations are in geographic contact.     

Single nucleotide polymorphism for population studies in the scallops Aequipecten opercularis and Mimachlamys varia

This study reports the first set of SNP markers for the queen scallop Aequipecten opercularis and the black scallop Mimachlamys varia. SNPs were identified in genomic regions amplified by PCR using primers designed for introns of conserved genes by direct sequencing of PCR products. The frequency of SNPs in the two species was around one every 100 bp. Four SNPs in A. opercularis and two in M. varia were genotyped in 317–327 queen scallops and 168 black scallops sampled in Spanish and Irish populations, using an approach based on a single-base extension reaction with mass-modified terminators. The SNPs characterized provided new estimates of genetic variation in the two scallops and constitute a complementary set of genetic markers for population studies in these exploited species.     

Starch phosphorylation associated SNPs found by genome-wide association studies in the potato (Solanum tuberosum L.)

Background

The natural variation of starch phosphate content in potatoes has been previously reported. It is known that, in contrast to raw starch, commercially phosphorylated starch is more stable at high temperatures and shear rates and has higher water capacity. The genetic improvement of phosphate content in potato starch by selection or engineering would allow the production of phosphorylated starch in a natural, environmentally friendly way without chemicals. The aim of the current research is to identify genomic SNPs associated with starch phosphorylation by carrying out a genome-wide association study in potatoes.

Results

A total of 90 S. tuberosum L. varieties were used for phenotyping and genotyping. The phosphorus content of starch in 90 potato cultivars was measured and then statistically analysed. Principal component analysis (PCA) revealed that the third and eighth principal components appeared to be sensitive to variation in phosphorus content (p = 0.0005 and p = 0.002, respectively). PC3 showed the correlation of starch phosphorus content with allelic variations responsible for higher phosphorylation levels, found in four varieties. Similarly, PC8 indicated that hybrid 785/8–5 carried an allele associated with high phosphorus content, while the Impala and Red Scarlet varieties carried alleles for low phosphorus content. Genotyping was carried out using an Illumina 22 K SNP potato array. A total of 15,214 scorable SNPs (71.7% success rate) was revealed. GWAS mapping plots were obtained using TASSEL based on several statistical models, including general linear models (GLMs), with and without accounting for population structure, as well as MLM. A total of 17 significant SNPs was identified for phosphorus content in potato starch, 14 of which are assigned to 8 genomic regions on chromosomes 1, 4, 5, 7, 8, 10, and 11. Most of the SNPs identified belong to protein coding regions; however, their allelic variation was not associated with changes in protein structure or function.

Conclusions

A total of 8 novel genomic regions possibly associated with starch phosphorylation on potato chromosomes 1, 4, 5, 7, 8, 10, and 11 was revealed. Further validation of the SNPs identified and the analysis of the surrounding genomic regions for candidate genes will allow better understanding of starch phosphorylation biochemistry. The most indicative SNPs may be useful for developing diagnostic markers to accelerate the breeding of potatoes with predetermined levels of starch phosphorylation.

     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.