Population-tailored mock genome enables genomic studies in species without a reference genome.

Based on molecular markers, genomic prediction enables us to speed up breeding schemes and increase the response to selection. There are several high-throughput genotyping platforms able to deliver thousands of molecular markers for genomic study purposes. However, even though its widely applied in plant breeding, species without a reference genome cannot fully benefit from genomic tools and modern breeding schemes. We used a method to assemble a population-tailored mock genome to call single-nucleotide polymorphism (SNP) markers without an available reference genome, and for the first time, we compared the results with standard genotyping platforms (array and genotyping-by-sequencing (GBS) using a reference genome) for performance in genomic prediction models. Our results indicate that using a population-tailored mock genome to call SNP delivers reliable estimates for the genomic relationship between genotypes. Furthermore, genomic prediction estimates were comparable to standard approaches, especially when considering only additive effects. However, mock genomes were slightly worse than arrays at predicting traits influenced by dominance effects, but still performed as well as standard GBS methods that use a reference genome. Nevertheless, the array-based SNP markers methods achieved the best predictive ability and reliability to estimate variance components. Overall, the mock genomes can be a worthy alternative for genomic selection studies, especially for those species where the reference genome is not available.

Consensus of experts from the Spanish Pharmacogenetics and Pharmacogenomics Society and the Spanish Society of Medical Oncology for the genotyping of DPYD in cancer patients who are candidates for treatment with fluoropyrimidines.

5-Fluorouracil (5-FU) and oral fluoropyrimidines, such as capecitabine, are widely used in the treatment of cancer, especially gastrointestinal tumors and breast cancer, but their administration can produce serious and even lethal toxicity. This toxicity is often related to the partial or complete deficiency of the dihydropyrimidine dehydrogenase (DPD) enzyme, which causes a reduction in clearance and a longer half-life of 5-FU. It is advisable to determine if a DPD deficiency exists before administering these drugs by genotyping DPYD gene polymorphisms. The objective of this consensus of experts, in which representatives from the Spanish Pharmacogenetics and Pharmacogenomics Society and the Spanish Society of Medical Oncology participated, is to establish clear recommendations for the implementation of genotype and/or phenotype testing for DPD deficiency in patients who are candidates to receive fluoropyrimidines. The genotyping of DPYD previous to treatment classifies individuals as normal, intermediate, or poor metabolizers. Normal metabolizers do not require changes in the initial dose, intermediate metabolizers should start treatment with fluoropyrimidines at doses reduced to 50%, and poor metabolizers are contraindicated for fluoropyrimidines.

Comparative analysis of different PCR-based strategies for HPV detection and genotyping from cervical samples.

BACKGROUND: Human papillomavirus (HPV) is the main cause of cervical cancer. Polymerase chain reaction (PCR)-based techniques are associated with accurate results with respect to HPV detection and genotyping, being able to identify viral DNA at low levels. However, differences in primer design influence their sensibility and specificity, depending on the HPV type assessed. OBJECTIVE: The aim of the study was to comparatively evaluate the effectiveness of three different PCR-based strategies for HPV detection and genotyping from cervical samples. STUDY DESIGN: The procedures were based on different primer design strategies, using MY09/MY11, EntroA, and type specific multiplex PCR primers. RESULTS: Out of 411 samples of cervical scrapings, 45 (10.9%), 50 (12.2%), and 117 (28.5%) were positive for MY09/MY11, EntroA, and multiplex PCR, respectively. For MY09/MY11 positive samples, 36 were negative for EntroA and 23 for multiplex PCR. For EntroA positive samples, 40 were negative for MY09/MY11 and 26 for multiplex PCR. For multiplex PCR positive samples, 96 were negative for MY09/MY11 and 94 for EntroA. MY09/MY11 identified 12 different HPV types, EntroA detected eight types and multiplex PCR detected 11 HPV types. EntroA primers were able to detect HPV in more samples than MY09/MY11, while multiplex PCR, despite the limited targeted HPV types, presented higher sensibility than the other methods. CONCLUSIONS: The three methods presented different advantages and disadvantages, and the present study reinforces the need to use more than one molecular strategy for HPV detection and genotyping, and the development of novel methods which could overcome the limitations of the existing tests.

Molecular screening of XY SRY-negative sex reversal cases in horses revealed anomalies in amelogenin testing.

Male-to-female sex reversal in horses is a developmental disorder in which phenotypic females have a male genetic constitution. Male-to-female sex reversal is the second most common genetic sex abnormality, after X chromosome monosomy. All male-to-female sex reversal cases studied to date have been found to be infertile. Therefore, a screening test is particularly useful in laboratories doing DNA genotyping in horses. Our laboratory has tested > 209,000 horses for parentage using a panel of microsatellite markers and the sex marker gene amelogenin (AMEL). Suspect XY sex reversal cases are reported females with a male profile by AMEL testing. After routine genotyping, 49 cases were detected and further tested using the sex-determining region Y (SRY) gene, confirming the XY SRY-negative genotype of suspect sex reversal cases. When some inconsistencies arose in the initial result, a molecular panel of X- and Y-linked markers was analyzed for these samples. Of the 49 cases, 33 were confirmed as XY SRY-negative. The remaining 16 cases were identified as false-positives as a result of anomalies of AMEL testing in horses.

Development of a Multilocus Sequence Typing Scheme for Giardia intestinalis.

Giardia intestinalis is an intestinal protozoan most commonly found in humans. It has been grouped into 8 assemblages (A-H). Markers such as the glutamate dehydrogenase gene, triose phosphate isomerase and beta-giardin (ß-giardin) have been widely used for genotyping. In addition, different genetic targets have been proposed as a valuable alternative to assess diversity and genetics of this microorganism. Thus, our objective was to evaluate new markers for the study of the diversity and intra-taxa genetic structure of G. intestinalis in silico and in DNA obtained from stool samples. We analysed nine constitutive genes in 80 complete genome sequences and in a group of 24 stool samples from Colombia. Allelic diversity was evaluated by locus and for the concatenated sequence of nine loci that could discriminate up to 53 alleles. Phylogenetic reconstructions allowed us to identify AI, AII and B assemblages. We found evidence of intra- and inter-assemblage recombination events. Population structure analysis showed genetic differentiation among the assemblages analysed.

Concordance of Giardia duodenalis assemblages determined by different PCR methodologies in three observational studies in Cuba.

Giardia duodenalis is one of the most important intestinal parasites globally, especially in children, and in Cuba is the leading cause of chronic paediatric diarrhoea in this population. G. duodenalis is composed of eight genetic groups (or assemblages), two of which (A and B) are apparently zoonotic, occurring in both humans and other animals. However, consensus on the most appropriate genotyping scheme for optimal characterization of G. duodenalis isolates is lacking. In this article we present the results of three descriptive observational studies conducted in Havana, Cuba between 2010 and 2013, with the aim of comparing the results from molecular (PCR) approaches targeting different genes in order to assign with confidence 224 isolates of G. duodenalis to the correct assemblages. In each sub-study, following DNA isolation by the phenol/chloroform/isoamyl alcohol extraction method, PCR targeting the triose phosphate isomerase (tpi) gene was used for molecular characterization, as well as one additional PCR-method targeting another gene or pair of genes. DNA amplification was obtained in 87%, 83%, and 80% in the three sub-studies. Although excellent agreement (kappa index = 1) was recorded between results from some pairs of genes, for other combinations only moderate or substantial agreement was achieved. These results highlight the importance of interpretation of genotyping data, especially when single genetic markers are used. From the results of our studies, PCR targeting a combination of the tpi gene and the intergenic spacer region of rDNA may be a useful approach for the molecular characterization of G. duodenalis isolates.

Evaluation of imputation accuracy using the combination of two high-density panels in Nelore beef cattle.

This study compared imputation from lower-density commercial and customized panels to high-density panels and a combined panel (Illumina and Affymetrix) in Nelore beef cattle. Additionally, linkage disequilibrium and haplotype block conformation were estimated in individual high-density panels and compared with corresponding values in the combined panel after imputation. Overall, 814 animals were genotyped using BovineHD BeadChip (IllumHD), and 93 of these animals were also genotyped using the Axion Genome-Wide BOS 1 Array Plate (AffyHD). In general, customization considering linkage disequilibrium and minor allele frequency had the highest accuracies. The IllumHD panel had higher values of linkage disequilibrium for short distances between SNPs than AffyHD and the combined panel. The combined panel had an increased number of small haplotype blocks. The use of a combined panel is recommended due to its increased density and number of haplotype blocks, which in turn increase the probability of a marker being close to a quantitative trait locus of interest. Considering common SNPs between IllumHD and AffyHD for the customization of a low-density panel increases the imputation accuracy for IllumHD, AffyHD and the combined panel.

Population structure and genetic relationships between Ethiopian and Brazilian Coffea arabica genotypes revealed by SSR markers.

Information about population structure and genetic relationships within and among wild and brazilian Coffea arabica L. genotypes is highly relevant to optimize the use of genetic resources for breeding purposes. In this study, we evaluated genetic diversity, clustering analysis based on Jaccard's coefficient and population structure in 33 genotypes of C. arabica and of three diploid Coffea species (C. canephora, C. eugenioides and C. racemosa) using 30 SSR markers. A total of 206 alleles were identified, with a mean of 6.9 over all loci. The set of SSR markers was able to discriminate all genotypes and revealed that Ethiopian accessions presented higher genetic diversity than commercial varieties. Population structure analysis indicated two genetic groups, one corresponding to Ethiopian accessions and another corresponding predominantly to commercial cultivars. Thirty-four private alleles were detected in the group of accessions collected from West side of Great Rift Valley. We observed a lower average genetic distance of the C. arabica genotypes in relation to C. eugenioides than C. canephora. Interestingly, commercial cultivars were genetically closer to C. eugenioides than C. canephora and C. racemosa. The great allelic richness observed in Ethiopian Arabica coffee, especially in Western group showed that these accessions can be potential source of new alleles to be explored by coffee breeding programs.

Cost-effective and fast KIR gene-content genotyping by multiplex melting curve analysis.

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 µM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time.

Identification of reference genes for RT-qPCR analysis in peach genotypes with contrasting chilling requirements.

Selecting and validating reference genes are the first steps in studying gene expression by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR). The present study aimed to evaluate the stability of five reference genes for the purpose of normalization when studying gene expression in various cultivars of Prunus persica with different chilling requirements. Flower bud tissues of nine peach genotypes from Embrapa's peach breeding program with different chilling requirements were used, and five candidate reference genes based on the RT-qPCR that were useful for studying the relative quantitative gene expression and stability were evaluated using geNorm, NormFinder, and bestKeeper software packages. The results indicated that among the genes tested, the most stable genes to be used as reference genes are Act and UBQ10. This study is the first survey of the stability of reference genes in peaches under chilling stress and provides guidelines for more accurate RT-qPCR results.

Genotyping of SNPs associated with meat tenderness: comparison of two PCR-based methods.

Single nucleotide polymorphisms (SNPs) carried in calpain (CAPN1), calpastatin (CAST), and leptin (LEP) genes are associated with meat tenderness. Due to the economic importance of this meat quality attribute, the development of fast, reliable, and affordable methods to identify bovine carriers of favorable alleles is of great importance for genetic improvement. Currently, PCR-RFLP is accepted as the standard gold method for genotyping SNPs associated with meat tenderness. But these SNPs can be detected by other techniques as high-resolution melting (HRM) analysis - a post-PCR method - that offers several advantages and has great application potential in the meat industry. In this study, we standardized, validated, and compared the performance of PCR-HRM to that of PCR-RFLP in genotyping bovine SNPs associated with meat tenderness: CAPN4751, CAPN316, CAST2959, CAST282, LEPE2FB, and LEPE2JW. We analyzed genotypes of a total of 380 bovines, 110 Bos taurus and 270 Bos indicus. Results obtained with PCR-HRM were consistent with those found by PCR-RLFP. Furthermore, HRM was found to be highly sensitive, and our results confirmed the repeatability (intra-assay precision) and reproducibility (inter-assay precision) of this assay. An internal control for endonuclease activity was created using site-directed mutagenesis to generate an additional enzymatic restriction point useful to discriminate SNP alleles. Our results show that PCR-HRM is an efficient method that produces reliable and rapid results. However, should be had in account that the method of DNA extraction, the quality and quantity of DNA, analyst-related variations, and primer design may generate challenges for allele discrimination.

Measurements of experimental precision for trials with cowpea (Vigna unguiculata L. Walp.) genotypes.

The aim of this study was to evaluate the suitability of statistics as experimental precision degree measures for trials with cowpea (Vigna unguiculata L. Walp.) genotypes. Cowpea genotype yields were evaluated in 29 trials conducted in Brazil between 2005 and 2012. The genotypes were evaluated with a randomized block design with four replications. Ten statistics that were estimated for each trial were compared using descriptive statistics, Pearson correlations, and path analysis. According to the class limits established, selective accuracy and F-test values for genotype, heritability, and the coefficient of determination adequately estimated the degree of experimental precision. Using these statistics, 86.21% of the trials had adequate experimental precision. Selective accuracy and the F-test values for genotype, heritability, and the coefficient of determination were directly related to each other, and were more suitable than the coefficient of variation and the least significant difference (by the Tukey test) to evaluate experimental precision in trials with cowpea genotypes.

An in-house real-time polymerase chain reaction: standardisation and comparison with the Cobas Amplicor HBV monitor and Cobas AmpliPrep/Cobas TaqMan HBV tests for the quantification of hepatitis B virus DNA.

This study aimed to standardise an in-house real-time polymerase chain reaction (rtPCR) to allow quantification of hepatitis B virus (HBV) DNA in serum or plasma samples, and to compare this method with two commercial assays, the Cobas Amplicor HBV monitor and the Cobas AmpliPrep/Cobas TaqMan HBV test. Samples from 397 patients from the state of São Paulo were analysed by all three methods. Fifty-two samples were from patients who were human immunodeficiency virus and hepatitis C virus positive, but HBV negative. Genotypes were characterised, and the viral load was measure in each sample. The in-house rtPCR showed an excellent success rate compared with commercial tests; inter-assay and intra-assay coefficients correlated with commercial tests (r = 0.96 and r = 0.913, p < 0.001) and the in-house test showed no genotype-dependent differences in detection and quantification rates. The in-house assay tested in this study could be used for screening and quantifying HBV DNA in order to monitor patients during therapy.

An in-house real-time polymerase chain reaction: standardisation and comparison with the Cobas Amplicor HBV monitor and Cobas AmpliPrep/Cobas TaqMan HBV tests for the quantification of hepatitis B virus DNA

This study aimed to standardise an in-house real-time polymerase chain reaction (rtPCR) to allow quantification of hepatitis B virus (HBV) DNA in serum or plasma samples, and to compare this method with two commercial assays, the Cobas Amplicor HBV monitor and the Cobas AmpliPrep/Cobas TaqMan HBV test. Samples from 397 patients from the state of São Paulo were analysed by all three methods. Fifty-two samples were from patients who were human immunodeficiency virus and hepatitis C virus positive, but HBV negative. Genotypes were characterised, and the viral load was measure in each sample. The in-house rtPCR showed an excellent success rate compared with commercial tests; inter-assay and intra-assay coefficients correlated with commercial tests (r = 0.96 and r = 0.913, p < 0.001) and the in-house test showed no genotype-dependent differences in detection and quantification rates. The in-house assay tested in this study could be used for screening and quantifying HBV DNA in order to monitor patients during therapy.

Quality control of genotypes using heritability estimates of gene content at the marker.

Quality control filtering of single-nucleotide polymorphisms (SNPs) is a key step when analyzing genomic data. Here we present a practical method to identify low-quality SNPs, meaning markers whose genotypes are wrongly assigned for a large proportion of individuals, by estimating the heritability of gene content at each marker, where gene content is the number of copies of a particular reference allele in a genotype of an animal (0, 1, or 2). If there is no mutation at the marker, gene content has an additive heritability of 1 by construction. The method uses restricted maximum likelihood (REML) to estimate heritability of gene content at each SNP and also builds a likelihood-ratio test statistic to test for zero error variance in genotyping. As a by-product, estimates of the allele frequencies of markers at the base population are obtained. Using simulated data with 10% permutation error (4% actual error) in genotyping, the method had a specificity of 0.96 (4% of correct markers are rejected) and a sensitivity of 0.99 (1% of wrong markers are accepted) if markers with heritability lower than 0.975 are discarded. Checking of Mendelian errors resulted in a lower sensitivity (0.84) for the same simulation. The proposed method is further illustrated with a real data set with genotypes from 3534 animals genotyped for 50,433 markers from the Illumina PorcineSNP60 chip and a pedigree of 6473 individuals; those markers underwent very little quality control. A total of 4099 markers with P-values lower than 0.01 were discarded based on our method, with associated estimates of heritability as low as 0.12. Contrary to other techniques, our method uses all information in the population simultaneously, can be used in any population with markers and pedigree recordings, and is simple to implement using standard software for REML estimation. Scripts for its use are provided.