Full-length Alleles Research Articles

Natural variation of GmFNSII-2 contributes to drought resistance by modulating enzyme activity in soybean

As an essential crop that provides vegetable oil and protein, soybean (Glycine max (L.) Merr.) is widely planted all over the world. However, the scarcity of water resources worldwide has seriously impacted on the quality and yield of soybean. To address this, exploring excellent genes for improving drought resistance in soybean is crucial. In this study, we identified natural variations of GmFNSII-2 (flavone synthase II) significantly affect the drought resistance of soybeans. Through sequence analysis of GmFNSII-2 in 632 cultivated and 44 wild soybeans nine haplotypes were identified. The full-length allele GmFNSII-2C, but not the truncated allele GmFNSII-2A possessing a nonsense nucleotide variation, increased enzyme activity. Further research found that GmDREB3, known to increase soybean drought resistance, bound to the promoter region of GmFNSII-2C. GmDREB3 positively regulated the expression of GmFNSII-2C, increased flavone synthase abundance and improved the drought resistance. Furthermore, a single-base mutation in the GmFNSII-2C promoter generated an additional drought response element (CCCCT), which had stronger interaction strength with GmDREB3 and increased its transcriptional activity under drought conditions. The frequency of drought-resistant soybean varieties with Hap 1 (Pro:GmFNSII-2C) has increased, suggesting that this haplotype may be selected during soybean breeding. In summary, GmFNSII-2C could be used for molecular breeding of drought-tolerant soybean.

MYB transcription factors encoded by diversified tandem gene clusters cause varied Morella rubra fruit color.

Chinese bayberry (Morella rubra) is a fruit tree with a remarkable variation in fruit color, ranging from white to dark red as determined by anthocyanin content. In dark red "Biqi" (BQ), red "Dongkui" (DK), pink "Fenhong" (FH), and white "Shuijing" (SJ), we identified an anthocyanin-related MYB transcription factor-encoding gene cluster of four members, i.e. MrMYB1.1, MrMYB1.2, MrMYB1.3, and MrMYB2. Collinear analysis revealed that the MYB tandem cluster may have occurred in a highly conserved region of many eudicot genomes. Two alleles of MrMYB1.1 were observed; MrMYB1.1-1 (MrMYB1.1n) was a full-length allele and homozygous in "BQ", MrMYB1.1-2 (MrMYB1.1d) was a nonfunctional allele with a single base deletion and homozygous in "SJ", and MrMYB1.1n/MrMYB1.1d were heterozygous in "DK" and "FH". In these four cultivars, expression of MrMYB1.1, MrMYB1.2, and MrMYB2 was enhanced during ripening. Both alleles were equally expressed in MrMYB1.1n/MrMYB1.1d heterozygous cultivars as revealed by a cleaved amplified polymorphic sequence marker. Expression of MrMYB1.3 was restricted to some dark red cultivars only. Functional characterization revealed that MrMYB1.1n and MrMYB1.3 can induce anthocyanin accumulation while MrMYB1.1d, MrMYB1.2, and MrMYB2 cannot. DNA-protein interaction assays indicated that MrMYB1.1n and MrMYB1.3 can directly bind to and activate the promoters of anthocyanin-related genes via interaction with a MYC-like basic helix-loop-helix protein MrbHLH1. We concluded that the specific genotype of MrMYB1.1 alleles, as well as the exclusive expression of MrMYB1.3 in some dark red cultivars, contributes to fruit color variation. The study provides insights into the mechanisms for regulation of plant anthocyanin accumulation by MYB tandem clusters.

Submitting Novel Full-Length HLA, MIC, and KIR Alleles with TypeLoader2.

The Immuno Polymorphism Database (IPD) plays a pivotal role for immunogenetics. Due to technical limitations, genotyping often focuses on specific key regions like the antigen recognition domain (ARD) for HLA genotyping, and the databases are populated accordingly. More recently, though, modern next generation sequencing (NGS) assays allow using larger gene segments or even complete genes for genotyping. It is therefore essential that the databases are updated with complete genetic reference sequences to fully serve current and future applications. However, the process of manually annotating and submitting full-length allele sequences to IPD is time-consuming and error-prone, which may discourage HLA-genotyping laboratories or researchers from submitting full-length sequences of novel alleles.Here, we detail the process of preparing and submitting novel HLA, MIC, and KIR alleles to ENA and IPD using TypeLoader2, a convenient software tool developed to streamline this process by automating the sequence annotation, the creation of all necessary files, as well as parts of the submission process itself. The software is freely available from GitHub ( https://github.com/DKMS-LSL/typeloader ).

Characterization and genetic diversity of MHC class II DRB genes in the Arabian camel (Camelus dromedarius).

This study investigated the MHC DRB genes in the Arabian camel (Camelus dromedarius). The results revealed the presence of - at least - two transcribed DRB-like genes in chromosome 20, designated MhcCadr-DRB1 and MhcCadr-DRB2. These genes are 155Kb apart, have similar gene structure, and are transcribed in opposite directions. Compared to DRB1, the DRB2 locus contains a deletion of 12 nucleotides in the second exon (270bp), exhibits lower transcript abundance, and is expressed as two splice variants differing by exon 2 skipping. This gene seems to be of minor functional relevance in the dromedary camel. Conversely, the DRB1 is thought to be the main gene in this species showing higher transcript abundance and polymorphism levels. A total of seven DRB1 exon 2 alleles were identified in the Tunisian dromedary camel resulting from 18 amino acid substitutions. Six full length alleles were characterized at the mRNA level. Although there is no clear evidence for balancing selection (i.e., heterozygote advantage), signals of weak historical positive selection acting on the DRB1 gene were detected, as indicated by the limited number of the sites being positively selected. This trend might be related to the low exposure to pathogens and to the demographic history of the species. Comparative analysis with Bactrian and wild camel genomes suggested occurrence of trans species polymorphism (TSP) in the Camelus genus. The results lay the foundation for the MHC DRB1 genetic diversity analysis in this genus since the developed genotyping protocols are fully applicable in the three Camelus species.

Allelic variation of TaWD40-4B.1 contributes to drought tolerance by modulating catalase activity in wheat

Drought drastically restricts wheat production, so to dissect allelic variations of drought tolerant genes without imposing trade-offs between tolerance and yield is essential to cope with the circumstance. Here, we identify a drought tolerant WD40 protein encoding gene TaWD40-4B.1 of wheat via the genome-wide association study. The full-length allele TaWD40-4B.1C but not the truncated allele TaWD40-4B.1T possessing a nonsense nucleotide variation enhances drought tolerance and grain yield of wheat under drought. TaWD40-4B.1C interacts with canonical catalases, promotes their oligomerization and activities, and reduces H2O2 levels under drought. The knock-down of catalase genes erases the role of TaWD40-4B.1C in drought tolerance. TaWD40-4B.1C proportion in wheat accessions is negatively correlative with the annual rainfall, suggesting this allele may be selected during wheat breeding. The introgression of TaWD40-4B.1C enhances drought tolerance of the cultivar harboring TaWD40-4B.1T. Therefore, TaWD40-4B.1C could be useful for molecular breeding of drought tolerant wheat.

Alternative splicing and genetic variation of mhc-e: implications for rhesus cytomegalovirus-based vaccines

Rhesus cytomegalovirus (RhCMV)-based vaccination against Simian Immunodeficiency virus (SIV) elicits MHC-E-restricted CD8+ T cells that stringently control SIV infection in ~55% of vaccinated rhesus macaques (RM). However, it is unclear how accurately the RM model reflects HLA-E immunobiology in humans. Using long-read sequencing, we identified 16 Mamu-E isoforms and all Mamu-E splicing junctions were detected among HLA-E isoforms in humans. We also obtained the complete Mamu-E genomic sequences covering the full coding regions of 59 RM from a RhCMV/SIV vaccine study. The Mamu-E gene was duplicated in 32 (54%) of 59 RM. Among four groups of Mamu-E alleles: three ~5% divergent full-length allele groups (G1, G2, G2_LTR) and a fourth monomorphic group (G3) with a deletion encompassing the canonical Mamu-E exon 6, the presence of G2_LTR alleles was significantly (p = 0.02) associated with the lack of RhCMV/SIV vaccine protection. These genomic resources will facilitate additional MHC-E targeted translational research.

HLA-Inception: A structure-based deep learning framework for MHC-I binding motif prediction

Abstract The ability to accurately identify peptide ligands for a given major histocompatibility complex class I (MHC-I) molecule has immense value for targeted anticancer and antiviral therapeutics. However, the highly polymorphic nature of the MHC-I protein makes universal prediction of peptide ligands challenging due to lack of experimental data describing most MHC-I variants, and the vast number of protein variants precludes comprehensive experimental determination. Therefore, there is a need for a framework to cluster MHC-I alleles to prioritize for experimental validation as well as identify alleles with potential disease associations. To address this challenge, we have developed a deep convolutional neural network, HLA-Inception, capable of predicting MHC-I peptide binding motif using data derived from the structure of the MHC-I binding pocket. By approaching this problem from a 3-dimensional perspective, we can fully consider the impact of sidechain arrangement and topology of the MHC-I binding pocket on peptide binding motif, which is not inherently captured by the popular protein sequence-based approaches. Through a combination of homology modeling and biophysical simulations, we created protein structure models for all full-length HLA-ABC alleles. The topology and interaction forces within the MHC-I binding pocket were accounted for by solving the 3-dimensional electrostatic potential near the surface of the protein. HLA-Inception was then trained on all MHC-I alleles with known MHC-I binding motifs and applied to the full set of MHC-I models. We found that predicted peptide binding motifs fell into distinct and well-defined clusters which maintained known peptide binding and disease associations.

Reliable Resolution of Full-Length Huntingtin Alleles by Quantitative Immunoblotting.

Therapeutics that lower mutant huntingtin (mHTT) have shown promise in preclinical studies and are in clinical development for the treatment of Huntington disease (HD). Multiple assays have been developed that either quantify mHTT or total HTT but may not accurately measure levels of wild type HTT (wtHTT) in biological samples. To optimize a method that can be used to resolve, quantify and directly compare levels of full length wtHTT and mHTT in HD samples. We provide a detailed quantitative immunoblotting protocol to reproducibly resolve full length wtHTT and mHTT in multiple HD mouse and patient samples. We show that this assay can be modified, depending on the sample, to resolve wtHTT and mHTT with a wide range of polyglutamine length differences (ΔQs 22-179). We also demonstrate that this method can be used to quantify allele-selective lowering of mHTT using an antisense oligonucleotide in HD patient-derived cells. This quantitative immunoblotting method can be used to reliably resolve full-length HTT alleles with ΔQs≥22 and allows for direct comparison of wtHTT and mHTT levels in HD samples.

Global distribution of IRC7 alleles in Saccharomyces cerevisiae populations: a genomic and phenotypic survey within the wine clade

The adaptation to the different biotic and abiotic factors of wine fermentation has led to the accumulation of numerous genomic hallmarks in Saccharomyces cerevisiae wine strains. IRC7, a gene encoding a cysteine-S-β-lyase enzyme related volatile thiols production in wines, has two alleles: a full-length allele (IRC7F ) and a mutated one (IRC7S ), harbouring a 38bp-deletion. Interestingly, IRC7S -encoding a less active enzyme- appears widespread amongst wine populations. Studying the global distribution of the IRC7S allele in different yeast lineages, we confirmed its high prevalence in the Wine clade and demonstrated a minority presence in other domesticated clades (Wine-PDM, Beer and Bread) while it is completely missing in wild clades. Here we show that IRC7S -homozygous (HS) strains exhibited both fitness and competitive advantages compared to IRC7F -homozygous (HF) strains. There are some pieces of evidence of the direct contribution of the IRC7S allele to the outstanding behaviour of HS strains (i.e. improved response to oxidative stress conditions and higher tolerance to high copper levels); however, we also identified a set of sequence variants with significant co-occurrence patterns with the IRC7S allele, which can be co-contributing to the fitness and competitive advantages of HS strains in wine fermentations. This article is protected by copyright. All rights reserved.

Estimation of Full-Length TprK Diversity in Treponema pallidum subsp. pallidum.

Immune evasion and disease progression of Treponema pallidum subsp. pallidum are associated with sequence diversity in the hypervariable outer membrane protein TprK. Previous attempts to study variation within TprK have sequenced at depths insufficient to fully appreciate the hypervariable nature of the protein, failed to establish linkage between the protein's seven variable regions, or were conducted on isolates passed through rabbits. As a consequence, a complete profile of tprK during infection in the human host is still lacking. Furthermore, prior studies examining how T. pallidum subsp. pallidum uses its repertoire of genomic donor sites to generate diversity within the variable regions of the tprK have yielded a partial understanding of this process due to the limited number of tprK alleles examined. In this study, we used short- and long-read deep sequencing to directly characterize full-length tprK alleles from T. pallidum subsp. pallidum collected from early lesions of patients attending two sexually transmitted infection clinics in Italy. We demonstrate that strains collected from cases of secondary syphilis contain significantly more unique variable region sequences and full-length TprK sequences than those from cases of primary syphilis. Our data, combined with recent data available on Chinese T. pallidum subsp. pallidum specimens, show the near-complete absence of overlap in TprK sequences among the 41 specimens profiled to date. We further estimate that the potential antigenic variability carried by TprK rivals that of current estimates of the human adaptive immune system. These data underscore the immunoevasive ability of TprK that allows T. pallidum subsp. pallidum to establish lifelong infection.IMPORTANCE Syphilis continues to be a significant public health issue in both low- and high-income countries, including the United States where the rate of syphilis infection has increased over the past 5 years. Treponema pallidum subsp. pallidum, the causative agent of syphilis, carries the outer membrane protein TprK that undergoes segmental gene conversion to constantly create new sequences. We performed full-length deep sequencing of TprK to examine TprK diversity in clinical T. pallidum subsp. pallidum strains. We then combined our results with data from all samples for which TprK deep sequencing results were available. We found almost no overlap in TprK sequences between different patients. Moreover, our data allowed us to estimate the total number of TprK variants that T. pallidum subsp. pallidum can potentially generate. Our results support how the T. pallidum subsp. pallidum TprK antigenic variation system is an equal adversary of the human immune system leading to pathogen persistence in the host.

Characterization of 17 full‐length MHC‐DQB1 alleles in Tibetan macaques (Macaca thibetana)

Seventeen full-length Math-DQB1 alleles were characterized in Tibetan macaques.

Functional analysis of an essential Ran-binding protein gene, CpRbp1, from the chestnut blight fungus Cryphonectria parasitica using heterokaryon rescue

A Ran binding protein (RanBP) homolog, CpRbp1, from Cryphonectria parasitica, has been identified as a protein that is affected by hypovirus infection or tannic acid supplementation. In this study, functional analyses of CpRbp1 were performed by constructing a knockout mutant and analyzing the resulting heterokaryon. Transformation-mediated gene replacement resulted in two putative CpRbp1-null mutants and genotype analyses identified these two mutants as heterokaryotic transformants consisting of two types of nuclei, one with the wild-type CpRbp1 allele and another with the CpRbp1-null mutant allele. Although stable mycelial growth of the heterokaryotic transformant was observed on selective medium containing hygromycin B, neither germination nor growth of the resulting conidia, which were single-cell monokaryotic progeny, was observed on the medium. In trans complementation of heterokaryons using a full-length wild-type allele of the CpRbp1 gene resulted in complemented transformants. These transformants sporulated single-cell monokaryotic conidia that were able to grow on media selective for replacing and/or complementing markers. These results clearly indicate that CpRbp1 is an essential gene, and heterokaryons allowed the fungus to maintain lethal CpRbp1-null mutant nuclei. Moreover, in trans complementation of heterokaryons using chimeric structures of the CpRbp1 gene allowed for analysis of its functional domains, which was previously hampered due to the lethality of the gene. In addition, in trans complementation using heterologous RanBP genes from Aspergillus nidulans was successful, suggesting that the function of RanBP is conserved during evolution. Furthermore, in trans complementation allowed for functional analyses of lethal orthologs. This study demonstrates that our fungal heterokaryon system can be applied effectively to determine whether a gene of interest is essential, perform functional analyses of a lethal gene, and analyze corresponding heterologous genes.

Analysis of Plasmodium falciparum Rh2b deletion polymorphism across different transmission areas

Despite significant progress in controlling malaria, the disease remains a global health burden. The intricate interactions the parasite Plasmodium falciparum has with its host allows it to grow and multiply in human erythrocytes. The mechanism by which P. falciparum merozoites invade human erythrocytes is complex, involving merozoite proteins as well as erythrocyte surface proteins. Members of the P. falciparum reticulocyte binding-like protein homolog (PfRh) family of proteins play a pivotal role in merozoite invasion and hence are important targets of immune responses. Domains within the PfRh2b protein have been implicated in its ability to stimulate natural protective antibodies in patients. More specifically, a 0.58 kbp deletion, at the C-terminus has been reported in high frequencies in Senegalese and Southeast Asian parasite populations, suggesting a possible role in immune evasion. We analysed 1218 P. falciparum clinical isolates, and the results show that this deletion is present in Ghanaian parasite populations (48.5% of all isolates), with Kintampo (hyper-endemic, 53.2%), followed by Accra (Hypo-endemic, 50.3%), Cape Coast (meso-endemic, 47.9%) and Sogakope (meso-endemic, 43.15%). Further analysis of parasite genomes stored in the MalariaGEN database revealed that the deletion variant was common across transmission areas globally, with an overall frequency of about 27.1%. Interestingly, some parasite isolates possessed mixed PfRh2b deletion and full-length alleles. We further showed that levels of antibodies to the domain of PfRh2 protein were similar to antibody levels of PfRh5, indicating it is less recognized by the immune system.

TypeLoader2: Automated submission of novel HLA and killer-cell immunoglobulin-like receptor alleles in full length.

The Immuno Polymorphism Database (IPD) databases provide global, curated repositories for information regarding polymorphisms of genes of the immune system, thereby generating immense value for the research and clinical communities. The advent of high‐throughput genotyping in immunogenetics has led to dramatically growing numbers of heretofore unknown HLA and lately also killer‐cell immunoglobulin‐like receptor (KIR) alleles, which are to be curated and deposited in the IPD‐IMGT/HLA and IPD‐KIR databases, respectively. It is highly desirable that these novel alleles are characterised and submitted in full length, and that known alleles are extended to cover the complete gene sequence.However, the manual annotation and submission of sequences to European Molecular Biology Laboratory's European Nucleotide Archive and the IPD‐IMGT/HLA and IPD‐KIR databases is time‐consuming and error‐prone. Here, we report the substantial extension of the HLA allele submission tool TypeLoader, which now also supports the annotation and submission of KIR alleles. To enable a more widespread use of this tool, we have made it available as a stand‐alone application that can easily be installed on standard Windows or Linux computers. Furthermore, an internal SQLite database was added to store a wide range of metadata about each allele. This allows TypeLoader2 to be used as a lab's central information platform for the annotation, curation and submission of full‐length HLA and KIR allele sequences. The software is freely available from GitHub (https://github.com/DKMS-LSL/typeloader).We hope that the increased convenience and scope of TypeLoader2 will foster the submission of more full‐length sequences to the IPD‐IMGT/HLA and IPD‐KIR databases, ultimately promoting the use of full‐length sequencing for genotyping both HLA and KIR.

Implications of HLA in allogeneic stem cell transplantation

Graft-versus-host disease (GVHD) is a major barrier to successful allogeneic stem cell transplantation. Therefore, the identification of an alternative donor with a lower risk of GVHD is important for patients lacking an HLA-identical sibling donor. To date, HLA studies of large cohorts of unrelated hematopoietic stem cell transplantation (UR-HSCT) have provided important and helpful information for donor selection. In UR-HSCT through Japan Marrow Donor Program, patient and donor HLA-A, -B, and -C, and HLA-DRB1 and -DQB1 double mismatches are significant risk factors for severe GVHD and mortality. HLA-DPB1 mismatch does not affect survival; however, it reduces the chances of leukemia relapse. In the analysis of a specific allele effect on transplant outcomes, HLA-C*14:02 was significantly associated with an increased risk of severe acute GVHD. The development of next-generation sequencing (NGS) has enabled full-length HLA allele typing. Evolutionary analysis of the entire HLA-DPB1 revealed that a highly conserved region from exon 3 to 3'UTR provoked acute GVHD that was different from a T-cell epitope mismatching algorithm, reflecting exon 2 polymorphisms. Furthermore, a recent study demonstrated the importance of full-length NGS HLA typing on UR-HSCT outcomes. The usage of NGS may provide important information on the implications of the HLA genes in allogeneic stem cell transplantation.

Challenging SNP impact on caspofungin resistance by full-length FKS1 allele replacement in Candida lusitaniae.

A strain of the opportunistic pathogenic yeast Candida lusitaniae was genetically engineered for full-length replacement of the FKS1 gene encoding the target of echinocandin antifungals in order to assess the impact of FKS mutations on echinocandin resistance and reduced echinocandin susceptibility (RES). FKS1 allelic exchange was achieved by transforming C. lusitaniae with two DNA fragments covering the entire FKS1 ORF. Both fragments overlap a 40 bp region where SNPs or small indels of interest were inserted. To target integration at the FKS1 locus, each DNA fragment was fused with split auxotrophic markers of which complementary truncated parts were previously inserted into the chromosomal regions flanking FKS1, allowing selection on minimal medium. Three SNPs described in the FKS1 hotspot (HS) regions HS1 or HS2 of clinical isolates of Candida albicans were expressed at an equivalent position in C. lusitaniae and were confirmed to confer either reduced susceptibility (F641V) or full resistance (S645P and R1361G) to caspofungin. The F659 deletion reported in an FKS2 allele of Candida glabrata and the naturally occurring P660A substitution in FKS1 of Candida parapsilosis were shown to confer a 256-fold and 6-fold increase in caspofungin MIC, respectively, when introduced into an FKS1 allele of C. lusitaniae. We have successfully developed a C. lusitaniae strain for the expression of full-length FKS1 alleles harbouring known mutations contributing to reduced susceptibility or resistance to caspofungin, thus opening the way for the screening of other FKS1/FKS2 mutations potentially involved in RES.

Reference Grade Characterization of Polymorphisms in Full-Length HLA Class I and II Genes With Short-Read Sequencing on the ION PGM System and Long-Reads Generated by Single Molecule, Real-Time Sequencing on the PacBio Platform.

Although NGS technologies fuel advances in high-throughput HLA genotyping methods for identification and classification of HLA genes to assist with precision medicine efforts in disease and transplantation, the efficiency of these methods are impeded by the absence of adequately-characterized high-frequency HLA allele reference sequence databases for the highly polymorphic HLA gene system. Here, we report on producing a comprehensive collection of full-length HLA allele sequences for eight classical HLA loci found in the Japanese population. We augmented the second-generation short read data generated by the Ion Torrent technology with long amplicon spanning consensus reads delivered by the third-generation SMRT sequencing method to create reference grade high-quality sequences of HLA class I and II gene alleles resolved at the genomic coding and non-coding level. Forty-six DNAs were obtained from a reference set used previously to establish the HLA allele frequency data in Japanese subjects. The samples included alleles with a collective allele frequency in the Japanese population of more than 99.2%. The HLA loci were independently amplified by long-range PCR using previously designed HLA-locus specific primers and subsequently sequenced using SMRT and Ion PGM sequencers. The mapped long and short-reads were used to produce a reference library of consensus HLA allelic sequences with the help of the reference-aware software tool LAA for SMRT Sequencing. A total of 253 distinct alleles were determined for 46 healthy subjects. Of them, 137 were novel alleles: 101 SNVs and/or indels and 36 extended alleles at a partial or full-length level. Comparing the HLA sequences from the perspective of nucleotide diversity revealed that HLA-DRB1 was the most divergent among the eight HLA genes, and that the HLA-DPB1 gene sequences diverged into two distinct groups, DP2 and DP5, with evidence of independent polymorphisms generated in exon 2. We also identified two specific intronic variations in HLA-DRB1 that might be involved in rheumatoid arthritis. In conclusion, full-length HLA allele sequencing by third-generation and second-generation technologies has provided polymorphic gene reference sequences at a genomic allelic resolution including allelic variations assigned up to the field-4 level for a stronger foundation in precision medicine and HLA-related disease and transplantation studies.

Patterns of non-ARD variation in more than 300 full-length HLA-DPB1 alleles

Our understanding of sequence variation in the HLA-DPB1 gene is largely restricted to the hypervariable antigen recognition domain (ARD) encoded by exon 2. Here, we employed a redundant sequencing strategy combining long-read and short-read data to accurately phase and characterise in full length the majority of common and well-documented (CWD) DPB1 alleles as well as alleles with an observed frequency of at least 0.0006% in our predominantly European sample set. We generated 664 DPB1 sequences, comprising 279 distinct allelic variants. This allows us to present the, to date, most comprehensive analysis of the nature and extent of DPB1 sequence variation.The full-length sequence analysis revealed the existence of two highly diverged allele clades. These clades correlate with the rs9277534 A → G variant, a known expression marker located in the 3′-UTR. The two clades are fully differentiated by 174 fixed polymorphisms throughout a 3.6 kb stretch at the 3′-end of DPB1. The region upstream of this differentiation zone is characterised by increasingly shared variation between the clades. The low-expression A clade comprises 59% of the distinct allelic sequences including the three by far most frequent DPB1 alleles, DPB1*04:01, DPB1*02:01 and DPB1*04:02. Alleles in the A clade show reduced nucleotide diversity with an excess of rare variants when compared to the high-expression G clade. This pattern is consistent with a scenario of recent proliferation of A-clade alleles.The full-length characterisation of all but the most rare DPB1 alleles will benefit the application of NGS for DPB1 genotyping and provides a helpful framework for a deeper understanding of high- and low-expression alleles and their implications in the context of unrelated haematopoietic stem-cell transplantation.

Characterization of 25 full-length S-RNase alleles, including flanking regions, from a pool of resequenced apple cultivars.

Data obtained from Illumina resequencing of 63 apple cultivars were used to obtain full-length S-RNase sequences using a strategy based on both alignment and de novo assembly of reads. The reproductive biology of apple is regulated by the S-RNase-based gametophytic self-incompatibility system, that is genetically controlled by the single, multi-genic and multi-allelic S locus. Resequencing of apple cultivars provided a huge amount of genetic data, that can be aligned to the reference genome in order to characterize variation to a genome-wide level. However, this approach is not immediately adaptable to the S-locus, due to some peculiar features such as the high degree of polymorphism, lack of colinearity between haplotypes and extensive presence of repetitive elements. In this study we describe a dedicated procedure aimed at characterizing S-RNase alleles from resequenced cultivars. The S-genotype of 63 apple accessions is reported; the full length coding sequence was determined for the 25S-RNase alleles present in the 63 resequenced cultivars; these included 10 previously incomplete sequences (S 5 , S 6a , S 6b , S 8 , S 11 , S 23 , S 39 , S 46 , S 50 and S 58 ). Moreover, sequence divergence clearly suggests that alleles S 6a and S 6b , proposed to be neutral variants of the same alleles, should be instead considered different specificities. The promoter sequences have also been analyzed, highlighting regions of homology conserved among all the alleles.

KIR2DL1 allele sequence extensions and discovery of 2DL1*0010102 and 2DL1*0010103 alleles by DNA sequencing.

Full-length KIR2DL1 allele sequence extensions characterised by single molecule real-time (SMRT) DNA sequencing.