Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

Design and Evaluation of a Robust CRISPR Kinetic Assay for Hot-Spot Genotyping.

Next-generation sequencing offers highly multiplexed and accurate detection of nucleic acid sequences but at the expense of complex workflows and high input requirements. The ease of use of CRISPR-Cas12 assays is attractive and may enable highly accurate detection of sequences implicated in, for example, cancer pathogenic variants. CRISPR assays often employ end-point measurements of Cas12 trans-cleavage activity after Cas12 activation by the target; however, end point-based methods can be limited in accuracy and robustness by arbitrary experimental choices. To overcome such limitations, we develop and demonstrate here an accurate assay targeting a mutation of the epidermal growth factor gene implicated in lung cancer (exon 19 deletion). The assay is based on characterizing the kinetics of Cas12 trans-cleavage to discriminate the mutant from wild-type targets. We performed extensive experiments (780 reactions) to calibrate key assay design parameters, including the guide RNA sequence, reporter sequence, reporter concentration, enzyme concentration, and DNA target type. Interestingly, we observed a competitive reaction between the target and reporter molecules that has important consequences for the design of CRISPR assays, which use preamplification to improve sensitivity. Finally, we demonstrate the assay on 18 tumor-extracted amplicons and 100 training iterations with 99% accuracy and discuss discrimination parameters and models to improve wild type versus mutant classification.

Continuous Global Improvement of Human Papillomavirus (HPV) Genotyping Services: The 2022 and 2023 HPV LabNet International Proficiency Studies.

The International Human Papillomavirus (HPV) Reference Center launches annual global HPV genotyping proficiency panels to enhance the precision and international standardization of HPV genotyping services. This study aims to assess the proficiency levels achieved in the global HPV genotyping proficiency panels conducted in 2022 and 2023, and to evaluate trends in performance over time. The proficiency panels comprised 44 blinded samples each, including 40 samples containing various purified plasmids corresponding to HPV types combined with human DNA, plus four control samples (one negative control and three extraction controls). Proficiency required a sensitivity of 50 International Units (IU)/5 µL for HPV 16 and HPV 18 500 IU/5 µL for HPVs 6, 11, 31, 33, 45, 52, and 58 and 500 genome equivalents (GE)/5 µL for other HPV types in both single and multiple infections, while avoiding false positivity. In 2022, 78 laboratories submitted a total of 154 data sets, and in 2023, 81 laboratories contributed 141 data sets. Most data sets (87%, 258/295) utilized commercially available HPV assays. Proficiency was common, with 77% of data sets meeting the proficiency criteria in 2022 and 79% in 2023. False positive results significantly decreased from 22% in 2022 to 13% in 2023. The high proficiency and increasing specificity in HPV genotyping services indicates progress toward more reliable HPV testing. High accuracy is crucial for supporting global efforts in HPV and cervical cancer elimination.

Ancestry-agnostic estimation of DNA sample contamination from sequence reads.

Detecting and estimating DNA sample contamination are important steps to ensure high-quality genotype calls and reliable downstream analysis. Existing methods rely on population allele frequency information for accurate estimation of contamination rates. Correctly specifying population allele frequencies for each individual in early stage of sequence analysis is impractical or even impossible for large-scale sequencing centers that simultaneously process samples from multiple studies across diverse populations. On the other hand, incorrectly specified allele frequencies may result in substantial bias in estimated contamination rates. For example, we observed that existing methods often fail to identify 10% contaminated samples at a typical 3% contamination exclusion threshold when genetic ancestry is misspecified. Such an incomplete screening of contaminated samples substantially inflates the estimated rate of genotyping errors even in deeply sequenced genomes and exomes. We propose a robust statistical method that accurately estimates DNA contamination and is agnostic to genetic ancestry of the intended or contaminating sample. Our method integrates the estimation of genetic ancestry and DNA contamination in a unified likelihood framework by leveraging individual-specific allele frequencies projected from reference genotypes onto principal component coordinates. Our method can also be used for estimating genetic ancestries, similar to LASER or TRACE, but simultaneously accounting for potential contamination. We demonstrate that our method robustly estimates contamination rates and genetic ancestries across populations and contamination scenarios. We further demonstrate that, in the presence of contamination, genetic ancestry inference can be substantially biased with existing methods that ignore contamination, while our method corrects for such biases.

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.

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.

Development and validation of a wart-associated human papilloma virus genotyping assay for detection of HPV in cutaneous warts.

Cutaneous warts are infectious disorders caused by human papillomavirus (HPV). A recent study revealed that the HPV genotype influences the natural course and response to treatment for plantar warts, suggesting that HPV genotyping could potentially be used to optimize wart treatment schemes. For this purpose, a wart-associated HPV genotyping assay was developed. The assay was subjected to an intensive validation process including, i.a., empiric determination of the annealing temperature, primer-probe optimization, evaluation of the analytical specificity and sensitivity, viral load quantification, and qualitative as well as quantitative analysis of intra-run repeatability and inter-run reproducibility. The newly developed assay was employed in a small-scale HPV genotyping study of wart biopsies (n = 50). The assay exhibited an analytical type-specific sensitivity and specificity of 100% (95% confidence interval [CI]: 83.9%-100%). The limit of quantification of the tested sequences corresponded to less than 17 viral copies/µl, while the limit of detection was less than 5 copies/µl. Very good to excellent agreements were gained between intra- and inter-run measurements (κ = 0.85-1.00) and coefficients of variation of the quantitative agreements were less then 3%. 22.5% (95% CI: 11%-39%) of the analyzed biopsies were negative for the tested HPV types, while 35% (95% CI: 21%-52%) contained multiple infections. The wart-associated HPV quantitative polymerase chain reaction assay was proven to be highly sensitive and specific. Multiple HPV infections were detected in 35% of lesions, contradicting the current literature claiming that in immunocompetent patients only 4%-16% of warts exhibit multiple HPV infections. This assay is qualified to be implemented in development of future genotype specific wart treatment strategies.

Standardization of a SNP panel for parentage verification and identification in the domestic cat (Felis silvestris catus).

The domestic cat (Felis silvestris catus) is a valued companion animal throughout the world. Over 60 different cat breeds are accepted for competition by the cat fancy registries in different countries. Genetic markers, including short tandem repeats and SNPs, are available to evaluate and manage levels of inbreeding and genetic diversity, population and breed structure relationships, and individual identification for forensic and registration purposes. The International Society of Animal Genetics (ISAG) hosts the Applied Genetics in Companion Animals Workshop, which supports the standardization of genetic marker panels and genotyping for the identification of cats via comparison testing. SNP panels have been in development for many species, including the domestic cat. An ISAG approved core panel of SNPs for use in cat identification and parentage analyses is presented. SNPs (n = 121) were evaluated by different university-based and commercial laboratories using 20 DNA samples as part of the ISAG comparison testing procedures. Different SNP genotyping technologies were examined, including DNA arrays, genotyping-by-sequencing and mass spectroscopy, to select a robust and efficient panel of 101 SNPs as the ISAG core panel for cats. The SNPs are distributed across all chromosomes including two on the X chromosome and an XY pseudo-autosomal sexing marker (zinc-finger XY; ZFXY). A population study demonstrated that the markers have an average polymorphic information content of 0.354 and a power of exclusion greater than 0.9999. The SNP panel should keep testing affordable while also allowing for the development of additional panels to monitor health, phenotypic traits, hybrid cats and highly inbred cats.

Transcriptome landscaping for gene mining and SSR marker development in Coriander (Coriandrum sativum L.).

Coriander (Coriandrum sativum L.) is an aromatic herb, widely used as a spice and is of great pharmaceutical interest. Despite high medicinal and economic value, there is a dearth of genomic information about profiling as well as the expressed sequence-based genic markers. In this study, transcriptome was sequenced from seeds, leaves, and flower for gene mining and identification of SSR markers. A total of 9746 SSR containing loci were identified, the most abundant type of SSR identified were the di-nucleotide repeat motifs (45.5%), followed by tri- (34.6%), tetra- (4.5%), penta- (1.5%) and hexanucleotide repeats (1%). A total of 3795 primers were designed, out of which 120 randomly selected were validated in 14 accessions of coriander cultivated in India. The current study provides useful information about preliminary transcriptome sketch and genic markers, which can be useful in breeding and genetic diversity estimation of coriander.

Development and characterization of non-coding RNA based simple sequence repeat markers in Capsicum species.

Plant growth and development are largely regulated by non-coding RNAs (ncRNA); thus ncRNA based markers would be rewarding in molecular breeding. In the present study, for the first time we developed total 623 ncRNA based SSRs including 119 microRNASSRs (miRNASSRs) and 504 long non-coding RNASSRs (lncRNASSRs) distributed across 12 Capsicum chromosomes. Out of 623 ncRNASSRs, 120 (including 60 each miRNASSRs and lncRNASSRs) were used for genotyping of 96 Capsicum accessions belonging to C. annuum, C. chinense and C. frutescens; and 75% SSRs were polymorphic. Model-based and distance-based cluster analyses identified three species specific clusters, i.e. cluster-I (C. annuum), cluster-II (C. frutescens) and cluster-III (C. chinense); therefore, these SSRs may have a potential role to play in interspecific Capsicum breeding. Tissue specific expression of SSR containing ncRNAs and versatile functions of their targets suggested the usefulness of SSRs for mapping of genes/QTLs and breeding of wide range of traits in Capsicum.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction.

BACKGROUND: Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes. RESULTS: We generated 96.7 × 109 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 109 paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females. CONCLUSIONS: The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Strategies for processing and quality control of Illumina genotyping arrays.

Illumina genotyping arrays have powered thousands of large-scale genome-wide association studies over the past decade. Yet, because of the tremendous volume and complicated genetic assumptions of Illumina genotyping data, processing and quality control (QC) of these data remain a challenge. Thorough QC ensures the accurate identification of single-nucleotide polymorphisms and is required for the correct interpretation of genetic association results. By processing genotyping data on > 100 000 subjects from >10 major Illumina genotyping arrays, we have accumulated extensive experience in handling some of the most peculiar scenarios related to the processing and QC of Illumina genotyping data. Here, we describe strategies for processing Illumina genotyping data from the raw data to an analysis ready format, and we elaborate on the necessary QC procedures required at each processing step. High-quality Illumina genotyping data sets can be obtained by following our detailed QC strategies.

Clinical performance of the HPV DNA Array genotyping assay in detection of CIN2+ lesions with BS GP5+/6+ MPG Luminex tested cervical samples.

Human papillomavirus (HPV) detection is used for screening of cervical cancer and genotype-specific persistence has shown to be mandatory for dysplasia development. Aim of this study was to evaluate the clinical performance of HPV DNA Array for cervical intraepithelial neoplasia 2+ (CIN2+) lesion detection. HPV DNA Array is a polymerase chain reaction-based assay that targets E1 sequences of 29 HPV types (6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 44, 45, 51, 52, 53, 54, 56, 58, 59, 66, 67, 68, 69, 70, 73, 82, 85, and 97). The clinical evaluation was performed against the reference assay, BS-GP5+/6+ multiplex genotyping (MPG)-Luminex, with 600 cervical smear samples of a referral population. HPV DNA Array detected CIN2+ lesions with a sensitivity of 90.2%, identical to that of MPG-Luminex. Detection of CIN3+ lesions was with a sensitivity of 90.3%, as compared with 88.7% of MPG-Luminex. It demonstrated very good agreement for HPV detection, irrespective of type, of 91.5% (κ = 0.832). HPV DNA Array is a simple and robust assay, with a short protocol of 4 hours hands-on time and automated readout by ELISpot AiDot software. It permits testing of up to 96 samples in one run and may be considered for use in organized screening programs and low resource settings.

Evaluation of the SureX HPV genotyping test for the detection of high-risk HPV in cervical cancer screening.

BACKGROUND: The SureX HPV genotyping test (SureX HPV test), which targets the human papillomavirus (HPV) E6/E7 genes was compared with the Cobas 4800 and Venus HPV tests for detecting 14 high-risk HPV (HR-HPV) types in clinical referral and follow-up patients to evaluate its value for cervical cancer screening. METHODS: Two different populations were enrolled in the study. The first population comprised 185 cases and was used for comparing the SureX HPV test (Health, China) with the Cobas 4800 test (Roche, USA). The second population comprised 290 cases and was used for comparing the SureX HPV test (Health, China) with the Venus HPV test (Zhijiang, China). Polymerase chain reaction (PCR) sequencing was performed for further confirmation of discordant results. RESULTS: In the first population, the overall agreement rate was 95.6% for 14 high-risk HPV types. Eight discordant cases were confirmed by PCR sequencing, which showed that the agreement rates were 75.0% between the SureX HPV test and PCR sequencing and 25.0% between the Cobas 4800 test and PCR sequencing (P < 0.01). In the second population, the overall agreement rate was 95.5%. Thirteen discordant cases were confirmed by PCR sequencing, which showed that the agreement rates were 76.9% between the SureX HPV test and PCR sequencing and 23.1% between the Venus HPV test and PCR sequencing (P < 0.01). With cervical intraepithelial neoplasia grade 2+ (CIN2+) as the reference standard, the sensitivity values of the SureX HPV test and the Venus HPV test were 93.5% and 92.0%, (P > 0.05), while the specificity values were 43.3% and 46.7%, respectively (P > 0.05). CONCLUSION: The SureX HPV test had good consistency with both the Cobas 4800 and Venus HPV tests for 14 HR-HPV types. In addition, it avoided some false negatives and false positives. Therefore, the SureX HPV test can be used for cervical cancer screening.

External quality assessment of noninvasive fetal RHD genotyping.

BACKGROUND AND OBJECTIVES: Fetal RHD genotyping of cell-free maternal plasma DNA from RhD negative pregnant women can be used to guide targeted antenatal and postnatal anti-D prophylaxis for the prevention of RhD immunization. To assure the quality of clinical testing, we conducted an external quality assessment workshop with the participation of 31 laboratories. MATERIALS AND METHODS: Aliquots of pooled maternal plasma from gestational week 25 were sent to each laboratory. One sample was fetal RHD positive, and a second sample was fetal RHD negative. A reporting scheme was supplied for data collection, including questions regarding the methodological setup, results and clinical recommendations. The samples were tested blindly. RESULTS: Different methodological approaches were used; 29 laboratories used qPCR and two laboratories used ddPCR, employing a total of eight different combinations of RHD exon targets. Fetal RHD genotyping was performed with no false-negative and no false-positive results. One inconclusive result was reported for the RHD positive sample. All clinical conclusions were satisfactory. CONCLUSION: This external quality assessment workshop demonstrates that despite the different approaches taken to perform the clinical assays, fetal RHD genotyping is a reliable laboratory assay to guide targeted use of Rh prophylaxis in a clinical setting.

Genetic diversity analysis for narrow-leafed lupin (Lupinus angustifolius L.) by SSR markers.

Narrow-leafed lupin (Lupinus angustifolius L.) is used as grain legumes, fodder for livestock and green manure in the world and has a great potential to be developed as a new crop in China. In this study, we assessed the genetic diversity among a set of 109 newly introduced accessions of narrow-leafed lupin using 76 genomic SSR markers. Data analysis suggested that the average gene diversity index and average polymorphism information content (PIC) were 0.4758 and 0.4328, respectively. The mean allele number per loci (Na) was 6.3816. The population structure analysis identified two subgroups based on delta K (ΔK) values. This result is in accordance with that of a PCA. The AMOVA analysis showed that most of molecular variance were within population. These results will be useful to guide the genetic improvement of the narrow-leafed lupin crop in China.

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.

Intraspecific mitochondrial genome comparison identified CYTB as a high-resolution population marker in a new pest Athetis lepigone.

A new outbreak pest, Athetis lepigone (Möschler) (Lepidoptera: Noctuidae), has caused severe economic loss in maize crops in China. In order to conduct population genetics study with a more polymorphic and scientific mitochondrial marker, we sequenced the complete mitochondrial genomes of 13 different A. lepigone individuals. Intraspecific comparison of all PCGs showed that the NADH dehydrogenase and cytochrome b genes had the highest nucleotide diversity. We also found evidence of episodic positive selection on two amino acids, which are encoded by NADH dehydrogenase genes (ND3 and ND4L), against a background of widespread neutral selection of all other mitochondrial PCGs. The genetic divergence observed in this study indicated that the cytochrome b gene (CYTB) is better than COI at recovering population structure. The preliminary population genetic analysis illustrated strong gene flow among A. lepigone populations in China. Our study provides basic information for further research on population genetics of A. lepigone.

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project.

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed "consensus calling," to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.

Quality and concordance of genotyping array data of 12,064 samples from 5840 cancer patients.

Genotyping arrays characterize genome-wide SNPs for a study cohort and were the primary technology behind genome wide association studies over the last decade. The Cancer Genome Atlas (TCGA) is one of the largest cancer consortium studies, and it collected genotyping data for all of its participants. Using TCGA SNP data genotyped using the Affymetrix 6.0 SNP array from 12,064 samples, we conducted a comprehensive comparisons across DNA sources (tumor tissue, normal tissue, and blood) and sample storage protocols (formalin-fixed paraffin-embedded (FFPE) vs. freshly frozen (FF)), examining genotypes, transition/transversion ratios, and mutation catalogues. During the analysis, we made important observations in relevance to the data quality issues. SNP concordance was excellent between blood and normal tissues, and slightly lower between blood and tumor tissue due to potential somatic mutations in the tumors. The observed poor SNP concordance between FFPE and FF samples suggested a batch effect. The transition/transversion ratio, a metric commonly used for quality control purpose in exome sequencing projects, appeared less applicable for genotyping array data due to the whole-genome coverage built into the array design. Moreover, there were substantially more loss of heterozygosity events than gain of heterozygosity when comparing tumors relative to normal tissues and blood. This might be a consequence of extensive copy number deletions in tumors. In summary, our thorough evaluation calls for more adequate quality control practices and provides guidelines for improved application of TCGA genotyping data.