Analysis of pooled genome sequences from Djallonke and Sahelian Sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment

Background: The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually. Results: Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis. Conclusions: Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases

Variants of HCMV UL18 sequenced directly from clinical specimens associate with antibody and T-cell responses to HCMV

Around 80% of adults worldwide carry human cytomegaloviris (HCMV). The HCMV gene UL18 is a homolog of HLA class I genes and encodes a protein with high affinity for the NK and T-cell cytotoxicity inhibitor LIR-1. UL18 was deep sequenced from blood, saliva or urine from Indonesian people with HIV (PWH) (n = 28), Australian renal transplant recipients (RTR) (n = 21), healthy adults (n = 7) and neonates (n = 4). 95% of samples contained more than one variant of HCMV UL18, as defined by carriage of nonsynonymous variations. When aligned with immunological markers of the host’s burden of HCMV, the S318N variation associated with high levels of antibody reactive with HCMV lysate in PWH over 12 months on antiretroviral therapy. The A107T variation associated with HCMV antibody levels and inflammatory biomarkers in PWH at early timepoints. Variants D32G, D248N, V250A and E252D aligned with elevated HCMV antibody levels in RTR, while M191K, E196Q and F165L were associated with HCMV-reactive T-cells and proportions of Vδ2− γδ T-cells—populations linked with high burdens of HCMV. We conclude that UL18 is a highly variable gene, where variation may alter the persistent burden of HCMV and/or the host response to that burden

Sequencing of the viral UL111a gene directly from clinical specimens reveals variants of HCMV-encoded IL-10 that are associated with altered immune responses to HCMV

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ~80% of adults worldwide. Acute infections are often asymptomatic in healthy individuals but generate diverse syndromes in neonates, renal transplant recipients (RTR), and people with HIV (PWH). The HCMV gene UL111a encodes a homolog of human interleukin-10 (IL-10) that interacts with the human IL-10 receptor. Deep sequencing technologies were used to sequence UL111a directly from 59 clinical samples from Indonesian PWH and Australian RTR, healthy adults, and neonates. Overall, 93% of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous variation. Carriage of these variants differed between neonates and adults, Australians and Indonesians, and between saliva and blood leukocytes. The variant alleles of N41D and S71Y occurred together in Australian RTR and were associated with higher T-cell responses to HCMV pp65. The variant P122S was associated with lower levels of antibodies reactive with a lysate of HCMV-infected fibroblasts. L174F was associated with increased levels of antibodies reactive with HCMV lysate, immediate-early 1 (IE-1), and glycoprotein B (gB) in Australian RTR and Indonesians PWH, suggesting a higher viral burden. We conclude that variants of UL111a are common in all populations and may influence systemic responses to HCMV

The mouse as a model for the effects of MHC genes on human disease

As mice are often used to model human major histocompatibility complex (MHC)-associated diseases, it is important to understand how their MHC regions differ at the DNA level. The sequencing of the mouse MHC (H2 region) has enabled a detailed map of this region to be assembled for comparison with the human MHC. Here, Richard Allcock and colleagues outline the similarities between the human and mouse MHC regions and discuss notable differences that might affect disease models

Investigation of NOTCH4 coding region polymorphisms in sporadic inclusion body myositis

The NOTCH4 gene, located within the MHC region, is involved in cellular differentiation and has varying effects dependent on tissue type. Coding region polymorphisms haplotypic of the sIBM-associated 8.1 ancestral haplotype were identified in NOTCH4 and genotyped in two different Caucasian sIBM cohorts. In both cohorts the frequency of the minor allele of rs422951 and the 12-repeat variation for rs72555375 was increased and was higher than the frequency of the sIBM-associated allele HLA-DRB1*0301. These NOTCH4 polymorphisms can be considered to be markers for sIBM susceptibility, but require further investigation to determine whether they are directly involved in the disease pathogenesis

Sporadic inclusion body myositis in Japanese is associated with the MHC ancestral haplotype 52.1

In Caucasians, sporadic inclusion body myositis has been associated with the MHC ancestral haplotypes; HLA-A1, B8, DR3 (8.1AH) and HLA-B35, DR1 (35.2AH). It is not known whether these haplotypes carry susceptibility for the disease in other ethnic groups. We report here the results of HLA-B and -DRB1 typing using a high-resolution sequence-based technique in a cohort of 31 Japanese patients with definite sIBM. Patient allele frequencies were 40.3% for HLA-B*5201 (10.7% in controls: p<0.001) and 37.1% for HLA-DRB1*1502 (10% in controls: p<0.001). Both alleles were found together as part of a conserved haplotype (52.1AH) at a frequency of 37.1% in patients (8.4% in controls: p<0.001). This is the first description of a haplotypic MHC association with sporadic inclusion body myositis in Japanese patients. These findings indicate that different MHC ancestral haplotypes are associated with sIBM in different ethnic groups and further emphasize the importance of genetic factors in this condition

Recombination mapping of the susceptibility region for sporadic inclusion body myositis within the major histocompatibility complex

Susceptibility to sporadic inclusion body myositis (sIBM) in Caucasians has been consistently associated with alleles of the major histocompatibility complex (MHC) 8.1 ancestral haplotype (AH) (defined by HLA-B*0801 and HLA-DRB1*0301). In this study recombination mapping was utilised to further refine the known 8.1AH susceptibility region near HLA-DRB1*0301. Caucasian sIBM patients carrying part of the 8.1AH were genotyped for a selection of 8.1AH-haplotypic polymorphisms. A common 8.1AH-specific susceptibility region was defined, spanning 172 kb and encompassing three genes — HLA-DRB3, HLA-DRA and BTNL2. It is thus likely that 8.1AH-derived susceptibility to sIBM originates from at least one of these genes

Microsatellite markers from the Ion Torrent: a multi-species contrast to 454 shotgun sequencing

The development and screening of microsatellite markers have been accelerated by next-generation sequencing (NGS) technology and in particular GS-FLX pyro-sequencing (454). More recent platforms such as the PGM semiconductor sequencer (Ion Torrent) offer potential benefits such as dramatic reductions in cost, but to date have not been well utilized. Here, we critically compare the advantages and disadvantages of microsatellite development using PGM semiconductor sequencing and GS-FLX pyro-sequencing for two gymnosperm (a conifer and a cycad) and one angiosperm species. We show that these NGS platforms differ in the quantity of returned sequence data, unique microsatellite data and primer design opportunities, mostly consistent with the differences in read length. The strength of the PGM lies in the large amount of data generated at a comparatively lower cost and time. The strength of GS-FLX lies in the return of longer average length sequences and therefore greater flexibility in producing markers with variable product length, due to longer flanking regions, which is ideal for capillary multiplexing. These differences need to be considered when choosing a NGS method for microsatellite discovery. However, the ongoing improvement in read lengths of the NGS platforms will reduce the disadvantage of the current short read lengths, particularly for the PGM platform, allowing greater flexibility in primer design coupled with the power of a larger number of sequences

Identification of a CD8+ T-cell response to a predicted neoantigen in malignant mesothelioma

Neoantigens present unique and specific targets for personalized cancer immunotherapy strategies. Given the low mutational burden yet immunotherapy responsiveness of malignant mesothelioma (MM) when compared to other carcinogen-induced malignancies, identifying candidate neoantigens and T cells that recognize them has been a challenge. We used pleural effusions to gain access to MM tumor cells as well as immune cells in order to characterize the tumor-immune interface in MM. We characterized the landscape of potential neoantigens from SNVs identified in 27 MM patients and performed whole transcriptome sequencing of cell populations from 18 patient-matched pleural effusions. IFNγ ELISpot was performed to detect a CD8+ T cell responses to predicted neoantigens in one patient. We detected a median of 68 (range 7–258) predicted neoantigens across the samples. Wild-type non-binding to mutant binding predicted neoantigens increased risk of death in a model adjusting for age, sex, smoking status, histology and treatment (HR: 33.22, CI: 2.55–433.02, p = .007). Gene expression analysis indicated a dynamic immune environment within the pleural effusions. TCR clonotypes increased with predicted neoantigen burden. A strong activated CD8+ T-cell response was identified for a predicted neoantigen produced by a spontaneous mutation in the ROBO3 gene. Despite the challenges associated with the identification of bonafide neoantigens, there is growing evidence that these molecular changes can provide an actionable target for personalized therapeutics in difficult to treat cancers. Our findings support the existence of candidate neoantigens in MM despite the low mutation burden of the tumor, and may present improved treatment opportunities for patients

Sequencing directly from clinical specimens reveals genetic variations in HCMV-Encoded Chemokine Receptor US28 that may influence antibody levels and interactions with human chemokines

Human cytomegalovirus (HCMV) is a beta-herpesvirus carried by ∼80% of the world’s population. Acute infections are asymptomatic in healthy individuals but generate diverse syndromes in neonates, solid organ transplant recipients, and HIV-infected individuals. The HCMV gene US28 encodes a homolog of a human chemokine receptor that is able to bind several chemokines and HIV gp120. Deep sequencing technologies were used to sequence US28 directly from 60 clinical samples from Indonesian HIV patients and Australian renal transplant recipients, healthy adults, and neonates. Molecular modeling approaches were used to predict whether nine nonsynonymous mutations in US28 may alter protein binding to a panel of six chemokines and two variants of HIV gp120. Ninety-two percent of samples contained more than one variant of HCMV, as defined by at least one nonsynonymous mutation. Carriage of these variants differed between neonates and adults, Australian and Indonesian samples, and saliva samples and blood leukocytes. Two nonsynonymous mutations (N170D and R267K) were associated with increased levels of immediate early protein 1 (IE-1) and glycoprotein B (gB) HCMV-reactive antibodies, suggesting a higher viral burden. Seven of the nine mutations were predicted to alter binding of at least one ligand. Overall, HCMV variants are common in all populations and have the potential to affect US28 interactions with human chemokines and/or gp120 and alter responses to the virus. The findings relied on deep sequencing technologies applied directly to clinical samples, so the variants exist in vivo