Common Ancestry-specific Ion Channel Variants Predispose to Drug-induced Arrhythmias

Background: Multiple reports associate the cardiac sodium channel gene ( SCN5A ) variants S1103Y and R1193Q with type 3 congenital long QT syndrome (LQTS) and drug-induced LQTS. These variants are, however, too common in ancestral populations to be highly arrhythmogenic at baseline: S1103Y allele frequency is 8.1% in Africans and R1193Q 6.1% in East Asians. R1193Q is known to increase late sodium current (I Na-L ) in cardiomyocytes derived from induced pluripotent stem cells (iPSC-CMs) but the role of these variants in modulating repolarization remains poorly-understood. Methods: We determined the effect of S1103Y on QT intervals among Africans in a large electronic health record. Using iPSC-CMs carrying naturally occurring or genome-edited variants, we studied action potential durations (APDs) at baseline and after challenge with the repolarizing potassium current (I Kr ) blocker dofetilide, and I Na-L and I Kr at baseline. Results: In 1479 African subjects with no confounding medications or diagnoses of heart disease, QT in S1103Y carriers was no different from that in non-carriers. Similarly, baseline APD was no different in cells expressing the Y allele (SY, YY cells) compared to isogenic cells with the reference allele (SS cells). However, I Na-L was increased in SY and YY cells and the I Na-L blocker GS967 shortened APD in SY/YY but not SS cells (p<0.001). I Kr was increased almost 2-fold in SY/YY cells compared to SS cells (tail current: 0.66±0.1 vs 1.2±0.1 pA/pF, p<0.001). Dofetilide challenge prolonged APD at much lower concentrations in SY (4.1 nM [IQR 1.5-9.3], n=11) and YY (4.2 nM [1.7- 5.0], n=5) than in SS cells (249 nM [22.3-2905], n=14, p<0.001 and p<0.01, respectively) and elicited afterdepolarizations in 8/16 SY/YY cells but only in 1/14 SS cells. R1193Q cells similarly displayed no difference in baseline APD but increased I Kr and increased dofetilide sensitivity. Conclusions: These common ancestry-specific variants do not affect baseline repolarization, despite generating increased I Na-L . We propose that increased I Kr serves to maintain normal repolarization but increases the risk of manifest QT prolongation with I Kr block in variant carriers. Our findings further emphasize the need for inclusion of diverse populations in the study of adverse drug reactions.

Constitutively active SARM1 variants that induce neuropathy are enriched in ALS patients

Background In response to injury, neurons activate a program of organized axon self-destruction initiated by the NAD+ hydrolase, SARM1. In healthy neurons SARM1 is autoinhibited, but single amino acid changes can abolish autoinhibition leading to constitutively active SARM1 enzymes that promote degeneration when expressed in cultured neurons. Methods To investigate whether naturally occurring human variants might disrupt SARM1 autoinhibition and potentially contribute to risk for neurodegenerative disease, we assayed the enzymatic activity of all 42 rare SARM1 alleles identified among 8507 amyotrophic lateral sclerosis (ALS) patients and 9671 controls. We then intrathecally injected mice with virus expressing SARM1 constructs to test the capacity of an ALS-associated constitutively active SARM1 variant to promote neurodegeneration in vivo. Results Twelve out of 42 SARM1 missense variants or small in-frame deletions assayed exhibit constitutive NADase activity, including more than half of those that are unique to the ALS patients or that occur in multiple patients. There is a > 5-fold enrichment of constitutively active variants among patients compared to controls. Expression of constitutively active ALS-associated SARM1 alleles in cultured dorsal root ganglion (DRG) neurons is pro-degenerative and cytotoxic. Intrathecal injection of an AAV expressing the common SARM1 reference allele is innocuous to mice, but a construct harboring SARM1V184G, the constitutively active variant found most frequently among the ALS patients, causes axon loss, motor dysfunction, and sustained neuroinflammation. Conclusions These results implicate rare hypermorphic SARM1 alleles as candidate genetic risk factors for ALS and other neurodegenerative conditions.

Assessing Bos taurus introgression in the UOA Bos indicus assembly

Background Reference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus—B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus—B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references. Results The alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity. Conclusions We report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.

Continuous Mapping Identifies Loci Associated With Weevil Resistance [Cosmopolites sordidus (Germar)] in a Triploid Banana Population

The first step toward marker-assisted selection is linking the phenotypes to molecular markers through quantitative trait loci (QTL) analysis. While the process is straightforward in self-pollinating diploid (2x) species, QTL analysis in polyploids requires unconventional methods. In this study, we have identified markers associated with weevil Cosmopolites sordidus (Germar) resistance in bananas using 138 triploid (2n = 3x) hybrids derived from a cross between a tetraploid “Monyet” (2n = 4x) and a 2x “Kokopo” (2n = 2x) banana genotypes. The population was genotyped by Diversity Arrays Technology Sequencing (DArTSeq), resulting in 18,009 polymorphic single nucleotide polymorphisms (SNPs) between the two parents. Marker–trait association was carried out by continuous mapping where the adjusted trait means for the corm peripheral damage (PD) and total cross-section damage (TXD), both on the logit scale, were regressed on the marker allele frequencies. Forty-four SNPs that were associated with corm PD were identified on the chromosomes 5, 6, and 8, with 41 of them located on chromosome 6 and segregated in “Kokopo.” Eleven SNPs associated with corm total TXD were identified on chromosome 6 and segregated in “Monyet.” The additive effect of replacing one reference allele with the alternative allele was determined at each marker position. The PD QTL was confirmed using conventional QTL linkage analysis in the simplex markers segregating in “Kokopo” (AAAA × RA). We also identified 43 putative genes in the vicinity of the markers significantly associated with the two traits. The identified loci associated with resistance to weevil damage will be used in the efforts of developing molecular tools for marker-assisted breeding in bananas.

THE CORRELATION OF SUBSTITUTION EFFECTS ACROSS POPULATIONS AND GENERATIONS IN THE PRESENCE OF NON-ADDITIVE FUNCTIONAL GENE ACTION

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of non-additive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a “focal” population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance and additive by additive variances; relationships (Nei’s minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the non-additive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.

Possible Protective Effect of LOXL1 Variant in the Cohort of Chernobyl Catastrophe Clean-Up Workers

Ionising radiation (IR) is an environmental factor known to alter genomes and therefore challenge organisms to adapt. Lithuanian clean-up workers of the Chernobyl nuclear disaster (LCWC) experienced high doses of IR, leading to different consequences. This study aims to characterise a unique protective genomic variation in a relatively healthy LCWC group. This variation influenced their individual reaction to IR and potentially protects against certain diseases such as exfoliation syndrome and glaucoma. Clinical and IR dosage data were collected using a questionnaire to characterise the cohort of 93 LCWC. Genome-wide genotyping using Illumina beadchip technology was performed. The control group included 466 unrelated, self-reported healthy individuals of Lithuanian descent. Genotypes were filtered out from the microarray dataset using a catalogue of SNPs. The data were used to perform association, linkage disequilibrium, and epistasis analysis. Phenotype data analysis showed the distribution of the most common disease groups among the LCWC. A genomic variant of statistical significance (Fishers’ exact test, p = 0.019), rs3825942, was identified in LOXL1 (NM_005576.4:c.458G>A). Linkage disequilibrium and epistasis analysis for this variant identified the genes LHFPL3, GALNT6, PIH1D1, ANKS1B, and METRNL as potentially involved in the etiopathogenesis of exfoliation syndrome and glaucoma, which were not previously associated with the disease. The LOXL1 variant is mostly considered a risk factor in the development of exfoliation syndrome and glaucoma. The influence of recent positive selection, the phenomenon of allele-flipping, and the fact that only individuals with the homozygous reference allele have glaucoma in the cohort of the LCWC suggest otherwise. The identification of rs3825942 and other potentially protective genomic variants may be useful for further analysis of the genetic architecture and etiopathogenetic mechanisms of other multifactorial diseases.

A locus conferring tolerance to Theileria infection in African cattle

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Boran cattle demonstrating heritable tolerance to infection by T. parva (h2 = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on Bos taurus chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained polymorphism in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa's most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

Human genetic diversity modifies therapeutic gene editing off-target potential

CRISPR gene editing holds great promise to modify somatic genomes to ameliorate disease. In silico prediction of homologous sites coupled with biochemical evaluation of possible genomic off-targets may predict genotoxicity risk of individual gene editing reagents. However, standard computational and biochemical methods focus on reference genomes and do not consider the impact of genetic diversity on off-target potential. Here we developed a web application called CRISPRme that explicitly and efficiently integrates human genetic variant datasets with orthogonal genomic annotations to predict and prioritize off-target sites at scale. The method considers both single-nucleotide variants (SNVs) and indels, accounts for bona fide haplotypes, accepts spacer:protospacer mismatches and bulges, and is suitable for personal genome analyses. We tested the tool with a guide RNA (gRNA) targeting the BCL11A erythroid enhancer that has shown therapeutic promise in clinical trials for sickle cell disease (SCD) and β-thalassemia. We find that the top predicted off-target site is produced by a non-reference allele common in African-ancestry populations (rs114518452, minor allele frequency (MAF)=4.5%) that introduces a protospacer adjacent motif (PAM) for SpCas9. We validate that SpCas9 generates indels (~9.6% frequency) and chr2 pericentric inversions in a strictly allele-specific manner in edited CD34+ hematopoietic stem/progenitor cells (HSPCs), although a high-fidelity Cas9 variant mitigates this off-target. This report illustrates how genetic variation may modify the genomic outcomes of therapeutic gene editing and provides a simple tool for comprehensive off-target assessment.

Novel functional sequences uncovered through a bovine multiassembly graph

Many genomic analyses start by aligning sequencing reads to a linear reference genome. However, linear reference genomes are imperfect, lacking millions of bases of unknown relevance and are unable to reflect the genetic diversity of populations. This makes reference-guided methods susceptible to reference-allele bias. To overcome such limitations, we build a pangenome from six reference-quality assemblies from taurine and indicine cattle as well as yak. The pangenome contains an additional 70,329,827 bases compared to the Bos taurus reference genome. Our multiassembly approach reveals 30 and 10.1 million bases private to yak and indicine cattle, respectively, and between 3.3 and 4.4 million bases unique to each taurine assembly. Utilizing transcriptomes from 56 cattle, we show that these nonreference sequences encode transcripts that hitherto remained undetected from the B. taurus reference genome. We uncover genes, primarily encoding proteins contributing to immune response and pathogen-mediated immunomodulation, differentially expressed between Mycobacterium bovis–infected and noninfected cattle that are also undetectable in the B. taurus reference genome. Using whole-genome sequencing data of cattle from five breeds, we show that reads which were previously misaligned against the Bos taurus reference genome now align accurately to the pangenome sequences. This enables us to discover 83,250 polymorphic sites that segregate within and between breeds of cattle and capture genetic differentiation across breeds. Our work makes a so-far unused source of variation amenable to genetic investigations and provides methods and a framework for establishing and exploiting a more diverse reference genome.

DR2S: an integrated algorithm providing reference-grade haplotype sequences from heterozygous samples

Background High resolution HLA genotyping of donors and recipients is a crucially important prerequisite for haematopoetic stem-cell transplantation and relies heavily on the quality and completeness of immunogenetic reference sequence databases of allelic variation. Results Here, we report on , an R package that leverages the strengths of two sequencing technologies—the accuracy of next-generation sequencing with the read length of third-generation sequencing technologies like PacBio’s SMRT sequencing or ONT’s nanopore sequencing—to reconstruct fully-phased high-quality full-length haplotype sequences. Although optimised for HLA and KIR genes, is applicable to all loci with known reference sequences provided that full-length sequencing data is available for analysis. In addition, integrates supporting tools for easy visualisation and quality control of the reconstructed haplotype to ensure suitability for submission to public allele databases. Conclusions is a largely automated workflow designed to create high-quality fully-phased reference allele sequences for highly polymorphic gene regions such as HLA or KIR. It has been used by biologists to successfully characterise and submit more than 500 HLA alleles and more than 500 KIR alleles to the IPD-IMGT/HLA and IPD-KIR databases.