Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions.

The human Y-chromosome does not recombine across its male-specific part and is therefore an excellent marker of human migrations. It also plays an important role in male fertility. However, its evolution is difficult to fully understand because of repetitive sequences, inverted repeats and the potentially large role of gene conversion. Here we perform an evolutionary analysis of 62 Y-chromosomes of Danish descent sequenced using a wide range of library insert sizes and high coverage, thus allowing large regions of these chromosomes to be well assembled. These include 17 father-son pairs, which we use to validate variation calling. Using a recent method that can integrate variants based on both mapping and de novo assembly, we genotype 10898 SNVs and 2903 indels (max length of 27241 bp) in our sample and show by father-son concordance and experimental validation that the non-recurrent SNP and indel variation on the Y chromosome tree is called very accurately. This includes variation called in a 0.9 Mb centromeric heterochromatic region, which is by far the most variable in the Y chromosome. Among the variation is also longer sequence-stretches not present in the reference genome but shared with the chimpanzee Y chromosome. We analyzed 2.7 Mb of large inverted repeats (palindromes) for variation patterns among the two palindrome arms and identified 603 mutation and 416 gene conversions events. We find clear evidence for GC-biased gene conversion in the palindromes (and a balancing AT mutation bias), but irrespective of this, also a strong bias towards gene conversion towards the ancestral state, suggesting that palindromic gene conversion may alleviate Muller's ratchet. Finally, we also find a large number of large-scale gene duplications and deletions in the palindromic regions (at least 24) and find that such events can consist of complex combinations of simultaneous insertions and deletions of long stretches of the Y chromosome.

Codon pair bias is a remarkably stable characteristic of a species. Although functionally uncharacterized, robust virus attenuation was achieved by recoding of viral proteins using underrepresented codon pairs. Because viruses replicate exclusively inside living cells, we posited that their codon pair preferences reflect those of their host(s). Analysis of many human viruses showed, however, that the encoding of viruses is influenced only marginally by host codon pair preferences. Furthermore, examination of codon pair preferences of vertebrate, insect, and arthropod-borne viruses revealed that the latter do not utilize codon pairs overrepresented in arthropods more frequently than other viruses. We found, however, that codon pair bias is a direct consequence of dinucleotide bias. We conclude that codon pair bias does not play a major role in the encoding of viral proteins and that virus attenuation by codon pair deoptimization has the same molecular underpinnings as attenuation based on an increase in CpG/TpA dinucleotides.

Chromatin topology is intricately linked to gene expression, yet its functional requirement remains unclear. Here, we comprehensively assessed the interplay between genome topology and gene expression using highly rearranged chromosomes (balancers) spanning ~75% of the Drosophila genome. Using transheterozyte (balancer/wild-type) embryos, we measured allele-specific changes in topology and gene expression in cis, while minimizing trans effects. Through genome sequencing, we resolved eight large nested inversions, smaller inversions, duplications and thousands of deletions. These extensive rearrangements caused many changes to chromatin topology, disrupting long-range loops, topologically associating domains (TADs) and promoter interactions, yet these are not predictive of changes in expression. Gene expression is generally not altered around inversion breakpoints, indicating that mis-appropriate enhancer-promoter activation is a rare event. Similarly, shuffling or fusing TADs, changing intra-TAD connections and disrupting long-range inter-TAD loops does not alter expression for the majority of genes. Our results suggest that properties other than chromatin topology ensure productive enhancer-promoter interactions.

B chromosomes constitute a heterogeneous mixture of genomic parasites that are sometimes derived intraspecifically from the standard genome of the host species, but result from interspecific hybridization in other cases. The mode of origin determines the DNA content, with the B chromosomes showing high similarity with the A genome in the first case, but presenting higher similarity with a different species in the second. The characid fish Moenkhausia sanctaefilomenae harbours highly invasive B chromosomes, which are present in all populations analyzed to date in the Parana and Tietê rivers. To investigate the origin of these B chromosomes, we analyzed two natural populations: one carrying B chromosomes and the other lacking them, using a combination of molecular cytogenetic techniques, nucleotide sequence analysis and high-throughput sequencing (Illumina HiSeq2000). Our results showed that i) B chromosomes have not yet reached the Paranapanema River basin; ii) B chromosomes are mitotically unstable; iii) there are two types of B chromosomes, the most frequent of which is lightly C-banded (similar to euchromatin in A chromosomes) (B1), while the other is darkly C-banded (heterochromatin-like) (B2); iv) the two B types contain the same tandem repeat DNA sequences (18S ribosomal DNA, H3 histone genes, MS3 and MS7 satellite DNA), with a higher content of 18S rDNA in the heterochromatic variant; v) all of these repetitive DNAs are present together only in the paracentromeric region of autosome pair no. 6, suggesting that the B chromosomes are derived from this A chromosome; vi) the two B chromosome variants show MS3 sequences that are highly divergent from each other and from the 0B genome, although the B2-derived sequences exhibit higher similarity with the 0B genome (this suggests an independent origin of the two B variants, with the less frequent, B2 type presumably being younger); and vii) the dN/dS ratio for the H3.2 histone gene is almost 4-6 times higher for B chromosomes than for A chromosome sequences, suggesting that purifying selection is relaxed for the DNA sequences located on the B chromosomes, presumably because they are mostly inactive.

The Labeoninae subfamily is a highly diversified but demonstrably monophyletic lineage of cyprinid fishes comprising five tribes and six incertae sedis genera. This widely distributed assemblage contains some 48 genera and around 480 recognized species distributed in freshwaters of Africa and Asia. In this study, the karyotypes and other chromosomal properties of five Labeoninae species found in Thailand Labeo chrysophekadion (Labeonini) and Epalzeorhynchos bicolor, Epalzeorhynchos munense, Henicorhynchus siamensis, Thynnichthys thynnoides (´Osteochilini´) were examined using conventional and molecular cytogenetic protocols. Our results confirmed a diploid chromosome number (2n) invariably 2n = 50, but the ratio of uni- and bi-armed chromosomes was highly variable among their karyotypes, indicating extensive structural chromosomal rearrangements. Karyotype of L. chrysophekadion contained 10m+6sm+20st+14a, 32m+10sm+8st for H. siamensis, 20m+12sm+10st+8a in E. bicolor, 20m+8sm+8st+14a in E. munense, and 18m+24sm+8st in T. thynnoides. Except for H. siamensis, which had four sites of 5S rDNA sites, other species under study had only one chromosome pair with those sites. In contrast, only one pair containing 18S rDNA sites were found in the karyotypes of three species, whereas two sites were found in that of E. bicolor. These cytogenetic patterns indicated that the cytogenomic divergence patterns of these labeonine species largely corresponded to the inferred phylogenetic tree. In spite of the 2n stability, diverse patterns of rDNA and microsatellite distribution as well as their various karyotype structures demonstrated significant evolutionary differentiation of Labeoninae genomes as exemplified in examined species. Labeoninae offers a traditional point of view on the evolutionary forces fostering biological diversity, and the recent findings add new pieces to comprehend the function of structural chromosomal rearrangements in adaption and speciation.

Wilms' tumor is the most common renal malignancy in children. In addition to staging, molecular risk stratification, such as loss of heterozygosity (LOH) in Chromosomes 1 and 16, is being increasingly used. Although genetic predisposition syndromes have been well-characterized in some Wilms' tumors, recent sequencing and biology efforts are expanding the classification of this malignancy. Here we present a case of siblings with remarkably similar presentations of bilateral Wilms' tumor at ∼12 mo of age. Thorough exam after the younger sibling was diagnosed did not reveal any signs to suggest one of the known Wilms' predisposition syndromes. Both were treated with standard therapies with good response and long-term sustained complete remission of 53 and 97 mo, respectively. Whole-exome sequencing was performed on a tumor sample from each patient and matched blood from one, revealing a shared truncation mutation of TRIM28 in all three samples with heterozygosity in the germline sample. TRIM28 loss has been recently implicated in early-stage Wilms' tumors with epithelioid morphology. These siblings expand the phenotype for presentation with multifocal disease with retained excellent response to standard therapy.

Sex chromosome divergence has been documented across phylogenetically diverse species, with amphibians typically having cytologically nondiverged ("homomorphic") sex chromosomes. With an aim of further characterizing sex chromosome divergence of an amphibian, we used "RAD-tags" and Sanger sequencing to examine sex specificity and heterozygosity in the Western clawed frog Silurana tropicalis (also known as Xenopus tropicalis). Our findings based on approximately 20 million genotype calls and approximately 200 polymerase chain reaction-amplified regions across multiple male and female genomes failed to identify a substantially sized genomic region with genotypic hallmarks of sex chromosome divergence, including in regions known to be tightly linked to the sex-determining region. We also found that expression and molecular evolution of genes linked to the sex-determining region did not differ substantially from genes in other parts of the genome. This suggests that the pseudoautosomal region, where recombination occurs, comprises a large portion of the sex chromosomes of S. tropicalis. These results may in part explain why African clawed frogs have such a high incidence of polyploidization, shed light on why amphibians have a high rate of sex chromosome turnover, and raise questions about why homomorphic sex chromosomes are so prevalent in amphibians.

Sex determination and the regulation of sexual dimorphism are among the most fascinating topics in modern biology. As the most species-rich group of sexually reproducing organisms on Earth, insects have multiple sex determination systems. Though sex chromosomes and sex-biased genes are well-studied in dozens of insects, their gene sequences are scattered in various databases. Moreover, a shortage of annotation hinders the deep mining of these data. Here, we collected the chromosome-level sex chromosome data of 49 insect species, including 34 X chromosomes, 15 Z chromosomes, 5 W chromosomes and 2 Y chromosomes. We also obtained Y-linked contigs of four insects species-Anopheles gambiae, Drosophila innubila, Drosophila yakuba and Tribolium castaneum. The unannotated chromosome-level sex chromosomes were annotated using a standard pipeline, yielding a total of 123 030 protein-coding genes, 2 159 427 repeat sequences, 894 miRNAs, 1574 rRNAs, 5105 tRNAs, 395 snoRNAs (small nucleolar RNA), 54 snRNAs (small nuclear RNA) and 5959 other ncRNAs (non-coding RNA). In addition, 36 781 sex-biased genes were identified by analyzing 62 RNA-seq (RNA sequencing) datasets. Together with 5707 sex-biased genes from the Drosophila genus collected from the Sex-Associated Gene Database, we obtained a total of 42 488 sex-biased genes from 13 insect species. All these data were deposited into InSexBase, a new user-friendly database of insect sex chromosomes and sex-biased genes. Database URL: http://www.insect-genome.com/Sexdb/.

The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12-15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79-116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

In genome-wide association studies (GWASs) of colorectal cancer, we have identified two genomic regions in which pairs of tagging-single nucleotide polymorphisms (tagSNPs) are associated with disease; these comprise chromosomes 1q41 (rs6691170, rs6687758) and 12q13.13 (rs7163702, rs11169552). We investigated these regions further, aiming to determine whether they contain more than one independent association signal and/or to identify the SNPs most strongly associated with disease. Genotyping of additional sample sets at the original tagSNPs showed that, for both regions, the two tagSNPs were unlikely to identify a single haplotype on which the functional variation lay. Conversely, one of the pair of SNPs did not fully capture the association signal in each region. We therefore undertook more detailed analyses, using imputation, logistic regression, genealogical analysis using the GENECLUSTER program and haplotype analysis. In the 1q41 region, the SNP rs11118883 emerged as a strong candidate based on all these analyses, sufficient to account for the signals at both rs6691170 and rs6687758. rs11118883 lies within a region with strong evidence of transcriptional regulatory activity and has been associated with expression of PDGFRB mRNA. For 12q13.13, a complex situation was found: SNP rs7972465 showed stronger association than either rs11169552 or rs7136702, and GENECLUSTER found no good evidence for a two-SNP model. However, logistic regression and haplotype analyses supported a two-SNP model, in which a signal at the SNP rs706793 was added to that at rs11169552. Post-GWAS fine-mapping studies are challenging, but the use of multiple tools can assist in identifying candidate functional variants in at least some cases.

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Antifungal drug tolerance is a response distinct from resistance, in which cells grow slowly above the MIC. Here, we found that the majority (69.2%) of 133 Candida albicans clinical isolates, including standard lab strain SC5314, exhibited temperature-enhanced tolerance at 37°C and 39°C, and were not tolerant at 30°C. Other isolates were either always tolerant (23.3%) or never tolerant (7.5%) at these three temperatures, suggesting that tolerance requires different physiological processes in different isolates. At supra-MIC fluconazole concentrations (8 to 128 μg/mL), tolerant colonies emerged rapidly at a frequency of ~10-3. In liquid passages over a broader range of fluconazole concentrations (0.25 to 128 μg/mL), tolerance emerged rapidly (within one passage) at supra-MICs. In contrast, resistance appeared at sub-MICs after 5 or more passages. Of 155 adaptors that evolved higher tolerance, all carried one of several recurrent aneuploid chromosomes, often including chromosome R, alone or in combination with other chromosomes. Furthermore, loss of these recurrent aneuploidies was associated with a loss of acquired tolerance, indicating that specific aneuploidies confer fluconazole tolerance. Thus, genetic background and physiology and the degree of drug stress (above or below the MIC) influence the evolutionary trajectories and dynamics with which antifungal drug resistance or tolerance emerges. IMPORTANCE Antifungal drug tolerance differs from drug resistance: tolerant cells grow slowly in drug, while resistant cells usually grow well, due to mutations in a few known genes. More than half of Candida albicans clinical isolates have higher tolerance at body temperature than they do at the lower temperatures used for most lab experiments. This implies that different isolates achieve drug tolerance via several cellular processes. When we evolved different strains at a range of high drug concentrations above inhibitory levels, tolerance emerged rapidly and at high frequency (one in 1,000 cells) while resistance appeared only later at very low drug concentrations. An extra copy of all or part of chromosome R was associated with tolerance, while point mutations or different aneuploidies were seen with resistance. Thus, genetic background and physiology, temperature, and drug concentration all influence how drug tolerance or resistance evolves.

Mammalian X and Y chromosomes evolved from a pair of autosomes. Although most ancestral genes have been lost from the Y chromosome, a small number of ancestral X-Y gene pairs are still present on the sex chromosomes. The KDM5C and KDM5D genes, which encode H3K4 histone demethylases, are a surviving ancestral gene pair located on the X and Y chromosomes, respectively. Mutations in KDM5C cause X-linked intellectual disability in human males, suggesting functional divergence between KDM5C and KDM5D in the nervous system. In this study, to explore the functional conservation and divergence between these two genes in other organs, we generated female mice lacking Kdm5c (homozygous X5c- X5c- females) and male mice lacking both Kdm5c and Kdm5d (compound hemizygous X5c- Y5d- males). Both X5c- X5c- females and X5c- Y5d- males showed lower body weights and postnatal lethality. Histological examination of the hearts showed prominent trabecular extension and a thin layer of compacted myocardium in the left and right ventricles, indicating noncompaction cardiomyopathy. However, hemizygous males lacking either Kdm5c or Kdm5d showed no signs of noncompaction cardiomyopathy. These results clearly demonstrate that the function of Kdm5c and Kdm5d in heart development is conserved.

It is widely known that castration has a significant effect on the accumulation of adipose tissue. microRNAs (miRNAs) are known to be involved in fat deposition and to be regulated by the androgen-induced androgen receptor (AR). However, there is little understanding of the relationship between miRNAs and fat deposition after castration. In this study, the high-throughput SOLiD sequencing approach was used to identify and characterize miRNA expression in backfat from intact and castrated full-sib male 23-week-old pigs. The patterns of adipogenesis and fat deposition were compared between castrated and intact male pigs.

Blood stream invasion by Escherichia coli is the commonest cause of bacteremia in the UK and elsewhere with an attributable mortality of about 15-20 %; antibiotic resistance to multiple agents is common in this microbe and is associated with worse outcomes. Genes conferring antimicrobial resistance, and their frequent location on horizontally transferred genetic elements is well-recognised, but the origin of these determinants, and their ability to be maintained and spread within clinically-relevant bacterial populations is unclear. Here, we set out to examine the distribution of antimicrobial resistance genes in chromosomes and plasmids of 16 bloodstream isolates of E. coli from patients within Scotland, and how these genes are maintained and spread. Using a combination of short and long-read whole genome sequencing methods, we were able to assemble complete sequences of 44 plasmids, with 16 Inc group F and 20 col plasmids; antibiotic resistance genes located almost exclusively within the F group. blaCTX-M15 genes had re-arranged in some strains into the chromosome alone (five strains), while others contained plasmid copies alone (two strains). Integrons containing multiple antibiotic genes were widespread in plasmids, notably many with a dfrA7 gene encoding resistance to trimethoprim, thus linking trimethoprim resistance to the other antibiotic resistance genes within the plasmids. This will allow even narrow spectrum antibiotics such as trimethoprim to act as a selective agent for plasmids containing antibiotic resistance genes mediating much broader resistance, including blaCTX-M15. To our knowledge, this is the first analysis to provide complete sequence data of chromosomes and plasmids in a collection of pathogenic human bloodstream isolates of E. coli. Our findings reveal the interplay between plasmids and integrative and conjugative elements in the maintenance and spread of antibiotic resistance genes within pathogenic E. coli.

One of the biggest challenges in chromosome biology is to understand the occurrence and complex genetics of the extra, non-essential karyotype elements, commonly known as supernumerary or B chromosomes (Bs). The non-Mendelian inheritance and non-pairing abilities of B chromosomes make them an interesting model for genomics studies, thus bringing to bear different questions about their genetic composition, evolutionary survival, maintenance and functional role inside the cell. This study uncovers these phenomena in multiple species that we considered as representative organisms of both vertebrate and invertebrate models for B chromosome analysis.

Sex chromosomes often carry large nonrecombining regions that can extend progressively over time, generating evolutionary strata of sequence divergence. However, some sex chromosomes display an incomplete suppression of recombination. Large genomic regions without recombination and evolutionary strata have also been documented around fungal mating-type loci, but have been studied in only a few fungal systems. In the model fungus Podospora anserina (Ascomycota, Sordariomycetes), the reference S strain lacks recombination across a 0.8-Mb region around the mating-type locus. The lack of recombination in this region ensures that nuclei of opposite mating types are packaged into a single ascospore (pseudohomothallic lifecycle). We found evidence for a lack of recombination around the mating-type locus in the genomes of ten P. anserina strains and six closely related pseudohomothallic Podospora species. Importantly, the size of the nonrecombining region differed between strains and species, as indicated by the heterozygosity levels around the mating-type locus and experimental selfing. The nonrecombining region is probably labile and polymorphic, differing in size and precise location within and between species, resulting in occasional, but infrequent, recombination at a given base pair. This view is also supported by the low divergence between mating types, and the lack of strong linkage disequilibrium, chromosomal rearrangements, transspecific polymorphism and genomic degeneration. We found a pattern suggestive of evolutionary strata in P. pseudocomata. The observed heterozygosity levels indicate low but nonnull outcrossing rates in nature in these pseudohomothallic fungi. This study adds to our understanding of mating-type chromosome evolution and its relationship to mating systems.

LC16m8 (m8), a highly attenuated vaccinia virus (VAC) strain, was developed as a smallpox vaccine, and its safety and immunogenicity have been confirmed. Here, we aimed to develop a system that recovers infectious m8 from a bacterial artificial chromosome (BAC) that retains the full-length viral genomic DNA (m8-BAC system). The infectious virus was successfully recovered from a VAC-BAC plasmid, named pLC16m8-BAC. Furthermore, the bacterial replicon-free virus was generated by intramolecular homologous recombination and was successfully recovered from a modified VAC-BAC plasmid, named pLC16m8.8S-BAC. Also, the growth of the recovered virus was indistinguishable from that of authentic m8. The full genome sequence of the plasmid, which harbors identical inverted terminal repeats (ITR) to that of authentic m8, was determined by long-read next-generation sequencing (NGS). The ITR contains x 18 to 32 of the 70 and x 30 to 45 of 54 base pair tandem repeats, and the number of tandem repeats was different between the ITR left and right. Since the virus recovered from pLC16m8.8S-BAC was expected to retain the identical viral genome to that of m8, including the ITR, a reference-based alignment following a short-read NGS was performed to validate the sequence of the recovered virus. Based on the pattern of coverage depth in the ITR, no remarkable differences were observed between the virus and m8, and the other region was confirmed to be identical as well. In summary, this new system can recover the virus, which is geno- and phenotypically indistinguishable from authentic m8.

Hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1) is a common housekeeping gene for sample normalization in the quantitative reverse transcriptase polymerase chain (qRT-PCR). However, co-amplification of HPRT1 pseudogenes may affect accurate results obtained in qRT-PCR. We designed a primer pair (HPSF) for pseudogene-free amplification of HPRT1 in qRT-PCR [1]. We showed specific amplification of HPRT1 mRNA in some common laboratory cell lines, including HeLa, NIH/3T3, CHO, BHK, COS-7 and VERO. This article provides data supporting the presence and location of HPRT1 pseudogenes within human and mouse genome, and the strategies used for designing primers that avoid the co-amplification of contaminating pseudogenes in qRT-PCR. In silico analysis of human genome showed three homologous sequences for HPRT1 on chromosomes 4, 5 and 11. The mRNA sequence of HPRT1 was aligned with the pseudogenes, and the primers were designed toward 5' end of HPRT1 mRNA that was only specific to HPRT1 mRNA not to the pseudogenes. The standard curve plot generated by HPSF primers showed the correlation coefficient of 0.999 and the reaction efficiency of 99.5%. Our findings suggest that HPSF primers can be recommended as a candidate primer pair for accurate and reproducible qRT-PCR assays.