Human Nucleotide Research Articles

Estimation of translational importance of mammalian mRNA nucleotide sequence characteristics based on ribosome profiling data

It is known that the 5’ untranslated region (5’ UTR) mRNA characteristics can influence translation initiation efficiency and specificity. Previous knowledge about 5’ UTR characteristics was obtained theoretically and in vitro for mRNA of individual genes. It did not allow systematic analysis of mRNA translationally important parameters. To identify the above mentioned 5’ UTR characteristics, it is necessary to analyze their relationships with the translational activity of the corresponding mRNAs. Until recently, there were no experimental data on translation efficiency. Thanks to ribosome profiling technology, genome-wide experimental data of translation efficiency have been obtained for many eukaryotic mRNAs. Now it seems to be possible to reveal translationally important mRNA parameters and predict translation efficiency based on their nucleotide sequences. The aim of this study was to determine the translational significance of individual 5’ UTR characteristics in accordance with experimental ribosome profiling data. A statistical analysis was carried out for revealing relationships between the human and mouse mRNA nucleotide sequence characteristics and ribosome profiling data. Some of the mRNA parameters influencing translation efficiency were most significant, and the same trends for all three samples analyzed were revealed: a purine at start codon context position –3, upstream AUG presence and G+C complementary nucleotide concentration reduce translation efficiency; whereas gexonucleotides CCGCCA (5’ UTR) and AAGAAA, AAGAAG, AAGCAG, AAAAAG (CDS) increase translation efficiency. A toolkit that allows analyzing the importance of 5’ UTR characteristics and a program for prediction of translation efficiency were developed on the base of the BioUML platform.

Comparison of the expression of human equilibrative nucleotide transporter 1 (hENT1) and ribonucleotide reductase subunit M1 (RRM1) genes in seven non-Hodgkin lymphoma cell lines.

We investigated the variability in the expression of human equilibrative nucleoside transporter 1 (hENT1) and ribonucleotide reductase subunit M1 (RRM1) in non-Hodgkin lymphoma cell lines. hENT1 and RRM1 mRNA expression levels in natural killer (NK) cells and seven non-Hodgkin lymphoma cell lines (YTS, SNK-6, Jeko-1, ly-1, Raji, Karpas, and Jurket) were studied using reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR) and the results were compared using the Student t-test. mRNA expression of hENT1 was detectable in YTS, SNK-6, Jeko-1, ly-1, Raji, Karpas, Jurket, and NK cells, which revealed variability in gene expression. There were significant differences in the mRNA expression values of hENT1 (P = 0.021) and RRM1 (P = 0.002) compared to those in NK cells. mRNA expression of both hENT1 and RRM1 was closely associated with non-Hodgkin lymphoma cell proliferation. Differential expression analysis of hENT1 and RRM1 in non-Hodgkin lymphoma cell lines may provide novel drug leads for precision medicine.

Genetic and epigenetic abnormalities in systemic sclerosis.

Systemic sclerosis (SSc) is a complicated autoimmune and connective tissue disease, of which the pathogenesis and treatment have not been fully elucidated. SSc patients may have a genetic predisposition, such as specific human leukocyte antigens and single nucleotide polymorphisms, but can also develop various clinical symptoms with individual differences. Epigenetics, such as DNA methylation, histone modification, long non-coding RNA and miRNA, could explain the crossroad between genetics and environmental factors. For instance, epigenetics is associated with environmental factors, which may cause the breakdown of immune tolerance and the development of SSc in patients with a particular genetic background. In the future, further investigations of the interplay between genetics and epigenetics will be beneficial to elucidate the complex molecular cross-talk and heterogeneity in the SSc pathogenesis. Furthermore, this research will lead to the discovery of new therapeutic approaches and biomarkers.

New one-pot synthesis of N-fused isoquinoline derivatives by palladium-catalyzed C-H arylation: potent inhibitors of nucleotide pyrophosphatase-1 and -3.

Various N-fused isoquinoline derivatives were synthesized using a new one-pot reaction of 1-bromo-2-(2,2-difluorovinyl)benzenes with N-H group containing heterocycles followed by intramolecular palladium-catalyzed C-H arylation. The method described gives convenient access to diverse structures of N-fused polycyclic isoquinolines. Sixteen of the synthesized compounds were screened as potential human nucleotide pyrophosphatase/phosphodiesterase 1 and 3 (h-NPP-1 and h-NPP-3) inhibitors. The most effective h-NPP-1 inhibitor showed an IC50 value as high as 0.36 ± 0.06 μM, whereas the most potent h-NPP-3 inhibitor posessed an inhibitory value of 0.48 ± 0.01 μM. Kinetic and molecular docking studies of both most effective inhibitors were carried out.

Brief Report: Common Genetic Variation in Chromosome 10 q22.1 Shows a Strong Sex Bias in Human Embryonic Stem Cell Lines and Directly Controls the Novel Alternative Splicing of Human NODAL which is Associated with XIST Expression in Female Cell Lines.

The potential use of pluripotent stem cells for personalized regenerative medicine necessitates an improved understanding of how germ-line genetic variation may affect pluripotency. Given previous reports of a female bias in established human embryonic stem cell (hESC) lines, sex-specific differences must also be considered. Herein we describe, for the first time, how genetic polymorphisms may affect the establishment of widely used hESC lines. We demonstrate that the minor allele of the human single nucleotide polymorphism (SNP) rs2231947 found within the NODAL gene locus is under-represented in male but not female hESC lines. We also show that this SNP is highly functional in hESC lines. The SNP rs2231947 directly controls the alternative splicing of a novel NODAL transcript isoform. Thus we demonstrate that genetic variation drastically affects the expression of a gene that plays a major role in the regulation of pluripotency and cell fate. Our work helps detail how genetic heterogeneity is manifested in hESC biology and highlights the need to identify how specific genetic variants can explain important differences between pluripotent cell line models both within and between species.

Human ISPD Is a Cytidyltransferase Required for Dystroglycan O-Mannosylation

A unique, unsolved O-mannosyl glycan on α-dystroglycan is essential for its interaction with proteinligands in the extracellular matrix. Defective O-mannosylation leads to a group of muscular dystrophies, called dystroglycanopathies. Mutations in isoprenoid synthase domain containing (ISPD) represent the second most common causeof these disorders, however, its molecular function remains uncharacterized. The human ISPD (hISPD) crystal structure showed a canonical N-terminal cytidyltransferase domain linked to a C-terminal domain that is absent in cytidyltransferase homologs. Functional studies demonstrated cytosolic localization of hISPD, and cytidyltransferase activity toward pentose phosphates, including ribulose5-phosphate, ribose 5-phosphate, and ribitol 5-phosphate. Identity of the CDP sugars was confirmed by liquid chromatography quadrupole time-of-flight mass spectrometry and two-dimensional nuclear magnetic resonance spectroscopy. Our combined results indicate that hISPDis a cytidyltransferase, suggesting the presence of a novel human nucleotide sugar essential for functional α-dystroglycan O-mannosylation in muscle and brain. Thereby, ISPD deficiency can be added to the growing list of tertiary dystroglycanopathies.

Analysis of Ribonucleotide Removal from DNA by Human Nucleotide Excision Repair

Ribonucleotides are incorporated into the genome during DNA replication. The enzyme RNase H2 plays a critical role in targeting the removal of these ribonucleotides from DNA, and defects in RNase H2 activity are associated with both genomic instability and the human autoimmune/inflammatory disorder Aicardi-Goutières syndrome. Whether additional general DNA repair mechanisms contribute to ribonucleotide removal from DNA in human cells is not known. Because of its ability to act on a wide variety of substrates, we examined a potential role for canonical nucleotide excision repair in the removal of ribonucleotides from DNA. However, using highly sensitive dual incision/excision assays, we find that ribonucleotides are not efficiently targeted by the human nucleotide excision repair system in vitro or in cultured human cells. These results suggest that nucleotide excision repair is unlikely to play a major role in the cellular response to ribonucleotide incorporation in genomic DNA in human cells.

Accelerated reproductive aging in females lacking a novel centromere protein SYCP2L.

Menopause results from loss of ovarian function and marks the end of a woman's reproductive life. Alleles of the human SYCP2L locus are associated with age at natural menopause (ANM). SYCP2L is a paralogue of the synaptonemal complex protein SYCP2 and is expressed exclusively in oocytes. Here we report that SYCP2L localizes to centromeres of dictyate stage oocytes, which represent the limited pool of primordial oocytes that are formed perinatally and remain arrested till ovulation. Centromere localization of SYCP2L requires its C-terminal portion, which is missing in truncated variants resulting from low-frequency nonsense mutations identified in humans. Female mice lacking SYCP2L undergo a significantly higher progressive loss of oocytes with age compared with wild-type females and are less fertile. Specifically, the pool of primordial oocytes becomes more rapidly depleted in SYCP2L-deficient than in wild-type females, such that with aging, fewer oocytes undergo maturation in developing follicles. We find that a human SYCP2L intronic single nucleotide polymorphism (SNP) rs2153157, which is associated with ANM, changes the splicing efficiency of U12-type minor introns and may therefore regulate the steady-state amount of SYCP2L transcript. Furthermore, the more efficiently spliced allele of this intronic SNP in SYCP2L is associated with increased ANM. Our results suggest that SYCP2L promotes the survival of primordial oocytes and thus provide functional evidence for its association with ANM in humans.

Computational and Functional Analysis of the Virus-Receptor Interface Reveals Host Range Trade-Offs in New World Arenaviruses.

Animal viruses frequently cause zoonotic disease in humans. As these viruses are highly diverse, evaluating the threat that they pose remains a major challenge, and efficient approaches are needed to rapidly predict virus-host compatibility. Here, we develop a combined computational and experimental approach to assess the compatibility of New World arenaviruses, endemic in rodents, with the host TfR1 entry receptors of different potential new host species. Using signatures of positive selection, we identify a small motif on rodent TfR1 that conveys species specificity to the entry of viruses into cells. However, we show that mutations in this region affect the entry of each arenavirus differently. For example, a human single nucleotide polymorphism (SNP) in this region, L212V, makes human TfR1 a weaker receptor for one arenavirus, Machupo virus, but a stronger receptor for two other arenaviruses, Junin and Sabia viruses. Collectively, these findings set the stage for potential evolutionary trade-offs, where natural selection for resistance to one virus may make humans or rodents susceptible to other arenavirus species. Given the complexity of this host-virus interplay, we propose a computational method to predict these interactions, based on homology modeling and computational docking of the virus-receptor protein-protein interaction. We demonstrate the utility of this model for Machupo virus, for which a suitable cocrystal structural template exists. Our model effectively predicts whether the TfR1 receptors of different species will be functional receptors for Machupo virus entry. Approaches such at this could provide a first step toward computationally predicting the "host jumping" potential of a virus into a new host species. We demonstrate how evolutionary trade-offs may exist in the dynamic evolutionary interplay between viruses and their hosts, where natural selection for resistance to one virus could make humans or rodents susceptible to other virus species. We present an algorithm that predicts which species have cell surface receptors that make them susceptible to Machupo virus, based on computational docking of protein structures. Few molecular models exist for predicting the risk of spillover of a particular animal virus into humans or new animal populations. Our results suggest that a combination of evolutionary analysis, structural modeling, and experimental verification may provide an efficient approach for screening and assessing the potential spillover risks of viruses circulating in animal populations.

Analysis of the correlation of CATSPER single nucleotide polymorphisms (SNPs) with idiopathic asthenospermia.

Idiopathic asthenospermia is the most common type of male infertility. Although the mechanisms causing asthenospermia are complex, recent studies have indicated an important role of cation channel of sperm (CATSPER) gene downregulation or abnormality in the etiology of idiopathic asthenospermia. In the present study, 192 patients with idiopathic asthenospermia and 288 healthy controls were enrolled, and a flight mass spectrometry using Sequenom's MassArray biochip system was applied for genotyping 16 CATSPER gene SNPs reported in the human single nucleotide polymorphism (SNP) database. Our results indicated a correlation between CATSPER1 SNPs and idiopathic asthenospermia. In particular, the exonal SNP rs1893316 in CATSPER1 significantly correlated with idiopathic asthenospermia risk and is a potential important factor in determining an individual's genetic susceptibility to idiopathic asthenospermia. These finding will help to further elucidate the role of CATSPER1 in idiopathic asthenospermia pathogenesis.

A model for integrating molecular-based testing in transfusion services.

Molecular-based laboratory tests can predict blood group antigens and supplement serological methods, adding a unique technology to assist in resolving discrepant or incomplete blood group typing or antibody identification. Hospital transfusion services have options for integrating molecular-based methods in their routine operations. We describe here the model of a hospital-reference laboratory partnership. Blood samples for compatibility testing were obtained from patients in a 609-bed hospital serving an urban multiethnic and multiracial population. When results of blood group phenotyping by serological methods were inconclusive, samples were referred for molecular-based testing. The reference laboratory used several methods for genotyping, including polymerase chain reaction followed by restriction enzyme-linked polymorphism analysis, sequence-specific primer polymerase chain reaction and array-based approaches. Human erythrocyte antigen, RHCE and RHD single nucleotide polymorphism arrays were integrated into the laboratory as they became commercially available. The hospital-reference laboratory model made it possible to integrate blood group genotyping promptly by current technology without the expense of new laboratory equipment or adding personnel with technical expertise. We describe ten cases that illustrate the categories of serological problems that were resolved by molecular methods. In-hospital molecular testing for transfusion services has logistical advantages, but is financially impractical for most hospitals. Our model demonstrates the advantages of a hospital-reference laboratory partnership. In conclusion, hospital transfusion services can integrate molecular-based testing in their routine services without delay by establishing a partnership with a molecular blood group reference laboratory. The hospital reference-laboratory model promotes genomic medicine without the expense of new equipment and skilled personnel, while supporting the economy of centralised large-scale laboratory operations.

Deep sequencing of the TCRβ repertoire reveals T cell clonal expansion is associated with HTLV-1 proviral load in HAM/TSP patients

HTLV-1 associated myelopathy/tropical spastic paraparesis (HAM/TSP) is associated with a chronic inflammatory central nervous system in which the host immune response against HTLV-1 infection is considered immunopathogenic. We have applied a new unbiased deep sequencing method of TCR repertoire to accurately measure the diversity and degree of clonal expansion of peripheral circulating T-cells in patients with HAM/TSP compared to age, gender, and ethnicity matched healthy controls. TCR libraries were generated using 5'rapid amplification of cDNA ends to amplify the V-D-J genes of unsorted T-cells obtained from peripheral blood. Over 100 million reads were analyzed using MiGec software (Shugay, et al., Nat Methods. 2014 Jun;11(6):653-5), that align and match the human TCR β-chain nucleotide sequences through IMGT database. T-cell clonotypes with high-quality CDR3 sequence were selected by a molecular identifier PCR strategy, which allowed for filtering PCR errors and correct for artificial sequences generated on a HiSeq 2000 Illumina system platform. By calculating the coefficient of variation of our dataset we considered clonal expansion as those unique TCR clonotypes showing reads >8. While we did not observe any significant statistical differences of the number of single clones (clonotypes with 1 read, P=0.62) between HAM/TSP patients and healthy controls groups, by contrast, a higher clonal expansion of the T-cell repertoire was observed in HAM/TSP patients when compared to healthy controls (unpaired parametric T-test with P=0.026). Interestingly, this clonal expansion of the HAM/TSP TCR repertoire correlated with the HTLV-I provirus load (as defined by digital PCR quantification of the HTLV-I tax gene) (r=0.64, P=0.008). Our findings strongly suggest that in this chronic neurological disorder, a strong dysregulated T-cell response may be driven by HTLV-I and that strategies to decrease virus may be of clinical benefit.

Prognostic Value of Excision Repair Cross-Complementing Gene 1, Dihydropyrimidine Dehydrogenase, and Human Equilibrative Nucleotide Transporter 1 Expression and Their Implications for Adjuvant Treatment in Patients With Ampullary Carcinoma.

The purpose of this study was to characterize the intratumoral expression profiles of excision repair cross-complementing gene 1 (ERCC1), dihydropyrimidine dehydrogenase (DPD), and human equilibrative nucleotide transporter 1 (hENT1) in ampullary carcinomas (ACs) to evaluate their prognostic values and better tailor adjuvant chemotherapy for individual patients with AC after surgery. This study included 49 patients with AC who underwent a curative pancreaticoduodenectomy. Various clinicopathological factors, including ERCC1, DPD, and hENT1, were analyzed in relation to postoperative disease recurrence and the patients' survival. The median recurrence-free survival and overall survival were 24.5 months and 32.4 months, respectively. Multivariate Cox regression analysis of recurrence-free survival identified a DPD expression (hazard ratio [HR], 8.18; 95% confidence interval [CI], 2.00-34.8; P = 0.003) and combined ERCC1/DPD expression (HR, 134.8; 95% CI, 11.8-1920; P < 0.001) as independent predictors of disease recurrence. Multivariate Cox regression analysis of overall survival also identified a DPD expression (HR, 8.48; 95% CI, 1.71-46.3; P = 0.008) and combined ERCC1/ DPD expression (HR, 135.6; 95% CI, 11.8-1940; P < 0.001) as independent predictors of survival. The DPD and ERCC1 expression profile could potentially serve as a useful prognostic biomarker and therapeutic target for surgically resected patients with AC.

Abstract 1417: Clonal succession in pancreatic cancer progression is not driven by genetic instability

Abstract Our group has recently shown that human pancreatic cancer (PDAC) progression is driven by a succession of transiently active tumor initiating cell (TIC) clones during serial xenotransplantation. Genetic labeling demonstrated that serial PDAC xenograft tumors and even tumors of parallel mice transplanted with cells from the same donor xenografts harbored very little to no overlap of active TIC clones, indicating substantial changes in the proliferative activity of individual TIC predominantly producing progeny without detectable tumor-initiating activity. We now asked whether observed clonal activation and inactivation is caused by acquisition of de novo mutations during evolution of genetic subclones or by functional plasticity of genetically stable TIC clones. Therefore, we monitored somatic non-synonymous mutations in culture and during PDAC progression in genetically marked serial xenografts of two patients. DNA was isolated from xenografts, primary TIC cultures and corresponding normal pancreas or primary tumor tissue. Following paired-end exome sequencing, reads were aligned to a concatenated hs37d5 human and mm10 mouse genome assembly and human specific single nucleotide variants (SNVs) and small insertions/deletions (indels) were identified. We found a total of 45 altered gene coding genomic loci (P1 = 10; P2 = 35) not present in control tissue. Strikingly, most SNVs detected were present in all samples, only very few SNV were acquired during serial transplantation. In P1, 4 novel SNVs not present in the original patient tumor sample were detected within coding regions of TTC13, OR4K15, SSPO and TPGS1. Allele frequencies ranged from 2-27% in serial xenografts. In xenografts of P2 we detected 35 SNVs not present in healthy tissue. Of these, one mutation in the gene C10orf12 aroused after serial transplantation with a maximum altered allele frequency of 17%. None of these mutations is a known cancer driver or was found as recurrent in large scale cancer sequencing approaches. To evaluate whether the clonal TIC dynamics within established tumors are recapitulated in vitro, we analysed individual TIC clone kinetics in serially passaged cultures and in cultures derived from transduced xenografts. Strikingly, the kinetics in vitro were similar to those observed within serially transplanted xenografts. Every culture passage was formed by a distinct set of actively proliferating cell clones without significant overlap between individual serial passages indicating that clonal succession of TIC activity in PDAC is not dependent on the cellular context in tumors in vivo. The remarkable genetic stability of xenografts during serial transplantation strongly indicates that changes in the functional state of PDAC cells and not genetic instability drive clonal succession of TIC activity in PDAC. Citation Format: Karl Roland Ehrenberg, Claudia R. Ball, Felix Oppel, Naveed Ishaque, Taronish D. Dubash, Sebastian M. Dieter, Christopher M. Hoffmann, Ulrich Abel, Moritz Koch, Jens Werner, Frank Bergmann, Manfred Schmidt, Christof von Kalle, Wilko Weichert, Jürgen Weitz, Benedikt Brors, Hanno Glimm. Clonal succession in pancreatic cancer progression is not driven by genetic instability. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1417. doi:10.1158/1538-7445.AM2015-1417

Novel mutations of the ATP7B gene in Han Chinese families with pre-symptomatic Wilson's disease.

Wilson's disease (WD) is an autosomal recessive genetic disorder of copper metabolism, caused by mutations in the ATP7B gene, resulting in copper accumulation in the liver, brain, kidney, and cornea and leading to significant disability or death if untreated. Early diagnosis and proper therapy usually predict a good prognosis, especially in pre-symptomatic WD. Genetic testing is the most accurate and effective diagnostic method for early diagnosis. The clinical and biochemical features of three unrelated Han Chinese families with pre-symptomatic WD were reported. The molecular defects in these families were investigated by polymerase chain reaction and DNA sequencing. Hundred healthy children with the same ethnic background served as controls. Bioinformatic tools (polymorphism phenotyping-2, sorting intolerant from tolerant, protein analysis through evolutionary relationships, and predictor of human deleterious single nucleotide polymorphisms) were combined and used to predict the functional effects of mutations. We identified 2 novel ATP7B mutations (p.Leu692Pro and p.Asn728Ser) and 3 known mutations (p.Met769fs, p.Arg778Leu and p.Val1216Met) in these Chinese WD families. These mutations were not observed in the 100 normal controls. The bioinformatic method showed that p.Leu692Pro and p.Asn728Ser mutations are pathogenic. Our research enriches the mutation spectrum of the ATP7B gene worldwide and provides valuable information for studying the mutation types and mode of inheritance of ATP7B in the Chinese population. Liver function analysis and genetic testing in young children with WD are necessary to shorten the time to the initiation of therapy, reduce damage to the liver and brain, and improve prognosis.

Selective Phosphorylation of South and North-Cytidine and Adenosine Methanocarba-Nucleosides by Human Nucleoside and Nucleotide Kinases Correlates with Their Growth Inhibitory Effects on Cultured Cells

Here bicyclo[3.1.0]hexane locked deoxycytidine (S-MCdC, N-MCdC), and deoxyadenosine analogs (S-MCdA and N-MCdA) were examined as substrates for purified preparations of human deoxynucleoside kinases: dCK, dGK, TK2, TK1, the ribonucleoside kinase UCK2, two NMP kinases (CMPK1, TMPK) and a NDP kinase.dCK can be important for the first step of phosphorylation of S-MCdC in cells, but S-MCdCMP was not a substrate for CMPK1, TMPK, or NDPK.dCK and dGK had a preference for the S-MCdA whereas N-MCdA was not a substrate for dCK, TK1, UCK2, TK2, dGK nucleoside kinases. The cell growth experiments suggested that N-MCdC and S-MCdA could be activated in cells by cellular kinases so that a triphosphate metabolite was formed.List of abbreviations: ddC, 2′, 3′-didioxycytosine, Zalcitabine; 3TC, β-L-(-)-2′,3′-dideoxy-3′-thiacytidine, Lamivudine; CdA, 2-cloro-2′-deoxyadenosine, Cladribine; AraA, 9-β-D-arabinofuranosyladenine; hCNT 1–3, human Concentrative Nucleoside Transporter type 1, 2 and 3; hENT 1–4, human Equilibrative Nucleoside Transporter type 1, 2, 3, and 4.

VarI-SIG 2014--From SNPs to variants: interpreting different types of genetic variants.

The decreasing cost of high-throughput sequencing is still rapidly increasing the number of known genetic [1,2]. For example, the size of the dbSNP database [3] has nearly doubled over the past years to 110 million human single nucleotide polymorphisms and short genomic variants. Unfortunately the characterization, annotation, and interpretation of these are still lagging. In particular, the interpretation of genetic and their implication in disease is one of the major challenges in personalized medicine [4-6]. The 4th edition of the Variant Interpretation Special Interest Group (VarI-SIG, formerly SNP-SIG) meeting [7-9] was held on July 12 at the ISMB 2014 in Boston (MA). The central meeting themes were Annotation and prediction of structural/functional impacts of coding variants and Genetic as effectors of change: disease and evolution. The VarI-SIG is organized as a venue for the development of a research network of scientists, necessary for facilitating the exchange of ideas and establishing new collaborations. This year's meeting attracted over 100 participants, with seven research talks and five presentations from the leading scientists of the field.

Blind study evaluation illustrates utility of the Ion PGM™ system for use in human identity DNA typing.

AimTo perform a blind study to assess the capability of the Ion Personal Genome Machine™ (PGM) system to sequence forensically relevant genetic marker panels and to characterize unknown individuals for ancestry and possible relatedness.MethodsTwelve genomic samples were provided by a third party for blinded genetic analysis. For these 12 samples, the mitochondrial genome and three PGM™ panels containing human identity single nucleotide polymorphisms (SNPs), ancestry informative SNPs, and short tandem repeats (STRs) were sequenced on the PGM™ system and analyzed.ResultsAll four genetic systems were run and analyzed on the PGM™ system in a reasonably quick time frame. Completeness of genetic profiles, depth of coverage, strand balance, and allele balance were informative metrics that illustrated the quality and reliability of the data produced. SNP genotypes allowed for identification of sex, paternal lineage, and population ancestry. STR genotypes were shown to be in complete concordance with genotypes generated by standard capillary electrophoresis-based technologies. Variants in the mitochondrial genome data provided information on population background and maternal relationships.ConclusionAll results from analysis of the 12 genomic samples were consistent with sample information provided by the sample providers at the end of the blinded study. The relatively easy identification of intra-STR allele SNPs offered the potential for increased discrimination power. The promising nature of these results warrants full validation studies of this massively parallel sequencing technology and its further development for forensic data analysis.

Genome-wide analysis of human global and transcription-coupled excision repair of UV damage at single-nucleotide resolution.

We developed a method for genome-wide mapping of DNA excision repair named XR-seq (excision repair sequencing). Human nucleotide excision repair generates two incisions surrounding the site of damage, creating an ∼30-mer. In XR-seq, this fragment is isolated and subjected to high-throughput sequencing. We used XR-seq to produce stranded, nucleotide-resolution maps of repair of two UV-induced DNA damages in human cells: cyclobutane pyrimidine dimers (CPDs) and (6-4) pyrimidine-pyrimidone photoproducts [(6-4)PPs]. In wild-type cells, CPD repair was highly associated with transcription, specifically with the template strand. Experiments in cells defective in either transcription-coupled excision repair or general excision repair isolated the contribution of each pathway to the overall repair pattern and showed that transcription-coupled repair of both photoproducts occurs exclusively on the template strand. XR-seq maps capture transcription-coupled repair at sites of divergent gene promoters and bidirectional enhancer RNA (eRNA) production at enhancers. XR-seq data also uncovered the repair characteristics and novel sequence preferences of CPDs and (6-4)PPs. XR-seq and the resulting repair maps will facilitate studies of the effects of genomic location, chromatin context, transcription, and replication on DNA repair in human cells.

Insight into neutral and disease-associated human genetic variants through interpretable predictors.

A variety of methods that predict human nonsynonymous single nucleotide polymorphisms (SNPs) to be neutral or disease-associated have been developed over the last decade. These methods are used for pinpointing disease-associated variants in the many variants obtained with next-generation sequencing technologies. The high performances of current sequence-based predictors indicate that sequence data contains valuable information about a variant being neutral or disease-associated. However, most predictors do not readily disclose this information, and so it remains unclear what sequence properties are most important. Here, we show how we can obtain insight into sequence characteristics of variants and their surroundings by interpreting predictors. We used an extensive range of features derived from the variant itself, its surrounding sequence, sequence conservation, and sequence annotation, and employed linear support vector machine classifiers to enable extracting feature importance from trained predictors. Our approach is useful for providing additional information about what features are most important for the predictions made. Furthermore, for large sets of known variants, it can provide insight into the mechanisms responsible for variants being disease-associated.