Footprints of recent selection and variability in breed composition in the Göttingen Minipig genome

The Göttingen Minipig (GMP) developed at the University of Göttingen is a synthetic breed that is widely used in medical research and toxicology. It combines the high fertility of the Vietnamese potbellied pig, the low body weight of the Minnesota Minipig and the white coat colour of the German Landrace pig. The aim of this study was to find genomic regions that may have undergone selection since the creation of the breed in the 1960s. Therefore, the whole genome was screened for footprints of recent selection based on single nucleotide polymorphism (SNP) genotypes from the Illumina Porcine SNP60 BeadChip using two methods: the extended haplotype homozygosity (EHH) test and the estimation of the genomic proportion of the three original breeds at each SNP using a Bayesian approach. Local deviations from the average genome‐wide breed composition were tested with a permutation‐based empirical test. Results for a comprehensive whole‐genome scan for both methods are presented. Several regions showing the highest P‐values in the EHH test are related to breeding goals relevant in the GMP, such as growth (SOCS2, TXN, DDR2 and GRB10 genes) and white colour (PRLR gene). Additionally, the calculated proportion of the founder breeds diverged significantly in many regions from the pedigree‐based expectations and the genome average. The results provide a genome‐wide map of selection signatures in the GMP, which leads to a better understanding of selection that took place over the last decades in GMP breed development.     

Runs of homozygosity and distribution of functional variants in the cattle genome

Background

Recent developments in sequencing technology have facilitated widespread investigations of genomic variants, including continuous stretches of homozygous genomic regions. For cattle, a large proportion of these runs of homozygosity (ROH) are likely the result of inbreeding due to the accumulation of elite alleles from long-term selective breeding programs. In the present study, ROH were characterized in four cattle breeds with whole genome sequence data and the distribution of predicted functional variants was detected in ROH regions and across different ROH length classes.

Results

On average, 19.5 % of the genome was located in ROH across four cattle breeds. There were an average of 715.5 ROH per genome with an average size of ~750 kbp, ranging from 10 (minimum size considered) to 49,290 kbp. There was a significant correlation between shared short ROH regions and regions putatively under selection (p < 0.001). By investigating the relationship between ROH and the predicted deleterious and non-deleterious variants, we gained insight into the distribution of functional variation in inbred (ROH) regions. Predicted deleterious variants were more enriched in ROH regions than predicted non-deleterious variants, which is consistent with observations in the human genome. We also found that increased enrichment of deleterious variants was significantly higher in short (<100 kbp) and medium (0.1 to 3 Mbp) ROH regions compared with long (>3 Mbp) ROH regions (P < 0.001), which is different than what has been observed in the human genome.

Conclusions

This study illustrates the distribution of ROH and functional variants within ROH in cattle populations. These patterns are different from those in the human genome but consistent with the natural history of cattle populations, which is confirmed by the significant correlation between shared short ROH regions and regions putatively under selection. These findings contribute to understanding the effects of inbreeding and probably selection in shaping the distribution of functional variants in the cattle genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1715-x) contains supplementary material, which is available to authorized users.     

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (F_ROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (F_HOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (F_PED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. F_HOM was highly correlated (0.671 to 0.985) with F_ROH measures in all breeds. Apulo-Calabrese and Casertana had the highest F_ROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (F_ROH1) to 16 Mbp (F_ROH16)), whereas the lowest F_ROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). F_ROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all F_ROH measures and F_PED (from 0.514 to 0.523) and between F_HOM and F_PED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between F_ROH and F_PED. This might be due to the imperfect measure of F_PED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.     

Microsatellite analysis of different ploidy offspring of artificial gynogenesis in Cyprinus carpio

Gynogenesis was induced by using UV-irradiated spermatozoa of blunt snout bream Megalobrama amblycephala to activate eggs of common carp Cyprinus carpio. The maternal genome was then duplicated by cold shock in 0 to 4° C cold water to retain the second polar body. Two kinds of fry, normal fry and abnormal tortuous fry, were hatched. Their DNA content was measured by flow cytometry. The normal fry were identified as diploid, representing the successful gynogenesis in C. carpio whereas the abnormal tortuous fry were haploid. Ten microsatellite loci were used to study the genetic diversity among C. carpio, diploid gynogenetic C. carpio and unduplicated haploid tortuous fry. The results indicated that the genetic homozygosity of gynogenetic C. carpio was significantly higher than that of C. carpio. The genetic homozygosity of the haploid C. carpio was intermediate between that of gynogenetic C. carpio and C. carpio. It might be easier for the allogenetic DNA fragments to be integrated into the haploid genome than into diploid gynogenetic genome.     

A genome scan for selection signatures in Taihu pig breeds using next-generation sequencing

Taihu pig breeds are the most prolific breeds of swine in the world, and they also have superior economic traits, including high resistance to disease, superior meat quality, high resistance to crude feed and a docile temperament. The formation of these phenotypic characteristics is largely a result of long-term artificial or natural selection. Therefore, exploring selection signatures in the genomes of the Taihu pigs will help us to identify porcine genes related to productivity traits, disease and behaviour. In this study, we used both intra-population (Relative Extend Haplotype Homozygosity Test (REHH)) and inter-population (the Cross-Population Extend Haplotype Homozygosity Test (XPEHH); F-STATISTICS, F_ST) methods to detect genomic regions that might be under selection process in Taihu pig breeds. As a result, we found 282 (REHH) and 112 (XPEHH) selection signature candidate regions corresponding to 159.78 Mb (6.15%) and 62.29 Mb (2.40%) genomic regions, respectively. Further investigations of the selection candidate regions revealed that many genes under these genomic regions were related to reproductive traits (such as the TLR9 gene), coat colour (such as the KIT gene) and fat metabolism (such as the CPT1A and MAML3 genes). Furthermore, gene enrichment analyses showed that genes under the selection candidate regions were significantly over-represented in pathways related to diseases, such as autoimmune thyroid and asthma diseases. In conclusion, several candidate genes potentially under positive selection were involved in characteristics of Taihu pig. These results will further allow us to better understand the mechanisms of selection in pig breeding.     

A genome‐wide scan for signatures of recent selection in Holstein cattle

The data from the newly available 50 K SNP chip was used for tagging the genome‐wide footprints of positive selection in Holstein–Friesian cattle. For this purpose, we employed the recently described Extended Haplotype Homozygosity test, which detects selection by measuring the characteristics of haplotypes within a single population. To assess formally the significance of these results, we compared the combination of frequency and the Relative Extended Haplotype Homozygosity value of each core haplotype with equally frequent haplotypes across the genome. A subset of the putative regions showing the highest significance in the genome‐wide EHH tests was mapped. We annotated genes to identify possible influence they have in beneficial traits by using the Gene Ontology database. A panel of genes, including FABP3, CLPN3, SPERT, HTR2A5, ABCE1, BMP4 and PTGER2, was detected, which overlapped with the most extreme P‐values. This panel comprises some interesting candidate genes and QTL, representing a broad range of economically important traits such as milk yield and composition, as well as reproductive and behavioural traits. We also report high values of linkage disequilibrium and a slower decay of haplotype homozygosity for some candidate regions harbouring major genes related to dairy quality. The results of this study provide a genome‐wide map of selection footprints in the Holstein genome, and can be used to better understand the mechanisms of selection in dairy cattle breeding.     

Dual role of FMN in flavodoxin function: Electron transfer cofactor and modulation of the protein-protein interaction surface

Flavodoxin (Fld) replaces Ferredoxin (Fd) as electron carrier from Photosystem I (PSI) to Ferredoxin-NADP⁺ reductase (FNR). A number of Anabaena Fld (AnFld) variants with replacements at the interaction surface with FNR and PSI indicated that neither polar nor hydrophobic residues resulted critical for the interactions, particularly with FNR. This suggests that the solvent exposed benzenoid surface of the Fld FMN cofactor might contribute to it. FMN has been replaced with analogues in which its 7- and/or 8-methyl groups have been replaced by chlorine and/or hydrogen. The oxidised Fld variants accept electrons from reduced FNR more efficiently than Fld, as expected from their less negative midpoint potential. However, processes with PSI (including reduction of Fld semiquinone by PSI, described here for the first time) are impeded at the steps that involve complex re-arrangement and electron transfer (ET). The groups introduced, particularly chlorine, have an electron withdrawal effect on the pyrazine and pyrimidine rings of FMN. These changes are reflected in the magnitude and orientation of the molecular dipole moment of the variants, both factors appearing critical for the re-arrangement of the finely tuned PSI:Fld complex. Processes with FNR are also slightly modulated. Despite the displacements observed, the negative end of the dipole moment points towards the surface that contains the FMN, still allowing formation of complexes competent for efficient ET. This agrees with several alternative binding modes in the FNR:Fld interaction. In conclusion, the FMN in Fld not only contributes to the redox process, but also to attain the competent interaction of Fld with FNR and PSI.     

Evaluation of an antisense RNA transgene for inhibiting growth hormone gene expression in transgenic rats

We compared the levels of growth hormone (GH) mRNA in the pituitary, plasma GH concentration, and altered phenotype in rats heterozygous and homozygous for an antisense RNA transgene targeted to the rat GH gene, with those in nontransgenic rats. We initially investigated whether the transgene promoter, which is connected to four copies of a thyroid hormone response element (TRE) that increases promoter activity, affected in vivo transgene expression in the pituitary of the transgenic rats. Plasma GH concentration correlated negatively with T, injection in surgically thyroidectomized heterozygous transgenic rats. There was a reduction of about ?35–40% in GH mRNA levels in the pituitary of homozygous animals compared with those in non-transgenic rats. Plasma GH concentration was significantly ?25–32 and ?29–41% lower in heterozygous and homozygous transgenic rats, respectively, compared with that in nontransgenic animals. Furthermore, the growth rates in homozygous transgenic rats were reduced by ?72–81 and ?51–70% compared with those of their heterozygous and nontransgenic littermates, respectively. The results of these studies suggested that the biological effect of GH in vivo is modulated dose-dependently by the antisense RNA transgene. The rat GH gene can therefore be targeted by antisense RNA produced from a transgene, as reflected in the protein and RNA levels. © 1995 Wiley-Liss, Inc.     

Induced Androgenesis in Fish: Obtaining Viable Nucleocytoplasmic Hybrids

A self-review of long-standing studies of induced androgenesis in fish. The data are provided on the induction of diploid androgenesis in teleost fish via blocking the first cleavage division and the development of the method of dispermic androgenesis in sturgeon fishes. Androgenetic nucleocytoplasmic hybrids were obtained using the method of dispermic androgenesis and eggs of previously produced true interspecific hybrids. Specific morphological features of androgenetic hybrids were described and dispermic androgenesis was realized using cryopreserved sperm. The mechanism of sex determination in sturgeon fishes is discussed with reference to possible use of androgenesis for restoration of rare and disappearing species.__________Translated from Ontogenez, Vol. 36, No. 4, 2005, pp. 254–264.Original Russian Text Copyright © 2005 by Grunina, Recoubratsky.     

Variational and genetic properties of developmental stability in Dalechampia scandens

Because low developmental stability may compromise the precision with which adaptations can be reached, the variability and genetic basis of developmental stability are important evolutionary parameters. Developmental stability is also an important clue to understanding how traits are regulated to achieve their phenotypic target value. However, developmental stability must be studied indirectly through proxy variables, such as fluctuating asymmetry, that are suggested to have noisy and often nonlinear relationships to the underlying variable of interest. In this paper we first show that mean-standardized measures of variance and covariance in fluctuating asymmetry, unlike heritabilities, repeatabilities, and correlations, are linearly related to corresponding measures of variation in underlying developmental stability. We then examine the variational properties of developmental stability in a population of the Neotropical vine, Dalechampia scandens (Euphorbiaceae). By studying fluctuating asymmetry in a large number of floral characters in both selfed and outcrossed individuals in a diallel design, we assemble strong evidence that both additive genetic and individual variation and covariation in developmental stability are virtually absent in this population.