A Bead-based Normalization for Uniform Sequencing depth (BeNUS) protocol for multi-samples sequencing exemplified by HLA-B

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

The Complete Genome Sequence of Bacillus thuringiensis Strain ATCC 13367 Features a cry-Containing Chromosome

We report the high-quality closed genome sequence of Bacillus thuringiensis ATCC 13367, which is used as a bioinsecticide. The genome features a 328-kbp plasmid and a 5.8-Mbp chromosome. The genome is atypical for the species since it lacks a plasmid pesticidal crystal protein (cry) gene and instead features a chromosomal cry gene like that of B. thuringiensis HER1410.

Characterisation of atypical Shiga toxin gene sequences and description of Stx2j, a new subtype

Shiga toxin (Stx) is the definitive virulence factor of Shiga toxin-producing Escherichia coli (STEC). Stx variants are currently organised into a taxonomic system of three Stx1 (a,c,d) and seven Stx2 (a,b,c,d,e,f,g) subtypes. In this study, seven STEC isolates from food and clinical samples possessing stx2 sequences that do not fit current Shiga toxin taxonomy were identified. Genome assemblies of the STEC strains was created from Oxford Nanopore and Illumina sequence data. The presence of atypical stx2 sequences were confirmed by Sanger sequencing, as were Stx2 expression and cytotoxicity. A strain of O157:H7 was found to possess stx1a and a truncated stx2a, which were originally misidentified as an atypical stx2. Two strains possessed unreported variants of Stx2a (O8:H28) and Stx2b (O146:H21). In four of the strains we found three Stx-subtypes that are not included in the current taxonomy. Stx2h (O170:H18) was identified in a Canadian sprout isolate; this subtype has only previously been reported in STEC from Tibetan Marmots. Stx2o (O85:H1) was identified in a clinical isolate. Finally, Stx2j (O158:H23 and O33:H14) was found in lettuce and clinical isolates. The results of this study expands the number of known Stx subtypes, the range of STEC serotypes, and isolation sources in which they may be found. The presence of the Stx2j and Stx2o in clinical isolates of STEC indicates that strains carrying these variants are potential human pathogens. Highlights Atypical Shiga toxin (stx) genes in Escherichia coli were sequenced. Two new variants of stx2a and stx2b are described. Two strains carried subtypes Stx2h and Stx2o, which have only one previous report. Two strains carried a previously undescribed subtype, Stx2j.

Delayed lactose utilization among Shiga toxin-producing Escherichia coli of serogroup O121

Two outbreaks of Shiga toxin-producing Escherichia coli O121:H19 associated with wheat flour, in the United States of America and Canada, involved strains with an unusual phenotype, delayed lactose utilization (DLU). These strains do not ferment lactose when initially cultured on MacConkey agar (MAC), but lactose fermentation occurs following subculture to a second plate of MAC. The prevalence of DLU was determined by examining the β-galactosidase activity of 49 strains of E. coli O121, and of 37 other strains of E. coli. Twenty four of forty three O121:H19 and one O121:NM displayed DLU. Two strains (O121:NM and O145:H34) did not have detectable β-galactosidase activity. β-glucuronidase activity of O121 strains was also determined. All but six DLU strains had normal β-glucuronidase activity. β-glucuronidase activity was suppressed on MAC for 17 of 23 O121 non-DLU strains. Genomic analysis found that DLU strains possessed an insertion sequence, IS600 (1267 bp), between lacZ (β-galactosidase) and lacY (β-galactoside permease), that was not present in strains exhibiting normal lactose utilization. The insert might reduce the expression of β-galactoside permease, delaying import of lactose, resulting in the DLU phenotype. The high probability of DLU should be considered when using lactose-containing media for the isolation of STEC O121.

Genomic analysis of human noroviruses using combined Illumina-Nanopore data

Whole-genome sequence analysis of noroviruses is routinely performed by employing a metagenomic approach. While this methodology has several advantages, such as allowing for the examination of co-infection, it has some limitations, such as the requirement of high viral load to achieve full-length or near full-length genomic sequences. In this study, we used a pre-amplification step to obtain full-length genomic amplicons from 39 Canadian GII isolates, followed by deep sequencing on Illumina and Oxford Nanopore platforms. This approach significantly reduced the required viral titre to obtain full-genome coverage. Herein, we compared the coverage and sequences obtained by both platforms and provided an in-depth genomic analysis of the obtained sequences, including the presence of single-nucleotide variants and recombination events.

Complete Genome Sequences of Priestia megaterium Type and Clinical Strains Feature Complex Plasmid Arrays

Here, we report the high-quality complete genome sequences and plasmid arrays of Priestia megaterium ATCC 14581^T and of two clinical strains (2008724129 and 2008724142) isolated from human samples in the United States.

Complete Genome Sequence of Clostridium botulinum CJ0611A1, a Type A(B) Isolate Associated with an International Outbreak of Botulism from Commercial Carrot Juice

We report the complete genome (3.9-Mb chromosome, 5.9-kb plasmid) of Clostridium botulinum CJ0611A1, a type A(B) strain isolated from carrot juice distributed in Canada and linked to an international 2006 foodborne botulism outbreak. This strain encodes a full-length bont/A1 gene and a truncated bont/B gene.

We report the complete genome (3.9-Mb chromosome, 5.9-kb plasmid) of Clostridium botulinum CJ0611A1, a type A(B) strain isolated from carrot juice distributed in Canada and linked to an international 2006 foodborne botulism outbreak. This strain encodes a full-length bont/A1 gene and a truncated bont/B gene.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.

Epimutations in Developmental Genes Underlie the Onset of Domestication in Farmed European Sea Bass

Domestication of wild animals induces a set of phenotypic characteristics collectively known as the domestication syndrome. However, how this syndrome emerges is still not clear. Recently, the neural crest cell deficit hypothesis proposed that it is generated by a mildly disrupted neural crest cell developmental program, but clear support is lacking due to the difficulties of distinguishing pure domestication effects from preexisting genetic differences between farmed and wild mammals and birds. Here, we use a farmed fish as model to investigate the role of persistent changes in DNA methylation (epimutations) in the process of domestication. We show that early domesticates of sea bass, with no genetic differences with wild counterparts, contain epimutations in tissues with different embryonic origins. About one fifth of epimutations that persist into adulthood are established by the time of gastrulation and affect genes involved in developmental processes that are expressed in embryonic structures, including the neural crest. Some of these genes are differentially expressed in sea bass with lower jaw malformations, a key feature of domestication syndrome. Interestingly, these epimutations significantly overlap with cytosine-to-thymine polymorphisms after 25 years of selective breeding. Furthermore, epimutated genes coincide with genes under positive selection in other domesticates. We argue that the initial stages of domestication include dynamic alterations in DNA methylation of developmental genes that affect the neural crest. Our results indicate a role for epimutations during the beginning of domestication that could be fixed as genetic variants and suggest a conserved molecular process to explain Darwin’s domestication syndrome across vertebrates.

Next-generation sequencing identifies contribution of both class I and II HLA genes on susceptibility of multiple sclerosis in Japanese

Background

The spectrum of classical and non-classical HLA genes related to the risk of multiple sclerosis (MS) and neuromyelitis optica spectrum disorder (NMOSD) in the Japanese population has not been studied in detail. We conducted a case-control analysis of classical and non-classical HLA genes.

Methods

We used next-generation sequencing (NGS)-based HLA genotyping methods for mapping risk for 45 MS patients, 31 NMOSD patients, and 429 healthy controls. We evaluated the association of the HLA variants with the risk of MS and NMOSD using logistic regression analysis and Fisher’s exact test.

Results

We confirmed that HLA-DRB1*15:01 showed the strongest association with MS (P = 2.1 × 10⁻⁵; odds ratio [OR] = 3.44, 95% confidence interval [95% CI] = 1.95–6.07). Stepwise conditional analysis identified HLA-DRB1*04:05, HLA-B*39:01, and HLA-B*15:01 as being associated with independent MS susceptibility (P_Conditional < 8.3 × 10⁻⁴). With respect to amino acid polymorphisms in HLA genes, we found that phenylalanine at HLA-DQβ1 position 9 had the strongest effect on MS susceptibility (P = 3.7 × 10⁻⁸, OR = 3.48, 95% CI = 2.23–5.43). MS risk at HLA-DQβ1 Phe9 was independent of HLA-DRB1*15:01 (P_Conditional = 1.5 × 10⁻⁵, OR = 2.91, 95% CI = 1.79–4.72), while HLA-DRB1*15:01 was just significant when conditioned on HLA-DQβ1 Phe9 (P_Conditional = 0.037). Regarding a case-control analysis for NMOSD, HLA-DQA1*05:03 had a significant association with NMOSD (P = 1.5 × 10⁻⁴, OR = 6.96, 95% CI = 2.55–19.0).

Conclusions

We identified HLA variants associated with the risk of MS and NMOSD. Our study contributes to the understanding of the genetic architecture of MS and NMOSD in the Japanese population.

Electronic supplementary material

The online version of this article (10.1186/s12974-019-1551-z) contains supplementary material, which is available to authorized users.

Genetic and phenotypic landscape of the major histocompatibilty complex region in the Japanese population

To perform detailed fine-mapping of the major-histocompatibility-complex region, we conducted next-generation sequencing (NGS)-based typing of the 33 human leukocyte antigen (HLA) genes in 1,120 individuals of Japanese ancestry, providing a high-resolution allele catalog and linkage-disequilibrium structure of both classical and nonclassical HLA genes. Together with population-specific deep-whole-genome-sequencing data (n = 1,276), we conducted NGS-based HLA, single-nucleotide-variant and indel imputation of large-scale genome-wide-association-study data from 166,190 Japanese individuals. A phenome-wide association study assessing 106 clinical phenotypes identified abundant, significant genotype-phenotype associations across 52 phenotypes. Fine-mapping highlighted multiple association patterns conferring independent risks from classical HLA genes. Region-wide heritability estimates and genetic-correlation network analysis elucidated the polygenic architecture shared across the phenotypes.

Dynamic epimarks in sex-related genes predict gonad phenotype in the European sea bass, a fish with mixed genetic and environmental sex determination

The integration of genomic and environmental influences into methylation patterns to bring about a phenotype is of central interest in developmental epigenetics, but many details are still unclear. The sex ratios of the species used here, the European sea bass, are determined by genetic and temperature influences. We created four families from parents known to produce offspring with different sex ratios, exposed larvae to masculinizing temperatures and examined, in juvenile gonads, the DNA methylation of seven genes related to sexual development by a targeted sequencing approach. The genes most affected by both genetics and environment were cyp19a1a and dmrt1, with contrasting sex-specific methylation and temperature responses. The relationship between cyp19a1a methylation and expression is relevant to the epigenetic regulation of vertebrate sex, and we report the evidence of such relationship only below a methylation threshold, ~ 80%, and that it was sex-specific: negatively correlated in females but positively correlated in males. From parents to offspring, the methylation in gonads was midway between oocytes and sperm, with bias towards oocytes for amh-r2, er-β2, fsh-r and cyp19a1a. In contrast, dmrt1 levels resembled those of sperm. The methylation of individual CpGs from foxl2, er-β2 and nr3c1 were conserved from parents to offspring, whereas those of cyp19a1a, dmrt1 and amh-r2 were affected by temperature. Utilizing a machine-learning procedure based on the methylation levels of a selected set of CpGs, we present the first, to our knowledge, system based on epigenetic marks capable of predicting sex in an animal with ~ 90% accuracy and discuss possible applications.

Quantification of massively parallel sequencing libraries - a comparative study of eight methods

Quantification of massively parallel sequencing libraries is important for acquisition of monoclonal beads or clusters prior to clonal amplification and to avoid large variations in library coverage when multiple samples are included in one sequencing analysis. No gold standard for quantification of libraries exists. We assessed eight methods of quantification of libraries by quantifying 54 amplicon, six capture, and six shotgun fragment libraries. Chemically synthesized double-stranded DNA was also quantified. Light spectrophotometry, i.e. NanoDrop, was found to give the highest concentration estimates followed by Qubit and electrophoresis-based instruments (Bioanalyzer, TapeStation, GX Touch, and Fragment Analyzer), while SYBR Green and TaqMan based qPCR assays gave the lowest estimates. qPCR gave more accurate predictions of sequencing coverage than Qubit and TapeStation did. Costs, time-consumption, workflow simplicity, and ability to quantify multiple samples are discussed. Technical specifications, advantages, and disadvantages of the various methods are pointed out.

Identification of an HLA class I allele closely involved in the autoantigen presentation in acquired aplastic anemia

Somatic mutations of HLA-B*40:02 are very frequently detected in granulocyte of patients with acquired aplastic anemia. Antigen presentation via HLA-B4002 may play a critical role in the pathophysiology of acquired aplastic anemia.

Use of the melting curve assay as a means for high-throughput quantification of Illumina sequencing libraries

Background. Multiplexed sequencing is commonly performed on massively parallel short-read sequencing platforms such as Illumina, and the efficiency of library normalisation can affect the quality of the output dataset. Although several library normalisation approaches have been established, none are ideal for highly multiplexed sequencing due to issues of cost and/or processing time.

Methods. An inexpensive and high-throughput library quantification method has been developed, based on an adaptation of the melting curve assay. Sequencing libraries were subjected to the assay using the Bio-Rad Laboratories CFX Connect^TM Real-Time PCR Detection System. The library quantity was calculated through summation of reduction of relative fluorescence units between 86 and 95 °C.

Results.PCR-enriched sequencing libraries are suitable for this quantification without pre-purification of DNA. Short DNA molecules, which ideally should be eliminated from the library for subsequent processing, were differentiated from the target DNA in a mixture on the basis of differences in melting temperature. Quantification results for long sequences targeted using the melting curve assay were correlated with those from existing methods (R² > 0.77), and that observed from MiSeq sequencing (R² = 0.82).

Discussion.The results of multiplexed sequencing suggested that the normalisation performance of the described method is equivalent to that of another recently reported high-throughput bead-based method, BeNUS. However, costs for the melting curve assay are considerably lower and processing times shorter than those of other existing methods, suggesting greater suitability for highly multiplexed sequencing applications.

The impact of next-generation sequencing technologies on HLA research

In the past decade, the development of next-generation sequencing (NGS) has paved the way for whole-genome analysis in individuals. Research on the human leukocyte antigen (HLA), an extensively studied molecule involved in immunity, has benefitted from NGS technologies. The HLA region, a 3.6-Mb segment of the human genome at 6p21, has been associated with more than 100 different diseases, primarily autoimmune diseases. Recently, the HLA region has received much attention because severe adverse effects of various drugs are associated with particular HLA alleles. Owing to the complex nature of the HLA genes, classical direct sequencing methods cannot comprehensively elucidate the genomic makeup of HLA genes. Thus far, several high-throughput HLA-typing methods using NGS have been developed. In HLA research, NGS facilitates complete HLA sequencing and is expected to improve our understanding of the mechanisms through which HLA genes are modulated, including transcription, regulation of gene expression and epigenetics. Most importantly, NGS may also permit the analysis of HLA-omics. In this review, we summarize the impact of NGS on HLA research, with a focus on the potential for clinical applications.