Consecutive Single Nucleotide Polymorphisms Research Articles

Single-step genome-wide association analyses for selected infrared-predicted cheese-making traits in Walloon Holstein cows

This study aimed to perform genome-wide association study to identify genomic regions associated with milk production and cheese-making properties (CMP) in Walloon Holstein cows. The studied traits were milk yield, fat percentage, protein percentage, casein percentage (CNP), calcium content, somatic cell score (SCS), coagulation time, curd firmness after 30 min from rennet addition, and titratable acidity. The used data have been collected from 2014 to 2020 on 78,073 first-parity (485,218 test-day records), 48,766 second-parity (284,942 test-day records), and 21,948 third-parity (105,112 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA) of 6,617 animals (1,712 males), were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 KB) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for positional candidate genes. Heritability estimates for the studied traits ranged from 0.10 (SCS) to 0.53 (CNP), 0.10 (SCS) to 0.50 (CNP), and 0.12 (SCS) to 0.49 (CNP) in the first, second, and third parity, respectively. Genome-wide association analyses identified 6 genomic regions (BTA1, BTA14 [4 regions], and BTA20) associated with the considered traits. Genes including the SLC37A1 (BTA1), SHARPIN, MROH1, DGAT1, FAM83H, TIGD5, MROH6, NAPRT, ADGRB1, GML, LYPD2, JRK (BTA14), and TRIO (BTA20) were identified as positional candidate genes for the studied CMP. The findings of this study help to unravel the genomic background of a cow's ability for cheese production and can be used for the future implementation and use of genomic evaluation to improve the cheese-making traits in Walloon Holstein cows.

Genome-wide association for milk urea concentration in Dual-Purpose Belgian Blue cows.

The objectives of this study were to estimate genetic parameters and identify genomic regions associated with milk urea concentration (MU) in Dual-Purpose Belgian Blue (DPBB) cows. The data were 29,693 test-day records of milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FP), protein percentage (PP) and MU collected between 2014 and 2020 on 2498 first parity cows (16,935 test-day records) and 1939 second-parity cows (12,758 test-day records) from 49 herds in the Walloon Region of Belgium. Data of 28,266 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), on 1699 animals (639 males and 1060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method using a single chain of 100,000 iterations after a burn-in period of 20,000. SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by windows of 25 consecutive SNPs (with an average size of ~2Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. The mean (SD) of MU was 22.89 (10.07) and 22.35 (10.07) mg/dl for first and second parity, respectively. The mean (SD) heritability estimates for daily MU were 0.18 (0.01) and 0.22 (0.02), for first and second parity, respectively. The mean (SD) genetic correlations between daily MU and MY, FY, PY, FP and PP were -0.05 (0.09), -0.07 (0.11), -0.03 (0.13), -0.05 (0.08) and -0.03 (0.11) for first parity, respectively. The corresponding values estimated for second parity were 0.02 (0.10), -0.02 (0.09), 0.02 (0.08), -0.08 (0.06) and -0.05 (0.05). The genome-wide association analyses identified three genomic regions (BTA2, BTA3 and BTA13) associated with MU. The identified regions showed contrasting results between parities and among different stages within each parity. This suggests that different groups of candidate genes underlie the phenotypic expression of MU between parities and among different lactation stages within a parity. The results of this study can be used for future implementation and use of genomic evaluation to reduce MU in DPBB cows.

Genomic regions and candidate genes linked with Phytophthora capsici root rot resistance in chile pepper (Capsicum annuum L.)

BackgroundPhytophthora root rot, caused by Phytophthora capsici, is a major disease affecting Capsicum production worldwide. A recombinant inbred line (RIL) population derived from the hybridization between ‘Criollo de Morellos-334’ (CM-334), a resistant landrace from Mexico, and ‘Early Jalapeno’, a susceptible cultivar was genotyped using genotyping-by-sequencing (GBS)-derived single nucleotide polymorphism (SNP) markers. A GBS-SNP based genetic linkage map for the RIL population was constructed. Quantitative trait loci (QTL) mapping dissected the genetic architecture of P. capsici resistance and candidate genes linked to resistance for this important disease were identified.ResultsDevelopment of a genetic linkage map using 1,973 GBS-derived polymorphic SNP markers identified 12 linkage groups corresponding to the 12 chromosomes of chile pepper, with a total length of 1,277.7 cM and a marker density of 1.5 SNP/cM. The maximum gaps between consecutive SNP markers ranged between 1.9 (LG7) and 13.5 cM (LG5). Collinearity between genetic and physical positions of markers reached a maximum of 0.92 for LG8. QTL mapping identified genomic regions associated with P. capsici resistance in chromosomes P5, P8, and P9 that explained between 19.7 and 30.4% of phenotypic variation for resistance. Additive interactions between QTL in chromosomes P5 and P8 were observed. The role of chromosome P5 as major genomic region containing P. capsici resistance QTL was established. Through candidate gene analysis, biological functions associated with response to pathogen infections, regulation of cyclin-dependent protein serine/threonine kinase activity, and epigenetic mechanisms such as DNA methylation were identified.ConclusionsResults support the genetic complexity of the P. capsici–Capsicum pathosystem and the possible role of epigenetics in conferring resistance to Phytophthora root rot. Significant genomic regions and candidate genes associated with disease response and gene regulatory activity were identified which allows for a deeper understanding of the genomic landscape of Phytophthora root rot resistance in chile pepper.

Sporadic Pseudohypoparathyroidism Type 1B in Monozygotic Twins: Insights Into the Pathogenesis of Methylation Defects

ContextSporadic pseudohypoparathyroidism type 1B (sporPHP1B) is an imprinting disease without a defined genetic cause, characterized by broad methylation changes in differentially methylated regions (DMRs) of the GNAS gene.ObjectiveThis work aims to provide insights into the causative event leading to the GNAS methylation defects through comprehensive molecular genetic analyses of a pair of female monozygotic twins concordant for sporPHP1B who were conceived naturally, that is, without assisted reproductive techniques.MethodsUsing the leukocyte genome of the twins and family members, we performed targeted bisulfite sequencing, methylation-sensitive restriction enzyme (MSRE)-quantitative polymerase chain reaction (qPCR), whole-genome sequencing (WGS), high-density single-nucleotide polymorphism (SNP) array, and Sanger sequencing.ResultsMethylation analyses by targeted bisulfite sequencing and MSRE-qPCR revealed almost complete losses of methylation at the GNAS AS, XL, and A/B DMRs and a gain of methylation at the NESP55 DMR in the twins, but not in other family members. Except for the GNAS locus, we did not find apparent methylation defects at other imprinted genome loci of the twins. WGS, SNP array, and Sanger sequencing did not detect the previously described genetic defects associated with familial PHP1B. Sanger sequencing also ruled out any novel genetic alterations in the entire NESP55/AS region. However, the analysis of 28 consecutive SNPs could not exclude the possibility of paternal heterodisomy in a span of 22 kb comprising exon NESP55 and AS exon 5.ConclusionOur comprehensive analysis of a pair of monozygotic twins with sporPHP1B ruled out all previously described genetic causes. Twin concordance indicates that the causative event was an imprinting error earlier than the timing of monozygotic twinning.

Genome-Wide Association Study Provides Insights into Important Genes for Reproductive Traits in Nelore Cattle.

Simple SummaryIn this study, we investigated the association between single nucleotide polymorphisms (SNPs) and reproductive traits in order to identify candidate genes and biological pathways associated with these traits in Nelore beef cattle. The genome-wide association analysis revealed genomic regions that could explain part of the genetic variance of the studied traits. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, being important for reproductive efficiency. The identification of candidate genes that were associated with the studied traits as well as genes enriched in the functional terms and pathways may be useful for exploring the genetic architecture underlying reproductive traits and may be used in Nelore breeding programs.The identification of genomic regions associated with reproductive traits as well as their biological processes allows a better understanding of the phenotypic variability of these traits. This information could be applied to animal breeding programs to accelerate genetic gain. The aim of this study was to evaluate the association between single nucleotide polymorphisms (SNP) with a scrotal circumference at 365 days of age (SC365) and at 450 days of age (SC450), gestation length (GL) as a calf trait, age at first calving (AFC), accumulated productivity (ACP), heifer early calving until 30 months (HC30), and stayability (STAY) traits, in order to identify candidate genes and biological pathways associated with reproductive traits in Nelore cattle. The data set consisted of pedigree, phenotypes, and genotypes of Nelore cattle from the “Associação Nacional de Criadores e Pesquisadores” (ANCP). The association analyses were performed using the Weighted Single-Step Genome-Wide Association method; the regions, consisting of 10 consecutive SNP, which explained more than 0.5% of additive genetic variance, were considered as a significant association. A total of 3, 6, 7, 5, 10, 25, and 12 windows were associated with SC355, SC450, GL, AFC, ACP, HC30, and STAY, respectively. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, among them CAMK1D, TASP1, ACOXL, RAB11FIP5, and SFXN5. Moreover, the genes were enriched in functional terms, like negative regulation of fat cell differentiation, fatty acid alpha-oxidation, and sphingolipids signaling pathway. The identification of the genes associated with the traits, as well as genes enriched in the terms and pathway mentioned above, should contribute to future biological validation studies and may be used as candidate genes in Nelore breeding programs.

Genomic Heritability and Genome-Wide Association Studies of Plasma Metabolites in Crossbred Beef Cattle.

Metabolites, substrates or products of metabolic processes, are involved in many biological functions, such as energy metabolism, signaling, stimulatory and inhibitory effects on enzymes and immunological defense. Metabolomic phenotypes are influenced by combination of genetic and environmental effects allowing for metabolome-genome-wide association studies (mGWAS) as a powerful tool to investigate the relationship between these phenotypes and genetic variants. The objectives of this study were to estimate genomic heritability and perform mGWAS and in silico functional enrichment analyses for a set of plasma metabolites in Canadian crossbred beef cattle. Thirty-three plasma metabolites and 45,266 single nucleotide polymorphisms (SNPs) were available for 475 animals. Genomic heritability for all metabolites was estimated using genomic best linear unbiased prediction (GBLUP) including genomic breed composition as covariates in the model. A single-step GBLUP implemented in BLUPF90 programs was used to determine SNP P values and the proportion of genetic variance explained by SNP windows containing 10 consecutive SNPs. The top 10 SNP windows that explained the largest genetic variation for each metabolite were identified and mapped to detect corresponding candidate genes. Functional enrichment analyses were performed on metabolites and their candidate genes using the Ingenuity Pathway Analysis software. Eleven metabolites showed low to moderate heritability that ranged from 0.09 ± 0.15 to 0.36 ± 0.15, while heritability estimates for 22 metabolites were zero or negligible. This result indicates that while variations in 11 metabolites were due to genetic variants, the majority are largely influenced by environment. Three significant SNP associations were detected for betaine (rs109862186), L-alanine (rs81117935), and L-lactic acid (rs42009425) based on Bonferroni correction for multiple testing (family wise error rate <0.05). The SNP rs81117935 was found to be located within the Catenin Alpha 2 gene (CTNNA2) showing a possible association with the regulation of L-alanine concentration. Other candidate genes were identified based on additive genetic variance explained by SNP windows of 10 consecutive SNPs. The observed heritability estimates and the candidate genes and networks identified in this study will serve as baseline information for research into the utilization of plasma metabolites for genetic improvement of crossbred beef cattle.

A Genome-Wide Association Study for Calving Interval in Holstein Dairy Cows Using Weighted Single-Step Genomic BLUP Approach.

Simple SummaryReproductive performance is an important factor, which determines productive life and drives culling decisions in dairy herds. There are strong motives for including reproductive performance in genetic selection programs of dairy cows; however, low heritability estimates reported for reproductive performance measures limit the genetic selection efficiency. More effective genetic selection could be achieved using genomic information. The aim of this study was to identify genomic region(s) associated with the length of the calving interval in Holstein cows. The accuracies of genomic estimated breeding values (GEBVs) with single-step genomic BLUP (ssGBLUP) and the pedigree-based BLUP were compared as well. The results showed that the accuracies of GEBVs using the single-step genomic BLUP were much higher than those estimated using the pedigree-based BLUP. We identified three genomic regions (BTA3, BTA6, and BTA7) associated with the length of the calving interval in Holstein dairy cows.The aim of the present study was to identify genomic region(s) associated with the length of the calving interval in primiparous (n = 6866) and multiparous (n = 5071) Holstein cows. The single nucleotide polymorphism (SNP) solutions were estimated using a weighted single-step genomic best linear unbiased prediction (WssGBLUP) approach and imputed high-density panel (777 k) genotypes. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. The results showed that the accuracies of GEBVs with WssGBLUP improved by +5.4 to +5.7, (primiparous cows) and +9.4 to +9.7 (multiparous cows) percent points over accuracies from the pedigree-based BLUP. The most accurate genomic evaluation was provided at the second iteration of WssGBLUP, which was used to identify associated genomic regions using a windows-based GWAS procedure. The proportion of additive genetic variance explained by windows of 50 consecutive SNPs (with an average of 165 Kb) was calculated and the region(s) that accounted for equal to or more than 0.20% of the total additive genetic variance were used to search for candidate genes. Three windows of 50 consecutive SNPs (BTA3, BTA6, and BTA7) were identified to be associated with the length of the calving interval in primi- and multiparous cows, while the window with the highest percentage of explained genetic variance was located on BTA3 position 49.42 to 49.52 Mb. There were five genes including ARHGAP29, SEC24D, METTL14, SLC36A2, and SLC36A3 inside the windows associated with the length of the calving interval. The biological process terms including alanine transport, L-alanine transport, proline transport, and glycine transport were identified as the most important terms enriched by the genes inside the identified windows.

A detailed characteristics of bias associated with long runs of homozygosity identification based on medium density SNP microarrays.

In the present study, runs of homozygosity (ROH) detected with the use of a standard bovine 54k single nucleotide polymorphism (SNP) genotyping assay and two different ROH detection approaches, based on 50 (M1) or 15 (M2) consecutive SNPs, were compared with results of whole genome sequencing. Both microarray-based methods accurately recognised medium-sized ROH, however, it was found that M2 method seemed to better than M1 identify short ROH, but highly overestimated their number, leading to numerous false positive calls. Moreover, long ROH identified with microarray data tended to break into shorter segments in sequencing data because of the presence of regions with high heterozygosity within the ROH sequences. This may indicate, that these long ROH are formed by closely positioned shorter homozygous segments that may be of older origin or may be created by two similar but not identical haplotypes, showing minor internal recombination signs. Such finding also suggests that at least some of the results of previous studies in regard to long ROH may be biased leading to inaccurate estimations of genomes autozygosity via ROH classification into length categories.

Selecting Closely-Linked SNPs Based on Local Epistatic Effects for Haplotype Construction Improves Power of Association Mapping.

Genome-wide association studies (GWAS) have gained central importance for the identification of candidate loci underlying complex traits. Single nucleotide polymorphism (SNP) markers are mostly used as genetic variants for the analysis of genotype-phenotype associations in populations, but closely linked SNPs that are grouped into haplotypes are also exploited. The benefit of haplotype-based GWAS approaches vs. SNP-based approaches is still under debate because SNPs in high linkage disequilibrium provide redundant information. To overcome some constraints of the commonly-used haplotype-based GWAS in which only consecutive SNPs are considered for haplotype construction, we propose a new method called functional haplotype-based GWAS (FH GWAS). FH GWAS is featured by combining SNPs into haplotypes based on the additive and epistatic effects among SNPs. Such haplotypes were termed functional haplotypes (FH). As shown by simulation studies, the FH GWAS approach clearly outperformed the SNP-based approach unless the minor allele frequency of the SNPs making up the haplotypes is low and the linkage disequilibrium between them is high. Applying FH GWAS for the trait flowering time in a large Arabidopsis thaliana population with whole-genome sequencing data revealed its potential empirically. FH GWAS identified all candidate regions which were detected in SNP-based and two other haplotype-based GWAS approaches. In addition, a novel region on chromosome 4 was solely detected by FH GWAS. Thus both the results of our simulation and empirical studies demonstrate that FH GWAS is a promising method and superior to the SNP-based approach even if almost complete genotype information is available.

Using Brachypodium distachyon natural populations to uncover genomic regions under selection.

Using Brachypodium distachyon natural populations to uncover genomic regions under selection.

Genome-wide association studies using a penalized moving-window regression.

Genome-wide association studies (GWAS) have played an important role in identifying genetic variants underlying human complex traits. However, its success is hindered by weak effect at causal variants and presence of noise at non-causal variants. In an effort to overcome these difficulties, a previous study proposed a regularized regression method that penalizes on the difference of signal strength between two consecutive single-nucleotide polymorphisms (SNPs). We provide a generalization to the afore-mentioned method so that more adjacent SNPs can be incorporated. The choice of optimal number of SNPs is studied. Simulation studies indicate that when consecutive SNPs have similar absolute coefficients our method performs better than using LASSO penalty. In other situations, our method is still comparable to using LASSO penalty. The practical utility of the proposed method is demonstrated by applying it to Genetic Analysis Workshop 16 rheumatoid arthritis GWAS data. An implementation of the proposed method is provided in R package MWLasso. kai-wang@uiowa.edu.

Prenatal and postnatal findings in small-for-gestational-age fetuses without structural ultrasound anomalies at 18-24 weeks.

To assess phenotypic and genotypic characteristics of small-for-gestational-age (SGA) fetuses without structural anomalies at 18-24 weeks' gestation. This retrospective study included structurally normal singleton fetuses with an abdominal circumference ≤ 5th percentile on detailed ultrasound examination between 18 and 24 weeks' gestation. Cases were stratified according to the absence or presence of other abnormal ultrasound findings, such as abnormal amniotic fluid or soft markers. All patients were offered invasive prenatal testing with rapid aneuploidy detection by qualitative fluorescence polymerase chain reaction (QF-PCR) and, if normal, consecutive single nucleotide polymorphism (SNP) array was also offered. Detailed postnatal follow-up (≥ 5 months) was performed. In cases in which a syndromic phenotype became apparent within 5 months after birth and SNP array had not been performed prenatally, it was performed postnatally. A total of 211 pregnancies were eligible for inclusion. Of the 158 cases with isolated SGA on ultrasound, 36 opted for invasive prenatal testing. One case of trisomy 21 and one case of a submicroscopic abnormality (a susceptibility locus for neurodevelopmental disease) were detected. Postnatal follow-up showed a postnatal apparent syndromic phenotype in 10 cases. In one case this was due to trisomy 21 and the other nine (5.8%; 95% CI, 2.8-10.0%) cases had normal SNP array results. In 32/53 cases with SGA and associated ultrasound abnormalities, parents opted for invasive testing. One case of trisomy 21 and one of triploidy were found. In 11 cases a syndromic phenotype became apparent after birth. One was due to trisomy 21 and in one case a submicroscopic anomaly (a susceptibility locus) was found. The remaining syndromic cases (17.3%; 95% CI, 8.7-29.0%) had normal SNP array results. Testing for chromosomal anomalies should be offered in cases of SGA between 18 and 24 weeks' gestation. Whole chromosome anomalies occur in 1.3% (95% CI, 0.2-3.9%) of isolated SGA and 5.8% (95% CI, 1.5-14.0%) of associated SGA. In 0.6% (95% CI, 0.1-2.8%) and 1.9% (95% CI, 0.2-8.2%), respectively, SNP array detected a susceptibility locus for neurodevelopmental disease that would not be detected by karyotyping, QF-PCR or non-invasive prenatal testing. Therefore, and because the genetic causes of SGA are diverse, we suggest SNP array testing in cases of SGA. Thorough postnatal examination and follow-up of infants that presented with reduced fetal growth is important because chromosomally normal syndromic phenotypes occur frequently in SGA fetuses. Copyright © 2016 ISUOG. Published by John Wiley & Sons Ltd.

In silico screening of the chicken genome for overlaps between genomic regions: microRNA genes, coding and non-coding transcriptional units, QTL, and genetic variations.

MicroRNAs (miRNAs) are a class of non-coding RNAs involved in posttranscriptional regulation of target genes. Regulation requires complementarity between target mRNA and the mature miRNA seed region, responsible for their recognition and binding. It has been estimated that each miRNA targets approximately 200 genes, and genetic variability of miRNA genes has been reported to affect phenotypic variability and disease susceptibility in humans, livestock species, and model organisms. Polymorphisms in miRNA genes could therefore represent biomarkers for phenotypic traits in livestock animals. In our previous study, we collected polymorphisms within miRNA genes in chicken. In the present study, we identified miRNA-related genomic overlaps to prioritize genomic regions of interest for further functional studies and biomarker discovery. Overlapping genomic regions in chicken were analyzed using the following bioinformatics tools and databases: miRNA SNiPer, Ensembl, miRBase, NCBI Blast, and QTLdb. Out of 740 known pre-miRNA genes, 263 (35.5%) contain polymorphisms; among them, 35 contain more than three polymorphisms The most polymorphic miRNA genes in chicken are gga-miR-6662, containing 23 single nucleotide polymorphisms (SNPs) within the pre-miRNA region, including five consecutive SNPs, and gga-miR-6688, containing ten polymorphisms including three consecutive polymorphisms. Several miRNA-related genomic hotspots have been revealed in chicken genome; polymorphic miRNA genes are located within protein-coding and/or non-coding transcription units and quantitative trait loci (QTL) associated with production traits. The present study includes the first description of an exonic miRNA in a chicken genome, an overlap between the miRNA gene and the exon of the protein-coding gene (gga-miR-6578/HADHB), and the first report of a missense polymorphism located within a mature miRNA seed region. Identified miRNA-related genomic hotspots in chicken can serve researchers as a starting point for further functional studies and association studies with poultry production and health traits and the basis for systematic screening of exonic miRNAs and missense/miRNA seed polymorphisms in other genomes.

The JR blood group system: genetic and molecular investigations

Background Jra is a high‐incidence erythrocyte antigen previously designated ISBT 901.005. The rare Jr(a−) phenotype appeared to be inherited as an autosomal recessive trait, but the gene controlling the expression of Jra was not known. We conducted genetic and molecular investigations in an attempt to resolve this issue.Aim This lecture and accompanying manuscript will detail the homozygosity by descent (HBD) gene mapping strategy we used to establish that ABCG2 null alleles define the Jr(a−) phenotype (Zelinski et al., Nat Genet 2012;44:131–132).Materials and Methods Genomic DNA was genotyped for single nucleotide polymorphisms (SNPs) using a Human GeneChip (Affymetrix, 250K NspI) Mapping array. The six Jr (a−) study subjects were two Caucasian sisters, two Asian siblings and two unrelated Asian individuals referred for serologic testing. Anti‐Jra was identified in the plasma from four of the six. Intronic primers flanking exons 2–16 of ABCG2 were designed and used for PCR amplification. Both the forward and reverse strands of each purified product were Sanger sequenced at The Centre for Applied Genomics (Toronto, Canada).Results Analysis of the SNP genotype data revealed that the Caucasian sisters were identically homozygous for large tracts of consecutive SNPs in seven different chromosomal regions; chromosome 1q (3·3 Mb), 3q (8·8 Mb), 4q (28·9 Mb), 7q (10·5 and 7·7 Mb), 9q (3·1 Mb) and 10q (2·0 Mb). DNA from the four Asian individuals did not consistently contain homozygous SNP haplotypes in any of these regions except for chromosome 4q22.1. In this region, the siblings and one unrelated Asian had blocks of homozygosity that overlapped for approximately 400 000 bp and were identical to each other, but, different from the Caucasians. The remaining unrelated Asian was homozygous (approximately 520 000 bp) for a third SNP haplotype. The 4q22.1 region of homozygosity identified in all six study subjects contains four validated genes, only one of which, ABCG2, is known to be highly expressed on erythrocytes. DNA sequencing revealed that homozygosity for each of the three SNP haplotypes corresponded with homozygosity for a different nonsense (c.736C&gt;T or c.376C&gt;T) or missense (c.34G&gt;A) mutation in ABCG2.Conclusion The genetic and molecular evidence revealed by this study proves that ABCG2 null alleles define the Jr(a−) blood group phenotype. Independently, our French colleagues established that Jra was carried on the ABCG2 protein (Saison et al., Nat Genet 2012;44:174–177). Based on these results the JR blood group system (ISBT 032) has been established.

Molecular Analysis of Mutant Alleles for Elevated Palmitate Concentration in Soybean

ABSTRACTSoybean [Glycine max (L.) Merr.] oil with increased palmitate concentration may be useful for some food and industrial applications. Chemical mutagenesis was used to develop mutant alleles for elevated palmitate concentration that were designated as fap2(A21), fap4(A24), fap5(A27), fap6(A25), and fap7(A30) based on classical genetic analysis. The objective of our study was to determine whether the elevated palmitate phenotypes in any of these lines were associated with mutations in either of the two known 3‐ketoacyl‐ACP synthase II (KAS II) genes of soybean, GmKAS IIA and GmKAS IIB. Deoxyribonucleic acid sequence analysis revealed single nucleotide polymorphisms (SNPs) that differentiated mutant and wild‐type alleles in one of the GmKAS genes in lines A21, A25, A27, and A30 but not in A24. The identified SNPs resulted in nonconservative amino acid substitutions that were predicted to impair enzyme function. The fap2(A21) allele was associated with two consecutive SNPs within the GmKAS IIA gene, confirming its allelic relationship to fap2(C1727) that was known to have a mutant SNP within the same gene. The allele designated fap5(A27) also had a SNP within the GmKAS IIA gene, indicating that its designation should be changed to fap2(A27). The fap6(A25) and fap7(A30) alleles each contained a single SNP at different locations within the GmKAS IIB gene. We propose that the designation of the mutant allele in A25 remain fap6(A25) and that the allele in A30 be changed from fap7(A30) to fap6(A30).

Identification of divergently selected regions between Japanese Black and Holstein cattle using bovine 50k SNP array

Differences between average allelic frequencies of genes that relate to traits suggest that it would be evidence of artificial selections. Sliding window approach is a useful method to identify genomic regions that have been differently selected between two breeds. The objective of this study was to identify the divergently selected regions between Japanese Black (JB) and Japanese Holstein (JH) cattle based on genotypic information obtained through a high-density single nucleotide polymorphism (SNP) panel. After genotyping of 54 001 SNP markers on 100 animals (50 JB and 50 JH), 40 635 SNPs were suitable for the analysis. For each of these SNPs, the absolute difference between allelic frequencies of JB and JH was calculated. In the current study, 10 consecutive SNPs were defined as components of a window. For each window, the average difference in allelic frequency was calculated. This was termed sliding window average difference (SWAD). Among 40 055 windows, we focused on 39 windows with the largest SWAD. This was equivalent to 0.1% of all windows and the SWAD was more than 0.435. Some of these windows overlapped and were distributed in 11 regions. These regions were in good agreement with reported quantitative trait locus, therefore would be selection signatures and good candidates that harbor the causative mutations.

Genomic breeding value prediction and QTL mapping of QTLMAS2010 data using Bayesian Methods

BackgroundBayesian methods allow prediction of genomic breeding values (GEBVs) using high-density single nucleotide polymorphisms (SNPs) covering the whole genome with effective shrinkage of SNP effects using appropriate priors. In this study we applied a modification of the well-known BayesA and BayesB methods to estimate the proportion of SNPs with zero effects (π) and a common variance for non-zero effects. The method, termed BayesCπ, was used to predict the GEBVs of the last generation of the QTLMAS2010 data. The accuracy of GEBVs from various methods was estimated by the correlation with phenotypes in the last generation. The methods were BayesCPi and BayesB with different π values, both with and without polygenic effects, and best linear unbiased prediction using an animal model with a genomic or numerator relationship matrix. Positions of quantitative trait loci (QTLs) were identified based on the variances of GEBVs for windows of 10 consecutive SNPs. We also proposed a novel approach to set significance thresholds for claiming QTL in this specific case by using pedigree-based simulation of genotypes. All analyses were focused on detecting and evaluating QTL with additive effects.ResultsThe accuracy of GEBVs was highest for BayesCπ, but the accuracy of BayesB with π equal to 0.99 was similar to that of BayesCπ. The accuracy of BayesB dropped with a decrease in π. Including polygenic effects into the model only had marginal effects on accuracy and bias of predictions. The number of QTL identified was 15 when based on a stringent 10% chromosome-wise threshold and increased to 21 when a 20% chromosome-wise threshold was used.ConclusionsThe BayesCπ method without polygenic effects was identified to be the best method for the QTLMAS2010 dataset, because it had highest accuracy and least bias. The significance criterion based on variance of 10-SNP windows allowed detection of more than half of the QTL, with few false positives.

Finding the mutation responsible for inherited rickets in Corriedale sheep

Inherited rickets of Corriedale sheep is of simple autosomal recessive inheritance. The gross and histological lesions include segmental thickening of physes, persistence of hypertrophic chondrocytes at sites of endochondral ossification, excessive osteoclastic resorption, microfractures and wide unmineralised osteoid seams lining trabeculae and filling secondary osteons. Initial testing revealed concentrations of 25-hydroxyvitamin D3 were normal, but those of 1,25-dihydroxyvitamin D3 were increased in serum of affected sheep, suggesting a defect in the vitamin D receptor. However, in vitro studies on cultured skin fibroblasts revealed normal binding and function of the vitamin D receptor. A genome-wide association study was conducted using the Illumina OvineSNP50 BeadChip. A homozygous region of 199 consecutive single-nucleotide polymorphism loci was identified in affected sheep, covering a region of 10 megabase pairs on ovine chromosome 6. Of the 91 candidate genes in this region, dentin matrix protein 1 gene (dmp1) was considered the most likely as it is known to play critical roles in osteocyte maturation and mineral metabolism. Sequencing of DNA from carriers showed a nonsense mutation 253T/C in dmp1. This T/C transition introduced a stop codon (R145X) that could truncate C-terminal amino acids. Using PCR restriction fragment length polymorphism for this mutation, genotyping showed that all affected sheep were “TT” genotypes, and phenotypically normal sheep were either “CT” or “CC”. These findings suggest that a mutation in dmp1 is responsible for inherited rickets of Corriedale sheep. A simple diagnostic test is currently being designed to identify carriers with the defective “T” alleles.

Genome-wide survey implicates the influence of copy number variants (CNVs) in the development of early-onset bipolar disorder

We used genome-wide single nucleotide polymorphism (SNP) data to search for the presence of copy number variants (CNVs) in 882 patients with bipolar disorder (BD) and 872 population-based controls. A total of 291 (33%) patients had an early age-at-onset < or =21 years (AO < or =21 years). We systematically filtered for CNVs that cover at least 30 consecutive SNPs and which directly affect at least one RefSeq gene. We tested whether (a) the genome-wide burden of these filtered CNVs differed between patients and controls and whether (b) the frequency of specific CNVs differed between patients and controls. Genome-wide burden analyses revealed that the frequency and size of CNVs did not differ substantially between the total samples of BD patients and controls. However, separate analysis of patients with AO < or =21 years and AO>21 years showed that the frequency of microduplications was significantly higher (P=0.0004) and the average size of singleton microdeletions was significantly larger (P=0.0056) in patients with AO < or =21 years compared with controls. A search for specific BD-associated CNVs identified two common CNVs: (a) a 160 kb microduplication on 10q11 was overrepresented in AO < or = 21 years patients (9.62%) compared with controls (3.67%, P=0.0005) and (b) a 248 kb microduplication on 6q27 was overrepresented in the AO< or = 21 years subgroup (5.84%) compared with controls (2.52%, P=0.0039). These data suggest that CNVs have an influence on the development of early-onset, but not later-onset BD. Our study provides further support for previous hypotheses of an etiological difference between early-onset and later-onset BD.

Associations of variants in CHRNA5/A3/B4 gene cluster with smoking behaviors in a Korean population.

Multiple genome-wide and targeted association studies reveal a significant association of variants in the CHRNA5-CHRNA3-CHRNB4 (CHRNA5/A3/B4) gene cluster on chromosome 15 with nicotine dependence. The subjects examined in most of these studies had a European origin. However, considering the distinct linkage disequilibrium patterns in European and other ethnic populations, it would be of tremendous interest to determine whether such associations could be replicated in populations of other ethnicities, such as Asians. In this study, we performed comprehensive association and interaction analyses for 32 single-nucleotide polymorphisms (SNPs) in CHRNA5/A3/B4 with smoking initiation (SI), smoking quantity (SQ), and smoking cessation (SC) in a Korean sample (N = 8,842). We found nominally significant associations of 7 SNPs with at least one smoking-related phenotype in the total sample (SI: P = 0.015∼0.023; SQ: P = 0.008∼0.028; SC: P = 0.018∼0.047) and the male sample (SI: P = 0.001∼0.023; SQ: P = 0.001∼0.046; SC: P = 0.01). A spectrum of haplotypes formed by three consecutive SNPs located between rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5′ end of CHRNB4 was associated with these three smoking-related phenotypes in both the total and the male sample. Notably, associations of these variants and haplotypes with SC appear to be much weaker than those with SI and SQ. In addition, we performed an interaction analysis of SNPs within the cluster using the generalized multifactor dimensionality reduction method and found a significant interaction of SNPs rs7163730 in LOC123688, rs6495308 in CHRNA3, and rs7166158, rs8043123, and rs11072793 in the intergenic region downstream from the 5′ end of CHRNB4 to be influencing SI in the male sample. Considering that fewer than 5% of the female participants were smokers, we did not perform any analysis on female subjects specifically. Together, our detected associations of variants in the CHRNA5/A3/B4 cluster with SI, SQ, and SC in the Korean smoker samples provide strong evidence for the contribution of this cluster to the etiology of SI, ND, and SC in this Asian population.