Centre d'etude du polymorphisme humain (CEPH): collaborative genetic mapping of the human genome

Experimental Data on PIRCHE and T-Cell Reactivity: HLA-DPB1-Derived Peptides Identified by PIRCHE-I Show Binding to HLA-A*02:01 in vitro and T-Cell Activation in vivo

Introduction

Human leukocyte antigen (HLA)-DPB1 mismatches during hematopoietic stem cell transplantation (HSCT) with an unrelated donor result in an increased risk for the development of graft-versus-host disease (GvHD). The number of CD8+ T-cell epitopes available for indirect allorecognition as predicted by the PIRCHE algorithm has been shown to be associated with GvHD development. As a proof of principle, PIRCHE-I predictions for HLA-DPB1 mismatches were validated in vitro and in vivo.

Methods

PIRCHE-I analysis was performed to identify HLA-DPB1-derived peptides that could theoretically bind to HLA-A*02:01. PIRCHE-I predictions for HLA-DPB1 mismatches were validated in vitro by investigating binding affinities of HLA-DPB1-derived peptides to the HLA-A*02:01 in a competition-based binding assay. To investigate the capacity of HLA-DPB1-derived peptides to elicit a T-cell response in vivo, mice were immunized with these peptides. T-cell alloreactivity was subsequently evaluated using an interferon-gamma ELISpot assay.

Results

The PIRCHE-I algorithm identified five HLA-DPB1-derived peptides (RMCRHNYEL, YIYNREEFV, YIYNREELV, YIYNREEYA, and YIYNRQEYA) to be presented by HLA-A*02:01. Binding of these peptides to HLA-A*02:01 was confirmed in a competition-based peptide binding assay, all showing an IC50 value of 21 μm or lower. The peptides elicited an interferon-gamma response in vivo.

Conclusion

Our results indicate that the PIRCHE-I algorithm can identify potential immunogenic HLA-DPB1-derived peptides present in recipients of an HLA-DPB1-mismatched donor. These combined in vitro and in vivo observations strengthen the validity of the PIRCHE-I algorithm to identify HLA-DPB1 mismatch-related GvHD development upon HSCT.

Shared genomic segments analysis identifies MHC class I and class III molecules as genetic risk factors for juvenile idiopathic arthritis

Juvenile idiopathic arthritis (JIA) is a complex rheumatic disease encompassing several clinically defined subtypes of varying severity. The etiology of JIA remains largely unknown, but genome-wide association studies (GWASs) have identified up to 22 genes associated with JIA susceptibility, including a well-established association with HLA-DRB1. Continued investigation of heritable risk factors has been hindered by disease heterogeneity and low disease prevalence. In this study, we utilized shared genomic segments (SGS) analysis on whole-genome sequencing of 40 cases from 12 multi-generational pedigrees significantly enriched for JIA. Subsets of cases are connected by a common ancestor in large extended pedigrees, increasing the power to identify disease-associated loci. SGS analysis identifies genomic segments shared among disease cases that are likely identical by descent and anchored by a disease locus. This approach revealed statistically significant signals for major histocompatibility complex (MHC) class I and class III alleles, particularly HLA-A∗02:01, which was observed at a high frequency among cases. Furthermore, we identified an additional risk locus at 12q23.2-23.3, containing genes primarily expressed by naive B cells, natural killer cells, and monocytes. The recognition of additional risk beyond HLA-DRB1 provides a new perspective on immune cell dynamics in JIA. These findings contribute to our understanding of JIA and may guide future research and therapeutic strategies.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Genome-Related Mechanisms Contributing to Differences in Alzheimer's Disease Incidence Between White and Black Older US Adults

In this manuscript, we leverage a modified GWAS algorithm adapted for use with multidimensional Cox models and data from the Health and Retirement Study to explore how genetic variation influences the size of the disparity in Alzheimer's disease (AD) incidence between older Black and White US adults. We identified four loci that were associated with higher AD incidence levels in older Black adults: (1) rs11077034 (hazard ratio (HR), 4.98) from the RBFOX1 gene; (2) rs7144494 (HR, 1.68) from the HISLA gene; (3) rs7660552 (HR, 3.07) from the SLC25A4 gene; and (4) rs12599485 (HR, 3.181) from the NIP30 gene. The RBFOX1, HISLA, SLC25A4, and NIP30 genes are known to be associated with AD (RBFOX1, NIP30) directly, and also influence the risk of AD risk-related morbidities such as hypertension (RBFOX1, SLC25A4), depression (SLC25A4), and certain cancers (HISLA, SLC25A4). A likely disparity-generating mechanism may lie in endocytosis and abnormal tissue growing mechanisms, depending on inherited gene mutations and the structure of proxies as well as gene-environment interactions with other risk factors such as lifestyle, education level, and access to adequate medical services.

Characterizing genetic variation in the regulation of the ER stress response through computational and cis-eQTL analyses

Misfolded proteins in the endoplasmic reticulum (ER) elicit the ER stress response, a large transcriptional response driven by 3 well-characterized transcription factors (TFs). This transcriptional response is variable across different genetic backgrounds. One mechanism in which genetic variation can lead to transcriptional variability in the ER stress response is through altered binding and activity of the 3 main TFs: XBP1, ATF6, and ATF4. This work attempts to better understand this mechanism by first creating a computational pipeline to identify potential binding sites throughout the human genome. We utilized GTEx data sets to identify cis-eQTLs that fall within predicted TF binding sites (TFBSs). We also utilized the ClinVar database to compare the number of pathogenic vs benign variants at different positions of the binding motifs. Finally, we performed a cis-eQTL analysis on human cell lines experiencing ER stress to identify cis-eQTLs that regulate the variable ER stress response. The majority of these cis-eQTLs are unique to a given condition: control or ER stress. Some of these stress-specific cis-eQTLs fall within putative binding sites of the 3 main ER stress response TFs, providing a potential mechanism by which these cis-eQTLs might be impacting gene expression under ER stress conditions through altered TF binding. This study represents the first cis-eQTL analysis on human samples experiencing ER stress and is a vital step toward identifying the genetic components responsible for the variable ER stress response.

Identification of common genetic polymorphisms associated with down-regulated gonadotropin levels in an exome-wide association study

OBJECTIVE

To investigate whether common genetic polymorphisms are associated with gonadotropin levels after down-regulation with daily gonadotropin-releasing hormone agonist and whether the polymorphisms of candidate variants influence the ovarian response to exogenous gonadotropins.

DESIGN

Genetic association study.

SETTING

University-affiliated in vitro fertilization center.

PATIENTS

Subjects enrolled in an exploratory exome-wide association study (n = 862), a replication exome-wide association study (n = 86), and a classifier validation study (n = 148) were recruited from September 2016 to October 2018, September 2019 to September 2020, and January 2021 to December 2021, respectively. The included patients were aged ≤40 years and had a basal follicle-stimulating hormone (FSH) ≤12 IU/L.

INTERVENTIONS

All participants received a luteal phase down-regulation long protocol. Genome DNA was extracted from the peripheral blood leukocytes. For the exploratory and replication cohorts, exome sequencing was conducted on a HiSeq 2500 sequencing platform. The multiplex polymerase chain reaction amplification technique and next-generation sequencing also were performed in the exploratory and replication cohorts. For the samples of the validation cohort, Sanger sequencing was performed.

MAIN OUTCOME MEASURES

The primary endpoint was the gonadotropin levels after down-regulation, and the secondary endpoints were hormone levels and follicle diameters during stimulation, the total dose of FSH, duration of FSH stimulation, number of oocytes retrieved, and clinical pregnancy rate.

RESULTS

In the exploratory cohort, we identified that FSHB rs6169 (P=2.71 × 10-24) and its single-nucleotide polymorphisms in high linkage disequilibrium were associated with the down-regulated FSH level. The same locus was confirmed in the replication cohort. Women carrying the C allele of FSHB rs6169 exhibited higher average estradiol level during stimulation (P=6.82 × 10-5), shorter duration of stimulation, and less amount of exogenous FSH (Pduration=0.0002; Pdose=0.0024). In the independent validation set, adding rs6169 genotypes into the prediction model for FSH level after down-regulation enhanced the area under the curve from 0.560 to 0.712 in a logistic regression model, and increased prediction accuracy by 41.05% when a support vector machine classifier was applied.

CONCLUSION

The C allele of FSHB rs6169 is a susceptibility site for the relatively high level of FSH after down-regulation, which may be associated with increased ovarian FSH sensitivity.

Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Genetic analysis of a large Han Chinese family line with schizoaffective psychosis

To locate the specific susceptibility genes of a high incidence of schizoaffective disease (SAD) with autonomic dominant inheritance, we recruited a family group from Henan Province with a high incidence of SAD, including 19 individuals sampled from five generations. We used a genome-wide high-density SNP chip to perform genotype detection. The LINKAGE package and MENDEL programs were used for. The two-point and multipoint analyses were calculated by Merlin and SimWalk2 software to obtain the nonparametric linkage (NPL) value, corresponding P value, and parameter linkage limit of detection (LOD) value. Genome-wide linkage analysis yielded a significant linkage signal located on the short arm of chromosome 19. In the dominant genetic model, the LOD of the multipoint parametric analysis was 2.5, and the nonparametric analysis was 19.4 (P < 0.00001). Further haploid genotype analysis localized the candidate region in the 19p13.3-13.2 region, beginning at rs178414 and ending at rs11668751 with a physical length of approximately 4.9 Mb. We believe that the genes responsible for SAD are in this region.

Annotation of structural variants with reported allele frequencies and related metrics from multiple datasets using SVAFotate

BACKGROUND

Identification of deleterious genetic variants using DNA sequencing data relies on increasingly detailed filtering strategies to isolate the small subset of variants that are more likely to underlie a disease phenotype. Datasets reflecting population allele frequencies of different types of variants serve as powerful filtering tools, especially in the context of rare disease analysis. While such population-scale allele frequency datasets now exist for structural variants (SVs), it remains a challenge to match SV calls between multiple datasets, thereby complicating estimates of a putative SV's population allele frequency.

RESULTS

We introduce SVAFotate, a software tool that enables the annotation of SVs with variant allele frequency and related information from existing SV datasets. As a result, VCF files annotated by SVAFotate offer a variety of metrics to aid in the stratification of SVs as common or rare in the broader human population.

CONCLUSIONS

Here we demonstrate the use of SVAFotate in the classification of SVs with regards to their population frequency and illustrate how SVAFotate's annotations can be used to filter and prioritize SVs. Lastly, we detail how best to utilize these SV annotations in the analysis of genetic variation in studies of rare disease.

A machine learning approach for missing persons cases with high genotyping errors

Estimating the relationships between individuals is one of the fundamental challenges in many fields. In particular, relationship.ip estimation could provide valuable information for missing persons cases. The recently developed investigative genetic genealogy approach uses high-density single nucleotide polymorphisms (SNPs) to determine close and more distant relationships, in which hundreds of thousands to tens of millions of SNPs are generated either by microarray genotyping or whole-genome sequencing. The current studies usually assume the SNP profiles were generated with minimum errors. However, in the missing person cases, the DNA samples can be highly degraded, and the SNP profiles generated from these samples usually contain lots of errors. In this study, a machine learning approach was developed for estimating the relationships with high error SNP profiles. In this approach, a hierarchical classification strategy was employed first to classify the relationships by degree and then the relationship types within each degree separately. As for each classification, feature selection was implemented to gain better performance. Both simulated and real data sets with various genotyping error rates were utilized in evaluating this approach, and the accuracies of this approach were higher than individual measures; namely, this approach was more accurate and robust than the individual measures for SNP profiles with genotyping errors. In addition, the highest accuracy could be obtained by providing the same genotyping error rates in train and test sets, and thus estimating genotyping errors of the SNP profiles is critical to obtaining high accuracy of relationship estimation.

A high-throughput real-time PCR tissue-of-origin test to distinguish blood from lymphoblastoid cell line DNA for (epi)genomic studies

Lymphoblastoid cell lines (LCLs) derive from blood infected in vitro by Epstein–Barr virus and were used in several genetic, transcriptomic and epigenomic studies. Although few changes were shown between LCL and blood genotypes (SNPs) validating their use in genetics, more were highlighted for other genomic features and/or in their transcriptome and epigenome. This could render them less appropriate for these studies, notably when blood DNA could still be available. Here we developed a simple, high-throughput and cost-effective real-time PCR approach allowing to distinguish blood from LCL DNA samples based on the presence of EBV relative load and rearranged T-cell receptors γ and β. Our approach was able to achieve 98.5% sensitivity and 100% specificity on DNA of known origin (458 blood and 316 LCL DNA). It was further applied to 1957 DNA samples from the CEPH Aging cohort comprising DNA of uncertain origin, identifying 784 blood and 1016 LCL DNA. A subset of these DNA was further analyzed with an epigenetic clock indicating that DNA extracted from blood should be preferred to LCL for DNA methylation-based age prediction analysis. Our approach could thereby be a powerful tool to ascertain the origin of DNA in old collections prior to (epi)genomic studies.

Microsatellite Markers in Biobanking: A New Multiplexed Assay

Microsatellites, or MSATs, offer a fast and cost-effective way for biobanks to establish a biospecimen genetic profile. Importantly, this genetic profile can be used to authenticate multiple submissions derived from the same individual as well as biospecimens derived from the same original sample submission over time. While the Certificate of Confidentiality provided by the National Institutes of Health offers some meaningful protection to prevent the disclosure of potentially identifiable information to entities within the United States, we consider, in this study, the potential to offer additional protection to participants who choose to donate to biobanks by minimizing the use of forensic Combined DNA Index System (CODIS) MSAT markers in biobanking. To this end, we report the design and validation of a new multiplexed MSAT assay that does not include CODIS markers for use in biobanking operations and quality control management.

Differences in the number of de novo mutations between individuals are due to small family-specific effects and stochasticity

The number of de novo mutations (DNMs) in the human germline is correlated with parental age at conception, but this explains only part of the observed variation. We investigated whether there is a family-specific contribution to the number of DNMs in offspring. The analysis of DNMs in 111 dizygotic twin pairs did not identify a substantial family-specific contribution. This result was corroborated by comparing DNMs of 1669 siblings to those of age-matched unrelated offspring following correction for parental age. In addition, by modeling DNM data from 1714 multi-offspring families, we estimated that the family-specific contribution explains ∼5.2% of the variation in DNM number. Furthermore, we found no substantial difference between the observed number of DNMs and those predicted by a stochastic Poisson process. We conclude that there is a small family-specific contribution to DNM number and that stochasticity explains a large proportion of variation in DNM counts.

History of the methodology of disease gene identification

The past 45 years have witnessed a triumph in the discovery of genes and genetic variation that cause Mendelian disorders due to high impact variants. Important discoveries and organized projects have provided the necessary tools and infrastructure for the identification of gene defects leading to thousands of monogenic phenotypes. This endeavor can be divided in three phases in which different laboratory strategies were employed for the discovery of disease-related genes: (i) the biochemical phase, (ii) the genetic linkage followed by positional cloning phase, and (iii) the sequence identification phase. However, much more work is needed to identify all the high impact genomic variation that substantially contributes to the phenotypic variation.

Unfazed: parent-of-origin detection for large and small de novo variants

SUMMARY

Unfazed is a command-line tool to determine the parental gamete of origin for de novo mutations from paired-end Illumina DNA sequencing reads. Unfazed uses variant information for a sequenced trio to identify the parental gamete of origin by linking phase-informative inherited variants to de novo mutations using read-based phasing. It achieves a high success rate by chaining reads into haplotype groups, thus increasing the search space for informative sites. Unfazed provides a simple command-line interface and scales well to large inputs, determining parent-of-origin for nearly 30,000 de novo variants in under 60 hours.

AVAILABILITY

Unfazed is available at https://github.com/jbelyeu/unfazed.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

Response and Toxicity to Cytarabine Therapy in Leukemia and Lymphoma: From Dose Puzzle to Pharmacogenomic Biomarkers

Simple Summary

In this review, the authors propose a crosswise examination of cytarabine-related issues ranging from the spectrum of clinical activity and severe toxicities, through updated cellular pharmacology and drug formulations, to the genetic variants associated with drug-induced phenotypes. Cytarabine (cytosine arabinoside; Ara-C) in multiagent chemotherapy regimens is often used for leukemia or lymphoma treatments, as well as neoplastic meningitis. Chemotherapy regimens can induce a suboptimal clinical outcome in a fraction of patients. The individual variability in clinical response to Leukemia & Lymphoma treatments among patients appears to be associated with intracellular accumulation of Ara-CTP due to genetic variants related to metabolic enzymes. The review provides exhaustive information on the effects of Ara-C-based therapies, the adverse drug reaction will also be provided including bone pain, ocular toxicity (corneal pain, keratoconjunctivitis, and blurred vision), maculopapular rash, and occasional chest pain. Evidence for predicting the response to cytarabine-based treatments will be highlighted, pointing at their significant impact on the routine management of blood cancers.

Abstract

Cytarabine is a pyrimidine nucleoside analog, commonly used in multiagent chemotherapy regimens for the treatment of leukemia and lymphoma, as well as for neoplastic meningitis. Ara-C-based chemotherapy regimens can induce a suboptimal clinical outcome in a fraction of patients. Several studies suggest that the individual variability in clinical response to Leukemia & Lymphoma treatments among patients, underlying either Ara-C mechanism resistance or toxicity, appears to be associated with the intracellular accumulation and retention of Ara-CTP due to genetic variants related to metabolic enzymes. Herein, we reported (a) the latest Pharmacogenomics biomarkers associated with the response to cytarabine and (b) the new drug formulations with optimized pharmacokinetics. The purpose of this review is to provide readers with detailed and comprehensive information on the effects of Ara-C-based therapies, from biological to clinical practice, maintaining high the interest of both researcher and clinical hematologist. This review could help clinicians in predicting the response to cytarabine-based treatments.

Neurexin 1 variants as risk factors for suicide death

Suicide is a significant public health concern with complex etiology. Although the genetic component of suicide is well established, the scope of gene networks and biological mechanisms underlying suicide has yet to be defined. Previously, we reported genome-wide evidence that neurexin 1 (NRXN1), a key synapse organizing molecule, is associated with familial suicide risk. Here we present new evidence for two non-synonymous variants (rs78540316; P469S and rs199784139; H885Y) associated with increased familial risk of suicide death. We tested the impact of these variants on binding interactions with known partners and assessed functionality in a hemi-synapse formation assay. Although the formation of hemi-synapses was not altered with the P469S variant relative to wild-type, both variants increased binding to the postsynaptic binding partner, leucine-rich repeat transmembrane neuronal 2 (LRRTM2) in vitro. Our findings indicate that variants in NRXN1 and related synaptic genes warrant further study as risk factors for suicide death.

ped_draw: pedigree drawing with ease

Background

Pedigree files are ubiquitously used within bioinformatics and genetics studies to convey critical information about relatedness, sex and affected status of study samples. While the text based format of ped files is efficient for computational methods, it is not immediately intuitive to a bioinformatician or geneticist trying to understand family structures, many of which encode the affected status of individuals across multiple generations. The visualization of pedigrees into connected nodes with descriptive shapes and shading provides a far more interpretable format to recognize visual patterns and intuit family structures. Despite these advantages of a visual pedigree, it remains difficult to quickly and accurately visualize a pedigree given a pedigree text file.

Results

Here we describe ped_draw a command line and web tool as a simple and easy solution to pedigree visualization. Ped_draw is capable of drawing complex multi-generational pedigrees and conforms to the accepted standards for depicting pedigrees visually. The command line tool can be used as a simple one liner command, utilizing graphviz to generate an image file. The web tool, https://peddraw.github.io, allows the user to either: paste a pedigree file, type to construct a pedigree file in the text box or upload a pedigree file. Users can save the generated image file in various formats.

Conclusions

We believe ped_draw is a useful pedigree drawing tool that improves on current methods due to its ease of use and approachability. Ped_draw allows users with various levels of expertise to quickly and easily visualize pedigrees.

Identification of a tumor-specific allo-HLA-restricted γδTCR

γδT cells are key players in cancer immune surveillance because of their ability to recognize malignant transformed cells, which makes them promising therapeutic tools in the treatment of cancer. However, the biological mechanisms of how γδT-cell receptors (TCRs) interact with their ligands are poorly understood. Within this context, we describe the novel allo-HLA-restricted and CD8α-dependent Vγ5Vδ1TCR. In contrast to the previous assumption of the general allo-HLA reactivity of a minor fraction of γδTCRs, we show that classic anti-HLA-directed, γδTCR-mediated reactivity can selectively act on hematological and solid tumor cells, while not harming healthy tissues in vitro and in vivo. We identified the molecular interface with proximity to the peptide-binding groove of HLA-A*24:02 as the essential determinant for recognition and describe the critical role of CD8 as a coreceptor. We conclude that alloreactive γδT-cell repertoires provide therapeutic opportunities, either within the context of haplotransplantation or as individual γδTCRs for genetic engineering of tumor-reactive T cells.

Germline mutation rates in young adults predict longevity and reproductive lifespan

Ageing may be due to mutation accumulation across the lifespan, leading to tissue dysfunction, disease, and death. We tested whether germline autosomal mutation rates in young adults predict their remaining survival, and, for women, their reproductive lifespans. Age-adjusted mutation rates (AAMRs) in 61 women and 61 men from the Utah CEPH (Centre d’Etude du Polymorphisme Humain) families were determined. Age at death, cause of death, all-site cancer incidence, and reproductive histories were provided by the Utah Population Database, Utah Cancer Registry, and Utah Genetic Reference Project. Higher AAMRs were significantly associated with higher all-cause mortality in both sexes combined. Subjects in the top quartile of AAMRs experienced more than twice the mortality of bottom quartile subjects (hazard ratio [HR], 2.07; 95% confidence interval [CI], 1.21–3.56; p = 0.008; median survival difference = 4.7 years). Fertility analyses were restricted to women whose age at last birth (ALB) was ≥ 30 years, the age when fertility begins to decline. Women with higher AAMRs had significantly fewer live births and a younger ALB. Adult germline mutation accumulation rates are established in adolescence, and later menarche in women is associated with delayed mutation accumulation. We conclude that germline mutation rates in healthy young adults may provide a measure of both reproductive and systemic ageing. Puberty may induce the establishment of adult mutation accumulation rates, just when DNA repair systems begin their lifelong decline.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Pedigree-based estimation of human mobile element retrotransposition rates

Germline mutation rates in humans have been estimated for a variety of mutation types, including single-nucleotide and large structural variants. Here, we directly measure the germline retrotransposition rate for the three active retrotransposon elements: L1, Alu, and SVA. We used three tools for calling mobile element insertions (MEIs) (MELT, RUFUS, and TranSurVeyor) on blood-derived whole-genome sequence (WGS) data from 599 CEPH individuals, comprising 33 three-generation pedigrees. We identified 26 de novo MEIs in 437 births. The retrotransposition rate estimates for Alu elements, one in 40 births, is roughly half the rate estimated using phylogenetic analyses, a difference in magnitude similar to that observed for single-nucleotide variants. The L1 retrotransposition rate is one in 63 births and is within range of previous estimates (1:20-1:200 births). The SVA retrotransposition rate, one in 63 births, is much higher than the previous estimate of one in 900 births. Our large, three-generation pedigrees allowed us to assess parent-of-origin effects and the timing of insertion events in either gametogenesis or early embryonic development. We find a statistically significant paternal bias in Alu retrotransposition. Our study represents the first in-depth analysis of the rate and dynamics of human retrotransposition from WGS data in three-generation human pedigrees.

Large, three-generation human families reveal post-zygotic mosaicism and variability in germline mutation accumulation

The number of de novo mutations (DNMs) found in an offspring's genome increases with both paternal and maternal age. But does the rate of mutation accumulation in human gametes differ across families? Using sequencing data from 33 large, three-generation CEPH families, we observed significant variability in parental age effects on DNM counts across families, ranging from 0.19 to 3.24 DNMs per year. Additionally, we found that ~3% of DNMs originated following primordial germ cell specification in a parent, and differed from non-mosaic germline DNMs in their mutational spectra. We also discovered that nearly 10% of candidate DNMs in the second generation were post-zygotic, and present in both somatic and germ cells; these gonosomal mutations occurred at equivalent frequencies on both parental haplotypes. Our results demonstrate that rates of germline mutation accumulation vary among families with similar ancestry, and confirm that post-zygotic mosaicism is a substantial source of human DNM.

High-efficiency genomic editing in Epstein-Barr virus-transformed lymphoblastoid B cells using a single-stranded donor oligonucleotide strategy

While human lymphoblastoid cell lines represent a valuable resource for population genetic studies, they have usually been regarded as difficult for CRISPR-mediated genomic editing because of very inefficient DNA transfection and retroviral or lentiviral transduction in these cells, which becomes a substantial problem when multiple constructs need to be co-expressed. Here we describe a protocol using a single-stranded donor oligonucleotide strategy for 'scarless' editing in lymphoblastoid cells, yielding 12/60 (20%) of clones with homology-directed recombination, when rates of <5-10% are frequently typical for many other cell types. The protocol does not require the use of lentiviruses or stable transfection, permitting lymphoblastoid cell lines to be used for CRISPR-mediated genomic targeting and screening in population genetic studies.

Is it time to change the reference genome?

The use of the human reference genome has shaped methods and data across modern genomics. This has offered many benefits while creating a few constraints. In the following opinion, we outline the history, properties, and pitfalls of the current human reference genome. In a few illustrative analyses, we focus on its use for variant-calling, highlighting its nearness to a 'type specimen'. We suggest that switching to a consensus reference would offer important advantages over the continued use of the current reference with few disadvantages.

Functional genetic variants can mediate their regulatory effects through alteration of transcription factor binding

Functional variants in the genome are usually identified by their association with local gene expression, DNA methylation or chromatin states. DNA sequence motif analysis and chromatin immunoprecipitation studies have provided indirect support for the hypothesis that functional variants alter transcription factor binding to exert their effects. In this study, we provide direct evidence that functional variants can alter transcription factor binding. We identify a multifunctional variant within the TBC1D4 gene encoding a canonical NFκB binding site, and edited it using CRISPR-Cas9 to remove this site. We show that this editing reduces TBC1D4 expression, local chromatin accessibility and binding of the p65 component of NFκB. We then used CRISPR without genomic editing to guide p65 back to the edited locus, demonstrating that this re-targeting, occurring ~182 kb from the gene promoter, is enough to restore the function of the locus, supporting the central role of transcription factors mediating the effects of functional variants.

Functional Analysis of Promoter Variants in Genes Involved in Sex Steroid Action, DNA Repair and Cell Cycle Control

Genetic variants affecting the regulation of gene expression are among the main causes of human diversity. The potential importance of regulatory polymorphisms is underscored by results from Genome Wide Association Studies, which have already implicated such polymorphisms in the susceptibility to complex diseases such as breast cancer. In this study, we re-sequenced the promoter regions of 24 genes involved in pathways related to breast cancer including sex steroid action, DNA repair, and cell cycle control in 60 unrelated Caucasian individuals. We constructed haplotypes and assessed the functional impact of promoter variants using gene reporter assays and electrophoretic mobility shift assays. We identified putative functional variants within the promoter regions of estrogen receptor 1 (ESR1), ESR2, forkhead box A1 (FOXA1), ubiquitin interaction motif containing 1 (UIMC1) and cell division cycle 7 (CDC7). The functional polymorphism on CDC7, rs13447455, influences CDC7 transcriptional activity in an allele-specific manner and alters DNA–protein complex formation in breast cancer cell lines. Moreover, results from the Breast Cancer Association Consortium show a marginal association between rs13447455 and breast cancer risk (p = 9.3 × 10⁻⁵), thus warranting further investigation. Furthermore, our study has helped provide methodological solutions to some technical difficulties that were encountered with gene reporter assays, particularly regarding inter-clone variability and statistical consistency.

The influence of Neanderthal alleles on cytotoxic response

Various studies have shown that people of Eurasian origin contain traces of DNA inherited from interbreeding with Neanderthals. Recent studies have demonstrated that these Neanderthal variants influence a range of clinically important traits and diseases. Thus, understanding the genetic factors responsible for the variability in individual response to drug or chemical exposure is a key goal of pharmacogenomics and toxicogenomics, as dose responses are clinically and epidemiologically important traits. It is well established that ethnic and racial differences are important in dose response traits, but to our knowledge the influence of Neanderthal ancestry on response to xenobiotics is unknown. Towards this aim, we examined if Neanderthal ancestry plays a role in cytotoxic response to anti-cancer drugs and toxic environmental chemicals. We identified common Neanderthal variants in lymphoblastoid cell lines (LCLs) derived from the globally diverse 1000 Genomes Project and Caucasian cell lines from the Children's Hospital of Oakland Research Institute. We analyzed the effects of these Neanderthal alleles on cytotoxic response to 29 anti-cancer drugs and 179 environmental chemicals at varying concentrations using genome-wide data. We identified and replicated single nucleotide polymorphisms (SNPs) from these association results, including a SNP in the SNORD-113 cluster. Our results also show that the Neanderthal alleles cumulatively lead to increased sensitivity to both the anti-cancer drugs and the environmental chemicals. Our results demonstrate the influence of Neanderthal ancestry-informative markers on cytotoxic response. These results could be important in identifying biomarkers for personalized medicine or in dissecting the underlying etiology of dose response traits.

The contribution of DNA microarray technology to gene expression profiling in Leishmania spp.: A retrospective view

The first completed genome project of any living organism, excluding viruses, was of the gammaproteobacteria Haemophilus influenzae in 1995. Until the last decade, genome sequencing was very tedious because genome survey sequences (GSS) and/or expressed sequence tags (ESTs) belonging to plasmid, cosmid, and artificial chromosome genome libraries had to be sequenced and assembled in silico. No genome is completely assembled because gaps and unassembled contigs are always remaining. However, most represent an organism's whole genome from a practical point of view. The first genome sequencing projects of trypanosomatid parasites Leishmania major, Trypanosoma cruzi, and T. brucei were completed in 2005 following those strategies. The functional genomics era developed on the basis of microarray technology and has been continuously evolving. In the case of the genus Leishmania, substantial information about differentiation in the digenetic life cycle of the parasite has been obtained. More recently, next generation sequencing has revolutionized genome sequencing and functional genomics, leading to more sensitive and accurate results by using much fewer resources. Though this new technology is more advantageous, it does not invalidate microarray results. In fact, promising vaccine candidates and drug targets have been found by means of microarray-based screening and preliminary proof-of-concept tests.

Genetic Ancestry and Natural Selection Drive Population Differences in Immune Responses to Pathogens

Individuals from different populations vary considerably in their susceptibility to immune-related diseases. To understand how genetic variation and natural selection contribute to these differences, we tested for the effects of African versus European ancestry on the transcriptional response of primary macrophages to live bacterial pathogens. A total of 9.3% of macrophage-expressed genes show ancestry-associated differences in the gene regulatory response to infection, and African ancestry specifically predicts a stronger inflammatory response and reduced intracellular bacterial growth. A large proportion of these differences are under genetic control: for 804 genes, more than 75% of ancestry effects on the immune response can be explained by a single cis- or trans-acting expression quantitative trait locus (eQTL). Finally, we show that genetic effects on the immune response are strongly enriched for recent, population-specific signatures of adaptation. Together, our results demonstrate how historical selective events continue to shape human phenotypic diversity today, including for traits that are key to controlling infection.

A reference data set of 5.4 million phased human variants validated by genetic inheritance from sequencing a three-generation 17-member pedigree

Improvement of variant calling in next-generation sequence data requires a comprehensive, genome-wide catalog of high-confidence variants called in a set of genomes for use as a benchmark. We generated deep, whole-genome sequence data of 17 individuals in a three-generation pedigree and called variants in each genome using a range of currently available algorithms. We used haplotype transmission information to create a phased “Platinum” variant catalog of 4.7 million single-nucleotide variants (SNVs) plus 0.7 million small (1–50 bp) insertions and deletions (indels) that are consistent with the pattern of inheritance in the parents and 11 children of this pedigree. Platinum genotypes are highly concordant with the current catalog of the National Institute of Standards and Technology for both SNVs (>99.99%) and indels (99.92%) and add a validated truth catalog that has 26% more SNVs and 45% more indels. Analysis of 334,652 SNVs that were consistent between informatics pipelines yet inconsistent with haplotype transmission (“nonplatinum”) revealed that the majority of these variants are de novo and cell-line mutations or reside within previously unidentified duplications and deletions. The reference materials from this study are a resource for objective assessment of the accuracy of variant calls throughout genomes.

Variable bleeding phenotype in an Amish pedigree with von Willebrand disease

Through a cross-sectional study design, the bleeding phenotype in the Amish in Indiana (IN) and Wisconsin (WI) was described using two different bleeding scores. von Willebrand factor (VWF) testing was performed and bleeding questionnaires from Centers for Disease Control and Prevention (CDC) and European MCMDM-1 (Tosetto bleeding score (BS)) were administered to the IN and WI cohort respectively. Seven hundred and seventy nine subjects were recruited, 17% were diagnosed with VWD based on Ristocetin cofactor, VWF:RCo < 30 IU/dl. Majority of the affected (AF), 67%, were tested and had a common mutation c.4120 C > T. The WI AF were much younger at a mean age 15 years vs 26 years in IN AF cohort. The AF subjects had a median VWF:RCo of 13IU/dl with a statistically significant higher median BS 1 versus 0 in the WI AF vs WI Unaffected (UA), 2 vs 1 in the IN AF vs IN UA, P < 0.01. Adults had a higher median BS compared to children in the WI and IN cohort, 2 vs 1 and 3 vs 1 respectively (P < 0.05) but there was no statistically significant difference in the BS between males and females in either cohort. The common symptoms reported were epistaxis and gingival oozing. BS ≥ 3 and BS ≥ 4 were observed in 46% of AF IN and 16.6% of AF WI, respectively. There was significant variability in the bleeding phenotype, with an overall low BS in the affected Amish with VWD, despite a unifying mutation. Am. J. Hematol. 91:E431-E435, 2016. © 2016 Wiley Periodicals, Inc.

Genomic variation within alpha satellite DNA influences centromere location on human chromosomes with metastable epialleles

Alpha satellite is a tandemly organized type of repetitive DNA that comprises 5% of the genome and is found at all human centromeres. A defined number of 171-bp monomers are organized into chromosome-specific higher-order repeats (HORs) that are reiterated thousands of times. At least half of all human chromosomes have two or more distinct HOR alpha satellite arrays within their centromere regions. We previously showed that the two alpha satellite arrays of Homo sapiens Chromosome 17 (HSA17), D17Z1 and D17Z1-B, behave as centromeric epialleles, that is, the centromere, defined by chromatin containing the centromeric histone variant CENPA and recruitment of other centromere proteins, can form at either D17Z1 or D17Z1-B. Some individuals in the human population are functional heterozygotes in that D17Z1 is the active centromere on one homolog and D17Z1-B is active on the other. In this study, we aimed to understand the molecular basis for how centromere location is determined on HSA17. Specifically, we focused on D17Z1 genomic variation as a driver of epiallele formation. We found that D17Z1 arrays that are predominantly composed of HOR size and sequence variants were functionally less competent. They either recruited decreased amounts of the centromere-specific histone variant CENPA and the HSA17 was mitotically unstable, or alternatively, the centromere was assembled at D17Z1-B and the HSA17 was stable. Our study demonstrates that genomic variation within highly repetitive, noncoding DNA of human centromere regions has a pronounced impact on genome stability and basic chromosomal function.

RhoB Mediates Phosphoantigen Recognition by Vγ9Vδ2 T Cell Receptor

Human Vγ9Vδ2 T cells respond to tumor cells by sensing elevated levels of phosphorylated intermediates of the dysregulated mevalonate pathway, which is translated into activating signals by the ubiquitously expressed butyrophilin A1 (BTN3A1) through yet unknown mechanisms. Here, we developed an unbiased, genome-wide screening method that identified RhoB as a critical mediator of Vγ9Vδ2 TCR activation in tumor cells. Our results show that Vγ9Vδ2 TCR activation is modulated by the GTPase activity of RhoB and its redistribution to BTN3A1. This is associated with cytoskeletal changes that directly stabilize BTN3A1 in the membrane, and the subsequent dissociation of RhoB from BTN3A1. Furthermore, phosphoantigen accumulation induces a conformational change in BTN3A1, rendering its extracellular domains recognizable by Vγ9Vδ2 TCRs. These complementary events provide further evidence for inside-out signaling as an essential step in the recognition of tumor cells by a Vγ9Vδ2 TCR.

Genetic variations in GPSM3 associated with protection from rheumatoid arthritis affect its transcript abundance

G protein signaling modulator 3 (GPSM3) is a regulator of G protein-coupled receptor signaling, with expression restricted to leukocytes and lymphoid organs. Previous genome-wide association studies have highlighted single-nucleotide polymorphisms (SNPs; rs204989 and rs204991) in a region upstream of the GPSM3 transcription start site as being inversely correlated to the prevalence of rheumatoid arthritis (RA)-this association is supported by the protection afforded to Gpsm3-deficient mice in models of inflammatory arthritis. Here, we assessed the functional consequences of these polymorphisms. We collected biospecimens from 50 volunteers with RA diagnoses, 50 RA-free volunteers matched to the aforementioned group and 100 unmatched healthy young volunteers. We genotyped these individuals for GPSM3 (rs204989, rs204991), CCL21 (rs2812378) and HLA gene region (rs6457620) polymorphisms, and found no significant differences in minor allele frequencies between the RA and disease-free cohorts. However, we identified that individuals homozygous for SNPs rs204989 and rs204991 had decreased GPSM3 transcript abundance relative to individuals homozygous for the major allele. In vitro promoter activity studies suggest that SNP rs204989 is the primary cause of this decrease in transcript levels. Knockdown of GPSM3 in THP-1 cells, a human monocytic cell line, was found to disrupt ex vivo migration to the chemokine MCP-1.

Minor histocompatibility Ags: identification strategies, clinical results and translational perspectives

Allogeneic stem cell transplantation (allo-SCT) and donor lymphocyte infusion are effective treatment modalities for various hematological malignancies. Their therapeutic effect, the graft-versus-tumor (GvT) effect, is based mainly on an alloimmune response of donor T cells directed at tumor cells, in which differences in the expression of minor histocompatibility Ags (mHags) on the cells of the patient and donor have a crucial role. However, these differences are also responsible for induction of sometimes detrimental GvHD. As relapse and development of GvHD pose major threats for a large proportion of allotransplanted patients, additional therapeutic strategies are required. To augment the GvT response without increasing the risk of GvHD, specific mHag-directed immunotherapeutic strategies have been developed. Over the past years, much effort has been put into the identification of therapeutically relevant mHags to enable these strategies for a substantial proportion of patients. Currently, the concept of mHag-directed immunotherapy is tested in clinical trials on feasibility, safety and efficacy. In this review, we will summarize the recent developments in mHag identification and the clinical data on mHag-specific immune responses and mHag-directed therapies in patients with hematological malignancies. Finally, we will outline the current challenges and future prospectives in the field.

Minor histocompatibility antigens: past, present, and future

Minor histocompatibility (H) antigens are key molecules driving allo-immune responses in both graft-versus-host-disease (GvHD) and in graft-versus-leukemia (GvL) reactivity in human leukocyte antigen (HLA)-matched hematopoietic stem-cell transplantation (HSCT). Dissection of the dual function of minor H antigens became evident through their different modes of tissue and cell expression, i.e. hematopoietic system-restricted or broad. Broadly expressed minor H antigens can cause both GvHD and GvL effects, while hematopoietic system-restricted minor H antigens are more prone to induce GvL responses. This phenomenon renders the latter group of minor H antigens as curative tools for HSCT-based immunotherapy of hematological malignancies and disorders, in which minor H antigen-specific responses are enhanced in order to eradicate the malignant cells. This article describes the immunogenetics of minor H antigens and methods that have been developed to identify them. Moreover, it summarizes the clinical relevance of minor H antigens in transplantation, with special regards to allogeneic HSCT and solid-organ transplantation.

Evaluating intra- and inter-individual variation in the human placental transcriptome

BACKGROUND

Gene expression variation is a phenotypic trait of particular interest as it represents the initial link between genotype and other phenotypes. Analyzing how such variation apportions among and within groups allows for the evaluation of how genetic and environmental factors influence such traits. It also provides opportunities to identify genes and pathways that may have been influenced by non-neutral processes. Here we use a population genetics framework and next generation sequencing to evaluate how gene expression variation is apportioned among four human groups in a natural biological tissue, the placenta.

RESULTS

We estimate that on average, 33.2%, 58.9%, and 7.8% of the placental transcriptome is explained by variation within individuals, among individuals, and among human groups, respectively. Additionally, when technical and biological traits are included in models of gene expression they each account for roughly 2% of total gene expression variation. Notably, the variation that is significantly different among groups is enriched in biological pathways associated with immune response, cell signaling, and metabolism. Many biological traits demonstrate correlated changes in expression in numerous pathways of potential interest to clinicians and evolutionary biologists. Finally, we estimate that the majority of the human placental transcriptome exhibits expression profiles consistent with neutrality; the remainder are consistent with stabilizing selection, directional selection, or diversifying selection.

CONCLUSIONS

We apportion placental gene expression variation into individual, population, and biological trait factors and identify how each influence the transcriptome. Additionally, we advance methods to associate expression profiles with different forms of selection.

Genome-wide association and pharmacological profiling of 29 anticancer agents using lymphoblastoid cell lines

AIM

Association mapping with lymphoblastoid cell lines (LCLs) is a promising approach in pharmacogenomics research, and in the current study we utilized LCLs to perform association mapping for 29 chemotherapy drugs.

MATERIALS & METHODS

Currently, we use LCLs to perform genome-wide association mapping of the cytotoxic response of 520 European-Americans to 29 different anticancer drugs; the largest LCL study to date. A novel association approach using a multivariate analysis of covariance design was employed with the software program MAGWAS, testing for differences in the dose-response profiles between genotypes without making assumptions about the response curve or the biologic mode of association. Additionally, by classifying 25 of the 29 drugs into eight families according to structural and mechanistic relationships, MAGWAS was used to test for associations that were shared across each drug family. Finally, a unique algorithm using multivariate responses and multiple linear regressions across pairs of response curves was used for unsupervised clustering of drugs.

RESULTS

Among the single-drug studies, suggestive associations were obtained for 18 loci, 12 within/near genes. Three of these, MED12L, CHN2 and MGMT, have been previously implicated in cancer pharmacogenomics. The drug family associations resulted in four additional suggestive loci (three contained within/near genes). One of these genes, HDAC4, associated with the DNA alkylating agents, shows possible clinical interactions with temozolomide. For the drug clustering analysis, 18 of 25 drugs clustered into the appropriate family.

CONCLUSION

This study demonstrates the utility of LCLs in identifying genes that have clinical importance in drug response and for assigning unclassified agents to specific drug families, and proposes new candidate genes for follow-up in a large number of chemotherapy drugs.

Analysis of genome-wide RNA-sequencing data suggests age of the CEPH/Utah (CEU) lymphoblastoid cell lines systematically biases gene expression profiles

In human, Lymphoblastoid cell lines (LCLs) from the CEPH/CEU (Centre d'Etude du Polymorphisme Humain – Utah) family resource have been extensively used for examining the genetics of gene expression levels. However, we noted that CEU/CEPH cell lines were collected and transformed approximately thirty years ago, much earlier than the other cell lines from the pertaining individuals, which we suspected could potentially affect gene expression, data analysis and results interpretation. In this study, by analyzing RNA sequencing data of CEU and the other three European populations as well as an African population, we systematically examined and evaluated the potential confounding effect of LCL age on gene expression levels and patterns. Our results indicated that gene expression profiles of CEU samples have been biased by the older age of CEU cell lines. Interestingly, most of CEU-specific expressions are associated with functions related to cell proliferation, which are more likely due to older age of cell lines than intrinsic characters of the population. We suggested the results be carefully explained when CEU LCLs are used for transcriptomic data analysis in future studies.

A genome-wide association study identifies a functional ERAP2 haplotype associated with birdshot chorioretinopathy

Birdshot chorioretinopathy (BSCR) is a rare form of autoimmune uveitis that can lead to severe visual impairment. Intriguingly, >95% of cases carry the HLA-A29 allele, which defines the strongest documented HLA association for a human disease. We have conducted a genome-wide association study in 96 Dutch and 27 Spanish cases, and 398 unrelated Dutch and 380 Spanish controls. Fine-mapping the primary MHC association through high-resolution imputation at classical HLA loci, identified HLA-A*29:02 as the principal MHC association (odds ratio (OR) = 157.5, 95% CI 91.6-272.6, P = 6.6 × 10(-74)). We also identified two novel susceptibility loci at 5q15 near ERAP2 (rs7705093; OR = 2.3, 95% CI 1.7-3.1, for the T allele, P = 8.6 × 10(-8)) and at 14q32.31 in the TECPR2 gene (rs150571175; OR = 6.1, 95% CI 3.2-11.7, for the A allele, P = 3.2 × 10(-8)). The association near ERAP2 was confirmed in an independent British case-control samples (combined meta-analysis P = 1.7 × 10(-9)). Functional analyses revealed that the risk allele of the polymorphism near ERAP2 is strongly associated with high mRNA and protein expression of ERAP2 in B cells. This study further defined an extremely strong MHC risk component in BSCR, and detected evidence for a novel disease mechanism that affects peptide processing in the endoplasmic reticulum.

Evidence for linkage and association of GABRB3 and GABRA5 to panic disorder

Panic disorder (PD) is a debilitating anxiety disorder characterized by episodes of intense fear with autonomic and psychological symptoms that lead to behavioral impairment. A convergence of genetic and biological evidence implicates gamma-aminobutyric acid type A receptor subunits on chromosome 15q12 as candidate genes for PD. This study investigated 120 Caucasian, multiplex PD pedigrees using regional microsatellites (chr15q11-13) and found support for linkage (logarithm of odds (LOD) ⩾2), with a prominent parent-of-origin effect. Genotyping with 10 single-nucleotide polymorphisms (SNPs) showed linkage to GABRB3 (rs11631421, LOD=4.6) and GABRA5 (rs2075716, LOD=2.2), and allelic association to GABRB3 (rs8024564, p=0.005; rs8025575, p=0.02) and GABRA5 (rs35399885, p=0.05). Genotyping of an independent Sardinian PD trio sample also supported association in the region, again with a parent-of-origin effect. These findings provide genetic evidence for the involvement of the genes GABRB3 and GABRA5 in the susceptibility to PD.

Application of next generation sequencing to CEPH cell lines to discover variants associated with FDA approved chemotherapeutics

BACKGROUND

The goal of this study was to perform candidate gene association with cytotoxicity of chemotherapeutics in cell line models through resequencing and discovery of rare and low frequency variants along with common variations. Here, an association study of cytotoxicity response to 30 FDA approved drugs was conducted and we applied next generation targeted sequencing technology to discover variants from 103 candidate genes in 95 lymphoblastoid cell lines from 14 CEPH pedigrees. In this article, we called variants across 95 cell lines and performed association analysis for cytotoxic response using the Family Based Association Testing method and software.

RESULTS

We called 2281 variable SNP genotypes across the 103 genes for these cell lines and identified three genes of significant association within this marker set. Specifically, ATP-binding cassette, sub-family C, member 5 (ABCC5), metallothionein 1A (MT1A) and NAD(P)H dehydrogenase quinone1 (NQO1) were significantly associated with oxaliplatin drug response. The significant SNP on NQO1 (rs1800566) has been linked with poor survival rates in patients with non-small cell lung cancer treated with cisplatin (which belongs to the same class of drugs as oxaliplatin). A SNP (rs1846692) near the 5' region of MT1A was associated with arsenic trioxide.

CONCLUSIONS

The results from this study are promising and this serves as a proof-of-principle demonstration of the use of sequencing data in the cytotoxicity models of human cell lines. With increased sample sizes, such studies will be a fast and powerful way to associate common and rare variants with drug response; while overcoming the cost and time limitations to recruit cohorts for association study.

"The race" to clone BRCA1

The existence of BRCA1 was proven in 1990 by mapping predisposition to young-onset breast cancer in families to chromosome 17q21. Knowing that such a gene existed and approximately where it lay triggered efforts by public and private groups to clone and sequence it. The press baptized the competition "the race" and reported on it in detail for the next 4 years. BRCA1 was positionally cloned in September 1994. Twenty years later, I reflect on "the race" and its consequences for breast cancer prevention and treatment.