Identification of Monoallelically Expressed Genes Associated with Economic Traits in Hanwoo (Korean Native Cattle)

Hanwoo, an indigenous Korean cattle breed, has been genetically improved by selecting superior sires called Korean-proven bulls. However, cows still contribute half of the genetic stock of their offspring, and allelic-specific expressed genes have potential, as selective targets of cows, to enhance genetic gain. The aim of this study is to identify genes that have MAEs based on both the genome and transcriptome and to estimate their effects on breeding values (BVs) for economically important traits in Hanwoo. We generated resequencing data for the parents and RNA-sequencing data for the muscle, fat, and brain tissues of the offspring. A total of 3801 heterozygous single nucleotide polymorphisms (SNPs) in offspring were identified and they were located in 1569 genes. Only 14 genes showed MAE (seven expressing maternal alleles and seven expressing paternal alleles). Tissue-specific MAE was observed, and LANCL1 showed maternal allele expression across all tissues. MAE genes were enriched for the biological process of cell death and angiogenesis, which included ACKR3 and PDCL3 genes, whose SNPs were significantly associated with BVs of lean meat production-related traits, such as weight at 12 months of age, carcass weight, and loin eye area. In the current study, monoallelically expressed genes were identified in various adult tissues and these genes were associated with genetic capacity in Hanwoo.

Pklr Is a Genetic Modifier of Sickle Cell Disease

Background: Acute pain, the most prominent complication of sickle cell disease (SCD), results from vasoocclusion triggered by sickling of deoxygenated red blood cells (RBCs). A key factor influencing HbS oxygenation is the intracellular concentration of 2,3- diphosphoglycerate (2,3-DPG). 2,3-DPG, an intermediate substrate in the glycolytic pathway, decreases oxygen binding and stabilizes the deoxygenated hemoglobin. Pyruvate kinase (gene PKLR, protein PKR) is a rate-limiting enzyme in glycolysis; variants in PKLR may affect PKR activity, 2,3-DPG levels in RBCs, subsequent frequency of sickling and acute pain episodes (APE). There is thus a strong biological basis for exploring PKLR as a candidate gene affecting acute pain in SCD. Methods: The study population for genetic association consists of 2 cohorts: 1) 242 adults with HbSS from King's College Hospital (KCH), London, UK, with complete hospitalisation records over 10 years (2004-2013 inclusive) as the "discovery" cohort; 2) 977 children with HbSS or HbSb 0 thalassemia from the Silent Infarct Transfusion (SIT) trial, with a 3-year history of severe vasoocclusive pain based on hospitalization, as the "validation" cohort. Both studies were approved by the local Institutional Review Boards at KCH and Vanderbilt University Medical Center, respectively. An independent cohort comprises 52 adults with SCD enrolled under 3 protocols - NCT00011648, NCT00081523, and NCT03685721 - approved by the NHLBI Review Board (NIH), for evaluation of imbalance in allele expression. Genome scan for the KCH cohort was performed using llumina's Infinium "MEGA" chip (1.7m markers). The SIT DNA samples were genotyped using Illumina HumanHap650Y array 5 (661K markers) or Illumina Infinium HumanOmni1-Quad array (1.1m markers). The results were quality controlled followed by genotype imputation based on the 1000 Genomes Project phase 3 data. An annualised "hospitalisation rate" as a measure of pain incidence rate, was calculated by dividing the number of hospital admissions for severe acute pain by the number of years of observation for KCH and SIT cohorts (Fig A). We performed association analysis with common SNPs at PKLR locus using data from our genome-wide SNP set and a linear mixed modelling approach incorporating a genetic relatedness matrix to take account of relatedness, plus sex and age as fixed covariates. We corrected for multiple testing after quantifying the linkage disequilibrium (LD) within PKLR and used this to calculate appropriate significance levels. For the PKLR region and hospitalisation rate, the modified significance level was p<0.001268 for the discovery (KCH) cohort. For the allele expression assays, a synonymous variant, rs1052176 (R596R), in exon 11 of PKLR acted as a marker of relative expression levels of the 2 alleles of the gene. Allele specific expression was carried using the Bio-Rad digital droplet PCR system. Results: 7 of 47 variants evaluated in PKLR were associated with hospitalization rate (LnLnHospRate) in the discovery cohort: intron 4 - rs071053, and intron 2 - rs8177970, rs116244351, rs114455416, rs12741350, rs3020781, and rs8177964). All 7 were validated in Fisher's meta-analyses of the KCH and the SIT cohorts using p<0.0071 as threshold to correct for multiple testing (Fig B). We examined the pairwise LD between PKLR variants, and found all the intron 2 variants in tight LD, while R596R belongs to another LD block (Fig C). 52 SCD individuals had the R596R variant, of which 29 were heterozygous and 23 homozygous for the intron 2 haplotype associated with APE in SCD. We performed a Wilcoxon rank sum test and compared the variation in PKLR expression between the 2 alleles in subjects homozygous and heterozygous for the wildtype intron 2 haplotype, using genomic DNA as internal control for each subject. The results reveal a significant deviation from the expected expression ratio in those heterozygous for the intron 2 haplotype (mean 0.2073, +/- SD 0.0135) when compared with to those without the variant (mean 0.1239, +/- SD 0.0682), p=0.0297 (Fig D). Conclusion: Intronic variants of PKLR are associated with hospitalization rate for acute pain episodes in adults and children with SCD. We show that the intronic variants are likely to influence acute pain by affecting expression of the PKLR gene using allele-specific expression analyses, although the causal variant is unclear. These results support PKLR as a genetic modifier of SCD. Figure 1 Figure 1. Disclosures No relevant conflicts of interest to declare.

A simple strategy for sample annotation error detection in cytometry datasets

Mislabeling samples or data with the wrong participant information can impact study integrity and lead investigators to draw inaccurate conclusions. Quality control to prevent these types of errors is commonly embedded into the analysis of genomic datasets, but a similar identification strategy is not standard for cytometric data. Here, we present a method for detecting sample identification errors in cytometric data using expression of HLA class I alleles. We measured HLA-A*02 and HLA-B*07 expression in 3 longitudinal samples from 41 participants using a 33-marker CyTOF panel designed to identify major immune cell types. 3/123 samples (2.4%) showed HLA allele expression that did not match their longitudinal pairs. Furthermore, these same three samples cytometric signature did not match qPCR HLA class I allele data, suggesting that they were accurately identified as mismatches. We conclude that this technique is useful for detecting sample labeling errors in cytometric analyses of longitudinal data. This technique could also be used in conjunction with another method, like GWAS or PCR, to detect errors in cross-sectional data. We suggest widespread adoption of this or similar techniques will improve the quality of clinical studies that utilize cytometry.

SCReadCounts: estimation of cell-level SNVs expression from scRNA-seq data

Background Recent studies have demonstrated the utility of scRNA-seq SNVs to distinguish tumor from normal cells, characterize intra-tumoral heterogeneity, and define mutation-associated expression signatures. In addition to cancer studies, SNVs from single cells have been useful in studies of transcriptional burst kinetics, allelic expression, chromosome X inactivation, ploidy estimations, and haplotype inference. Results To aid these types of studies, we have developed a tool, SCReadCounts, for cell-level tabulation of the sequencing read counts bearing SNV reference and variant alleles from barcoded scRNA-seq alignments. Provided genomic loci and expected alleles, SCReadCounts generates cell-SNV matrices with the absolute variant- and reference-harboring read counts, as well as cell-SNV matrices of expressed Variant Allele Fraction (VAFRNA) suitable for a variety of downstream applications. We demonstrate three different SCReadCounts applications on 59,884 cells from seven neuroblastoma samples: (1) estimation of cell-level expression of known somatic mutations and RNA-editing sites, (2) estimation of cell- level allele expression of biallelic SNVs, and (3) a discovery mode assessment of the reference and each of the three alternative nucleotides at genomic positions of interest that does not require prior SNV information. For the later, we applied SCReadCounts on the coding regions of KRAS, where it identified known and novel somatic mutations in a low-to-moderate proportion of cells. The SCReadCounts read counts module is benchmarked against the analogous modules of GATK and Samtools. SCReadCounts is freely available (https://github.com/HorvathLab/NGS) as 64-bit self-contained binary distributions for Linux and MacOS, in addition to Python source. Conclusions SCReadCounts supplies a fast and efficient solution for estimation of cell-level SNV expression from scRNA-seq data. SCReadCounts enables distinguishing cells with monoallelic reference expression from those with no gene expression and is applicable to assess SNVs present in only a small proportion of the cells, such as somatic mutations in cancer.

HLA-Cw6 Allele Expression is Associated with Good Narrowband Ultraviolet B Response in Javanese-Indonesian Psoriasis Subjects

BACKGROUND: Psoriasis is an autoimmune disease involving genetic-environmental factors. Human Leukocyte Antigen (HLA)-Cw6 allele is the main genetic risk factor of psoriasis, but its prevalence varies widely. HLA-Cw6 allele correlates with psoriasis clinical type and treatment responses. Narrowband ultraviolet B (NB-UVB) and methotrexate (MTX) are effective psoriasis treatments but their association with HLA-Cw6 allele has not been identified. The study aims to determine the association between HLA-Cw6 allele expression, NB-UVB, and MTX treatment responses in Javanese-Indonesian psoriasis subjects.METHODS: Ninety Javanese-Indonesian psoriasis subjects were recruited in this study, 45 subjects were treated using NB UVB, while the other 45 subjects were treated using MTX, respectively, for 12 weeks. The psoriasis diagnosis and treatment responses were evaluated by dermatologists using Psoriasis Area Severity Index (PASI) score. The HLA-Cw6 allele was examined using the single-specific-primer polymerase chain reaction method. Fisher's Exact and Chi-square tests were employed, where p-value<0.05 was considered of significant association.RESULTS: The HLA-Cw6 allele positivity was identified in 23 psoriasis subjects (25.56%), while the other 67 subjects expressed HLA-Cw6 allele negative (74.44%). Female with HLA-Cw6 allele positivity who did not have comorbid disease show good response to NB-UVB than MTX. Meanwhile, subjects who were treated with MTX showed no association between therapeutic response and HLA-Cw6. HLA-Cw6 status was not correlated with the onset of psoriasis, family history, and comorbid diseases in all subjects.CONCLUSION: HLA-Cw6 allele expression is more associated with good NB-UVB treatment response than with MTX treatment response in female Javanese-Indonesian psoriasis subjects without comorbid disease.KEYWORDS: HLA-Cw6, narrowband ultraviolet B, methotrexate, Javanese, psoriasis

Allele expression biases in mixed-ploid sugarcane accessions

Allele-specific expression (ASE) represents differences in the magnitude of expression between alleles of the same gene. This is not straightforward for polyploids, especially autopolyploids, as knowledge about the dose of each allele is required for accurate estimation of ASE. This is the case for the genomically complex Saccharum species, characterized by high levels of ploidy and aneuploidy. We used a Beta-Binomial model to test for allelic imbalance in Saccharum, with adaptations for mixed-ploid organisms. The hierarchical Beta-Binomial model was used to test if allele expression followed the expectation based on genomic allele dosage. The highest frequencies of ASE occurred in sugarcane hybrids, suggesting a possible influence of interspecific hybridization in these genotypes. For all accessions, ASEGs were less frequent than those with balanced allelic expression. These genes were related to a broad range of processes, mostly associated with general metabolism, organelles, responses to stress and responses to stimuli. In addition, the frequency of ASEGs in high-level functional terms was similar among the genotypes, with a few genes associated with more specific biological processes. We hypothesize that ASE in Saccharum is largely a genotype-specific phenomenon, as a large number of ASEGs were exclusive to individual accessions.

Diversity of MHC IIB genes and parasitism in hybrids of evolutionarily divergent cyprinoid species indicate heterosis advantage

The genes of the major histocompatibility complex (MHC) are an essential component of the vertebrate immune system and MHC genotypes may determine individual susceptibility to parasite infection. In the wild, selection that favors MHC variability can create situations in which interspecies hybrids experience a survival advantage. In a wild system of two naturally hybridizing leuciscid fish, we assessed MHC IIB genetic variability and its potential relationships to hosts’ ectoparasite communities. High proportions of MHC alleles and parasites were species-specific. Strong positive selection at specific MHC codons was detected in both species and hybrids. MHC allele expression in hybrids was slightly biased towards the maternal species. Controlling for a strong seasonal effect on parasite communities, we found no clear associations between host-specific parasites and MHC alleles or MHC supertypes. Hybrids shared more MHC alleles with the more MHC-diverse parental species, but expressed intermediate numbers of MHC alleles and positively selected sites. Hybrids carried significantly fewer ectoparasites than either parent species, suggesting a hybrid advantage via potential heterosis.

Dynamic Spatial-Temporal Expression Ratio of X Chromosome to Autosomes but Stable Dosage Compensation in Mammals

In the evolutionary model of dosage compensation, per-allele expression level of the X chromosome was proposed to have two-fold upregulation, compensating for its dose reduction in males (XY) compared to females (XX). However, the upregulation of X chromosome is still in dispute, and comprehensive evaluations are still lacking. By integrating multi-omics datasets in mammals, we investigated the expression ratios and underlying pattern of X to autosomes (X:AA ratio) and X to orthologs (X:XX ratio) at the transcriptome, translatome, and proteome layers. The results indicated a dynamic spatial-temporal X:AA ratio during development in human and mouse. Meanwhile, by tracing the evolution of orthologous gene expressions in chicken, platypus, and opossum, we found a constant expression ratio between X-linked genes in human and their autosomal orthologs in other species (X:XX ~1) across tissues and developmental stages, demonstrating stable dosage compensation in mammals. We also revealed that different epigenetic regulations could shape the higher tissue- and stage-specificity of X-linked gene expression, and affect X:AA ratios. We conclude that the dynamics of X:AA ratios are attributed to the different gene contents and expression preferences of the X chromosome, instead of the stable dosage compensation.

Evidence for a shift in placental HLA-G allelic dominance and the HLA-G isoform profile during a healthy pregnancy and pre-eclampsia

Human leukocyte antigen (HLA)-G is a non-classical class Ib major expressed by placental trophoblast cells plays a central role in establishing tolerance to the semi-allogeneic fetus and in placentation. HLA-G exists in different soluble or membrane-bound isoforms. Pre-eclampsia, a major cause of fetal and maternal morbidity and mortality, has been linked to insufficient placentation and an altered immune response in pregnancy, including altered HLA-G expression. The 14 bp insertion/deletion polymorphism in the 3′ untranslated region of the gene and the isoform profile may affect HLA-G expression. The aim of the current pilot study was to characterize the expression patterns of HLAG mRNA, protein and isoform profile in uncomplicated term pregnancies and in cases of pre-eclampsia. Maternal sHLA-G mRNA and protein levels was slightly reduced in pre-eclampsia. No difference was found for placental blood, and no correlation between peripheral and placental sHLA-G levels was found. We observed no association between neither fetal nor maternal HLA-G 14 bp insertion/deletion genotypes and pre-eclampsia, nor a significant difference in isoform profiles. However, in HLA-G 14 bp insertion/deletion heterozygous placental samples, we observed abundant HLA-G1 14 bp insertion allele expression in the term placentae, which is contrary to previous findings in first trimester trophoblast. Increased HLA-G1 14 bp insertion allele expression in the placenta was associated with reduced levels of placental sHLA-G and an altered isoform profile with increased relative levels of HLA-G1 and -G5 and reduced levels of HLA-G3. The results indicate that an allelic shift in heterozygous individuals could represent a novel regulatory pathway.

Clinical Evidence for the Importance of the Wild-Type PRPF31 Allele in the Phenotypic Expression of RP11

PRPF31-associated retinopathy (RP11) is a common form of autosomal dominant retinitis pigmentosa (adRP) that exhibits wide variation in phenotype ranging from non-penetrance to early-onset RP. Herein, we report inter-familial and intra-familial variation in the natural history of RP11 using multimodal imaging and microperimetry. Patients were recruited prospectively. The age of symptom onset, best-corrected visual acuity, microperimetry mean sensitivity (MS), residual ellipsoid zone span and hyperautofluorescent ring area were recorded. Genotyping was performed using targeted next-generation and Sanger sequencing and copy number variant analysis. PRPF31 mutations were found in 14 individuals from seven unrelated families. Four disease patterns were observed: (A) childhood onset with rapid progression (N = 4), (B) adult-onset with rapid progression (N = 4), (C) adult-onset with slow progression (N = 4) and (D) non-penetrance (N = 2). Four different patterns were observed in a family harbouring c.267del; patterns B, C and D were observed in a family with c.772_773delins16 and patterns A, B and C were observed in 3 unrelated individuals with large deletions. Our findings suggest that the RP11 phenotype may be related to the wild-type PRPF31 allele rather than the type of mutation. Further studies that correlate in vitro wild-type PRPF31 allele expression level with the disease patterns are required to investigate this association.