Saliva DNA quality and genotyping efficiency in a predominantly elderly population

Validated WGS and WES protocols proved saliva-derived gDNA as an equivalent to blood-derived gDNA for clinical and population genomic analyses

BACKGROUND

Whole exome sequencing (WES) and whole genome sequencing (WGS) have become standard methods in human clinical diagnostics as well as in population genomics (POPGEN). Blood-derived genomic DNA (gDNA) is routinely used in the clinical environment. Conversely, many POPGEN studies and commercial tests benefit from easy saliva sampling. Here, we evaluated the quality of variant call sets and the level of genotype concordance of single nucleotide variants (SNVs) and small insertions and deletions (indels) for WES and WGS using paired blood- and saliva-derived gDNA isolates employing genomic reference-based validated protocols.

METHODS

The genomic reference standard Coriell NA12878 was repeatedly analyzed using optimized WES and WGS protocols, and data calls were compared with the truth dataset published by the Genome in a Bottle Consortium. gDNA was extracted from the paired blood and saliva samples of 10 participants and processed using the same protocols. A comparison of paired blood-saliva call sets was performed in the context of WGS and WES genomic reference-based technical validation results.

RESULTS

The quality pattern of called variants obtained from genomic-reference-based technical replicates correlates with data calls of paired blood-saliva-derived samples in all levels of tested examinations despite a higher rate of non-human contamination found in the saliva samples. The F1 score of 10 blood-to-saliva-derived comparisons ranged between 0.8030-0.9998 for SNVs and between 0.8883-0.9991 for small-indels in the case of the WGS protocol, and between 0.8643-0.999 for SNVs and between 0.7781-1.000 for small-indels in the case of the WES protocol.

CONCLUSION

Saliva may be considered an equivalent material to blood for genetic analysis for both WGS and WES under strict protocol conditions. The accuracy of sequencing metrics and variant-detection accuracy is not affected by choosing saliva as the gDNA source instead of blood but much more significantly by the genomic context, variant types, and the sequencing technology used.

Pre-Analytical Factors Affecting Extracellular DNA in Saliva

Salivary DNA is widely used for genetic analyses because of its easy collection. However, its extracellular fraction in particular, similar to the extracellular DNA (ecDNA) in plasma, could be a promising biomarker for oral or systemic diseases. In contrast to genetics, the quantity of salivary ecDNA is of importance and can be affected by the pre-analytical processing of samples, but the details are not known. The aim of our study was to analyze the effects of centrifugation and freezing of saliva on the concentration of ecDNA in saliva. Fifteen healthy volunteers, free of any known systemic or oral diseases, were asked to collect unstimulated saliva samples. Aliquots were centrifuged at 1600× g and frozen or directly processed. The fresh or thawed cell-free saliva samples underwent subsequent centrifugation at 16,000× g. The supernatants were used for DNA isolation and quantification using fluorometry and real-time PCR. While freezing had minimal effects on the salivary ecDNA concentration, another centrifugation step decreased ecDNA considerably in both fresh and frozen samples (by 97.8% and 98.4%, respectively). This was mirrored in the quantitative PCR targeting a nuclear (decrease by 93.5%) and mitochondrial (decrease by 97.7%) ecDNA sequence. In conclusion, in this first study focusing on the technical aspects of salivary ecDNA quantitation, we show that, regardless of its subcellular origin, the concentration of ecDNA in saliva is mainly affected by additional centrifugation and not by the freezing of centrifuged cell-free saliva samples. This suggests that most salivary ecDNA likely is associated with cell debris and apoptotic bodies. Which fraction is affected by a particular disease should be the focus of further targeted studies.

A Genetic Lab-on-Phone Test for Point-of-Care Diagnostic of Lactose Intolerance near Patient and in less than 90 Minutes

BACKGROUND

The -13910 C/T single nucleotide polymorphism located within the MCM6 gene, an enhancer region located upstream of the lactase-phlorizin hydrolase gene, is associated with lactase persistence/non-persistence traits among the Caucasian population. The performance of a new point-of-care CE-IVD (In Vitro Diagnostic) marked isothermal lab-on-phone lactose intolerance assay, using crude samples, was assessed in comparison with Sanger sequencing using purified DNA, as reference method.

METHODS

The study was conducted following a non-probability sampling using direct buccal swab (n = 63) and capillary blood (n = 43) clinical samples from a total of 63 volunteers. A 3 × 3 confusion matrix/contingency table was used to evaluate the performance of the isothermal lab-on-phone lactose intolerance assay.

RESULTS

The isothermal lab-on-phone lactose intolerance assay successfully detected the -13910 C/T variant with a limit of detection of 5 cells/assay and demonstrated an overall accuracy of 98.41% (95% CI, 91.47%-99.96%) for buccal swab samples and 100% (95% CI, 91.19%-100%) for capillary blood, taking just 90 min from sample to result, with only 2 min hands-on.

CONCLUSIONS

The lab-on-phone pocket-sized assay displayed good performance when using direct buccal swab and capillary blood samples, enabling a low-cost, real-time, and accurate genotyping of the -13910 C/T region for the rapid diagnosis of primary lactose intolerance at point-of-care, which enables a prompt implementation of appropriate diet habits and/or intolerance therapies. To our knowledge, this is the first point-of-care genetic test for lactose intolerance to be made available on the market.

Prevalence and factors associated with optic disc grey crescent in the Primary Open-Angle African Ancestry Glaucoma Genetics (POAAGG) Study

Aim

To investigate the prevalence and factors associated with optic disc grey crescent (GC) in African Americans with glaucoma.

Methods

Stereo optic disc image features from subjects with glaucoma in the Primary Open-Angle African Ancestry Glaucoma Genetics Study were evaluated independently by non-physician graders and discrepancies adjudicated by an ophthalmologist. Risk factors for GC were evaluated by logistic regression models with intereye correlation accounted for by generalised estimating equations. Adjusted ORs (aORs) were generated.

Results

GC was present in 227 (15%) of 1491 glaucoma cases, with 57 (3.82%) bilateral and 170 (11.4%) unilateral. In multivariable analysis, factors associated with GC were younger age (aOR 1.27, 95% CI 1.11 to 1.43 for every decade younger in age, p=0.001), diabetes (aOR 1.46, 95% CI 1.09 to 1.96, p=0.01), optic disc tilt (aOR 1.84, 95% CI 1.36 to 2.48, p<0.0001), a sloping retinal region adjacent to the outer disc margin (aOR 2.37, 95% CI 1.74 to 3.32, p<0.0001) and beta peripapillary atrophy (aOR 2.32, 95% CI 1.60 to 3.37, p<0.0001). Subjects with GC had a lower mean (SD) value of the ancestral component q0 than those without GC (0.22 (0.15) vs 0.27 (0.20), p=0.001), consistent with higher degrees of African ancestry.

Conclusions

More than 1 in 10 glaucoma cases with African ancestry have GC, occurring more frequently in younger subjects, higher degrees of African ancestry and diabetes. GC was associated with several ocular features, including optic disc tilt and beta peripapillary atrophy. These associations should be considered when evaluating black patients with primary open-angle glaucoma.

The significance of amelogenin loci from toothpicks as forensic evidence for sex determination

Objective

DNA analysis has become the gold standard of all identification methods in forensic science. There are several primary and secondary sources of DNA samples in the field of forensics. Primary sources can be obtained directly from individuals, whereas secondary sources of DNA samples include items used by someone such as clothes, toothbrushes, and toothpicks. In Asian countries such as Korea, Japan, and Indonesia, toothpicks are often used, as the tip of toothpicks used to clear food particles in the interdental regions can be utilized to capture DNA samples from saliva. The use of saliva samples for DNA analysis in sex identification is based on short tandem repeats with the amelogenin (AMG) locus. The purpose of this study was to examine the concentration and purity of DNA on toothpicks as evidence in criminal case investigations, as well as its potential as sex determination material.

Methods

Eight subjects were instructed to clean the interproximal side of their posterior teeth using a toothpick. Each toothpick sample was kept for 0 days (as a control), 14 days, and 20 days. The purity and DNA concentration of each sample were determined through DNA examination. After determining the concentration and purity of DNA from each sample, electrophoresis of the AMG loci was performed for sex determination.

Results

This study showed that the average concentration of DNA on toothpicks ranged from 425.25 to 796.25 μg/ml, and the average purity of DNA ranged from 1.09 to 1.13 μg/ml. The AMG gene produces 112 and 116 base pair amplicons from the X and Y chromosomes.

Conclusion

Sex determination using DNA can be done using AMG loci, a protein found on the sex chromosomes (X and Y). The value of DNA concentration on toothpicks could be used to support forensic identification after 20 days at room temperature.

Establishing analytical validity of BeadChip array genotype data by comparison to whole-genome sequence and standard benchmark datasets

BACKGROUND

Clinical use of genotype data requires high positive predictive value (PPV) and thorough understanding of the genotyping platform characteristics. BeadChip arrays, such as the Global Screening Array (GSA), potentially offer a high-throughput, low-cost clinical screen for known variants. We hypothesize that quality assessment and comparison to whole-genome sequence and benchmark data establish the analytical validity of GSA genotyping.

METHODS

To test this hypothesis, we selected 263 samples from Coriell, generated GSA genotypes in triplicate, generated whole genome sequence (rWGS) genotypes, assessed the quality of each set of genotypes, and compared each set of genotypes to each other and to the 1000 Genomes Phase 3 (1KG) genotypes, a performance benchmark. For 59 genes (MAP59), we also performed theoretical and empirical evaluation of variants deemed medically actionable predispositions.

RESULTS

Quality analyses detected sample contamination and increased assay failure along the chip margins. Comparison to benchmark data demonstrated that > 82% of the GSA assays had a PPV of 1. GSA assays targeting transitions, genomic regions of high complexity, and common variants performed better than those targeting transversions, regions of low complexity, and rare variants. Comparison of GSA data to rWGS and 1KG data showed > 99% performance across all measured parameters. Consistent with predictions from prior studies, the GSA detection of variation within the MAP59 genes was 3/261.

CONCLUSION

We establish the analytical validity of GSA assays using quality analytics and comparison to benchmark and rWGS data. GSA assays meet the standards of a clinical screen although assays interrogating rare variants, transversions, and variants within low-complexity regions require careful evaluation.

Recruitment strategies and lessons learned from a large genetic study of African Americans

Genetic studies must enroll large numbers of participants to obtain adequate statistical power. Data are needed on how researchers can best use limited financial and practical resources to achieve these targets, especially in under-represented populations. This paper provides a retrospective analysis of the recruitment strategies for a large glaucoma genetics study in African Americans. The Primary Open-Angle African American Glaucoma Genetics study enrolled 10,192 African American subjects from the Philadelphia region. Major recruitment approaches included clinic enrollment from University of Pennsylvania (UPenn) sites, clinic enrollment from external sites, sampling of Penn Medicine Biobank (PMBB), and community outreach. We calculated the enrollment yield, cost per subject, and seasonal trends of these approaches. The majority (65%) of subject were enrolled from UPenn sites with an average cost of $133/subject. Over time, monthly case enrollment declined as the pool of eligible subjects was depleted. Expanding to external sites boosted case numbers ($129/subject) and the biobank provided additional controls at low cost ($5/subject), in large part due to the generosity of PMBB providing samples free of cost. Community outreach was costly with low return on enrollment ($978/subject for 220 subjects). Summer months (Jun-Aug) produced the highest recruitment yields (p<0.001). Genetic studies will benefit from a multi-pronged and culturally sensitive recruitment approach. In our experience, the biobank was most cost-effective for control enrollment, while recruitment from clinics (including expansion to new sites) was necessary to recruit fully phenotyped cases.

Detection of Spinal Muscular Atrophy Patients Using Dried Saliva Spots

Spinal muscular atrophy (SMA) is a lower motor neuron disease, once considered incurable. The main symptoms are muscle weakness and muscular atrophy. More than 90% of cases of SMA are caused by homozygous deletion of survival motor neuron 1 (SMN1). Emerging treatments, such as splicing modulation of SMN2 and SMN gene replacement therapy, have improved the prognoses and motor functions of patients. However, confirmed diagnosis by SMN1 testing is often delayed, suggesting the presence of diagnosis-delayed or undiagnosed cases. To enable patients to access the right treatments, a screening system for SMA is essential. Even so, the current newborn screening system using dried blood spots is still invasive and cumbersome. Here, we developed a completely non-invasive screening system using dried saliva spots (DSS) as an alternative DNA source to detect SMN1 deletion. In this study, 60 DSS (40 SMA patients and 20 controls) were tested. The combination of modified competitive oligonucleotide priming-polymerase chain reaction and melting peak analysis clearly distinguished DSS samples with and without SMN1. In conclusion, these results suggest that our system with DSS is applicable to SMA patient detection in the real world.

The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine, the hihi

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.

Collection of genetic data at scale for a nationally representative population: the UK Millennium Cohort Study

A DNA bank has been created from the Millennium Cohort Study (MCS) saliva samples. A total of 23,336 samples are available, from 9,259 cohort members (4,630 males and 4,629 females), 8,898 mothers and 5,179 fathers. There are 4,533 mother, child, father 'triads'. This paper describes the collection of the saliva samples from cohort members and their biological parents in the MCS. It analyses response rates and predictors of response, and details the DNA extraction, genotyping and imputation procedures performed on the data.

The Role of Genetic Ancestry as a Risk Factor for Primary Open-angle Glaucoma in African Americans

Purpose

POAG is the leading cause of irreversible blindness in African Americans. In this study, we quantitatively assess the association of autosomal ancestry with POAG risk in a large cohort of self-identified African Americans.

Methods

Subjects recruited to the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study were classified as glaucoma cases or controls by fellowship-trained glaucoma specialists. POAAGG subjects were genotyped using the MEGA Ex array (discovery cohort, n = 3830; replication cohort, n = 2135). Population structure was interrogated using principal component analysis in the context of the 1000 Genomes Project superpopulations.

Results

The majority of POAAGG samples lie on an axis between African and European superpopulations, with great variation in admixture. Cases had a significantly lower mean value of the ancestral component q0 than controls for both cohorts (P = 6.14-4; P = 3-6), consistent with higher degree of African ancestry. Among POAG cases, higher African ancestry was also associated with thinner central corneal thickness (P = 2-4). Admixture mapping showed that local genetic ancestry was not a significant risk factor for POAG. A polygenic risk score, comprised of 23 glaucoma-associated single nucleotide polymorphisms from the NHGRI-EBI genome-wide association study catalog, was significant in both cohorts (P < 0.001), suggesting that both known POAG single nucleotide polymorphisms and an omnigenic ancestry effect influence POAG risk.

Conclusions

In sum, the POAAGG study population is very admixed, with a higher degree of African ancestry associated with an increased POAG risk. Further analyses should consider social and environmental factors as possible confounding factors for disease predisposition.

Spectral fingerprinting to evaluate effects of storage conditions on biomolecular structure of filter-dried saliva samples and recovered DNA

Saliva has been widely recognized as a non-invasive, painless and easy-to-collect bodily fluid, which contains biomarkers that can be used for diagnosis of both oral and systemic diseases. Under ambient conditions, salivary biomarkers are subject to degradation. Therefore, in order to minimize degradation during transport and storage, saliva specimens need to be stabilized. The aim of this study was to investigate the feasibility of preserving saliva samples by drying to provide a shelf-stable source of DNA. Human saliva was dried on filters under ambient conditions using sucrose as lyoprotective agent. Samples were stored under different conditions, i.e. varying relative humidity (RH) and temperature. In addition to assessment of different cell types in saliva and their DNA contents, Fourier transform infrared spectroscopy (FTIR) was used to evaluate the effects of storage on biomolecular structure characteristics of saliva. FTIR analysis showed that saliva dried without a lyoprotectant exhibits a higher content of extended β-sheet protein secondary structures compared to samples that were dried with sucrose. In order to evaluate differences in characteristic bands arising from the DNA backbone among differently stored samples, principal component analysis (PCA) was performed, allowing a clear discrimination between groups with/without sucrose as well as storage durations and conditions. Our results indicated that saliva dried on filters in the presence of sucrose exhibits higher biomolecular stability during storage.

Contaminating DNA in human saliva alters the detection of variants from whole genome sequencing

Cells obtained from human saliva are commonly used as an alternative DNA source when blood is difficult or less convenient to collect. Although DNA extracted from saliva is considered to be of comparable quality to that derived from blood, recent studies have shown that non-human contaminating DNA derived from saliva can confound whole genome sequencing results. The most concerning complication is that non-human reads align to the human reference genome using standard methodology, which can critically affect the resulting variant genotypes identified in a genome. We identified clusters of anomalous variants in saliva DNA derived reads which aligned in an atypical manner. These reads had only short regions of identity to the human reference sequence, flanked by soft clipped sequence. Sequence comparisons of atypically aligning reads from eight human saliva-derived samples to RefSeq genomes revealed the majority to be of bacterial origin (63.46%). To partition the non-human reads during the alignment step, a decoy of the most prevalent bacterial genome sequences was designed and utilised. This reduced the number of atypically aligning reads when trialled on the eight saliva-derived samples by 44% and most importantly prevented the associated anomalous genotype calls. Saliva derived DNA is often contaminated by DNA from other species. This can lead to non-human reads aligning to the human reference genome using current alignment best-practices, impacting variant identification. This problem can be diminished by using a bacterial decoy in the alignment process.

Molecular Genetics and Functional Analysis Implicate CDKN2BAS1-CDKN2B Involvement in POAG Pathogenesis

The genes in the 9p21 locus (CDKN2B-AS1 & CDKN2B) are widely associated with Primary open-angle glaucoma (POAG). However, the functional importance of this locus in POAG pathogenesis is still unexplored. This study investigated the role of CDKN2BAS1-CDKN2B axis in POAG. We observed significant association of CDKN2B-AS1 SNP rs4977756 with POAG and its endophenotypic traits (vertical cup-disc ratio (p = 0.033) and central corneal thickness (p = 0.008)) by screening African American POAG cases (n = 1567) and controls (n = 1600). A luciferase reporter assay in Human embryonic kidney 293T (HEK293T) cells revealed that the region surrounding rs4977756 likely serves as a transcriptional repressor. siRNA-mediated knockdown of CDKN2B-AS1 in HEK293T cells and trabecular meshwork (TM) cells resulted in significantly increased expression of CDKN2B, which was also observed in human POAG ocular tissues. Pathway focused qRT-PCR gene expression analysis showed increased cellular senescence, TGFβ signaling and ECM deposition in TM cells after CDKN2B-AS1 suppression. In conclusion, we report that CDKN2B-AS1 may act as a regulator, and it could function by modulating the expression of CDKN2B. In addition, increase in CDKN2B levels due to CDKN2B-AS1 suppression may result in the senescence of trabecular meshwork cells leading to POAG pathogenesis.

The Identification of Admixture Patterns Could Refine Pharmacogenetic Counseling: Analysis of a Population-Based Sample in Mexico

Pharmacogenetic analysis has generated translational data that could be applied to guide treatments according to individual genetic variations. However, pharmacogenetic counseling in some mestizo (admixed) populations may require tailoring to different patterns of admixture. The identification and clustering of individuals with related admixture patterns in such populations could help to refine the practice of pharmacogenetic counseling. This study identifies related groups in a highly admixed population-based sample from Mexico, and analyzes the differential distribution of actionable pharmacogenetic variants. A subsample of 1728 individuals from the Mexican Genomic Database for Addiction Research (MxGDAR/Encodat) was analyzed. Genotyping was performed with the commercial PsychArray BeadChip, genome-wide ancestry was estimated using EIGENSOFT, and model-based clustering was applied to defined admixture groups. Actionable pharmacogenetic variants were identified with a query to the Pharmacogenomics Knowledge Base (PharmGKB) database, and functional prediction using the Variant Effect Predictor (VEP). Allele frequencies were compared with chi-square tests and differentiation was estimated by F_ST. Seven admixture groups were identified in Mexico. Some, like Group 1, Group 4, and Group 5, were found exclusively in certain geographic areas. More than 90% of the individuals, in some groups (Group 1, Group 4 and Group 5) were found in the Central-East and Southeast region of the country. MTRR p.I49M, ABCG2 p.Q141K, CHRNA5 p.D398N, SLCO2B1 rs2851069 show a low degree of differentiation between admixture groups. ANKK1 p.G318R and p.H90R, had the lowest allele frequency of Group 1. The reduction in these alleles reduces the risk of toxicity from anticancer and antihypercholesterolemic drugs. Our analysis identified different admixture patterns and described how they could be used to refine the practice of pharmacogenetic counseling for this admixed population.

A comparison of genetic sampling methodologies for candidate-gene analyses

Much of the recent gains in knowledge regarding the influence of patient genetics on medication pharmacokinetics (drug absorption, distribution, metabolism and elimination) how patients process medications) and pharmacodynamics (drug response) have been attributed to the technologic advances in genetic testing methodologies and the involvement of large clinical data sets and biobanks. Indeed, Genome Wide Association Studies (GWAS) and Phenome Wide Association Studies (PWAS) along with ever-evolving biomedical informatics techniques and the expansion of the -omics sciences (e.g., transcriptomics, metabolomics, proteomics) have brought about unprecedented advances in precision medicine. Although the simpler candidate-gene analysis technique is not considered cutting-edge, it is reliable and important to the translation of pharmacogenomic research and the advancement of precision medicine. Leveraging the knowledge of biological plausibility (i.e., genetic mutation → altered function of protein product → altered drug pharmacokinetics/dynamics) to appropriately select genes for inclusion, the candidate-gene analysis technique does not necessitate large patient cohorts nor extensive multi-gene genetic analysis arrays. It is often the ideal method for clinicians to begin evaluating whether genetic information might improve their pharmacologic treatment strategies for their patients. Having access to specific patient populations and expertise regarding their medical subspecialty, physician scientists can implement a candidate-gene analysis in small cohorts. Even with less than 100 patients, results can often be used to determine whether further investigation is warranted and to inform future studies. Herein, we present a comparison of select contemporary methodologies regarding collection, processing and genotype testing applicable to the efficient implementation of candidate-gene studies.

SNP located in an AluJb repeat downstream of TMCO1, rs4657473, is protective for POAG in African Americans

AIMS

To determine the association of single nucleotide polymorphisms (SNPs) downstream from the TMCO1 gene with primary open-angle glaucoma (POAG) in African Americans (AA).

METHODS

AA subjects were recruited for the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study from the Scheie Eye Institute and its satellite sites in Philadelphia. A region containing an AluJb repeat and seven SNPs, including rs4656461 near the TMCO1 gene, were PCR-Sanger sequenced from POAAGG cases (n=1537) and controls (n=1570). Association between POAG and SNPs near TMCO1 was investigated by logistic regression analysis. Phenotypic trait associations with these SNPs were assessed by analysis of variance. Electrophoretic mobility shift assay (EMSA) was performed to assess the affinity of human T-box 5 (TBX5) protein for a predicted binding motif in the TMCO1 region. Dual Luciferase assays were performed by transfecting recombinant plasmids containing the region surrounding the above SNPs in HEK293T and trabecular meshwork cells.

RESULTS

The SNP rs4657473 (C>T) was associated with POAG; the TT genotype was protective (OR 0.20, 95% CI 0.09 to 0.42; p<0.001). No significant associations were found between the TMCO1 variants and phenotypic traits. EMSA confirmed the affinity of TBX5 for a predicted binding motif containing TMCO1 SNP rs4657475. Luciferase assays demonstrated a regulatory function for the genomic region around SNP rs4656561, located within AluJb repeat.

CONCLUSION

Our results demonstrate that a SNP downstream of TMCO1, rs4657473, is associated with POAG in an AA population. Our studies suggest a regulatory role for the previously POAG-associated locus near the TMCO1 gene that may affect gene expression.

Impact of DNA source on genetic variant detection from human whole-genome sequencing data

BACKGROUND

Whole blood is currently the most common DNA source for whole-genome sequencing (WGS), but for studies requiring non-invasive collection, self-collection, greater sample stability or additional tissue references, saliva or buccal samples may be preferred. However, the relative quality of sequencing data and accuracy of genetic variant detection from blood-derived, saliva-derived and buccal-derived DNA need to be thoroughly investigated.

METHODS

Matched blood, saliva and buccal samples from four unrelated individuals were used to compare sequencing metrics and variant-detection accuracy among these DNA sources.

RESULTS

We observed significant differences among DNA sources for sequencing quality metrics such as percentage of reads aligned and mean read depth (p<0.05). Differences were negligible in the accuracy of detecting short insertions and deletions; however, the false positive rate for single nucleotide variation detection was slightly higher in some saliva and buccal samples. The sensitivity of copy number variant (CNV) detection was up to 25% higher in blood samples, depending on CNV size and type, and appeared to be worse in saliva and buccal samples with high bacterial concentration. We also show that methylation-based enrichment for eukaryotic DNA in saliva and buccal samples increased alignment rates but also reduced read-depth uniformity, hampering CNV detection.

CONCLUSION

For WGS, we recommend using DNA extracted from blood rather than saliva or buccal swabs; if saliva or buccal samples are used, we recommend against using methylation-based eukaryotic DNA enrichment. All data used in this study are available for further open-science investigation.

The association of mitochondrial DNA haplogroups with POAG in African Americans

Mitochondrial dysfunction has been implicated in the pathogenesis of primary open-angle glaucoma (POAG). However, the potential significance of mitochondrial DNA (mtDNA) haplogroups to POAG has not been evaluated in the overaffected African American population. To investigate the association of mtDNA haplogroups with POAG and its phenotypic characteristics, genotyping data from 4081 African American subjects (1919 cases and 2162 controls) was analyzed using 1293 positions on mtDNA. The overall frequency of mtDNA haplogroups in the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study cohort was 37% L3, 29% L2, 21% L1, 4% L0, and 10% non-African haplogroups (non-L). When all haplogroups (L0, L1, L2, and non-L) were compared against theL3 reference group, after adjusting by age and principal component of ancestry, the non-L3 haplogroups showed higher risk of POAG (OR-1.19, p = 0.02), with a particularly strong association among males (OR = 1.41, p = 0.003). More specifically the non-L group was associated with higher POAG risk than the L3 haplogroup (OR = 1.77, p = 0.007, Bonferroni adjusted p = 0.027) and to the L3e (n = 256, OR = 1.92, p = 0.007, Bonferroni adjusted p = 0.029). No significant association was found when genders were analyzed together or in female only analysis. There were no significant differences in various POAG endophenotypes across mtDNA haplogroups. This study expands our knowledge of mitochondrial genetics and mtDNA haplogroup associations in African American POAG. Further work is needed to better understand the functional role of mtDNA polymorphisms and their interactions with nuclear genes that affect POAG.

The utility of DNA extracted from saliva for genome-wide molecular research platforms

OBJECTIVE

The study aimed to investigate the suitability of DNA extracted from saliva for high throughput molecular genotyping and DNA methylation platforms by comparing its performance with that of DNA extracted from blood. The genome-wide methylation profile, using the Infinium HumanMethylation450 Beadchip array^® (Illumina, San Diego, CA), was measured for 20 DNA samples. Common genetic variation was measured, using the Infinium HumanCore Beadchip^® (Illumina, San Diego, CA) for 4 samples (matching samples from 2 people).

RESULTS

DNA from blood and saliva returned genotyping call rates and reproducibility frequencies of > 99%. High-quality DNA methylation data was obtained from both saliva and blood DNA, with average detection p-values for each sample ranging from 0.001 to 0.006. Slightly higher global DNA methylation levels were observed in whole blood DNA than saliva DNA. Correlations between individuals for each sample type were generally greater than correlations between two sample types from the same individual (Pearson's correlation, r = 0.9696 in 10 pairs of matched blood and saliva derived DNA, r = 0.9702 between saliva samples, and r = 0.9769 between blood derived DNA). Saliva yields DNA of sufficient quantity and quality to compare favourably with blood as a source of DNA for genetic and epigenetic research purposes.

Assessment of the stability of DNA in specimens collected under conditions for drug testing-A pilot study

For forensic biological sample collections, the specimen donor is linked solidly to his or her specimen through a chain of custody (CoC) sometimes referenced as a chain of evidence. Rarely, a donor may deny that a urine or oral fluid (OF) specimen is his or her specimen even with a patent CoC. The goal of this pilot study was to determine the potential effects of short-term storage on the quality and quantity of DNA in both types of specimen under conditions that may be encountered with employment-related drug testing specimens. Fresh urine and freshly collected oral fluid all produced complete STR profiles. For the "pad" type OF collectors, acceptable DNA was extractable both from the buffer/preservative and the pad. Although fresh urine and OF produced complete STR profiles, partial profiles were obtained after storage for most samples. An exception was the DNA in the Quantisal OF collector, from which a complete profile was obtained for both freshly collected OF and stored OF.

A Preliminary Genome-Wide Association Study of Pain-Related Fear: Implications for Orofacial Pain

Background

Acute and chronic orofacial pain can significantly impact overall health and functioning. Associations between fear of pain and the experience of orofacial pain are well-documented, and environmental, behavioral, and cognitive components of fear of pain have been elucidated. Little is known, however, regarding the specific genes contributing to fear of pain.

Methods

A genome-wide association study (GWAS; N = 990) was performed to identify plausible genes that may predispose individuals to various levels of fear of pain. The total score and three subscales (fear of minor, severe, and medical/dental pain) of the Fear of Pain Questionnaire-9 (FPQ-9) were modeled in a variance components modeling framework to test for genetic association with 8.5 M genetic variants across the genome, while adjusting for sex, age, education, and income.

Results

Three genetic loci were significantly associated with fear of minor pain (8q24.13, 8p21.2, and 6q26; p < 5 × 10⁻⁸ for all) near the genes TMEM65, NEFM, NEFL, AGPAT4, and PARK2. Other suggestive loci were found for the fear of pain total score and each of the FPQ-9 subscales.

Conclusions

Multiple genes were identified as possible candidates contributing to fear of pain. The findings may have implications for understanding and treating chronic orofacial pain.