A reliable, non-invasive PCR method for giant panda (Ailuropoda melanoleuca) sex identification

We developed an inexpensive, fast and reliable PCR method for sex identification of giant panda (Ailuropoda melanoleuca) by using one pair of primers to co-amplify homologous fragments with size polymorphism that located at amelogenin (AMEL) exon 5. In giant panda, a 63 bp deletion in exon 5 of Y-linked allele provides a significant discrimination between AMELX and AMELY, thus the amplification products can be distinguished simply by agarose gel electrophoresis, exhibiting sex-specific banding patterns (male: 237 bp, 174 bp; female: 237 bp). Both blood and feces samples from known-sex giant pandas were successfully amplified. Cross species test also revealed that this method could be applied to other Ursidae species. These authors contributed equally to this work.     

Female‐only sex‐linked amplified fragment length polymorphism markers support ZW/ZZ sex determination in the giant freshwater prawn Macrobrachium rosenbergii

Sex determination mechanisms in many crustacean species are complex and poorly documented. In the giant freshwater prawn, Macrobrachium rosenbergii, a ZW/ZZ sex determination system was previously proposed based on sex ratio data obtained by crosses of sex‐reversed females (neomales). To provide molecular evidence for the proposed system, novel sex‐linked molecular markers were isolated in this species. Amplified fragment length polymorphism (AFLP) using 64 primer combinations was employed to screen prawn genomes for DNA markers linked with sex loci. Approximately 8400 legible fragments were produced, 13 of which were uniquely identified in female prawns with no indication of corresponding male‐specific markers. These AFLP fragments were reamplified, cloned and sequenced, producing two reliable female‐specific sequence characterized amplified region (SCAR) markers. Additional individuals from two unrelated geographic populations were used to verify these findings, confirming female‐specific amplification of single bands. Detection of internal polymorphic sites was conducted by designing new primer pairs based on these internal fragments. The internal SCAR fragments also displayed specificity in females, indicating high levels of variation between female and male specimens. The distinctive feature of female‐linked SCAR markers can be applied for rapid detection of prawn gender. These sex‐specific SCAR markers and sex‐associated AFLP candidates unique to female specimens support a sex determination system consistent with female heterogamety (ZW) and male homogamety (ZZ).     

A simple and improved PCR-based technique for white-tailed deer (<Emphasis Type="Italic">Odocoileus virginianus</Emphasis>) sex identification

We describe a simple single-reaction technique for identifying the sex of white-tailed deer (Odocoileus virginianus) based on the PCR amplification of a zinc-finger intron using one pair of primers. Although Sry-coamplification confirmed sex identities, use of the Sry marker was unnecessary due to dimorphic alleles on the X and Y chromosomes at the zinc-finger locus. Insertions in intron 7 of the Y-linked allele (417 bp) make it nearly twice as long as the X-linked allele (236 bp) and thus the amplification products are easily discernable by simple agarose gel electrophoresis. The relatively short size of these products makes them useful for DNA-based sex identification from potentially low-yield tissue samples (e.g., hair, feces). This technique will provide ecologists, conservation geneticists and wildlife managers with a mechanism to readily and reliably identify the sex of unknown white-tailed deer tissue samples, and likely similar samples from other cervid species.     

Molecular sexing of African rhinoceros

We report the development of a fast and reliable PCR-based method for sex identification of African Rhinoceros (Ceratotherium simum and Diceros bicornis) that could easily be incorporated into fluorescent short tandem repeat (STR) profiling. A single primer pair, consisting of a fluorescently labelled forward primer and an unlabelled reverse primer, is used to co-amplify homologous fragments of a zinc finger (ZF) protein intron which exhibits size polymorphism between the X and Y chromosomes. In both species, the amplified ZFX and ZFY amplicons differ in size by 7 bp and can thus be differentiated by capillary electrophoresis. Blood, tissue, horn, and faecal samples were correctly sexed using this method. Cross species testing also demonstrated that this method could be used to sex Indian rhinoceros (Rhinoceros unicornis) samples.     

A Comparative Analysis of ISSR and RAPD Markers for Study of Genetic Diversity in Shisham (<Emphasis Type="Italic">Dalbergia sissoo</Emphasis>)

Shisham (Dalbergia sissoo) is one of the most preferred timber tree species of South Asia. Two DNA-based molecular marker techniques, intersimple sequence repeat (ISSR) and random amplified polymorphism DNA (RAPD), were compared to study the genetic diversity in this species. A total of 30 polymorphic primers (15 ISSR and 15 random) were used. Amplification of genomic DNA of 22 genotypes, using ISSR analysis, yielded 117 fragments, of which 64 were polymorphic. Number of amplified fragments with ISSR primers ranged from five to ten and varied in size from 180 to 1,900 bp. Percentage polymorphism ranged from 0 to 87.5. The 15 RAPD primers produced 144 bands across 22 genotypes, of which 84 were polymorphic. The number of amplified bands varied from five to 13, with size range from 180 to 2,400 bp. Percentage polymorphism ranged from 0 to 100, with an average of 58.3 across. RAPD markers were relatively more efficient than the ISSR assay. The mental test between two Jaccard’s similarity matrices gave r ≥ 0.90, showing very good fit correlation in between ISSR- and RAPD-based similarities. Clustering of isolates remained more or less the same in RAPD and combined data of RAPD and ISSR. The similarity coefficient ranged from 0.734 to 0.939, 0.563 to 0.946, and 0.648 to 0.920 with ISSR, RAPD, and combined dendrogram, respectively.     

A sensitive and specific multiplex PCR approach for sex identification of ursine and tremarctine bears suitable for non‐invasive samples

We report a new approach for molecular sex identification of extant Ursinae and Tremarctinae bears. Two Y‐specific fragments (SMCY and 318.2) and one X‐specific fragment (ZFX) are amplified in a multiplex PCR, yielding a double test for male‐specific amplification and an internal positive control. The primers were designed and tested to be bear‐specific, thereby minimizing the risk of cross‐amplification in other species including humans. The high sensitivity and small amplicon sizes (100, 124, 160 base pairs) facilitate analysis of non‐invasively obtained DNA material. DNA from tissue and blood as well as from 30 non‐invasively collected hair and faeces yielded clear and easily interpretable results. The fragments were detected both by standard gel electrophoresis and automated capillary electrophoresis.     

Simultaneous identification of Grapholita molesta Busck and Grapholita dimorpha Komai by PCR‐RFLP

Although several molecular diagnostic techniques are available for the identification of the apple‐feeding pests Grapholita molesta Busck and Grapholita dimorpha Komai, these pests are severely affecting apple orchards in Korea. These two pests may be misidentified or the available molecular diagnostic techniques may not facilitate the simultaneous identification of the morphological features of both species. In this study, we developed a multiplex assay for these two species using the polymerase chain reaction – restriction fragment length polymorphism (PCR‐RFLP) method. Sixty‐two specimens were collected from apples presumed infested with moth larvae and from pheromone traps from 2013 to 2014. Both species were identified morphologically, and a partial region of the cytochrome b gene was sequenced to design primers for PCR‐RFLP. Digestion profiles of G. molesta and G. dimorpha, using the Sau3A1 restriction enzyme, were characterized using three DNA fragments each for G. molesta (363 bp, 91 bp and 31 bp) and G. dimorpha (220 bp, 234 bp and 31 bp). The RFLP assay developed for both species in this study was more efficient and accurate than other currently used diagnostic assays and would be helpful to identify field‐collected specimens for pest control research.     

Molecular sexing and sources of CHD1‐Z/W sequence variation in Hawaiian birds

Sequence information from 28 CHD1 gene fragments reveals that a primary source of variability in CHD1‐W genes is a variable intron microsatellite; a single‐codon deletion was found in the 3′ exon in one species. Sequence variation of CHD1‐Z genes was detected in males that altered polymerase chain reaction (PCR) fragment length. Three sets of CHD1‐based primers were evaluated for sex determination in 12 endemic and 8 alien Hawaiian species, including one of the last po’o‐uli. Combined, these primers provide a reliable means of sex determination in most species (including the po’o‐uli), and have produced a valuable reference database for future expanded population‐level studies.     

Heteroduplex molecules cause sexing errors in a standard molecular protocol for avian sexing

Molecular methods are a necessary tool for sexing monomorphic birds. These molecular approaches are usually reliable, but sexing protocols should be evaluated carefully because biochemical interactions may lead to errors. We optimized laboratory protocols for genetic sexing of a monomorphic shorebird, the upland sandpiper (Bartramia longicauda), using two independent sets of primers, P2/P8 and 2550F/2718R, to amplify regions of the sex‐linked CHD‐Z and CHD‐W genes. We discovered polymorphisms in the region of the CHD‐Z intron amplified by the primers P2/P8 which caused four males to be misidentified as females (n = 90 mated pairs). We cloned and sequenced one CHD‐W allele (370 bp) and three CHD‐Z alleles in our population: Z° (335 bp), Z′ (331 bp) and Z″ (330 bp). Normal (Z°Z°) males showed one band in agarose gel analysis and were easily differentiated from females (Z°W), which showed two bands. However, males heterozygous for CHD‐Z alleles (Z′Z″) unexpectedly showed two bands in a pattern similar to females. While the Z′ and Z″ fragments contained only short deletions, they annealed together during the polymerase chain reaction (PCR) process and formed heteroduplex molecules that were similar in size to the W fragment. Errors previously reported for molecular sex‐assignment have usually been due to allelic dropout, causing females to be misidentified as males. Here, we report evidence that events in PCRs can lead to the opposite error, with males misidentified as females. We recommend use of multiple primer sets and large samples of known‐sex birds for validation when designing protocols for molecular sex analysis.     

High‐resolution melting analysis for bird sexing: a successful approach to molecular sex identification using different biological samples

High‐resolution melting (HRM) analysis is a very attractive and flexible advanced post‐PCR method with high sensitivity/specificity for simple, fast and cost‐effective genotyping based on the detection of specific melting profiles of PCR products. Next generation real‐time PCR systems, along with improved saturating DNA‐binding dyes, enable the direct acquisition of HRM data after quantitative PCR. Melting behaviour is particularly influenced by the length, nucleotide sequence and GC content of the amplicons. This method is expanding rapidly in several research areas such as human genetics, reproductive biology, microbiology and ecology/conservation of wild populations. Here we have developed a successful HRM protocol for avian sex identification based on the amplification of sex‐specific CHD1 fragments. The melting curve patterns allowed efficient sexual differentiation of 111 samples analysed (plucked feathers, muscle tissues, blood and oral cavity epithelial cells) of 14 bird species. In addition, we sequenced the amplified regions of the CHD1 gene and demonstrated the usefulness of this strategy for the genotype discrimination of various amplicons (CHD1Z and CHD1W), which have small size differences, ranging from 2 bp to 44 bp. The established methodology clearly revealed the advantages (e.g. closed‐tube system, high sensitivity and rapidity) of a simple HRM assay for accurate sex differentiation of the species under study. The requirements, strengths and limitations of the method are addressed to provide a simple guide for its application in the field of molecular sexing of birds. The high sensitivity and resolution relative to previous real‐time PCR methods makes HRM analysis an excellent approach for improving advanced molecular methods for bird sexing.     

Sex identification of Eurasian otter (<Emphasis Type="Italic">Lutra lutra</Emphasis>) non-invasive DNA samples using ZFX/ZFY sequences

We developed a new PCR/RFLP system targeted to amplify and cut a segment of the ZFX/ZFY gene in non-invasive otter (Lutra lutra) samples. This assay produced one sex-specific fragment in females (XX genotypes) and two fragments in males (XY genotypes), and is intrinsically more reliable than alternative systems (e.g., SRY genes) that are based on the amplification of a single Y-linked fragment. The new primer pair correctly identified the sex of 23 DNA extracted from sexed otter tissues. This procedure was used to sex unknown DNA extracted from otter scats. A multi-tube approach led to identify the sex of 72/91 (79%) samples, using a minimum of two PCR replicates. This procedure is currently used in non-invasive genetic monitoring of Italian otter populations.     

Screening and identification of male‐specific DNA fragments in common carps Cyprinus carpio using suppression subtractive hybridization

In this study, a sex subtractive genomic DNA library was constructed using suppression subtractive hybridization (SSH) between male and female Cyprinus carpio. Twenty‐two clones with distinguishable hybridization signals were selected and sequenced. The specific primers were designed based on the sequence data. Those primers were then used to amplify the sex‐specific fragments from the genomic DNA of male and female carp. The amplified fragments from two clones showed specificity to males but not to females, which were named as Ccmf2 [387 base pairs (bp)] and Ccmf3 (183 bp), respectively. The sex‐specific pattern was analysed in a total of 40 individuals from three other different C. carpio. stocks and grass carp Ctenopharyngodon idella using Ccmf2 and Ccmf3 as dot‐blotting probes. The results revealed that the molecular diversity exists on the Y chromosome of C. carpio. No hybridization signals, however, were detected from individuals of C. idella, suggesting that the two sequences are specific to C. carpio. No significant homologous sequences of Ccmf2 and Ccmf3 were found in GenBank. Therefore, it was interpreted that the results as that Ccmf2 and Ccmf3 are two novel male‐specific sequences; and both fragments could be used as markers to rapidly and accurately identify the genetic sex of part of C. carpio. This may provide a very efficient selective tool for practically breeding monosex female populations in aquacultural production.     

A comparative analysis of genetic diversity in blackgram genotypes using RAPD and ISSR markers

Random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) markers were used to study the DNA polymorphism in elite blackgram genotypes. A total of 25 random and 16 ISSR primers were used. Amplification of genomic DNA of the 18 genotypes, using RAPD analysis, yielded 104 fragments that could be scored, of which 44 were polymorphic, with an average of 1.8 polymorphic fragments per primer. Number of amplified fragments with random primers ranged from two (OPA-13) to nine (OPK-4) and varied in size from 200 bp to 2,500 bp. Percentage polymorphism ranged from 16.6% (OPK-7) to a maximum of 66.6% (OPE-5, OPH-2, and OPK-8), with an average of 42.7%. The 16 ISSR primers used in the study produced 101 bands across 18 genotypes, of which 55 were polymorphic. The number of amplified bands varied from two (ISSR 858) to ten (ISSR 810), with a size range of 200–2,200 bp. The average numbers of bands per primer and polymorphic bands per primer were 6.3 and 3.4, respectively. Percentage polymorphism ranged from 25% (ISSR 885) to 100% (ISSR 858), with an average percentage polymorphism of 57.5% across all the genotypes. The 3-anchored primers based on poly(GA) and poly(AG) motifs produced high average polymorphisms of 54.98% and 58.32%, respectively. ISSR markers were more efficient than the RAPD assay, as they detected 57.4% polymorphic DNA markers in Vigna mungo as compared to 42.7% for RAPD markers. The Mantel test between the two Jaccards similarity matrices gave r =0.32, showing low correlation between RAPD- and ISSR-based similarities. Clustering of genotypes within groups was not similar when RAPD and ISSR derived dendrogram were compared, whereas the pattern of clustering of the genotypes remained more or less the same in ISSR and combined data of RAPD and ISSR.     

Design of Mini‐barcode for Catfishes for assessment of archival biodiversity

Recovery of DNA barcode sequences is often challenging from the archived specimens. However, short fragments of DNA may be recovered, which would significantly improve many unresolved taxonomic conflicts. Here, we designed a mini‐barcode for catfishes comprising several species and many cryptic taxa. We analysed a data set of 3048 publicly available COI barcode sequences representing 547 worldwide catfish species and performed 152 628 interspecies comparisons. A significantly more positively correlated interspecies distance was detected with transversion (0.78, P < 0.001) than with transition (0.70, P < 0.001). This suggested that transversions were better diagnostics for species identification. In the aligned data set, two transversion‐rich fragments (53 bp and 119 bp) were identified. Transition/transversion bias value was 1.04 in 53‐bp fragment, 1.23 in 119‐bp fragment and 1.50 in full‐length barcode. The interspecies distance with full‐length barcode was 0.212 ± 0.037, while that with 53‐bp and 119‐bp fragments was 0.325 ± 0.039 and 0.218 ± 0.045, respectively. Survey of 53‐bp fragment showed a possibility of only 1144 barcodes, while that of 119‐bp fragment showed >4 million barcodes. Thus, the 119‐bp fragment is a viable mini‐barcode for catfishes comprising >3000 extant species. Experiment with 82 archived catfishes showed successful recovery of this mini‐barcode using the designed primer. The mini‐barcode sequences showed species‐specific similarity in the range of 98‐100% with the global database. Therefore, survey of a transversion‐rich fragment within the full‐length barcode would be an ideal approach of mini‐barcode design for biodiversity assessment.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

CONSPECIFICITY AND INDO-PACIFIC DISTRIBUTION OF SYMBIODINIUM GENOTYPES (DINOPHYCEAE) FROM GIANT CLAMS

We previously reported the occurrence of genetically‐diverse symbiotic dinoflagellates (zooxanthellae) within and between 7 giant clam species (Tridacnidae) from the Philippines based on the algal isolates' allozyme and random amplified polymorphic DNA (RAPD) patterns. We also reported that these isolates all belong to clade A of the Symbiodinium phylogeny with identical 18S rDNA sequences. Here we extend the genetic characterization of Symbiodinium isolates from giant clams and propose that they are conspecific. We used the combined DNA sequences of the internal transcribed spacer (ITS)1, 5.8S rDNA, and ITS2 regions (rDNA‐ITS region) because the ITS1 and ITS2 regions evolve faster than 18S rDNA and have been shown to be useful in distinguishing strains of other dinoflagellates. DGGE of the most variable segment of the rDNA‐ITS region, ITS1, from clonal representatives of clades A, B, and C showed minimal intragenomic variation. The rDNA‐ITS region shows similar phylogenetic relationships between Symbiodinium isolates from symbiotic bivalves and some cnidarians as does 18S rDNA, and that there are not many different clade A species or strains among cultured zooxanthellae (CZ) from giant clams. The CZ from giant clams had virtually identical sequences, with only a single nucleotide difference in the ITS2 region separating two groups of isolates. These data suggest that there is one CZ species and perhaps two CZ strains, each CZ strain containing individuals that have diverse allozyme and RAPD genotypes. The CZ isolated from giant clams from different areas in the Philippines (21 isolates, 7 clam species), the Australian Great Barrier Reef (1 isolate, 1 clam species), Palau (8 isolates, 7 clam species), and Okinawa, Japan (1 isolate, 1 clam species) shared the same rDNA‐ITS sequences. Furthermore, analysis of fresh isolates from giant clams collected from these geographical areas shows that these bivalves also host indistinguishable clade C symbionts. These data demonstrate that conspecific Symbiodinium genotypes, particularly clade A symbionts, are distributed in giant clams throughout the Indo‐Pacific.     

Development of a simple PCR-based assay for the identification of triazine resistance in the noxious plant common ragweed (Ambrosia artemisiifolia) and its applicability in higher plants

Bidirectional allele-specific PCR (Bi-PASA) was used to detect a point mutation causing triazine resistance in common ragweed (Ambrosia artemisiifolia). The external primers amplified a 278 bp standard DNA fragment in all genotypes. In the susceptible S264S genotypes, a 208 bp fragment was expected while in resistant S264G common ragweed genotypes a 109 bp band was expected. In resistant plants, both the wild and mutant type fragments were detected, indicating that the original triazine sensitive cpDNA is maintained in a heteroplasmic state in the resistant S264G genotypes. Additionally, in silico analysis confirmed the potential applicability of our diagnostic assay for other plant species. In 24 out of 74 taxa (32%), the assay could be used without any change, while in the others some of the primers should be redesigned before use.     

CHARACTERIZATION OF GENETIC MARKERS LINKED TO SEX DETERMINATION IN THE HAPLOID‐DIPLOID RED ALGA GRACILARIA CHILENSIS1

Bulk segregant analysis, random amplified polymorphic DNA (RAPD), and sequence characterized amplified region (SCAR) methods were used to identify sex‐linked molecular markers in the haploid‐diploid rhodophyte Gracilaria chilensis C. J. Bird, McLachlan et E. C. Oliveira. One hundred and eighty 10 bp primers were tested on three bulks of DNA: haploid males, haploid females, and diploid tetrasporophytes. Three RAPD primers (OPD15, OPG16, and OPN20) produced male‐specific bands; and one RAPD primer (OPD12), a female‐specific band. The sequences of the cloned putative sex‐specific PCR fragments were used to design specific primers for the female marker SCAR‐D12‐386 and the male marker SCAR‐G16‐486. Both SCAR markers gave unequivocal band patterns that allowed sex and phase to be determined in G. chilensis. Thus, all the females presented only the female band, and all the males only the male band, while all the tetrasporophytes amplified both male and female bands. Despite this sex‐specific association, we were able to amplify SCAR‐D12‐386 and SCAR‐G16‐486 in both sexes at low melting temperature. The differences between male and female sequences were of 8%–9% nucleotide divergence for SCAR‐D12‐386 and SCAR‐G16‐486, respectively. SCAR‐D12‐386 and SCAR‐G16‐486 could represent degenerated or diverged sequences located in the nonrecombining region of incipient sex chromosomes or heteromorphic sex chromosomes with sequence differences at the DNA level such that PCR primers amplify only one allele and not the other in highly specific PCR conditions. Seven gametic progenies composed of 19 males, 19 females, and the seven parental tetrasporophytes were analyzed. In all of them, the two SCAR markers segregated perfectly with sexual phenotypes.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.     

Indels within promoter and intron 1 of bovine prion protein gene modulate the gene expression levels in the medulla oblongata of two Japanese cattle breeds

Genetic differences which exist in the prion protein gene (PRNP) have been reported to influence susceptibility of humans, sheep and goats to prion diseases. In cattle, however, none of the known coding polymorphisms has a direct effect on bovine spongiform encephalopathy (BSE). It has been reported that 23‐bp insertion/deletion (indel) polymorphisms within the promoter region have a tentative association to BSE susceptibility in German cattle, and a lower number of 24‐bp repeat units in the open reading frame (ORF) was reported to reduce BSE susceptibility in transgenic mice. In this study, because of the hypothesis that bovine PRNP promoter polymorphisms cause changes in PRNP expression, we genotyped PRNP polymorphisms in the promoter and intron 1 using 218 genomic DNA samples from two Japanese cattle breeds. We also analysed the expression levels of prion in 40 animals by quantification of real‐time PCR using mRNAs extracted from the medulla oblongata to study the relationship between PRNP genotypes and PRNP expression. We found a significant correlation between promoter indel polymorphisms and PRNP‐mRNA expression (P_0.0413) and therefore hypothesize that differences in polymorphisms could be one of the causes of differences in PRNP expression levels. We also report a novel difference in PRNP expression (P < 0.0001) between Japanese Black and Japanese Brown cattle breeds. There was no significant difference based on age and sex of the animals.