Estimation of the autosomal contribution to total additive genetic variability of carcass traits in Japanese Black cattle.

We attempted to estimate the additive genetic variance explained by each autosome, using genotype data of 33,657 single nucleotide polymorphism (SNP) markers in 2271 Japanese Black fattened steers. Traits were cold carcass weight, ribeye area, rib thickness, subcutaneous fat thickness, estimated yield percentage, and marbling score. Two mixed linear models were used: One is that (model 1) incorporating a genomic relationship matrix (G matrix) constructed by using all available SNPs, and another (model 2), incorporating two G matrices constructed by using the SNPs on one autosome and using those on the remaining autosomes. Genomic heritabilities estimated using model 1 were moderate to high. The sums of the proportions of the additive genetic variance explained by each autosome to the total genetic variance estimated by using model 2 were >90%. For carcass weight, the proportions explained by Bos taurus autosomes 6, 8, and 14 were higher than those explained by the remaining autosomes. In some cases, the estimated proportion was close to 0. The results obtained from model 2 could provide a novel insight into the genetic architecture, such as heritability per chromosome, of carcass traits in Japanese Black cattle, although further careful investigation would be required.

Fitting Genomic Prediction Models with Different Marker Effects among Prefectures to Carcass Traits in Japanese Black Cattle.

We fitted statistical models, which assumed single-nucleotide polymorphism (SNP) marker effects differing across the fattened steers marketed into different prefectures, to the records for cold carcass weight (CW) and marbling score (MS) of 1036, 733, and 279 Japanese Black fattened steers marketed into Tottori, Hiroshima, and Hyogo prefectures in Japan, respectively. Genotype data on 33,059 SNPs was used. Five models that assume only common SNP effects to all the steers (model 1), common effects plus SNP effects differing between the steers marketed into Hyogo prefecture and others (model 2), only the SNP effects differing between Hyogo steers and others (model 3), common effects plus SNP effects specific to each prefecture (model 4), and only the effects specific to each prefecture (model 5) were exploited. For both traits, slightly lower values of residual variance than that of model 1 were estimated when fitting all other models. Estimated genetic correlation among the prefectures in models 2 and 4 ranged to 0.53 to 0.71, all <0.8. These results might support that the SNP effects differ among the prefectures to some degree, although we discussed the necessity of careful consideration to interpret the current results.

An intersection network based on combining SNP coassociation and RNA coexpression networks for feed utilization traits in Japanese Black cattle.

Genome-wide association studies (GWAS) of quantitative traits have detected numerous genetic associations, but they encounter difficulties in pinpointing prominent candidate genes and inferring gene networks. The present study used a systems genetics approach integrating GWAS results with external RNA-expression data to detect candidate gene networks in feed utilization and growth traits of Japanese Black cattle, which are matters of concern. A SNP coassociation network was derived from significant correlations between SNPs with effects estimated by GWAS across 7 phenotypic traits. The resulting network genes contained significant numbers of annotations related to the traits. Using bovine transcriptome data from a public database, an RNA coexpression network was inferred based on the similarity of expression patterns across different tissues. An intersection network was then generated by superimposing the SNP and RNA networks and extracting shared interactions. This intersection network contained 4 tissue-specific modules: nervous system, reproductive system, muscular system, and glands. To characterize the structure (topographical properties) of the 3 networks, their scale-free properties were evaluated, which revealed that the intersection network was the most scale-free. In the subnetwork containing the most connected transcription factors (URI1, ROCK2, and ETV6), most genes were widely expressed across tissues, and genes previously shown to be involved in the traits were found. Results indicated that the current approach might be used to construct a gene network that better reflects biological information, providing encouragement for the genetic dissection of economically important quantitative traits.

Structure and polymorphisms of the major histocompatibility complex in the Oriental stork, Ciconia boyciana.

The major histocompatibility complex (MHC) is highly polymorphic and plays a central role in the vertebrate immune system. Despite its functional consistency, the MHC genomic structure differs substantially among organisms. In birds, the MHCs of Galliformes and the Japanese crested ibis (Pelecaniformes) are well-characterized, but information about other avian MHCs remains scarce. The Oriental stork (Ciconia boyciana, order Ciconiiformes) is a large endangered migrant. The current Japanese population of this bird originates from a few founders; thus, understanding the genetic diversity among them is critical for effective population management. We report the structure and polymorphisms in C. boyciana MHC. One contig (approximately 128 kb) was assembled by screening of lambda phage genomic library and its complete sequence was determined, revealing a gene order of COL11A2, two copies of MHC-IIA/IIB pairs, BRD2, DMA/B1/B2, MHC-I, TAP1/2, and two copies each of pseudo MHC-I and TNXB. This structure was highly similar to that of the Japanese crested ibis, but largely different from that of Galliformes, at both the terminal regions. Genotyping of the MHC-II region detected 10 haplotypes among the six founders. These results provide valuable insights for future studies on the evolution of the avian MHCs and for conservation of C. boyciana.

Change of fatty acid composition of the lumbar longissimus during the final stage of fattening in the Japanese Black cattle.

Consideration of the shortened fattening period seems to be worthwhile for the realization of profitable beef production. In this study, change of fatty acid composition of the lumbar longissimus during the final stage of fattening was investigated in Japanese Black cattle. Each of 110 fattening animals was sampled three times: the initial two samples were taken by biopsy (25.7 months and 27.5 months on average) and the final one was from carcasses (29.9 months on average). Preliminary analysis indicated that removing muscle tissues from the constant body position of the living animals should be essential for sampling. Average monounsaturated fatty acids (MUFA) at three sampling points were 58.1%, 58.5% and 60.5%, and the differences of the third sampling with the first and second samplings were significant. Both in steers and heifers, MUFA also increased as the fattening stage proceeded, and MUFA of the heifers at all the sampling points were significantly higher than those of the steers. The increasing rate of MUFA rose from 0.21 percentage points (pp)/month at period 1 (from the first sampling to the second sampling) to 0.84 pp/month at period 2 (from the second sampling to the slaughter).

Estimation of variance and genomic prediction using genotypes imputed from low-density marker subsets for carcass traits in Japanese black cattle.

The influence of genotype imputation using low-density single nucleotide polymorphism (SNP) marker subsets on the genomic relationship matrix (G matrix), genetic variance explained, and genomic prediction (GP) was investigated for carcass weight and marbling score in Japanese Black fattened steers, using genotype data of approximately 40,000 SNPs. Genotypes were imputed using equally spaced SNP subsets of different densities. Two different linear models were used. The first (model 1) incorporated one G matrix, while the second (model 2) used two different G matrices constructed using the selected and remaining SNPs. When using model 1, the estimated additive genetic variance was always larger when using all SNPs obtained via genotype imputation than when using only equally spaced SNP subsets. The correlations between the genomic estimated breeding values obtained using genotype imputation with at least 3,000 SNPs and those using all available SNPs without imputation were higher than 0.99 for both traits. While additive genetic variance was likely to be partitioned with model 2, it did not enhance the accuracy of GP compared with model 1. These results indicate that genotype imputation using an equally spaced low-density panel of an appropriate size can be used to produce a cost-effective, valid GP.

Accuracy of imputation of single nucleotide polymorphism marker genotypes from low-density panels in Japanese Black cattle.

Using target and reference fattened steer populations, the performance of genotype imputation using lower-density marker panels in Japanese Black cattle was evaluated. Population imputation was performed using BEAGLE software. Genotype information for approximately 40,000 single nucleotide polymorphism (SNP) markers by Illumina BovineSNP50 BeadChip was available, and imputation accuracy was assessed based on the average concordance rates of the genotypes, varying equally spaced SNP densities, and the number of individuals in the reference population. Two additional statistics were also calculated as indicators of imputation performance. The concordance rates tended to be lower for SNPs with greater minor allele frequencies, or those located near the ends of the chromosomes. Longer autosomes yielded greater imputation accuracies than shorter ones. When SNPs were selected based on linkage disequilibrium information, relative imputation accuracy was slightly improved. When 3000 and 10,000 equally spaced SNPs were used, the imputation accuracies were greater than 90% and approximately 97%, respectively. These results indicate that combining genotyping using a lower-density SNP chip with genotype imputation based on a population of individuals genotyped using a higher-density SNP chip is a cost-effective and valid approach for genomic prediction.

Genetic Analyses in the Japanese Captive Population of Japanese Crested Ibis (Nipponia nippon) using Pedigree Information.

The Japanese captive population of Japanese crested ibis (Nipponia nippon) was established using 5 founders derived from the Chinese captive population. Its size has increased rapidly, and the maintenance phase is about to start. Thus, this study was designed to perform genetic analyses in this population with pedigree information, considering the adoption of mean kinship strategy as the breeding strategy suited to the maintenance phase. Because the relationships among the 5 founders were unknown, different assumptions were set up ranging from 0 to 0.25 of kinship coefficients between the 5 founders. Assuming that the 5 founders were non-inbred in all the assumptions, the results showed that the gene diversity and the mean inbreeding coefficient would fluctuate largely from ∼65% to ∼82% and from ∼0.07 to ∼0.29, respectively. Moreover, the genetic importance of individuals based on mean kinship shifted largely. This study suggested that the Japanese captive population had low gene diversity and high mean inbreeding coefficient even under the assumption that the 5 founders were unrelated and non-inbred. In addition, the study also suggested that it became more effective to analyze the genetic status and to introduce mean kinship strategy into this population with more credible molecular evaluation of the relationships among founders.

Predicting the phenotypic values of physiological traits using SNP genotype and gene expression data in mice.

Predicting phenotypes using genome-wide genetic variation and gene expression data is useful in several fields, such as human biology and medicine, as well as in crop and livestock breeding. However, for phenotype prediction using gene expression data for mammals, studies remain scarce, as the available data on gene expression profiling are currently limited. By integrating a few sources of relevant data that are available in mice, this study investigated the accuracy of phenotype prediction for several physiological traits. Gene expression data from two tissues as well as single nucleotide polymorphisms (SNPs) were used. For the studied traits, the variance of the effects of the expression levels was more likely to differ among the genes than were the effects of SNPs. For the glucose concentration, the total cholesterol amount, and the total tidal volume, the accuracy by cross validation tended to be higher when the gene expression data rather than the SNP genotype data were used, and a statistically significant increase in the accuracy was obtained when the gene expression data from the liver were used alone or jointly with the SNP genotype data. For these traits, there were no additional gains in accuracy from using the gene expression data of both the liver and lung compared to that of individual use. The accuracy of prediction using genes that were selected differently was examined; the use of genes with a higher tissue specificity tended to result in an accuracy that was similar to or greater than that associated with the use of all of the available genes for traits such as the glucose concentration and total cholesterol amount. Although relatively few animals were evaluated, the current results suggest that gene expression levels could be used as explanatory variables. However, further studies are essential to confirm our findings using additional animal samples.

Structure and polymorphism of the major histocompatibility complex class II region in the Japanese Crested Ibis, Nipponia nippon.

The major histocompatibility complex (MHC) is a highly polymorphic genomic region that plays a central role in the immune system. Despite its functional consistency, the genomic structure of the MHC differs substantially among organisms. In birds, the MHC-B structures of Galliformes, including chickens, have been well characterized, but information about other avian MHCs remains sparse. The Japanese Crested Ibis (Nipponia nippon, Pelecaniformes) is an internationally conserved, critically threatened species. The current Japanese population of N. nippon originates from only five founders; thus, understanding the genetic diversity among these founders is critical for effective population management. Because of its high polymorphism and importance for disease resistance and other functions, the MHC has been an important focus in the conservation of endangered species. Here, we report the structure and polymorphism of the Japanese Crested Ibis MHC class II region. Screening of genomic libraries allowed the construction of three contigs representing different haplotypes of MHC class II regions. Characterization of genomic clones revealed that the MHC class II genomic structure of N. nippon was largely different from that of chicken. A pair of MHC-IIA and -IIB genes was arranged head-to-head between the COL11A2 and BRD2 genes. Gene order in N. nippon was more similar to that in humans than to that in chicken. The three haplotypes contained one to three copies of MHC-IIA/IIB gene pairs. Genotyping of the MHC class II region detected only three haplotypes among the five founders, suggesting that the genetic diversity of the current Japanese Crested Ibis population is extremely low. The structure of the MHC class II region presented here provides valuable insight for future studies on the evolution of the avian MHC and for conservation of the Japanese Crested Ibis.

Effects of single nucleotide polymorphism marker density on degree of genetic variance explained and genomic evaluation for carcass traits in Japanese Black beef cattle.

BACKGROUND: Japanese Black cattle are a beef breed whose meat is well known to excel in meat quality, especially in marbling, and whose effective population size is relatively low in Japan. Unlike dairy cattle, the accuracy of genomic evaluation (GE) for carcass traits in beef cattle, including this breed, has been poorly studied. For carcass weight and marbling score in the breed, as well as the extent of whole genome linkage disequilibrium (LD), the effects of equally-spaced single nucleotide polymorphisms (SNPs) density on genomic relationship matrix (G matrix), genetic variance explained and GE were investigated using the genotype data of about 40,000 SNPs and two statistical models. RESULTS: Using all pairs of two adjacent SNPs in the whole SNP set, the means of LD (r2) at ranges 0-0.1, 0.1-0.2, 0.2-0.5 and 0.5-1 Mb were 0.22, 0.13, 0.10 and 0.08, respectively, and 25.7, 13.9, 10.4 and 6.4% of the r2 values exceeded 0.3, respectively. While about 90% of the genetic variance for carcass weight estimated using all available SNPs was explained using 4,000-6,000 SNPs, the corresponding percentage for marbling score was consistently lower. With the conventional linear model incorporating the G matrix, correlation between the genomic estimated breeding values (GEBVs) obtained using 4,000 SNPs and all available SNPs was 0.99 for carcass weight and 0.98 for marbling score, with an underestimation of the former GEBVs, especially for marbling score. CONCLUSIONS: The Japanese Black is likely to be in a breed group with a relatively high extent of whole genome LD. The results indicated that the degree of marbling is controlled by only QTLs with relatively small effects, compared with carcass weight, and that using at least 4,000 equally-spaced SNPs, there is a possibility of ranking animals genetically for these carcass traits in this breed.

Estimation of variance components for carcass traits in Japanese Black cattle using 50K SNP genotype data.

Genomic selection using high-density single nucleotide polymorphism (SNP) genotype data may accelerate genetic improvements in livestock animals. In this study, we attempted to estimate the variance components of six carcass traits in fattened Japanese Black steers using SNP genotype data. Six hundred and seventy-three steers were genotyped using an Illumina Bovine SNP50 BeadChip and phenotyped for cold carcass weight, ribeye area, rib thickness, subcutaneous fat thickness, estimated yield percent and marbling score. Additive polygenic variance and the variance attributable to a set of SNPs that had statistically significant effects on the trait were estimated via Gibbs sampling with two models: (i) a model with the chosen SNPs and the additive polygenic effects; and (ii) a model with the polygenic effects alone. The proportion of the estimated variance attributable to the SNPs became higher as the number of SNP effects that fit increased. High correlations between breeding values estimated with the model containing the polygenic effect alone and those estimated by chosen SNPs were obtained. No fraction of the total genetic variance was explained by SNPs associated with the trait at P ≥ 0.1. Our results suggest that for the carcass traits of Japanese Black cattle, a maximum of half of the total additive genetic variance may be explained by SNPs between 100 several tens to several 100s.

Identification of novel genetic markers and evaluation of genetic structure in a population of Japanese crested ibis.

Japanese population of the Japanese crested ibis Nipponia nippon was founded by five individuals gifted from the People's Republic of China. In order to exactly evaluate genetic structure, we first performed development of novel genetic makers using 89 microsatellite primer pairs of related species for cross-amplification. Of these, only three primer pairs were useful for the genetic markers. Additionally, we sequenced allelic PCR products of these three markers together with 10 markers previously identified. Most markers showed typical microsatellite repeat units, but two markers were not simple microsatellites. Moreover, over half of the markers did not have the same repeat units as those of the original species. These results suggested that development of novel genetic markers in this population by cross-amplification is not efficient, partly because of low genetic diversity. Furthermore, the cluster analysis by STRUCTURE program using 17 markers showed that the five founders were divided into two clusters. However, the genetic relationships among the founders indicated by the clustering seemed to be questionable, because the analysis relied largely on a small number of triallelic markers, in spite of the addition of the three useful markers. Therefore, more efficient methods for identifying large numbers of single nucleotide polymorphisms are desirable.

Genome-wide SNP and STR discovery in the Japanese crested ibis and genetic diversity among founders of the Japanese population.

The Japanese crested ibis is an internationally conserved, critically threatened bird. Captive-breeding programs have been established to conserve this species in Japan. Since the current Japanese population of crested ibis originates only from 5 founders donated by the Chinese government, understanding the genetic diversity between them is critical for an effective population management. To discover genome-wide single nucleotide polymorphisms (SNPs) and short tandem repeats (STRs) while obtaining genotype data of these polymorphic markers in each founder, reduced representation libraries were independently prepared from each of the founder genomes and sequenced on an Illumina HiSeq2000. This yielded 316 million 101-bp reads. Consensus sequences were created by clustering sequence reads, and then sequence reads from each founder were mapped to the consensus sequences, resulting in the detection of 52,512 putative SNPs and 162 putative STRs. The numbers of haplotypes and STR alleles and the investigation of genetic similarities suggested that the total genetic diversity between the founders was lower, although we could not identify a pair with closely related genome sequences. This study provided important insight into protocols for genetic management of the captive breeding population of Japanese crested ibis in Japan and towards the national project for reintroduction of captive-bred individuals into the wild. We proposed a simple, efficient, and cost-effective approach for simultaneous detection of genome-wide polymorphic markers and their genotypes for species currently lacking a reference genome sequence.

Genetic diversity and structure in the Sado captive population of the Japanese crested ibis.

The Japanese crested ibis Nipponia nippon is a critically threatened bird. We assessed genetic diversity and structure in the Sado captive population of the Japanese crested ibis based on 24 and 50 microsatellite markers developed respectively for the same and related species. Of a total of 74 loci, 19 showed polymorphisms in the five founder birds of the population, and therefore were useful for the analysis of genetic diversity and structure. Genetic diversity measures, A, ne, He, Hoand PIC, obtained by genotyping of the 138 descendants were similar to those of other species with population bottlenecks, and thus considerably low. The low level of genetic diversity resulting from such bottlenecks was consistent with the results of lower genetic diversity measures for the Sado captive relative to the Chinese population that is the source population for the Sado group as determined using previously reported data and heterozygosity excess by Hardy-Weinberg equilibrium tests. Further, individual clustering based on the allele-sharing distance and Bayesian model-based clustering revealed that the founder genomes were equally at population in total, and with various admixture patterns at individual levels inherited by the descendants. The clustering results, together with the result of inheritance of all alleles of the microsatellites from the founders to descendants, suggest that planned mating in captive-breeding programs for the population has succeeded in maintaining genetic diversity and minimizing kinship. In addition, the Bayesian model-based clustering assumed two different components of genomes in the Sado captive Japanese crested ibis, supporting a considerably low level of genetic diversity.

The W- and Z-linked EE0.6 sequences used for molecular sexing of captive Japanese crested ibis on Sado Island.

The Japanese crested ibis Nipponia nippon is a critically threatened bird. Accurate sexing is necessary to perform effective management of captive breeding toward a national project for a tentative release of the Japanese crested ibis on Sado Island. A PCR-based sexing method targeting a 0.6 kb EcoRI fragment (EE0.6) sequence on W chromosome with AWS03 and USP3 primers has been developed for the Japanese crested ibis. However, the primers were selected from the EE0.6 sequences from bird species other than the Japanese crested ibis. In this study, we determined the W- and Z-linked EE0.6 sequences in the Japanese crested ibis, and clarified Japanese crested ibis sequence mismatch in the binding sites of the primers. Further, we found no polymorphism in the primer binding sites among five founder birds for the Sado captive Japanese crested ibis population. These findings validated the PCR-based sexing method with the AWS03 and USP3 as accurate molecular sexing methods of captive Japanese crested ibis on the Sado Island. Additionally, we designed a primer set for a novel PCR-based sexing, based on the EE0.6 sequences obtained in this study. This novel sexing method may be useful for future ecological research following the release of Japanese crested ibis on Sado Island. This is the first report to show the EE0.6 sequences in Japanese crested ibis.

QTL/microarray approach using pathway information.

BACKGROUND: A combined quantitative trait loci (QTL) and microarray-based approach is commonly used to find differentially expressed genes which are then identified based on the known function of a gene in the biological process governing the trait of interest. However, a low cutoff value in individual gene analyses may result in many genes with moderate but meaningful changes in expression being missed. RESULTS: We modified a gene set analysis to identify intersection sets with significantly affected expression for which the changes in the individual gene sets are less significant. The gene expression profiles in liver tissues of four strains of mice from publicly available microarray sources were analyzed to detect trait-associated pathways using information on the QTL regions of blood concentrations of high density lipoproteins (HDL) cholesterol and insulin-like growth factor 1 (IGF-1). Several metabolic pathways related to HDL levels, including lipid metabolism, ABC transporters and cytochrome P450 pathways were detected for HDL QTL regions. Most of the pathways identified for the IGF-1 phenotype were signal transduction pathways associated with biological processes for IGF-1's regulation. CONCLUSION: We have developed a method of identifying pathways associated with a quantitative trait using information on QTL. Our approach provides insights into genotype-phenotype relations at the level of biological pathways which may help to elucidate the genetic architecture underlying variation in phenotypic traits.

Non-invasive sampling technique for DNA extraction from captive Japanese Crested Ibis on Sado Island.

The Japanese Crested Ibis Nipponia nippon is a critically threatened bird. The post-hatch eggs of the current captive population of this species on Sado Island have been stored at room temperature for the long-term. In this study, we investigated the suitability of the vascularized chorioallantois membrane from the eggs as a non-invasive DNA source. Using microsatellite loci developed for the Japanese Crested Ibis, we performed three experiments for comparison of genotypes obtained among DNA. First, DNA from five different sites of the identical membrane showed the same genotypes at either of two loci examined. Second, DNA from the membrane of each full-sibling birds and blood of their parents showed the genotypes that were consistent with Mendelian parent-offspring relationships at any of eight loci examined. Third, DNA from the membrane and blood of the same bird showed the matched genotypes at any of eight loci examined. These results indicate that the vascularized chorioallantois membrane from post-hatch eggs stored at room temperature for the long- term can be used as a reliable DNA source of offspring that had hatched from the egg. This study will promote a molecular genetics study on genetic diversity of the current captive Japanese Crested Ibis population on Sado Island.

Investigation of Gibbs sampling conditions to estimate variance components from Japanese Black carcass field data.

The genetic evaluation using the carcass field data in Japanese Black cattle has been carried out employing an animal model, implementing the restricted maximum likelihood (REML) estimation of additive genetic and residual variances. Because of rapidly increasing volumes of the official data sets and therefore larger memory spaces required, an alternative approach like the REML estimation could be useful. The purpose of this study was to investigate Gibbs sampling conditions for the single-trait variance component estimation using the carcass field data. As prior distributions, uniform and normal distributions and independent scaled inverted chi-square distributions were used for macro-environmental effects, breeding values, and the variance components, respectively. Using the data sets of different sizes, the influences of Gibbs chain length and thinning interval were investigated, after the burn-in period was determined using the coupling method. As would be expected, the chain lengths had obviously larger effects on the posterior means than those of thinning intervals. The posterior means calculated using every 10th sample from 90,000 of samples after 10,000 samples discarded as burn-in period were all considered to be reasonably comparable to the corresponding estimates by REML.

Estimation of breeding values from large-sized routine carcass data in Japanese Black cattle using Bayesian analysis.

Volumes of official data sets have been increasing rapidly in the genetic evaluation using the Japanese Black routine carcass field data. Therefore, an alternative approach with smaller memory requirement to the current one using the restricted maximum likelihood (REML) and the empirical best linear unbiased prediction (EBLUP) is desired. This study applied a Bayesian analysis using Gibbs sampling (GS) to a large data set of the routine carcass field data and practically verified its validity in the estimation of breeding values. A Bayesian analysis like REML-EBLUP was implemented, and the posterior means were calculated using every 10th sample from 90,000 of samples after 10,000 samples discarded. Moment and rank correlations between breeding values estimated by GS and REML-EBLUP were very close to one, and the linear regression coefficients and the intercepts of the GS on the REML-EBLUP estimates were substantially one and zero, respectively, showing a very good agreement between breeding value estimation by the current GS and the REML-EBLUP. The current GS required only one-sixth of the memory space with REML-EBLUP. It is confirmed that the current GS approach with relatively small memory requirement is valid as a genetic evaluation procedure using large routine carcass data.