Genomic selection for marker-assisted improvement in line crosses

Efficiency of genomic selection with low-cost genotyping in a composite line from a cross between inbred lines was evaluated for a trait with heritability 0.10 or 0.25 using a low-density marker map. With genomic selection, selection was on the sum of estimates of effects of all marker intervals across the genome, fitted either as fixed (fixed GS) or random (random GS) effects. Reponses to selection over 10 generations, starting from the F₂, were compared with standard BLUP selection. Estimates of variance for each interval were assumed independent and equal. Both GS strategies outperformed BLUP selection, especially in initial generations. Random GS outperformed fixed GS in early generations and performed slightly better than fixed GS in later generations. Random GS gave higher genetic gain when the number of marker intervals was greater (180 or 10 cM intervals), whereas fixed GS gave higher genetic gain when the number of marker intervals was low (90 or 20 cM). Including interactions between generation and marker scores in the model resulted in lower genetic gains than models without interactions. When phenotypes were available only in the F₂ for GS, treating marker scores as fixed effects led to considerably lower genetic gain than random GS. Benefits of GS over standard BLUP were lower with high heritability. Genomic selection resulted in greater response than MAS based on only significant marker intervals (standard MAS) by increasing the frequency of QTL with both large and small effects. The efficiency of genomic selection over standard MAS depends on stringency of the threshold used for QTL detection. In conclusion, genomic selection can be effective in composite lines using a sparse marker map.     

Genomic selection and complex trait prediction using a fast EM algorithm applied to genome-wide markers

Background  

The information provided by dense genome-wide markers using high throughput technology is of considerable potential in human disease studies and livestock breeding programs. Genome-wide association studies relate individual single nucleotide polymorphisms (SNP) from dense SNP panels to individual measurements of complex traits, with the underlying assumption being that any association is caused by linkage disequilibrium (LD) between SNP and quantitative trait loci (QTL) affecting the trait. Often SNP are in genomic regions of no trait variation. Whole genome Bayesian models are an effective way of incorporating this and other important prior information into modelling. However a full Bayesian analysis is often not feasible due to the large computational time involved.     

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs

The Bama Xiang pig(BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat.To uncover the genetic basis of the BMX phenotype,we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs,respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2 NL, FRS2,KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.     

Imputation of genotypes in Danish purebred and two-way crossbred pigs using low-density panels

Background

Genotype imputation is commonly used as an initial step in genomic selection since the accuracy of genomic selection does not decline if accurately imputed genotypes are used instead of actual genotypes but for a lower cost. Performance of imputation has rarely been investigated in crossbred animals and, in particular, in pigs. The extent and pattern of linkage disequilibrium differ in crossbred versus purebred animals, which may impact the performance of imputation. In this study, first we compared different scenarios of imputation from 5 K to 8 K single nucleotide polymorphisms (SNPs) in genotyped Danish Landrace and Yorkshire and crossbred Landrace-Yorkshire datasets and, second, we compared imputation from 8 K to 60 K SNPs in genotyped purebred and simulated crossbred datasets. All imputations were done using software Beagle version 3.3.2. Then, we investigated the reasons that could explain the differences observed.

Results

Genotype imputation performs as well in crossbred animals as in purebred animals when both parental breeds are included in the reference population. When the size of the reference population is very large, it is not necessary to use a reference population that combines the two breeds to impute the genotypes of purebred animals because a within-breed reference population can provide a very high level of imputation accuracy (correct rate ≥ 0.99, correlation ≥ 0.95). However, to ensure that similar imputation accuracies are obtained for crossbred animals, a reference population that combines both parental purebred animals is required. Imputation accuracies are higher when a larger proportion of haplotypes are shared between the reference population and the validation (imputed) populations.

Conclusions

The results from both real data and pedigree-based simulated data demonstrate that genotype imputation from low-density panels to medium-density panels is highly accurate in both purebred and crossbred pigs. In crossbred pigs, combining the parental purebred animals in the reference population is necessary to obtain high imputation accuracy.

Electronic supplementary material

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

Maize transformation using xylose isomerase gene as a selection marker]

The xylA gene, encoding xylose isomerase, was cloned as a 1342-bp BamHI/SacI fragment from the E. coli. As a selection marker, the xylA gene was fused between the enhanced CaMV 35S promoter (E35S) and terminator (35St) in pBAC413 (Fig.2). pBAC413 was constructed to prevent the expression of sbeIIb in maize. PDS1000/He was used to bombard maize calli, which were induced to form by the elite inbred lines. The selection was carried out on the media containing concentrations of xylose from 0 to 100%. The results showed that the media containing 50% to 100% D-xylose were better, but differed with the genotype of maize (Tables 1 and 2). Successful integration of xylA gene into the maize genome was confirmed by DNA dot blotting, PCR and PCR-Southern hybridization (Figs.4 to 6). A method was established in which transformed maize cells were successively screened on a medium containing xylose instead of antibiotic and herbicide for bio-safety.     

Genomic regions involved in response to grain yield selection at high and low nitrogen fertilization in maize

In order to validate the role of genomic regions involved in nitrogen use efficiency and detected in a population of recombinant inbred lines (RIL), we have applied from the same population a recurrent selection for adaptation to low N-input (N0) and to high N-input (N1). Variation of allele frequency at neutral marker during the two cycles of recurrent selection may provide information about markers linked to QTLs. Significant temporal variation of allele frequency was investigated using the test of Waples, which tests the hypothesis of genetic drift versus selection. Most genomic regions (12/19) responding to selection were detected for selection at high N-input and only two were common to selection at high and low N-inputs. This was consistent with the greater grain yield response to selection observed for the population selected under high N-input compared with the population selected under low N-input, when they were evaluated at high N-fertilization. In contrast, when they were evaluated at low N-input both types of selection gave similar yield. As was expected, in the first cycle we observed selection of markers linked to grain yield QTLs. In the course of the second cycle three situations were observed: the confirmation of most regions already selected in C1 including all C1 regions overlapping with grain yield QTLs; the non-confirmation of some C1 regions (2/9); and the identification of new genomic zones (10/17). The detected marker–QTL associations revealed the consistency of the involvement of some traits, such as root architecture and glutamine synthetase activity, which would be of major importance for grain yield setting whatever the nitrogen fertilization.     

Cardiomyocyte marker expression in a human lymphocyte cell line using mouse cardiomyocyte extract

Cell transplantation shows potential for the treatment of cardiac diseases. Embryonic stem cells, cord blood and mesenchymal stem cells have been suggested as sources for transplantation therapy. Because of some technical limitations with the use of stem cells, transdifferentiation of fully differentiated cells is a potentially useful alternative. We investigated whether human peripheral blood cells could transdifferentiate into cardiomyocyte. Transdifferentiation was induced in a human B lymphocyte cell line (Raji). Cardiomyocyte extract was prepared from adult mouse cardiomyocytes. The cells were treated with 5-aza-2-deoxycytidine and trichostatin A, permeabilized with streptolysin O, and exposed to the mouse cardiomyocyte extract. They were cultured for 10 days, 3 weeks and 4 weeks. Cardiomyocyte markers were detected with immunohistochemistry and flow cytometry. Immunocytochemistry revealed that some cells expressed myosin heavy chain, α-actinin and cardiac troponin T after 3 and 4 weeks. Flow cytometry confirmed these data. In cells exposed to trichostatin A and 5-aza-2-deoxycytidine and permeabilized in the presence of the cardiomyocyte extract, troponin T expression was seen in 3.53% of the cells and 3.11% of them expressed α-actinin. After exposure to the cardiomyocyte extract, some permeabilized cells adhered to the plate loosely; however, the morphology did not change significantly, and they continued to show a rounded shape after 4 weeks. Our treated lymphocytes expressed cardiomyocyte markers. Our results suggest that lymphocytes may be useful in future research as a source of cells for reprogramming procedures.     

Economic aspects of implementing genomic evaluations in a pig sire line breeding scheme

Background

Replacing pedigree-based BLUP evaluations by genomic evaluations in pig breeding schemes can result in greater selection accuracy and genetic gains, especially for traits with limited phenotypes. However, this methodological change would generate additional costs. The objective of this study was to determine whether additional expenditures would be more profitably devoted to implementing genomic evaluations or to increasing phenotyping capacity while retaining traditional evaluations.

Methods

Stochastic simulation was used to simulate a population with 1050 breeding females and 50 boars that was selected for 10 years for a breeding goal with two uncorrelated traits with heritabilities of 0.4. The reference breeding scheme was based on phenotyping 13 770 candidates per year for trait 1 and 270 sibs of candidates per year for trait 2, with selection based on pedigree-based BLUP estimated breeding values. Increased expenditures were allocated to either increasing the phenotyping capacity for trait 2 while maintaining traditional evaluations, or to implementing genomic selection. The genomic scheme was based on two training populations: one for trait 2, consisting of phenotyped sibs of the candidates whose number increased from 1000 to 3430 over time, and one for trait 1, consisting of the selection candidates. Several genomic scenarios were tested, where the size of the training population for trait 1, and the number of genotyped candidates pre-selected based on their parental estimated breeding value, varied.

Results

Both approaches resulted in higher genetic trends for the population breeding goal and lower rates of inbreeding compared to the reference scheme. However, even a very marked increase in phenotyping capacity for trait 2 could not match improvements achieved with genomic selection when the number of genotyped candidates was large. Genotyping just a limited number of pre-selected candidates significantly reduced the extra costs, while preserving most of the benefits in terms of genetic trends and inbreeding. Implementing genomic evaluations was the most efficient approach when major expenditure was possible, whereas increasing phenotypes was preferable when limited resources were available.

Conclusions

Economic decisions on implementing genomic evaluations in a pig nucleus population must take account of population characteristics, phenotyping and genotyping costs, and available funds.     

Bayesian analysis of response to selection: a case study using litter size in Danish Yorkshire pigs

Implementation of a Bayesian analysis of a selection experiment is illustrated using litter size [total number of piglets born (TNB)] in Danish Yorkshire pigs. Other traits studied include average litter weight at birth (WTAB) and proportion of piglets born dead (PRBD). Response to selection for TNB was analyzed with a number of models, which differed in their level of hierarchy, in their prior distributions, and in the parametric form of the likelihoods. A model assessment study favored a particular form of an additive genetic model. With this model, the Monte Carlo estimate of the 95% probability interval of response to selection was (0.23; 0.60), with a posterior mean of 0.43 piglets. WTAB showed a correlated response of -7.2 g, with a 95% probability interval equal to (-33.1; 18.9). The posterior mean of the genetic correlation between TNB and WTAB was -0.23 with a 95% probability interval equal to (-0.46; -0.01). PRBD was studied informally; it increases with larger litters, when litter size is >7 piglets born. A number of methodological issues related to the Bayesian model assessment study are discussed, as well as the genetic consequences of inferring response to selection using additive genetic models.     

Genomic assisted selection for enhancing line breeding: merging genomic and phenotypic selection in winter wheat breeding programs with preliminary yield trials

Key message

Early generation genomic selection is superior to conventional phenotypic selection in line breeding and can be strongly improved by including additional information from preliminary yield trials.

Abstract

The selection of lines that enter resource-demanding multi-environment trials is a crucial decision in every line breeding program as a large amount of resources are allocated for thoroughly testing these potential varietal candidates. We compared conventional phenotypic selection with various genomic selection approaches across multiple years as well as the merit of integrating phenotypic information from preliminary yield trials into the genomic selection framework. The prediction accuracy using only phenotypic data was rather low (r = 0.21) for grain yield but could be improved by modeling genetic relationships in unreplicated preliminary yield trials (r = 0.33). Genomic selection models were nevertheless found to be superior to conventional phenotypic selection for predicting grain yield performance of lines across years (r = 0.39). We subsequently simplified the problem of predicting untested lines in untested years to predicting tested lines in untested years by combining breeding values from preliminary yield trials and predictions from genomic selection models by a heritability index. This genomic assisted selection led to a 20% increase in prediction accuracy, which could be further enhanced by an appropriate marker selection for both grain yield (r = 0.48) and protein content (r = 0.63). The easy to implement and robust genomic assisted selection gave thus a higher prediction accuracy than either conventional phenotypic or genomic selection alone. The proposed method took the complex inheritance of both low and high heritable traits into account and appears capable to support breeders in their selection decisions to develop enhanced varieties more efficiently.

     

Impaired binding of low density lipoprotein to hepatic membranes from uremic guinea pigs.

Chronic renal failure is associated with abnormalities in lipoprotein metabolism that may contribute to premature atherosclerosis and early mortality in patients on dialysis. In previous studies, we found that plasma clearance of radiolabelled low density lipoprotein (LDL) was retarded in nephrectomized guinea pigs left with one-sixth of normal functioning renal mass. To elucidate potential mechanisms of delayed LDL clearance, we compared binding of LDL to hepatic membranes from both normal and uremic guinea pigs. One hundred micrograms of the 8000-100,000 X g hepatic microsomal protein was incubated with 125I-labelled normal guinea pig LDL (10-150 micrograms/mL) for 1 h at 37 degrees C, and the membrane washed and pelleted by centrifugation in a Beckman Ti 42.2 rotor. Parallel incubations with excess unlabelled LDL were done to determine specific binding. LDL specific binding to uremic hepatic membranes was significantly impaired compared with normal ones. The major abnormality, as determined by Scatchard transformation of the binding data, was a reduction of the apparent maximal binding of LDL to uremic membranes, with an average Bmax of 4.1 micrograms/mg protein compared with 6.6 micrograms/mg protein for normal hepatic microsomes. The affinity of LDL for uremic liver membranes was only slightly diminished with a mean apparent Kd of 35.2 micrograms/mL in comparison with 21.8 micrograms/mL for normal liver membranes. These results provide a biochemical explanation for the diminished LDL clearance in uremia and may account for the dyslipidemia of renal failure.     

Efficient selection of transgenic mouse embryos using EGFP as a marker gene

We have established a reliable method that uses the EGFP (Enhanced Green Fluorescent Protein) gene as a marker for selecting transgenic embryos from preimplantation embryos. Embryos that were subjected to the pronuclear microinjection of the CMV/β‐actin/EGFP fusion gene were cultured in vitro until they developed into the morulae‐ or blastocyst‐stage. The expression of EGFP was easily observed by a fluorescent microscopy. There appeared to be no damage to the in vivo developmental ability of the embryos in response to the EGFP excitation light, which utilized an IB filter for a period of 30 min. Modified PCR analysis using Dpn I and Bal 31 digestion of the embryonic DNA showed that all of the embryos expressing EGFP in all their cells were transgenic, while more than half with mosaic expression of EGFP were not transgenic. Approximately 77% of pups born from the embryos that uniformly expressed the EGFP gene were transgenic, while 21.4% of pups from the embryos with mosaic expression were transgenics. The results showed that the use of EGFP as a marker is very useful and reliable for selecting transgenic embryos, and that it is important to transfer the embryos expressing EGFP in all their cells to obtain truly transgenic animals. Mol. Reprod. Dev. 54:43–48, 1999. © 1999 Wiley‐Liss, Inc.     

Systems for the removal of a selection marker and their combination with a positive marker

Many systems have been developed for the removal of a selection marker in order to generate marker-free transgenic plants. These systems consist of (1) a site-specific recombination system (Cre/lox) or a phage-attachment region (attP) to remove the selectable marker gene and (2) a transposable element system (Ac) or a co-transformation system to segregate the gene of interest from the selectable marker gene. Overall, the process is more time-consuming than conventional transformation methods because two rounds of transformation - two steps of regeneration or sexual crossings - are required to obtain the desired transgenic plants. Recently, removal systems combined with a positive marker, denoted as MAT vectors, have been developed to save time and effort by generating marker-free transgenic plants through a single-step transformation. We summarize here the transformation procedures using these systems and discuss their feasibility for practical use.     

Generation of double gene disruptions in Dictyostelium discoideum using a single antibiotic marker selection

Gene targeting is a powerful molecular genetic technique that has been widely used to understand specific gene function in vivo. This technique allows the ablation of an endogenous gene by recombination between an introduced DNA fragment and the homologous target gene. However, when multiple gene disruptions are needed, the availability of only a limited number of marker genes becomes a complication. Here we describe a new approach to perform double gene disruptions in Dictyostelium discoideum by simultaneous transfection of two gene targeting cassettes followed by performing clonal selection against only one marker gene. The subsequent PCR-based screens of blasticidin-resistant clones revealed the integration of both the selected and the nonselected targeting cassettes at their original respective loci creating complete gene disruptions. For the genes we have tested in these studies (myosin heavy chain kinases B and C), the efficiency of the double gene targeting event is found in the range of 2%-5% of all blasticidin-resistant colonies following the transfection step. This approach for the simultaneous disruptions of multiple genes should prove to be a valuable tool for other laboratories interested in creating multiple gene disruptants in Dictyostelium or other organisms where a limited number of selectable markers are available.     

Direct transformation of a clinical isolate of Candida parapsilosis using a dominant selection marker

Candida parapsilosis is a human pathogenic fungus with increasing importance, particularly in nosocomial infections. For detailed molecular genetic explorations of prototrophic clinical isolates of C. parapsilosis, we developed an efficient transformation system based on a dominant selectable marker. The gene encoding resistance to mycophenolic acid (MPA) was used for selection in yeast transformation. C. parapsilosis cells were transformed with a plasmid vector containing the Candida albicans inosine monophosphate dehydrogenase gene (IMH3) responsible for mycophenolic acid resistance. Transformation was carried out both by electroporation and by the lithium acetate (LiAc) method. The LiAc method resulted in very poor transformation efficiency, while the modified electroporation method yielded a high number of mitotically stable transformants exhibiting unambiguous MPA resistance. Two hundred transformants were analysed for the presence of the C. albicans IMH3(r) gene by polymerase chain reaction. Integration of single or multiple plasmid copies into the genomic DNA of C. parapsilosis was determined by Southern hybridization. To our knowledge, the present study is the first report about a method based on a dominant selectable marker for the transformation of a prototrophic, clinical isolate of C. parapsilosis. The described technique may prove to be an efficient tool for the examination of the biology and virulence of this pathogenic yeast.     

Toward genetic transformation of mitochondria in mammalian cells using a recoded drug-resistant selection marker

Due to technical difficulties, the genetic transformation of mitochondria in mammalian cells is still a challenge. In this report, we described our attempts to transform mammalian mitochondria with an engineered mitochondrial genome based on selection using a drug resistance gene. Because the standard drug-resistant neomycin phosphotransferase confers resistance to high concentrations of G418 when targeted to the mitochondria, we generated a recoded neomycin resistance gene that uses the mammalian mitochondrial genetic code to direct the synthesis of this protein in the mitochondria, but not in the nucleus (mitochondrial version). We also generated a universal version of the recoded neomycin resistance gene that allows synthesis of the drug-resistant proteins both in the mitochondria and nucleus. When we transfected these recoded neomycin resistance genes that were incorporated into the mouse mitochondrial genome clones into mouse tissue culture cells by electroporation, no DNA constructs were delivered into the mitochondria. We found that the universal version of the recoded neomycin resistance gene was expressed in the nucleus and thus conferred drug resistance to G418 selection, while the synthetic mitochondrial version of the gene produced no background drug-resistant cells from nuclear transformation. These recoded synthetic drug-resistant genes could be a useful tool for selecting mitochondrial genetic transformants as a precise technology for mitochondrial transformation is developed.     

Genomic selection in soybean: accuracy and time gain in relation to phenotypic selection

Genomic selection (GS) can potentially accelerate genetic improvement of soybean [Glycine max L. (Merrill)] by reducing the time to complete breeding cycles. The objectives of this study were to (1) explore the accuracy of GS in soybean, (2) evaluate the contribution of intrapopulational structure to the accuracy of GS, and (3) compare the efficiencies of phenotypic selection and GS in soybean. For this, phenotypic and genotypic data were collected from 324 soybean genotypes (243 recombinant inbred lines and 81 cultivars) and GS was performed for five yield related traits. BayesB methodology with a 10-fold cross-validation was used to compute accuracies. The GS accuracies were evaluated for grain yield, plant height, insertion of first pod, days to maturity, and 1000-grain weight at eight locations. We found that GS can reduce the time required to complete a selection cycle in soybean, which can lead to increased production of this commercially important crop. Furthermore, genotypic accuracy was similar regardless of population structure correction.     

Genomic islands as a marker to differentiate between clinical and environmental Burkholderia pseudomallei

Burkholderia pseudomallei, as a saprophytic bacterium that can cause a severe sepsis disease named melioidosis, has preserved several extra genes in its genome for survival. The sequenced genome of the organism showed high diversity contributed mainly from genomic islands (GIs). Comparative genome hybridization (CGH) of 3 clinical and 2 environmental isolates, using whole genome microarrays based on B. pseudomallei K96243 genes, revealed a difference in the presence of genomic islands between clinical and environmental isolates. The largest GI, GI8, of B. pseudomallei was observed as a 2 sub-GI named GIs8.1 and 8.2 with distinguishable %GC content and unequal presence in the genome. GIs8.1, 8.2 and 15 were found to be more common in clinical isolates. A new GI, GI16c, was detected on chromosome 2. Presences of GIs8.1, 8.2, 15 and 16c were evaluated in 70 environmental and 64 clinical isolates using PCR assays. A combination of GIs8.1 and 16c (positivity of either GI) was detected in 70% of clinical isolates and 11.4% of environmental isolates (P<0.001). Using BALB/c mice model, no significant difference of time to mortality was observed between K96243 isolate and three isolates without GIs under evaluation (P>0.05). Some virulence genes located in the absent GIs and the difference of GIs seems to contribute less to bacterial virulence. The PCR detection of 2 GIs could be used as a cost effective and rapid tool to detect potentially virulent isolates that were contaminated in soil.