Family-based association of HLA-DRB1 and DQB1 alleles and haplotypes in a group of Iranian Type 1 diabetes children

This family-based study was conducted in a group of Iranians with Type 1 diabetes (T1D) to investigate the transmission from parents of risk and non-risk HLA alleles and haplotypes, and to estimate the genetic risk score for this disease within this population. A total of 240 T1D subjects including 111 parent-child trios (111 children with T1D, 133 siblings, and 222 parents) and 330 ethnically matched healthy individuals were recruited. High-resolution HLA typing for DRB1/DQB1 loci was performed for all study subjects (n = 925) using polymerase chain reaction-sequence-specific oligonucleotide probe method. The highest predisposing effect on developing T1D was conferred by the following haplotypes both in all subjects and in probands compared to controls: DRB1*04:05-DQB1*03:02 (Pc = 2.97e-06 and Pc = 6.04e-10, respectively), DRB1*04:02-DQB1*03:02 (Pc = 5.94e-17 and Pc = 3.86e-09, respectively), and DRB1*03:01-DQB1*02:01 (Pc = 8.26e-29 and Pc = 6.56e-16, respectively). Conversely, the major protective haplotypes included DRB1*13:01-DQB1*06:03 (Pc = 6.99e-08), DRB1*15:01-DQB1*06:02 (Pc = 2.97e-06) in the cases versus controls. Also, DRB1*03:01-DQB1*02:01/DRB1*04:02|05-DQB1*03:02 and DRB1*03:01-DQB1*02:01/DRB1*03:01-DQB1*02:01 diplotypes conferred the highest predisposing effect in the cases (Pc = 8.65e-17 and Pc = 6.26e-08, respectively) and in probands (Pc = 5.4e-15 and Pc = 0.001, respectively) compared to controls. Transmission disequilibrium test showed that the highest risk was conferred by DRB1*04:02-DQB1*03:02 (Pc = 3.26e-05) and DRB1*03:01-DQB1*02:01 (Pc = 1.78e-12) haplotypes and the highest protection by DRB1*14:01-DQB1*05:03 (Pc = 8.66e-05), DRB1*15:01-DQB1*06:02 (Pc = 0.002), and DRB1*11:01-DQB1*03:01 (Pc = 0.0003) haplotypes. Based on logistic regression analysis, carriage of risk haplotypes increased the risk of T1D development 24.5 times in the Iranian population (p = 5.61e-13). Also, receiver operating characteristic curve analysis revealed a high predictive power of those risk haplotypes in discrimination of susceptible from healthy individuals (area under curve: 0.88, p = 5.5e-32). Our study highlights the potential utility of genetic risk assessment based on HLA diplotypes for predicting T1D risk in individuals, particularly among family members of affected children in our population.

The novel HLA-DQA1*05:05:17:03 allele, identified in a potential organ donor

HLA-DQA1*05:05:17:03 differs from HLA-DQA1*05:05:01:02 by a single base substitution in exon 1 and HLA-DQA1*05:05:17:01 within introns 1 and 2.

Clinical application of a real-time polymerase chain reaction test for CYP2C19 genotyping based on genotype distribution in a healthy Korean population

OBJECTIVE

With the recent reports of additional alleles of the CYP2C19 gene with decreased or no function, the clinical utility of real-time polymerase chain reaction (PCR)-based testing that detects only a small number of variant targets needs to be evaluated.

METHOD

We retrospectively reviewed 7-year data for real-time PCR test records from a single hospital and analyzed CYP2C19 genotypes from publicly available whole-genome or whole-exome data from a healthy Korean population.

RESULTS

Of the 2327 test results in this hospital, the *1 allele was most common (60.5%), followed by *2 (28.0%), *3 (10.1%), and *17 (1.4%). Among 5305 healthy Korean individuals, the frequencies of the *2, *3, and *17 alleles were 28.6%, 9.9%, and 1.0%, respectively, which were not statistically different from those of the hospital data (P = .4439, P = .6025, and P = .1142, respectively). Interestingly, the total frequency of additional nonfunctional alleles (*4, *6, *22, and *24) that could not be detected using real-time PCR was only 0.1%, with a total allele count of 8.

CONCLUSION

Our study shows that the clinical utility of real-time PCR for CYP2C19 genotyping remains satisfactory. However, caution should be exercised because the test can miss patients with decreased CYP2C19 function.

The novel HLA-C*06:372 allele, identified in a stem cell donor of an old order mennonite ethnicity

HLA-C*06:372 differs from HLA-C*06:02:01:01 by a single substitution in exon 4.

Genetic Association of ICAM-1 (rs5498) Gene Polymorphism With Susceptibility to Stage II Grade B Periodontitis: A Case-Control Study in South Indian Population

INTRODUCTION

In the contemporary perspective, periodontitis is considered a complex issue triggered and perpetuated by bacteria but strongly influenced by the way the body reacts to bacterial plaque. Recent research has indicated that variations in genes might have an impact on the development of periodontitis. This study was conducted to explore a probable link between the genetic variations in intercellular adhesion molecule-1 (ICAM-1) represented by rs5498 and the occurrence of periodontitis.  Methods: A total of 100 participants, 50 with periodontitis and 50 with periodontally healthy or mild gingivitis, were recruited for this study. Whole blood drawn from the participants was used to obtain genomic DNA. The ICAM-1 gene polymorphism (rs5498) was determined using polymerase chain reaction (PCR) amplification and digestion. The ICAM-1 gene's flanking primers were used to amp up the DNA. For statistical analysis, the genotype that was analyzed using the pattern of restriction fragment length polymorphism was recorded. The Chi-square test compared genotype and allele frequency distributions between both groups. The odds ratio with 95% confidence intervals with each individual allele or genotype was used to compute the risk. Statistical significance was established in all tests when the p-value was less than 0.05.

RESULTS

 There was no discernible difference between the genotype frequencies of patients and controls χ2df (P = 0.6065). The findings demonstrated that no significant difference was present between the two groups for homozygous or heterozygous mutant genotypes (AA vs. AG+GG; P = 0.6854). There was no discernible difference in the detected frequencies of the A allele (58% vs. 61%), G allele (42% vs. 39%), TT (16% vs. 24%), AG (40% vs. 36%), and TT genotypes in the studied groups.

CONCLUSION

According to the results of the current investigation, the ICAM-1 (rs5498) gene polymorphism is not associated with periodontitis in the population investigated.

A novel HLA-C*02:10 variant, HLA-C*02:10:09, identified in a healthy individual from Brazil

Identification of the novel HLA-C*02:10:09 allele that differs from HLA-C*02:10:01:01 at one position in exon 1.

The novel HLA-DRB1*08:127 allele, identified in a deceased organ donor of Scandinavian descent

HLA-DRB1*08:127 differs from HLA-DRB1*08:01:01:01 by a single substitution in exon 4.

Addressing Ancestry and Sex Bias in Pharmacogenomics

The association of an individual's genetic makeup with their response to drugs is referred to as pharmacogenomics. By understanding the relationship between genetic variants and drug efficacy or toxicity, we are able to optimize pharmacological therapy according to an individual's genotype. Pharmacogenomics research has historically suffered from bias and underrepresentation of people from certain ancestry groups and of the female sex. These biases can arise from factors such as drugs and indications studied, selection of study participants, and methods used to collect and analyze data. To examine the representation of biogeographical populations in pharmacogenomic data sets, we describe individuals involved in gene-drug response studies from PharmGKB, a leading repository of drug-gene annotations, and showcaseCYP2D6, a gene that metabolizes approximately 25% of all prescribed drugs. We also show how the historical underrepresentation of females in clinical trials has led to significantly more adverse drug reactions in females than in males.

Natural and artificial selection of multiple alleles revealed through genomic analyses

Genome-to-phenome research in agriculture aims to improve crops through in silico predictions. Genome-wide association study (GWAS) is potent in identifying genomic loci that underlie important traits. As a statistical method, increasing the sample quantity, data quality, or diversity of the GWAS dataset positively impacts GWAS power. For more precise breeding, concrete candidate genes with exact functional variants must be discovered. Many post-GWAS methods have been developed to narrow down the associated genomic regions and, ideally, to predict candidate genes and causative mutations (CMs). Historical natural selection and breeding-related artificial selection both act to change the frequencies of different alleles of genes that control phenotypes. With higher diversity and more extensive GWAS datasets, there is an increased chance of multiple alleles with independent CMs in a single causal gene. This can be caused by the presence of samples from geographically isolated regions that arose during natural or artificial selection. This simple fact is a complicating factor in GWAS-driven discoveries. Currently, none of the existing association methods address this issue and need to identify multiple alleles and, more specifically, the actual CMs. Therefore, we developed a tool that computes a score for a combination of variant positions in a single candidate gene and, based on the highest score, identifies the best number and combination of CMs. The tool is publicly available as a Python package on GitHub, and we further created a web-based Multiple Alleles discovery (MADis) tool that supports soybean and is hosted in SoyKB (https://soykb.org/SoybeanMADisTool/). We tested and validated the algorithm and presented the utilization of MADis in a pod pigmentation L1 gene case study with multiple CMs from natural or artificial selection. Finally, we identified a candidate gene for the pod color L2 locus and predicted the existence of multiple alleles that potentially cause loss of pod pigmentation. In this work, we show how a genomic analysis can be employed to explore the natural and artificial selection of multiple alleles and, thus, improve and accelerate crop breeding in agriculture.

An intermediate-sized donor registry experience: HLA barriers in matching procedures

The data enabling the estimation of the possibility of finding a matched unrelated donor (MUD) within a relatively short time is important for the success of hematopoietic stem cell transplantation (HSCT). In the present study, 738 unrelated Croatian patients in the program of unrelated HSCT were retrospectively analyzed for gender matching, donor origin (national or international), the distribution of HLA alleles and haplotypes, as well as for the probability of finding a 9-10/10 MUD. Almost 70% of the patients in our study group had a 10/10 MUD, while among the patients with a 9/10 MUD, a 1st field resolution level mismatched donor was selected for 55.0% of patients. The majority of pairs were HLA-A mismatched (33.8%). A comparison of HLA allele frequencies between two subgroups of patients revealed significant differences for 13 alleles. However, after p value correction, the difference in frequency remained significant only for four alleles; three HLA alleles (B*08:01, C*07:01, and DRB1*03:01) demonstrated a significantly higher frequency among patients with a 10/10 MUD (Pcorr < 0.0001, Pcorr = 0.0096, and Pcorr < 0.0001, respectively), while the B*35:08 allele was significantly more present among patients with a 9/10 MUD (Pcorr = 0.0328). The comparison of the distribution of HLA haplotypes between patients with a 10/10 MUD and patients with a 9/10 MUD showed significant differences for a number of two-locus and three-locus haplotypes, as well as for one five-locus haplotype (HLA-A*01:01~B*08:01~C*07:01~DRB1*03:01~DQB1*02:01), which was significantly more present in the group of patients with a 10/10 MUD. At least one HLA haplotype from the group of non-frequent HLA haplotypes (positions >1000) was carried by patients with a 9/10 MUD. The data obtained by the present study will contribute to a better estimation of the probability of finding a suitable 9-10/10 MUD for Croatian patients in need of HSCT.

Hypertension Trends and White Matter Brain Injury in the Offspring Framingham Heart Study Cohort

BACKGROUND

Hypertension is the most potent stroke risk factor and is also related to cerebral small vessel disease. We studied the relation between mid-to-late-life hypertension trends and cerebral white matter injury in community-dwelling individuals from the FHS (Framingham Heart Study).

METHODS

FHS Offspring cohort participants with available mid-life and late-life blood pressure measurements and brain magnetic resonance imaging were included. Multiple regression analyses were used to relate hypertension trends (normotension-normotension [reference], normotension-hypertension, and hypertension-hypertension) to white matter injury metrics on diffusion tensor imaging (free water, fractional anisotropy, and peak skeletonized mean diffusivity) and Fluid Attenuated Inversion Recovery (white matter hyperintensity volume) by different blood pressure cutoffs (130/80, 140/90, and 150/90 mm Hg).

RESULTS

We included 1018 participants (mean age 47.3±7.4 years at mid-life and 73.2±7.3 at late-life). At the 140/90 mm Hg cutoff, the hypertension-hypertension trend was associated with higher free water (β, 0.16 [95% CI, 0.03-0.30]; P=0.021) and peak skeletonized mean diffusivity (β, 0.15 [95% CI, 0.01-0.29]; P=0.033). At a 130/80 mm Hg cutoff, the hypertension-hypertension trend had significantly higher free water (β, 0.16 [95% CI, 0.01-0.30]; P=0.035); and the normotension-hypertension (β, 0.24 [95% CI, 0.03-0.44]; P=0.027) and hypertension-hypertension (β, 0.22 [95% CI, 0.04-0.41]; P=0.022) trends had significantly increased white matter hyperintensity volume. Exploratory stratified analysis showed effect modifications by APOE ɛ4 allele and age.

CONCLUSIONS

Mid-to-late-life hypertension exposure is significantly associated with microstructural and to a lesser extent, visible white matter injury; the effects are observed at both conventional and lower blood pressure cutoffs and are associated with longer duration of hypertension.

Nonsense Variant PRDM16-Q187X Causes Impaired Myocardial Development and TGF-β Signaling Resulting in Noncompaction Cardiomyopathy in Humans and Mice

BACKGROUND

PRDM16 plays a role in myocardial development through TGF-β (transforming growth factor-beta) signaling. Recent evidence suggests that loss of PRDM16 expression is associated with cardiomyopathy development in mice, although its role in human cardiomyopathy development is unclear. This study aims to determine the impact of PRDM16 loss-of-function variants on cardiomyopathy in humans.

METHODS

Individuals with PRDM16 variants were identified and consented. Induced pluripotent stem cell-derived cardiomyocytes were generated from a proband hosting a Q187X nonsense variant as an in vitro model and underwent proliferative and transcriptional analyses. CRISPR (clustered regularly interspaced short palindromic repeats)-mediated knock-in mouse model hosting the Prdm16Q187X allele was generated and subjected to ECG, histological, and transcriptional analysis.

RESULTS

We report 2 probands with loss-of-function PRDM16 variants and pediatric left ventricular noncompaction cardiomyopathy. One proband hosts a PRDM16-Q187X variant with left ventricular noncompaction cardiomyopathy and demonstrated infant-onset heart failure, which was selected for further study. Induced pluripotent stem cell-derived cardiomyocytes prepared from the PRDM16-Q187X proband demonstrated a statistically significant impairment in myocyte proliferation and increased apoptosis associated with transcriptional dysregulation of genes implicated in cardiac maturation, including TGF-β-associated transcripts. Homozygous Prdm16Q187X/Q187X mice demonstrated an underdeveloped compact myocardium and were embryonically lethal. Heterozygous Prdm16Q187X/WT mice demonstrated significantly smaller ventricular dimensions, heightened fibrosis, and age-dependent loss of TGF-β expression. Mechanistic studies were undertaken in H9c2 cardiomyoblasts to show that PRDM16 binds TGFB3 promoter and represses its transcription.

CONCLUSIONS

Novel loss-of-function PRDM16 variant impairs myocardial development resulting in noncompaction cardiomyopathy in humans and mice associated with altered TGF-β signaling.

Association of matrix metalloproteinase-2 gene variants with diabetic nephropathy risk

BACKGROUND

Diabetic nephropathy is a highly destructive microvascular complication of diabetes. Genetic predisposition is involved in the pathogenesis of diabetic nephropathy, with multiple allelic polymorphisms associated with the development and progression of the disease, thereby increasing the overall risk. To date, no study is available that shows the association of matrix metalloproteinase-2 (MMP-2) gene polymorphisms with diabetic nephropathy risk. Thus, we investigated the potential genetic influence of MMP-2 promoter variants in the development of diabetic nephropathy in type 2 diabetic patients.

METHODS

In total, 726 type 2 diabetic patients and 310 healthy controls were included in the study and genotyped for MMP-2, -1306C/T, -790T/G, -1575G/T and -735C/T by real-time PCR. The analysis of the outcomes was performed assuming three genetic models. The threshold for statistical significance was set at 0.05.

RESULTS

The results showed that the minor allele frequency of the -790T/G variant was significantly higher in patients with and without nephropathy compared to controls. Furthermore, the distribution analysis revealed a significant association of the -790T/G variant, in all genetic models, with increased risk of diabetic nephropathy that persisted after adjusting for key covariates. No significant associations between MMP-2, -1306C/T, -1575G/T, -735C/T and the risk of diabetic nephropathy were detected. Haplotype analysis identified two risk haplotypes GCGC and GTAC associated with diabetic nephropathy.

CONCLUSIONS

The present study is the first to demonstrate the allelic and genotypic association of the MMP-2-790T/G variant and two haplotypes with an increased risk of diabetic nephropathy in a Tunisian population with type 2 diabetes.

Dual Antiplatelet Therapies and Causes in Minor Stroke or Transient Ischemic Attack: A Prespecified Analysis in the CHANCE-2 Trial

BACKGROUND

It is unclear whether patients with different stroke/transient ischemic attack etiologies benefit differently from gene-directed dual antiplatelet therapy. This study explored the efficacy and safety of ticagrelor-aspirin versus clopidogrel-aspirin in transient ischemic attack or minor stroke with different causes in the CHANCE-2 trial (Clopidogrel in High-Risk Patients With Acute Nondisabling Cerebrovascular Events-II).

METHODS

This was a prespecified analysis of the CHANCE-2 trial, which enrolled 6412 patients with minor stroke or transient ischemic attack who carried CYP2C19 loss-of-function alleles. Patients with centralized evaluation of TOAST (Trial of ORG 10172 in Acute Stroke Treatment) classification of large-artery atherosclerosis, small-vessel occlusion, and stroke of undetermined cause were included. The primary efficacy outcome was new stroke, and the primary safety outcome was severe or moderate bleeding, both within 90 days. Cox proportional hazards models were used to assess the interaction of TOAST classification with the effects of dual antiplatelet therapy with ticagrelor-aspirin versus clopidogrel-aspirin.

RESULTS

A total of 6336 patients were included in this study. In patients administered ticagrelor-aspirin and clopidogrel-aspirin, respectively, stroke recurred in 85 (9.8%) and 88 (10.7%) patients with large-artery atherosclerosis (hazard ratio, 0.86 [95% CI, 0.63-1.18]; P=0.34); 32 (3.6%) and 61 (7.0%) patients with small-vessel occlusion (hazard ratio, 0.51 [95% CI, 0.33-0.79]; P=0.002); and 68 (4.8%) and 87 (5.9%) patients with stroke of undetermined cause (hazard ratio, 0.80 [95% CI, 0.58-1.10]; P=0.17), with P=0.08 for the treatment×cause subtype interaction effect. There were no significant differences in severe or moderate bleeding events in patients with different cause and different treatment.

CONCLUSIONS

In this prespecified analysis of the CHANCE-2 trial, the efficacy and safety of ticagrelor-aspirin versus clopidogrel-aspirin in preventing new stroke were consistent in patients with different causes. The influence of stroke cause on benefit of gene-guided antiplatelet therapy should be explored by further trials.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT04078737.

Conceptualizing Epigenetics and the Environmental Landscape of Autism Spectrum Disorders

Complex interactions between gene variants and environmental risk factors underlie the pathophysiological pathways in major psychiatric disorders. Autism Spectrum Disorder is a neuropsychiatric condition in which susceptible alleles along with epigenetic states contribute to the mutational landscape of the ailing brain. The present work reviews recent evolutionary, molecular, and epigenetic mechanisms potentially linked to the etiology of autism. First, we present a clinical vignette to describe clusters of maladaptive behaviors frequently diagnosed in autistic patients. Next, we microdissect brain regions pertinent to the nosology of autism, as well as cell networks from the bilateral body plan. Lastly, we catalog a number of pathogenic environments associated with disease risk factors. This set of perspectives provides emerging insights into the dynamic interplay between epigenetic and environmental variation in the development of Autism Spectrum Disorders.

Relation of resistin gene variants to resistin plasma levels and altered susceptibility to polycystic ovary syndrome: A case control study

BACKGROUND

A role for resistin in the pathogenesis of polycystic ovarian syndrome (PCOS) and related features were described for various ethnicities. As its expression is partly inherited, a role for RETN polymorphisms in regulating resistin levels and PCOS risk was shown, but with varied results.

AIM

To investigate the association of rs34124816 (-537A>C), rs1862513 (-420C>G), rs3219175 (-358G>A), rs3745367 (+299G>A), rs3745369 (+1263G>C), and rs1423096 (+4965C>T) RETN SNPs with PCOS.

METHODS

Study subjects included 583 women with PCOS, and 713 eumenorrheic women serving as controls. Genotyping was done by real-time PCR.

RESULTS

Higher minor allele frequency (MAF) of rs34124816, rs3219175, and rs3745369, and lower MAF of rs1862513 and rs1423096 were seen in PCOS cases. Reduced PCOS risk was found with rs3745367 minor-allele homozygotes and rs1423096 minor-allele homozygotes, while increased risk was linked with rs3745367 heterozygotes, and with rs3745369 heterozygotes and minor-allele homozygotes. While it did not reach statistical significance, serum resistin levels were elevated in PCOS cases than in control women and major-allele homozygotes of rs34124816 and rs1862513, and in rs1423096 minor-allele-containing carriers. Carriage of rs34124816 correlated positively with age and LH, whereas rs1862513 positively and rs3745367 negatively correlated with fasting glucose. Six-locus (rs34124816-rs1862513-rs3219175-rs3745367-rs3745369-rs1423096) haplotype analysis demonstrated a significant reduction in AGGGGG and a marked increase in AGGGCG haplotypes between cases and controls, thus assigning PCOS protective and susceptible nature to these haplotypes, respectively.

CONCLUSIONS

This study is the first to document the contribution of rs34124816 and rs1423096 RETN variants to the risk of PCOS. The varied association of RETN gene variants with PCOS suggests an ethnic contribution of RETN association with PCOS.

Psychological Stress Is Associated with Increased Cancer Risk in Dogs

Although there is evidence that psychological stress may be associated with increased cancer risk, the effect of stress on cancer risk is difficult to study, both in humans, due to socioeconomic factors, and in animal models, due to questionable biological relevance. Here, we test whether heritable canine temperament that increases psychological stress is associated with cancer risk. The study data are breed-specific averages of incidences of multiple cancer types and of temperament classes. The latter are derived from a latent class analysis of behavioral questionnaires completed by owners (C-BARQ). We thus classified the dogs according to whether they are calm vs. reactive within and across breeds. Using meta-analysis approaches, we modeled the risk of multiple cancer types in calm vs. reactive dogs. We adjusted for breed averages of body mass and lifespan, which are common confounders that impact cancer. Our study confirms that body size has a significant effect of on risk of multiple types of cancers in dogs and shows for the first time that temperament also has a moderate effect. These findings suggest dog models of heritable psychological stress are suitable for molecular epidemiological and translational studies on its effects on cancer risk.

Contribution of Previously Unrecognized RNA Splice-Altering Variants to Congenital Heart Disease

BACKGROUND

Known genetic causes of congenital heart disease (CHD) explain <40% of CHD cases, and interpreting the clinical significance of variants with uncertain functional impact remains challenging. We aim to improve diagnostic classification of variants in patients with CHD by assessing the impact of noncanonical splice region variants on RNA splicing.

METHODS

We tested de novo variants from trio studies of 2649 CHD probands and their parents, as well as rare (allele frequency, <2×10-6) variants from 4472 CHD probands in the Pediatric Cardiac Genetics Consortium through a combined computational and in vitro approach.

RESULTS

We identified 53 de novo and 74 rare variants in CHD cases that alter splicing and thus are loss of function. Of these, 77 variants are in known dominant, recessive, and candidate CHD genes, including KMT2D and RBFOX2. In 1 case, we confirmed the variant's predicted impact on RNA splicing in RNA transcripts from the proband's cardiac tissue. Two probands were found to have 2 loss-of-function variants for recessive CHD genes HECTD1 and DYNC2H1. In addition, SpliceAI-a predictive algorithm for altered RNA splicing-has a positive predictive value of ≈93% in our cohort.

CONCLUSIONS

Through assessment of RNA splicing, we identified a new loss-of-function variant within a CHD gene in 78 probands, of whom 69 (1.5%; n=4472) did not have a previously established genetic explanation for CHD. Identification of splice-altering variants improves diagnostic classification and genetic diagnoses for CHD.

REGISTRATION

URL: https://clinicaltrials.gov; Unique identifier: NCT01196182.

Mutations in the Second Alternative Oxidase Gene: A New Approach to Group Aspergillus niger Strains

Alternative oxidase is a terminal oxidase in the branched mitochondrial electron transport chain of most fungi including Aspergillus niger (subgenus Circumdati, section Nigri). A second, paralogous aox gene (aoxB) is extant in some A. niger isolates but also present in two divergent species of the subgenus Nidulantes-A. calidoustus and A. implicatus-as well as in Penicillium swiecickii. Black aspergilli are cosmopolitan opportunistic fungi that can cause diverse mycoses and acute aspergillosis in immunocompromised individuals. Amongst the approximately 75 genome-sequenced A. niger strains, aoxB features considerable sequence variation. Five mutations were identified that rationally affect transcription or function or terminally modify the gene product. One mutant allele that occurs in CBS 513.88 and A. niger neotype strain CBS 554.65 involves a chromosomal deletion that removes exon 1 and intron 1 from aoxB. Another aoxB allele results from retrotransposon integration. Three other alleles result from point mutations: a missense mutation of the start codon, a frameshift, and a nonsense mutation. A. niger strain ATCC 1015 has a full-length aoxB gene. The A. niger sensu stricto complex can thus be subdivided into six taxa according to extant aoxB allele, which may facilitate rapid and accurate identification of individual species.

Allelic resolution of insect and spider silk genes reveals hidden genetic diversity

Arthropod silk is vital to the evolutionary success of hundreds of thousands of species. The primary proteins in silks are often encoded by long, repetitive gene sequences. Until recently, sequencing and assembling these complex gene sequences has proven intractable given their repetitive structure. Here, using high-quality long-read sequencing, we show that there is extensive variation-both in terms of length and repeat motif order-between alleles of silk genes within individual arthropods. Further, this variation exists across two deep, independent origins of silk which diverged more than 500 Mya: the insect clade containing caddisflies and butterflies and spiders. This remarkable convergence in previously overlooked patterns of allelic variation across multiple origins of silk suggests common mechanisms for the generation and maintenance of structural protein-coding genes. Future genomic efforts to connect genotypes to phenotypes should account for such allelic variation.

Patients with Non-Alcoholic Fatty Liver Disease and Alcohol Dehydrogenase 1B/Aldehyde Dehydrogenase 2 Mutant Gene Have Higher Values of Serum Alanine Transaminase

Patients with non-alcoholic fatty liver disease (NAFLD) share similar pathophysiologies to those of patients with alcohol liver disease. Alcoholic metabolic enzyme-related genes (alcohol dehydrogenase 1B (ADH1B) and aldehyde dehydrogenase 2 (ALDH2)) may be associated with pathophysiology in NAFLD patients. In this study, the association between ADH1B/ALDH2 gene polymorphism and serum metabolic factors, body statures, and hepatic steatosis/fibrosis status was evaluated in patients with NAFLD. Using biochemistry data, abdominal ultrasonography, fibrosis evaluation (Kpa), and steatosis evaluation (CAP), ADH1B gene SNP rs1229984 and ALDH2 gene SNP rs671 polymorphism were analyzed in sixty-six patients from 1 January 2022 to 31 December 2022. The percentage of the mutant type (GA + AA) was 87.9% (58/66) in the ADH1B allele and 45.5% (30/66) in the ALDH2 allele. Patients with the mutant-type ADH1B/ALDH2 allele had higher values of alanine aminotransferase (ALT) than the wild type (β = 0.273, p = 0.04). No association was observed between body mass index, serum metabolic factors (sugar and lipid profile), CAP, kPa, and ADH1B/ALDH2. A high proportion of the mutant-type ADH1B allele (87.9%) and ALDH2 allele (45.5%) was observed in patients with NAFLD. No association was observed between ADH1B/ALDH2 allele, BMI, and hepatic steatosis/fibrosis. Patients with the mutant-type ADH1B/ALDH2 allele had higher values of ALT than those with the wild type.

Oligogenic Architecture of Rare Noncoding Variants Distinguishes 4 Congenital Heart Disease Phenotypes

BACKGROUND

Congenital heart disease (CHD) is highly heritable, but the power to identify inherited risk has been limited to analyses of common variants in small cohorts.

METHODS

We performed reimputation of 4 CHD cohorts (n=55 342) to the TOPMed reference panel (freeze 5), permitting meta-analysis of 14 784 017 variants including 6 035 962 rare variants of high imputation quality as validated by whole genome sequencing.

RESULTS

Meta-analysis identified 16 novel loci, including 12 rare variants, which displayed moderate or large effect sizes (median odds ratio, 3.02) for 4 separate CHD categories. Analyses of chromatin structure link 13 of the genome-wide significant loci to key genes in cardiac development; rs373447426 (minor allele frequency, 0.003 [odds ratio, 3.37 for Conotruncal heart disease]; P=1.49×10^-8) is predicted to disrupt chromatin structure for 2 nearby genes BDH1 and DLG1 involved in Conotruncal development. A lead variant rs189203952 (minor allele frequency, 0.01 [odds ratio, 2.4 for left ventricular outflow tract obstruction]; P=1.46×10^-8) is predicted to disrupt the binding sites of 4 transcription factors known to participate in cardiac development in the promoter of SPAG9. A tissue-specific model of chromatin conformation suggests that common variant rs78256848 (minor allele frequency, 0.11 [odds ratio, 1.4 for Conotruncal heart disease]; P=2.6×10^-8) physically interacts with NCAM1 (P_FDR=1.86×10^-27), a neural adhesion molecule acting in cardiac development. Importantly, while each individual malformation displayed substantial heritability (observed h2 ranging from 0.26 for complex malformations to 0.37 for left ventricular outflow tract obstructive disease) the risk for different CHD malformations appeared to be separate, without genetic correlation measured by linkage disequilibrium score regression or regional colocalization.

CONCLUSIONS

We describe a set of rare noncoding variants conferring significant risk for individual heart malformations which are linked to genes governing cardiac development. These results illustrate that the oligogenic basis of CHD and significant heritability may be linked to rare variants outside protein-coding regions conferring substantial risk for individual categories of cardiac malformation.

A report on a modified protocol for flow cytometry-based assessment of blood group erythrocyte antigens potentially suitable for analysis of weak ABO subgroups

BACKGROUND

Flow cytometry (FC) has proven its utility in scrutinizing AB antigen expression in red blood cells (RBCs), cooperating with serological tests for accurate blood group typing. However, technical difficulties may impair the characterization of weak ABO subtypes when background noises appear at non-negligible levels.

STUDY DESIGN AND METHODS

We sought to establish an FC method that could prevent antibody-induced hemagglutination and an increase in cellular autofluorescence, two major issues inherent to RBC-FC analysis of AB expression. We optimized fixatives, multicolor-staining protocols, and sequential gating strategies. Blood samples from weak ABO subtype cases, B_m and A_el , were analyzed with the established protocol.

RESULTS

The optimized mixture of glutaraldehyde and formaldehyde successfully generated fixed RBCs resistant to agglutination while maintaining low autofluorescence. These features allowed co-staining of leukocyte- and erythrocyte-markers, which enabled sequential gating strategies facilitating the precise AB antigen analysis in purely single RBCs with minimum background noises. By the established FC analysis, we could detect in the B_m sample a small RBC population exhibiting weak B antigen expression. The assay also proved it feasible to identify a small population (0.04%) of RBCs weakly expressing the A antigen in the A_el sample confirmed as harboring a rare c.816dupG ABO variant allele.

CONCLUSION

The RBC-FC analysis described here allows the detection of AB antigens weakly expressed in RBCs while achieving minimum background noise levels in negative control samples. Overall, the modified protocol provides a quick and reliable assay valuable in transfusion medicine and is potentially applicable to the characterization of rare weak ABO variants.

HLA allele and haplotype frequencies in Kazakhstani Russians and their relationship with other populations

HLA class I and class II genotypes from 947 Kazakhstani individuals of Russian origin were analyzed for investigating their most likely origin. The results were compared with similar data from other Russians (East and West), and also Worldwide populations, using standard genetic distances, neighbor-joining dendrograms, correspondence and haplotype analysis. Of the five HLA loci analyzed (HLA-A, HLA-C, HLA-B, HLA-DRB1, and HLA-DQB1) genotyped, 216 HLA alleles were identified. The most frequent alleles were A*02:01 (26.5%), B*07:02 (11.1%), C*04:01 (13.5%) and C*06:02 (12.1%), DRB1*07:01 (13.8%) and DRB1*15:01 (12.2%), and DQB1*03:01 (19.7%). Significant linkage disequilibrium was noted between all HLA pairs. DRB1*15:01 ~ DQB1*06:02 (10.5%), B*07:02 ~ C*07:02 (10.0%), B*07:02 ~ DRB1*15:01 (6.3%), and A*01:01 ~ B*08:01 (4.5%) were the most frequent two-locus haplotypes identified. Subsequent analyses showed that Kazakhstani Russians were closely related to West Russia-residing populations (Northwest Slavic, Vologda, Chelyabinsk, Moscow), East Europeans (Belarus Brest, Ukraine, Poland) and Scandinavians (Swedish, Finns), but distinct from East Russia-residing populations (Tuvians, Siberians from Chukotka, Kamchatka, and Ulchi) and East Mediterraneans (Levantines, Turks, North Macedonians, Albanians), and East Asians (Koreans, Japanese, Taiwanese, Mongolians). These results are in accordance with historical data indicating that the Russians of central Asia originate mainly from European Russia during the migratory flow of 18th and 19th centuries.

Hypomorphic PKD1 Alleles Impact Disease Variability in Autosomal Dominant Polycystic Kidney Disease

Autosomal dominant polycystic kidney disease (ADPKD) manifesting earlier than expected on the basis of family history can identify clinically tolerant PKD1 alleles with reduced expression. Hypomorphic PKD1 alleles can cause mild kidney disease or liver cysts in the absence of clinically manifest kidney involvement. The presented data highlight pleiotropic ADPKD clinical presentations and varying severity of kidney disease from PKD1 allele combinations.

Association of the risk factor UNC13A with survival and upper motor neuron involvement in amyotrophic lateral sclerosis

Background

The UNC13A gene is an established susceptibility locus for amyotrophic lateral sclerosis (ALS) and a determinant of shorter survival after disease onset, with up to 33.0 months difference in life expectancy for carriers of the rs12608932 risk genotype. However, its overall effect on other clinical features and ALS phenotypic variability is controversial.

Methods

Genotype data of the UNC13A rs12608932 SNP (A-major allele; C-minor allele) was obtained from a cohort of 972 ALS patients. Demographic and clinical variables were collected, including cognitive and behavioral profiles, evaluated through the Edinburgh Cognitive and Behavioral ALS Screen (ECAS) - Italian version and the Frontal Behavioral Inventory (FBI); upper and lower motor neuron involvement, assessed by the Penn Upper Motor Neuron Score (PUMNS) and the Lower Motor Neuron Score (LMNS)/Medical Research Council (MRC) scores, respectively; the ALS Functional Rating Scale Revised (ALSFRS-R) score at evaluation and progression rate; age and site of onset; survival. The comparison between the three rs12608932 genotypes (AA, AC, and CC) was performed using the additive, dominant, and recessive genetic models.

Results

The rs12608932 minor allele frequency was 0.31 in our ALS cohort, in comparison to 0.33-0.41 reported in other Caucasian ALS populations. Carriers of at least one minor C allele (AC + CC genotypes) had a shorter median survival than patients with the wild-type AA genotype (-11.7 months, p = 0.013), even after adjusting for age and site of onset, C9orf72 mutational status and gender. Patients harboring at least one major A allele (AA + AC genotypes) and particularly those with the wild-type AA genotype showed a significantly higher PUMNS compared to CC carriers (p = 0.015 and p_adj = 0.037, respectively), thus indicating a more severe upper motor neuron involvement. Our analysis did not detect significant associations with all the other clinical parameters considered.

Conclusion

Overall, our findings confirm the role of UNC13A as a determinant of survival in ALS patients and show the association of this locus also with upper motor neuron involvement.

Diversity and Spread of Acetolactate Synthase Allelic Variants at Position 574 Endowing Resistance in Amaranthus hybridus in Italy

Poor control of Amaranthus spp. with herbicides inhibiting acetolactate synthase (ALS) has been observed for several years in soybean fields in north-eastern Italy, but to date only a few ALS-resistant populations have been confirmed. An extensive sampling of putatively resistant Amaranthus accessions was completed in the Friuli Venezia Giulia region, across an arable land area of about 3000 km². In total, 58 accessions were tested to confirm their resistance status, recognize the Amaranthus species, identify the mutant ALS alleles endowing the resistance and determine the efficacy of 3 pre-emergence herbicides. Most accessions resulted in cross-resistance to thifensulfuron-methyl and imazamox. Genomic DNA were extracted from single seeds with a newly developed protocol; an allele-specific PCR assay revealed the presence of the 574-leucine in 20 accessions, of the 574-methionine in 22, and of both alleles in 9 accessions. The two variants showed a different spatial distribution. All resistant populations were ascribed to A. hybridus. A. hybridus resistant to ALS herbicides is well-established in this Italian region and its resistance is due to two ALS mutant alleles. Metribuzin, clomazone and metobromuron can be used as alternative herbicides to be applied in pre-emergence and they should be integrated into the management strategies to limit the spread of resistance.

The Effect of Interleukin-2 -330 T/G Polymorphism in People with Schizophrenia Among The Batak and Javanese Ethnic Groups

Background

Schizophrenia is a complex disorder involving multiple genes with mild to moderate effects and non-genetic risk factors such as environmental and psychological influences that alter brain chemistry. Significant reduction in interleukin-2 production by peripheral lymphocytes is an immunological finding replicated in schizophrenia across various countries. Investigations on the interleukin-2 -330 T/G polymorphism in people with schizophrenia (PWS) are still minimal, and the study location involves only a few countries with different results. Therefore, this study aims to examine the effect of interleukin-2 -330 T/G polymorphisms in people with schizophrenia among the Batak and Javanese ethnic groups in Indonesia, particularly in North Sumatra.

Method

This study used purposive non-probability sampling to recruit people with schizophrenia with 120 Batak and 120 Javanese subjects who were hospitalized at Prof. M. Ildrem Mental Hospital, Medan, Indonesia. The interleukin-2 -330 T/G polymorphism was examined by the PCR method.

Result

The results showed that the genotype frequency of the Batak people with schizophrenia is as follows: GG 11.7%, TG 53.3% and TT 35%. Furthermore, the group of Javanese people with schizophrenia is as follows: GG 23.3%, TG 44.2% and TT 32.5%. The OR for the genotypic comparison of GG was found to be 2.154 with 95% CI 0.992-4.678, p=0.053, while that of the TG genotype was 0.892 with 95% CI 0.505-1.574 and p=0.693. The T allele was higher than the G allele in Batak and Javanese ethnic groups, as demonstrated by p=0.713, OR=0.919 with 95% CI 0.641-1.318.

Conclusion

There is no statistically significant difference between the occurrence frequency of alleles. In addition, there is no significant relationship between the GG and the TG genotype of the interleukin-2 - 330T/G polymorphism in people with schizophrenia among the Batak and the Javanese ethnic groups.

Designing of cytotoxic T lymphocyte-based multi-epitope vaccine against SARS-CoV2: a reverse vaccinology approach

SARS-CoV2 is a single-stranded RNA virus, gaining much attention after it out broke in China in December 2019. The virus rapidly spread to several countries around the world and caused severe respiratory illness to humans. Since the outbreak, researchers around the world have devoted maximum resources and effort to develop a potent vaccine that would offer protection to uninfected individuals against SARS-CoV2. Reverse vaccinology is a relatively new approach that thrives faster in vaccine research. In this study, we constructed Cytotoxic T Lymphocytes (CTL)-based multi-epitope vaccine using hybrid epitope prediction methods. A total of 121 immunogenic CTL epitopes were screened by various sequence-based prediction methods and docked with their respective HLA alleles using the AutoDock Vina v1.1.2. In all, 17 epitopes were selected based on their binding affinity, followed by the construction of multi-epitope vaccine by placing the appropriate linkers between the epitopes and tuberculosis heparin-binding hemagglutinin (HBHA) adjuvant. The final vaccine construct was modeled by the I-TASSER server and the best model was further validated by ERRAT, ProSA, and PROCHECK servers. Furthermore, the molecular interaction of the constructed vaccine with TLR4 was assessed by ClusPro 2.0 and PROtein binDIng enerGY prediction (PRODIGY) server. The immune simulation analysis confirms that the constructed vaccine was capable of inducing long-lasting memory T helper (Th) and CTL responses. Finally, the nucleotide sequence was codon-optimized by the JCAT tool and cloned into the pET21a (+) vector. The current results reveal that the candidate vaccine is capable of provoking robust CTL response against the SARS-CoV2.Communicated by Ramaswamy H. Sarma.

Identification of MYB Transcription Factors Involving in Fruit Quality Regulation of Fragaria × ananassa Duch

The cultivated strawberry (Fragaria × ananassa Duch.) is an important horticultural crop. The economic values of strawberry cultivars are decided by their fruit qualities including taste, color and aroma. The important role of MYB transcription factors in fruit quality regulation is recognized increasingly with the identification of MYB genes involved in metabolism. A total of 407 MYB genes of F. × ananassa (FaMYBs) were identified in the genome-wide scale and named according to subgenome locations. The 407 FaMYBs were clustered into 36 groups based on phylogenetic analysis. According to synteny analysis, whole genome duplication and segmental duplication contributed over 90% of the expansion of the FaMYBs family. A total of 101 FaMYB loci with 1-6 alleles were identified by the homologous gene groups on homologous chromosomes. The differentially expressed FaMYB profiles of three cultivars with different fruit quality and fruit ripe processes provided the 8 candidate loci involved in fruit quality regulation. In this experiment, 7, 5, and 4 FaMYBs were screeded as candidate genes involved in the regulation of metabolism/transportation of anthocyanins, sugars or organic acids and 4-hydroxy-2, 5-dimethyl-3(2H)-furanone, respectively. These results pointed out the key FaMYBs for further functional analysis of gene regulation of strawberry fruit quality and would be helpful in the clarification on ofe roles of MYBs in the metabolism of fruit crops.

Current status of 4-aminoquinoline resistance markers 18 years after cessation of chloroquine use for the treatment of uncomplicated falciparum malaria in the littoral coastline region of Cameroon

The onset and rapid spread of chloroquine resistance and the introduction of amodiaquine for the treatment of uncomplicated malaria in Cameroon have influenced the proportion of Plasmodium falciparum sensitive and resistant alleles related to 4-aminoquinoline drugs. This study was undertaken to determine the prevalence of resistance markers to antimalarial 4-aminoquinolines in Douala in the Littoral Region, and Buea in the South West Region in June 2020. Dry blood spots were prepared from malaria microscopy positive cases and used for parasite DNA extraction by chelex-100 method. Plasmodium species identification was carried out by PCR amplification/agarose gel electrophoresis of 18srRNA. The Pfcrt and Pfmdr1 genes were amplified by PCR followed by restriction digestion. The prevalence of single nucleotide polymorphisms (SNPs) was compared between study sites and with previous studies carried out between 2003-2005 and 2009-2011 using the Chi square test. The results showed that Plasmodium falciparum was the dominant species occurring as mono-infections (84.6%). The wild type K76 allele of the Pfcrt gene was found in 74.9% of isolates while the wild N86, Y184 and D1246 alleles of the Pfmdr1 gene were found respectively in 87.2%, 89.6% and 100% of field isolates. The results showed a significant reduction in the mutant alleles compared to results obtained in 2003-2005 and 2009-2013. The KNYD haplotype was observed to be the most prevalent. The results indicated that there is a gradual erosion of the mutant Pfcrt and Pfmdr1 genotype and a gradual return to the sensitive P. falciparum genotype in Cameroon.

Diversity of NKG2C genotypes in a European population: Conserved and recombinant haplotypes in the coding, promoter, and 3'-untranslated regions

NK cells monitor altered molecular patterns in tumors and infected cells through an ample array of receptors. Two families of evolutionarily distant receptors have converged to enable human NK cells to sense levels of HLA class I ligands, frequently abnormal in altered cells. Whilst different forms of polymorphism are a hallmark of killer-cell immunoglobulin-like receptors and their classic HLA-A, B, and C ligands, genetic diversity of killer-cell lectin-like receptors for the non-classical HLA-E (CD94/NKG2 heterodimers) is less conspicuous and has attracted less attention. A common pattern of diversification in both receptor families is evolution of pairs of inhibitory and activating homologs for a common ligand, the genes encoding activating receptors being more frequently affected by copy number variation (CNV). This is exemplified by the gene encoding the activating NKG2C subunit (KLRC2 or NKG2C), which marks an NK-cell subpopulation that differentiates or expands in response to cytomegalovirus. We have studied NKG2C diversity in 240 South European individuals, using polymerase chain reaction and sequencing methods to assess both gene CNV and single-nucleotide polymorphisms (SNPs) affecting its promoter, coding and 3'-untranslated (3'UT) regions. Sequence analysis revealed eight common SNPs-one in the promoter, two in the coding sequence, and five in the 3'UT region. These SNPs associate strongly with each other, forming three conserved extended haplotypes (frequencies: 0.456, 0.221, and 0.117). Homo- and heterozygous combination of these, together with complete gene deletion (0.175) and additional haplotypes with frequencies lower than 0.015, generate a diversity of NKG2C genotypes of potential immunological importance.

Association of Human Leucocyte Antigen Polymorphism with Coronavirus Disease 19 in Renal Transplant Recipients

Human leucocyte antigens (HLAs) are highly polymorphic glycoproteins expressed at the surface of all nucleated cells. It is required for the SARS-CoV-2 peptide antigen presentation to immune cells for their effector response. However, polymorphism in HLA significantly impacts the binding of SARS-CoV-2 antigenic peptide to the HLA pocket and regulates immune activation. In this study, 514 renal transplant recipients (RTRs) were recruited from the outpatient department and categorized either into symptomatic (n = 173) or asymptomatic groups (n = 341) based on Coronavirus disease-19 (COVID-19) symptoms. The anti-SARS-CoV-2 spike protein-specific IgG antibody titer was measured by chemiluminescent microparticle immune-assay methods in 310 RTRs. The HLA details of 514 patients were retrieved from the electronic medical records and analyzed retrospectively. We found that HLA antigen allele A*24 was significantly associated with asymptomatic infection in 22.78%, HLA C*02 in 4.51%, DRB1*12 in 10.85%, and HLA DQA1*02 in 27.74% of RTRs. Whereas HLA A*29 in 3.46%, A*33 in 26.01%, B*13 in 10.40%, DRB1*10 in 4.62%, DRB1*15 in 39.30%, DRB1*30 in 1.15%, and DQA1*60 in 3.57% of RTRs were associated with symptomatic infection. HLA DRB1*13 and DRB1*15 were associated with moderate to severe degrees of COVID-19 disease. The seroconversion rate in asymptomatic patients was 118/137 (86.13%), had a median titer of 647.80 au/ml, compared to symptomatic patients 148/173 (85.54%) with a median titer of 400.00 au/ml, which was not significant between the two groups (P = 0.88 and 0.13). In conclusion, HLA alleles A*24, C*02, DRB1*12, and DQA1*02 were significantly associated with asymptomatic infection, and A*29, A*33, B*13, DRB1*10, DRB*15, and DRB1*30 were significantly associated with symptomatic infection. HLA DRB1*13 and DRB1*15 were associated with moderate to severe degrees of COVID-19 disease.

A first complete phylogenomic hypothesis for diploid blueberries (Vaccinium section Cyanococcus)

PREMISE

The true blueberries (Vaccinium sect. Cyanococcus; Ericaceae), endemic to North America, have been intensively studied for over a century. However, with species estimates ranging from nine to 24 and much confusion regarding species boundaries, this ecologically and economically valuable group remains inadequately understood at a basic evolutionary and taxonomic level. As a first step toward understanding the evolutionary history and taxonomy of this species complex, we present the first phylogenomic hypothesis of the known diploid blueberries.

METHODS

We used flow cytometry to verify the ploidy of putative diploid taxa and a target-enrichment approach to obtain a genomic data set for phylogenetic analyses.

RESULTS

Despite evidence of gene flow, we found that a primary phylogenetic signal is present. Monophyly for all morphospecies was recovered, with two notable exceptions: one sample of V. boreale was consistently nested in the V. myrtilloides clade and V. caesariense was nested in the V. fuscatum clade. One diploid taxon, Vaccinium pallidum, is implicated as having a homoploid hybrid origin.

CONCLUSIONS

This foundational study represents the first attempt to elucidate evolutionary relationships of the true blueberries of North America with a phylogenomic approach and sets the stage for multiple avenues of future study such as a taxonomic revision of the group, the verification of a homoploid hybrid taxon, and the study of polyploid lineages within the context of a diploid phylogeny.

Immunodominant MHC-II (Major Histocompatibility Complex II) Restricted Epitopes in Human Apolipoprotein B

BACKGROUND

CD (cluster of differentiation) 4+ T-cell responses to APOB (apolipoprotein B) are well characterized in atherosclerotic mice and detectable in humans. CD4+ T cells recognize antigenic peptides displayed on highly polymorphic HLA (human leukocyte antigen)-II. Immunogenicity of individual APOB peptides is largely unknown in humans. Only 1 HLA-II-restricted epitope was validated using the DRB1*07:01-APOB3036-3050 tetramer. We hypothesized that human APOB may contain discrete immunodominant CD4+ T-cell epitopes that trigger atherosclerosis-related autoimmune responses in donors with diverse HLA alleles.

METHODS

We selected 20 APOB-derived peptides (APOB20) from an in silico screen and experimentally validated binding to the most commonly occurring human HLA-II alleles. We optimized a restimulation-based workflow to evaluate antigenicity of multiple candidate peptides in HLA-typed donors. This included activation-induced marker assay, intracellular cytokine staining, IFNγ (interferon gamma) enzyme-linked immunospot and cytometric bead array. High-throughput sequencing revealed TCR (T-cell receptor) clonalities of APOB-reactive CD4+ T cells.

RESULTS

Using stringent positive, negative, and crossover stimulation controls, we confirmed specificity of expansion-based protocols to detect CD4+ T cytokine responses to the APOB20 pool. Ex vivo assessment of AIM+CD4+ T cells revealed a statistically significant autoimmune response to APOB20 but not to a ubiquitously expressed negative control protein, actin. Resolution of CD4+ T responses to the level of individual peptides using IFNγ enzyme-linked immunospot led to the discovery of 6 immunodominant epitopes (APOB6) that triggered robust CD4+ T activation in most donors. APOB6-specific responding CD4+ T cells were enriched in unique expanded TCR clonotypes and preferentially expressed memory markers. Cytometric bead array analysis detected APOB6-induced secretion of both proinflammatory and regulatory cytokines. In clinical samples from patients with angiographically verified coronary artery disease, APOB6 stimulation induced higher activation and memory phenotypes and augmented secretion of proinflammatory cytokines TNF (tumor necrosis factor) and IFNγ, compared with patients with low coronary artery disease.

CONCLUSIONS

Using 3 cohorts, each with ≈20 donors, we discovered and validated 6 immunodominant, HLA-II-restricted APOB epitopes. The immune response to these APOB epitopes correlated with coronary artery disease severity.

ENPP2 (Endothelial Ectonucleotide Pyrophosphatase/Phosphodiesterase 2) Increases Atherosclerosis in Female and Male Mice

BACKGROUND

Maladapted endothelial cells (ECs) secrete ENPP2 (ectonucleotide pyrophosphatase/phosphodiesterase 2; autotaxin)-a lysophospholipase D that generates lysophosphatidic acids (LPAs). ENPP2 derived from the arterial wall promotes atherogenic monocyte adhesion induced by generating LPAs, such as arachidonoyl-LPA (LPA20:4), from oxidized lipoproteins. Here, we aimed to determine the role of endothelial ENPP2 in the production of LPAs and atherosclerosis.

METHODS

We quantified atherosclerosis in mice harboring loxP-flanked Enpp2 alleles crossed with Apoe^-/- mice expressing tamoxifen-inducible Cre recombinase under the control of the EC-specific bone marrow X kinase promoter after 12 weeks of high-fat diet feeding.

RESULTS

A tamoxifen-induced EC-specific Enpp2 knockout decreased atherosclerosis, accumulation of lesional macrophages, monocyte adhesion, and expression of endothelial CXCL (C-X-C motif chemokine ligand) 1 in male and female Apoe^-/- mice. In vitro, ENPP2 mediated the mildly oxidized LDL (low-density lipoprotein)-induced expression of CXCL1 in aortic ECs by generating LPA20:4, palmitoyl-LPA (LPA16:0), and oleoyl-LPA (LPA18:1). ENPP2 and its activity were detected on the endothelial surface by confocal imaging. The expression of endothelial Enpp2 established a strong correlation between the plasma levels of LPA16:0, stearoyl-LPA (LPA18:0), and LPA18:1 and plaque size and a strong negative correlation between the LPA levels and ENPP2 activity in the plasma. Moreover, endothelial Enpp2 knockout increased the weight of high-fat diet-fed male Apoe^-/- mice.

CONCLUSIONS

We demonstrated that the expression of ENPP2 in ECs promotes atherosclerosis and endothelial inflammation in a sex-independent manner. This might be due to the generation of LPA20:4, LPA16:0, and LPA18:1 from mildly oxidized lipoproteins on the endothelial surface.

Single-Nucleotide Polymorphism on Spermatogenesis Associated 16 Gene-Coding Region Affecting Bovine Leukemia Virus Proviral Load

Bovine leukemia virus (BLV) is an etiological agent of malignant lymphoma in cattle and is endemic in many cattle-breeding countries. Thus, the development of cattle genetically resistant to BLV is desirable. The purpose of this study was to identify novel single-nucleotide polymorphisms (SNPs) related to resistance to BLV. A total of 146 DNA samples from cattle with high BLV proviral loads (PVLs) and 142 samples from cattle with low PVLs were used for a genome-wide association study (GWAS). For the verification of the GWAS results, an additional 1342 and 456 DNA samples from BLV-infected Japanese Black and Holstein cattle, respectively, were used for an SNP genotyping PCR to compare the genotypes for the identified SNPs and PVLs. An SNP located on the spermatogenesis associated 16 (SPATA16)-coding region on bovine chromosome 1 was found to exceed the moderate threshold (p < 1.0 × 10−5) in the Additive and Dominant models of the GWAS. The SNP genotyping PCR revealed that the median values of the PVL were 1278 copies/50 ng of genomic DNA for the major homozygous, 843 for the heterozygous, and 621 for the minor homozygous genotypes in the Japanese Black cattle (p < 0.0001). A similar tendency was also observed in the Holstein cattle. We found that cattle with the minor allele for this SNP showed 20−25% lower PVLs. Although the mechanisms through which this SNP impacts the PVL remain unknown, we found a novel SNP related to BLV resistance located on the SPATA16 gene-coding region on bovine chromosome 1.

Whole genome analysis of rotavirus strains circulating in Benin before vaccine introduction, 2016-2018

Species A Rotaviruses (RVA) still play a major role in causing acute diarrhea in children under five years old worldwide. Currently, an 11-gene classification system is used to designate the full genotypic constellations of circulating strains. Viral proteins and non-structural proteins in the order VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6 are represented by the genotypes Gx-P[x]-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx, respectively. In Benin, ROTAVAC® vaccine was introduced into the Expanded Programme on Immunization in December 2019. To monitor circulating RVA strains for changes that may affect vaccine performance, in-depth analysis of strains prior to vaccine introduction are needed. Here we report, the whole-gene characterization (11 ORFs) for 72 randomly selected RVA strains of common and unusual genotypes collected in Benin from the 2016-2018 seasons. The sequenced strains were 15 G1P[8], 20 G2P[4], 5 G9P[8], 14 G12P[8], 9 G3P[6], 2 G1P[6], 3 G2P[6], 2 G9P[4], 1 G12P[6], and 1 G1G9P[8]/P[4]. The study strains exhibited two genetic constellations designed as Wa-like G1/G9/G12-P[6]/P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1 and DS-1-like G2/G3/G12-P[4]/P[6]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Genotype G9P[4] strains possessed a DS-1-like genetic constellation with an E6 NSP4 gene, G9-P[4]-I2-R2-C2-M2-A2-N2-T2-E6-H2. The mixed genotype showed both Wa-like and DS-1-like profiles with a T6 NSP3 gene G1/G9P[8]/[4]-I1/I2-R1/R2-C1/C2-M1/M2-A1/A2-N1/N2-T1/T6-E1/E6-H1/H2. At the allelic level, the analysis of the Benin strains, reference strains (with known alleles), vaccine strains (with known alleles) identified 2-13 and 1-17 alleles for DS-1-like and Wa-like strains, respectively. Most of the study strains clustered into previously defined alleles, but we defined 3 new alleles for the VP7 (G3=1 new allele and G12=2 new alleles) and VP4 (P[4]=1 new allele and P[6]=2 new alleles) genes which formed the basis of the VP7 and VP4 gene clusters, respectively. For the remaining 9 genes, 0-6 new alleles were identified for both Wa-like and DS-1-like strains. This analysis of whole genome sequences of RVA strains circulating in Benin described genetic point mutations and reassortment events as well as novel alleles. Further detailed studies on these new alleles are needed and these data can also provide a baseline for studies on RVA in the post-vaccination period.

SNP-based Phylogenomic Inference in Holarctic Ground Squirrels (Urocitellus)

Resolution of rapid evolutionary radiations requires harvesting maximal signal from phylogenomic datasets. However, studies of non-model clades often target conserved loci that are characterized by reduced information content, which can negatively affect gene tree precision and species tree accuracy. Single nucleotide polymorphism (SNP)-based methods are an underutilized but potentially valuable tool for estimating phylogeny and divergence times because they do not rely on resolved gene trees, allowing information from many or all variant loci to be leveraged in species tree reconstruction. We evaluated the utility of SNP-based methods in resolving phylogeny of Holarctic ground squirrels (Urocitellus), a radiation that has been difficult to disentangle, even in prior phylogenomic studies. We inferred phylogeny from a dataset of >3,000 ultraconserved element loci (UCEs) using two methods (SNAPP, SVDquartets) and compared our results with a new mitogenome phylogeny. We also systematically evaluated how phasing of UCEs improves per-locus information content, and inference of topology and other parameters within each of these SNP-based methods. Phasing improved topological resolution and branch length estimation at shallow levels (within species complexes), but less so at deeper levels, likely reflecting true uncertainty due to ancestral polymorphisms segregating in these rapidly diverging lineages. We resolved several key clades in Urocitellus and present targeted opportunities for future phylogenomic inquiry. Our results extend the roadmap for use of SNPs to address vertebrate radiations and support comparative analyses at multiple temporal scales.

BACH1 Binding Links the Genetic Risk for Severe Periodontitis with ST8SIA1

Genome-wide association studies identified various loci associated with periodontal diseases, but assigning causal alleles remains difficult. Likewise, the generation of biological meaning underlying a statistical association has been challenging. Here, we characterized the genetic association at the gene ST8SIA1 that increases the risk for severe periodontitis in smokers. We used CRISPR/dCas9 activation and RNA-sequencing to identify genetic interaction partners of ST8SIA1 and to determine its function in the cell. We used reporter gene assays to identify regulatory elements at the associated single-nucleotide polymorphisms (SNPs) and to determine effect directions and allele-specific changes of enhancer activity. Antibody electrophoretic mobility shift assays proved allele-specific transcription factor binding at the putative causal SNPs. We found the reported periodontitis risk gene ABCA1 as the top upregulated gene following ST8SIA1 activation. Gene set enrichment analysis showed highest effects on integrin cell surface interactions (area under the curve [AUC] = 0.85; q = 4.9 × 10-6) and cell cycle regulation (AUC = 0.89; q = 1.6 × 10-5). We identified 2 associated repressor elements in the introns of ST8SIA1 that bind the transcriptional repressor BACH1. The putative causative variant rs2012722 decreased BACH1 binding by 40%. We also pinpointed ST8SIA1 as the target gene of the association. ST8SIA1 inhibits cell adhesion with extracellular matrix proteins, integrins, and cell cycle, as well as enhances apoptosis. Likewise, tobacco smoke reportedly results in an inhibition of cell adhesion and a decrease in integrin-positive cells and cell growth. We conclude that impaired ST8SIA1 repression, independently caused by reduced BACH1 binding at the effect T allele, as well as by tobacco smoke, contributes to higher ST8SIA1 levels, and in smokers who carry the effect T allele, both factors would be additive with damaging effects on the gingival barrier integrity. The activity of ST8SIA1 is also linked with the periodontitis risk gene ABCA1.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

The novel HLA-DQB1*05:240 allele was likely generated by recombination between DQB1*05:01 and DQB1*03:02

Identification of the novel HLA-DQB1*05:240 allele that differs from HLA-DQB1*05:01:09 at three positions in exon 2.

Chronic wasting disease in Norway-A survey of prion protein gene variation among cervids

Susceptibility of cervids to Chronic Wasting Disease (CWD), a prion disease, can be modulated by variations in the prion protein gene (PRNP), encoding the cellular prion protein (PrP^C ). In prion diseases, PrP^C is conformationally converted to pathogenic conformers (PrP^Sc ), aggregates of which comprise infectious prions. CWD has recently been observed in its contagious form in Norwegian reindeer (Rangifer tarandus) and in novel, potentially sporadic forms, here called 'atypical CWD', in moose (Alces alces) and red deer (Cervus elaphus). To estimate relative susceptibility of different Norwegian cervid species to CWD, their non-synonymous PRNP variants were analyzed. In reindeer, seven PRNP alleles were observed and in red deer and moose two alleles were present, whereas roe deer (Capreolus capreolus) PRNP was monomorphic. One 'archetypal' PRNP allele associated with susceptibility was common to all four cervid species. The distribution of PRNP alleles differed between wild and semi-domesticated reindeer, with alleles associated with a high susceptibility occurring, on average, above 55% in wild reindeer and below 20% in semi-domesticated reindeer. This difference may reflect the diverse origins of the populations and/or selection processes during domestication and breeding. Overall, PRNP genetic data indicate considerable susceptibility to CWD among Norwegian cervids and suggest that PRNP homozygosity may be a risk factor for the atypical CWD observed in moose. The CWD isolates found in the Norwegian cervid species differ from those previously found in Canada and USA. Our study provides an overview of the PRNP genetics in populations exposed to these emerging strains that will provide a basis for understanding these strains' dynamics in relation to PRNP variability.

A simple genotyping method for CD247 3'-untranslated region polymorphism rs1052231 and characterization of a reference cell panel

CD247 (or CD3‐ζ chain) is an essential adaptor and signal‐transducing molecule of the T‐cell antigen receptor (TCR) complex, and it also couples to NK‐cell activating receptors such as NKp46, NKp30 and CD16A (FcγRIII). Noncoding sequence polymorphisms and variations in CD247 expression, a tightly regulated process, have been related with an altered immune response in multiple health conditions. A single nucleotide polymorphism (T > A) at nucleotide 844 of the CD247 3′‐untranslated region, rs1052231, has been related with lower CD247 gene expression and it has been investigated as a potential biomarker of autoimmune disease. We present here a simple, accurate, reliable, time‐efficient, and cost‐effective method for CD247‐rs1052231 genotyping. Using this method, based on polymerase chain reaction with confronting two‐pair primers (PCR‐CTPP), we have also characterized the CD247‐rs1052231 genotypes in a panel of worldwide available cell lines, which should facilitate study of the role of this polymorphism in immunity and human health.

HybPhaser: A workflow for the detection and phasing of hybrids in target capture data sets

PREMISE

Hybrids contain divergent alleles that can confound phylogenetic analyses but can provide insights into reticulated evolution when identified and phased. We developed a workflow to detect hybrids in target capture data sets and phase reads into parental lineages using a similarity and phylogenetic framework.

METHODS

We used Angiosperms353 target capture data for Nepenthes, including known hybrids to test the novel workflow. Reference mapping was used to assess heterozygous sites across the data set and to detect hybrid accessions and paralogous genes. Hybrid samples were phased by mapping reads to multiple references and sorting reads according to similarity. Phased accessions were included in the phylogenetic framework.

RESULTS

All known Nepenthes hybrids and nine additional samples had high levels of heterozygous sites, had reads associated with multiple divergent clades, and were phased into accessions resembling divergent haplotypes. Phylogenetic analysis including phased accessions increased clade support and confirmed parental lineages of hybrids.

DISCUSSION

HybPhaser provides a novel approach to detect and phase hybrids in target capture data sets, which can provide insights into reticulations by revealing origins of hybrids and reduce conflicting signal, leading to more robust phylogenetic analyses.

Kidney-Specific Klotho Gene Deletion Causes Aortic Aneurysm via Hyperphosphatemia

[Figure: see text].