Assessment of a novel NRAS in-frame tandem duplication causing a myelodysplastic/myeloproliferative neoplasm.

Myelodysplastic/myeloproliferative diseases of childhood cause a relevant disease burden, and many of these diseases may have a fatal course. The use of next-generation sequencing (NGS) has led to the identification of novel genetic variants in patients with these diseases, advancing our understanding of the underlying pathophysiology. However, novel mutations can often only be interpreted as variants of unknown significance (VUS), hindering adequate diagnosis and the use of a targeted therapy. To improve variant interpretation and test targeted therapies in a preclinical setting, we are using a rapid zebrafish embryo model that allows functional evaluation of the novel variant and possible therapeutic approaches within days. Thereby, we accelerate the translation from genetic findings to treatment options. Here, we establish this workflow on a novel in-frame tandem duplication in NRAS (c.192_227dup; p.G75_E76insDS65_G75) identified by Sanger sequencing in a 2.5-year-old patient with an unclassifiable myelodysplastic/myeloproliferative neoplasm (MDS/MPN-U). We show that this variant results in a myeloproliferative phenotype in zebrafish embryos with expansion of immature myeloid cells in the caudal hematopoietic tissue, which can be reversed by MEK inhibition. Thus, we could reclassify the variant from likely pathogenic to pathogenic using the American College of Medical Genetics (ACMG) criteria.

Assessment of BCOR Internal Tandem Duplications in Pediatric Cancers by Targeted RNA Sequencing.

Alterations in the BCOR gene, including internal tandem duplications (ITDs) of exon 15 have emerged as important oncogenic changes that define several diagnostic entities. In pediatric cancers, BCOR ITDs have recurrently been described in clear cell sarcoma of kidney (CCSK), primitive myxoid mesenchymal tumor of infancy (PMMTI), and central nervous system high-grade neuroepithelial tumor with BCOR ITD in exon 15 (HGNET-BCOR ITDex15). In adults, BCOR ITDs are also reported in endometrial and other sarcomas. The utility of multiplex targeted RNA sequencing for the identification of BCOR ITD in pediatric cancers was investigated. All available archival cases of CCSK, PMMTI, and HGNET-BCOR ITDex15 were collected. Each case underwent anchored multiplex PCR library preparation with a custom-designed panel, with BCOR targeted for both fusions and ITDs. BCOR ITD was detected in all cases across three histologic subtypes using the RNA panel, with no other fusions identified in any of the cases. All BCOR ITDs occurred in the final exon, within 16 codons from the stop sequence. Multiplex targeted RNA sequencing from formalin-fixed, paraffin-embedded tissue is successful at identifying BCOR internal tandem duplications. This analysis supports the use of anchored multiplex PCR targeted RNA next-generation sequencing panels for identification of BCOR ITDs in pediatric tumors. The use of post-analytic algorithms to improve the detection of BCOR ITD using DNA panels was also explored.

A simple and low cost tetra-primer ARMS-PCR method for detection triazole-resistant Aspergillus fumigatus.

The mutation at codon L98 accompanied by a tandem repeat of 34 base pairs (TR34/L98H) in the 5´upstream region of cyp51A is the principal mechanism of triazole resistance of Aspergillus fumigatus. We aimed to evaluate a simple and low-cost tetra-primer amplification refractory mutation system (ARMS)-PCR technique for detection of TR34/L98H mutations in the cyp51A gene of azole-resistant A. fumigatus. The tetra-primer ARMS-PCR assay optimized by four primers in one reaction consists of external primers for detection of tandem repeats in the promoter region and internal primers for detection of a point mutation in codon 98 (L98H) in the cyp51A gene of azole-resistant A. fumigatus. The specificity of TR34/L98H mutation detection was assessed by testing 36 clinical and environmental A. fumigatus strains. The tetra-primer ARMS-PCR assay from A. fumigatus, containing wild-type sequence (T allele) and L98H mutation at cyp51A (A allele), yielded two DNA fragments of 908 bp and 740 bp and two of 942 bp and 212 bp, respectively. None of the A. fumigatus isolates without the TR34/L98H mutation yielded false-positive results. The ARMS-PCR assay was 100% concordant with DNA sequencing results. Prevalence and screening of the TR34/L98H mutation in the cyp51A gene in A. fumigatus isolates may now be determined by a fast, low-cost, and simple method in resource-poor settings.

Estimation of duplication history under a stochastic model for tandem repeats.

BACKGROUND: Tandem repeat sequences are common in the genomes of many organisms and are known to cause important phenomena such as gene silencing and rapid morphological changes. Due to the presence of multiple copies of the same pattern in tandem repeats and their high variability, they contain a wealth of information about the mutations that have led to their formation. The ability to extract this information can enhance our understanding of evolutionary mechanisms. RESULTS: We present a stochastic model for the formation of tandem repeats via tandem duplication and substitution mutations. Based on the analysis of this model, we develop a method for estimating the relative mutation rates of duplications and substitutions, as well as the total number of mutations, in the history of a tandem repeat sequence. We validate our estimation method via Monte Carlo simulation and show that it outperforms the state-of-the-art algorithm for discovering the duplication history. We also apply our method to tandem repeat sequences in the human genome, where it demonstrates the different behaviors of micro- and mini-satellites and can be used to compare mutation rates across chromosomes. It is observed that chromosomes that exhibit the highest mutation activity in tandem repeat regions are the same as those thought to have the highest overall mutation rates. However, unlike previous works that rely on comparing human and chimpanzee genomes to measure mutation rates, the proposed method allows us to find chromosomes with the highest mutation activity based on a single genome, in essence by comparing (approximate) copies of the pattern in tandem repeats. CONCLUSION: The prevalence of tandem repeats in most organisms and the efficiency of the proposed method enable studying various aspects of the formation of tandem repeats and the surrounding sequences in a wide range of settings. AVAILABILITY: The implementation of the estimation method is available at http://ips.lab.virginia.edu/smtr .

Chimerism as an important marker in post-transplant monitoring chimerism monitoring.

Cell chimerism determination is important for the monitoring of engraftment dynamics and for relapse prediction. Our cohort of 474 patients was divided into two groups according to the determination methods used over time, and by their chimerism status. A significant difference in survival was observed between mixed vs complete chimerism (P < 0.0001 vs P < 0.0002) in both patient groups, and also vs microchimerism (P = 0.0201) in the second group. Detection of mixed chimerism is thus a high-risk factor, and microchimerism is potentially a risk factor in the post-transplantation course. Methods with a high sensitivity for monitoring cell chimerism significantly improve the assessment of patients post-transplant, and they enable the identification of patients with high relapse risk. Supported by MH CZ-DRO (00023736, UHKT).

PacBio-Based Mitochondrial Genome Assembly of Leucaena trichandra (Leguminosae) and an Intrageneric Assessment of Mitochondrial RNA Editing.

Reconstructions of vascular plant mitochondrial genomes (mt-genomes) are notoriously complicated by rampant recombination that has resulted in comparatively few plant mt-genomes being available. The dearth of plant mitochondrial resources has limited our understanding of mt-genome structural diversity, complex patterns of RNA editing, and the origins of novel mt-genome elements. Here, we use an efficient long read (PacBio) iterative assembly pipeline to generate mt-genome assemblies for Leucaena trichandra (Leguminosae: Caesalpinioideae: mimosoid clade), providing the first assessment of non-papilionoid legume mt-genome content and structure to date. The efficiency of the assembly approach facilitated the exploration of alternative structures that are common place among plant mitochondrial genomes. A compact version (729 kbp) of the recovered assemblies was used to investigate sources of mt-genome size variation among legumes and mt-genome sequence similarity to the legume associated root holoparasite Lophophytum. The genome and an associated suite of transcriptome data from select species of Leucaena permitted an in-depth exploration of RNA editing in a diverse clade of closely related species that includes hybrid lineages. RNA editing in the allotetraploid, Leucaena leucocephala, is consistent with co-option of nearly equal maternal and paternal C-to-U edit components, generating novel combinations of RNA edited sites. A preliminary investigation of L. leucocephala C-to-U edit frequencies identified the potential for a hybrid to generate unique pools of alleles from parental variation through edit frequencies shared with one parental lineage, those intermediate between parents, and transgressive patterns.

Annotation resource of tandem repeat-containing secretory proteins in sixty fungi.

Fungal secretory proteins that interact with host plants are regarded as effectors. Because fungal effectors rarely contain conserved sequence features, identification and annotation of fungal effectors from predicted secretory proteins are difficult using outward comparison methods such as BLAST or hidden Markov model. In desire of more sequence features to prioritize research interests of fungal secretory proteins, this study developed a pipeline to identify tandem repeat (TR) domain within putative secretory proteins and tested a hypothesis that at least one type of TR domain in non-orthologous secretory proteins has emerged from convergent evolution for plant pathogenicity. There were 2804 types of TR domains and a total of 2925 TR-containing secretory proteins found from 60 fungi. There was no conserved type of TR domain shared only by plant pathogens, indicating functional divergence for different types of TR domain and TR-containing secretory proteins. The annotation resource of putative fungal TR-containing secretory proteins provides new sequence features that will be useful for the community interested in fungal effector biology.

A next-generation sequencing-based assay for minimal residual disease assessment in AML patients with FLT3-ITD mutations.

Internal tandem duplications in fms-like tyrosine kinase 3 (FLT3-ITDs) are common in acute myeloid leukemia (AML) and confer a poor prognosis. A sensitive and specific assay for the detection of minimal residual disease (MRD) in FLT3-ITD mutated AML could guide therapy decisions. Existing assays for MRD in FLT3-ITD AML have not been particularly useful because of limited sensitivity. We developed a sensitive and specific MRD assay for FLT3-ITD mutations using next-generation sequencing. The initial validation of this assay was performed by spiking fixed amounts of mutant DNA into wild-type DNA to establish a sensitivity of detection equivalent to ≥1 FLT3-ITD-containing cell in 10 000, with a minimum input of 100 000 cell equivalents of DNA. We subsequently validated the assay in bone marrow samples from patients with FLT3-ITD AML in remission. Finally, we analyzed bone marrow samples from 80 patients with FLT3-ITD relapsed/refractory AML participating in a trial of a novel FLT3 inhibitor, gilteritinib, and demonstrated a relationship between the mutation burden, as detected by the assay, and overall survival. This novel MRD assay is specific and 2 orders of magnitude more sensitive than currently available polymerase chain reaction- or next-generation sequencing-based FLT3-ITD assays. The assay is being prospectively validated in ongoing randomized clinical trials.

Evaluation and assessment of read-mapping by multiple next-generation sequencing aligners based on genome-wide characteristics.

Massive data produced due to the advent of next-generation sequencing (NGS) technology is widely used for biological researches and medical diagnosis. The crucial step in NGS analysis is read alignment or mapping which is computationally intensive and complex. The mapping bias tends to affect the downstream analysis, including detection of polymorphisms. In order to provide guidelines to the biologist for suitable selection of aligners; we have evaluated and benchmarked 5 different aligners (BWA, Bowtie2, NovoAlign, Smalt and Stampy) and their mapping bias based on characteristics of 5 microbial genomes. Two million simulated read pairs of various sizes (36bp, 50bp, 72bp, 100bp, 125bp, 150bp, 200bp, 250bp and 300bp) were aligned. Specific alignment features such as sensitivity of mapping, percentage of properly paired reads, alignment time and effect of tandem repeats on incorrectly mapped reads were evaluated. BWA showed faster alignment followed by Bowtie2 and Smalt. NovoAlign and Stampy were comparatively slower. Most of the aligners showed high sensitivity towards long reads (>100bp) mapping. On the other hand NovoAlign showed higher sensitivity towards both short reads (36bp, 50bp, 72bp) and long reads (>100bp) mappings; It also showed higher sensitivity towards mapping a complex genome like Plasmodium falciparum. The percentage of properly paired reads aligned by NovoAlign, BWA and Stampy were markedly higher. None of the aligners outperforms the others in the benchmark, however the aligners perform differently with genome characteristics. We expect that the results from this study will be useful for the end user to choose aligner, thus enhance the accuracy of read mapping.

DNA Phenotyping: Snapshot of a Criminal.

Recent progresses in genomic technologies provide both opportunities and challenges to forensics.

Deep conservation of human protein tandem repeats within the eukaryotes.

Tandem repeats (TRs) are a major element of protein sequences in all domains of life. They are particularly abundant in mammals, where by conservative estimates one in three proteins contain a TR. High generation-scale duplication and deletion rates were reported for nucleic TR units. However, it is not known whether protein TR units can also be frequently lost or gained providing a source of variation for rapid adaptation of protein function, or alternatively, tend to have conserved TR unit configurations over long evolutionary times. To obtain a systematic picture, we performed a proteome-wide analysis of the mode of evolution for human protein TRs. For this purpose, we propose a novel method for the detection of orthologous TRs based on circular profile hidden Markov models. For all detected TRs, we reconstructed bispecies TR unit phylogenies across 61 eukaryotes ranging from human to yeast. Moreover, we performed additional analyses to correlate functional and structural annotations of human TRs with their mode of evolution. Surprisingly, we find that the vast majority of human TRs are ancient, with TR unit number and order preserved intact since distant speciation events. For example, ≥ 61% of all human TRs have been strongly conserved at least since the root of all mammals, approximately 300 Ma. Further, we find no human protein TR that shows evidence for strong recent duplications and deletions. The results are in contrast to the high generation-scale mutability of nucleic TRs. Presumably, most protein TRs fold into stable and conserved structures that are indispensable for the function of the TR-containing protein. All of our data and results are available for download from http://www.atgc-montpellier.fr/TRE.

Quantitative assessment of ratiometric bimolecular beacons as a tool for imaging single engineered RNA transcripts and measuring gene expression in living cells.

Recently, we developed an oligonucleotide-based probe, ratiometric bimolecular beacon (RBMB), which generates a detectable fluorescent signal in living cells that express the target RNA. Here, we show that RBMBs can also be used to image single RNA transcripts in living cells, when the target RNA is engineered to contain as few as four hybridization sites. Moreover, comparison with single-molecule fluorescence in situ hybridization confirmed that RBMBs could be used to accurately quantify the number of RNA transcripts within individual cells. Measurements of gene expression could be acquired within 30 min and using a wide range of RBMB concentrations. The ability to acquire accurate measurements of RNA copy number in both HT-1080 cells and CHO cells also suggests that RBMBs can be used to image and quantify single RNA transcripts in a wide range of cell lines. Overall, these findings highlight the robustness and versatility of RBMBs as a tool for imaging RNA in live cells. We envision that the unique capabilities of RBMBs will open up new avenues for RNA research.

An assessment of karyotype restructuring in the neoallotetraploid Tragopogon miscellus (Asteraceae).

Tragopogon miscellus and Tragopogon mirus are two rare examples of allopolyploids that have formed recently in nature. Molecular cytogenetic studies have revealed chromosome copy number variation and intergenomic translocations in both allotetraploids. Due to a lack of interstitial chromosome markers, there remained the possibility of additional karyotype restructuring in these neopolyploids, via intrachromosomal and intragenomic rearrangements. To address this issue, we searched for additional high-copy tandem repeats in genomic sequences of the diploid progenitor species-Tragopogon dubius, Tragopogon pratensis and Tragopogon porrifolius-for application to the chromosomes of the allotetraploids. Eight novel repeats were localised by fluorescence in situ hybridisation (FISH) in the diploids; one of these repeats, TTR3, provided interstitial coverage. TTR3 was included in a cocktail with other previously characterised probes, producing better-resolved karyotypes for the three diploids. The cocktail was then used to test a hypothesis of karyotype restructuring in the recent allotetraploid T. miscellus by comparing repeat distributions to its diploid progenitors, T. dubius and T. pratensis. Five individuals of T. miscellus were selected from across the range of karyotypic variation previously observed in natural populations. FISH signal distributions mostly matched those observed in the diploid progenitors, with the exception of several losses or gains of signal at chromosomal subtermini and previously noted intergenomic translocations. Thus, in T. miscellus, we find most changes restricted to the subterminal regions, and although some were recurrent, none of the changes were common to all individuals analysed. We consider these findings in relation to the gene loss reported previously for T. miscellus.

Management practices associated with the bulk tank milk prevalence of Mycoplasma spp. in dairy herds in Northwestern Portugal.

The objective of this study was to evaluate the effect of some management practices on the prevalence of Mycoplasma spp. in Northwestern Portuguese dairy farms from bulk tank milk (BTM) samples. Additionally, the within-herd prevalence of Mycoplasma spp. was also determined, but only in BTM positive herds. From May 2007 to November 2008, 492 BTM samples from 164 dairies randomly chosen in a population of 1234 dairy farms were analyzed. Five herds (3.0%) had positive mycoplasmal culture results, from which 4 out of 164 (2.4%) were Mycoplasma bovis, with simultaneous presence of Mycoplasma bovigenitalium or Mycoplasma canadense in two of those samples. In one out of 164 (0.6%) herds Mycoplasma capricolum subsp. capricolum was also found. In BTM positive Mycoplasma spp. herds, the apparent intra-herd prevalence was low and varied between 2.5% and 4.5%. Multiple locus variable-number of tandem-repeat analysis was conducted in order to compare the genetic relationship between the isolates. Mycoplasma spp. was found to be present in cows with subclinical mastitis with or without California Mastitis Test positive results, hence all cows should be tested when the agent is isolated from bulk tank rather than selecting suspected cows. A multivariable logistic regression using the Firth's penalized maximum likelihood estimation was performed showing that increasing number of lactating cows (OR=1.05; P<0.01) was associated with a higher probability of isolating Mycoplasma spp. On the other hand, identifying problem cows was associated with a lower probability (OR=0.06; P<0.05). Particular importance was given to the prevalence of M. bovis, and the results obtained highlight the need to include this agent in mastitis control protocols in national dairies and in sanitary controls of transitioned animals between European countries.

Assessment of fidelity and utility of the whole-genome amplification for the clinical tests offered in a histocompatibility and immunogenetics laboratory.

Increasing emphasis on the use of molecular tests in a histocompatibility and immunogenetics laboratory (HIL) poses a potential problem of lack of sufficient DNA to perform multiple genetic analyses. In this study, we report the feasibility, fidelity and utility of multiple displacement amplification (MDA) method to perform whole-genome amplification (WGA) to generate DNA specimens that can be analyzed by multiple molecular techniques and can be used for different clinical tests offered by an HIL. The MDA-generated DNA when compared with the native DNA showed 100% congruency in genotyping of 37 genes/loci using multiple downstream molecular techniques: sequence-based typing and sequence-specific primer-based typing for 5 human leukocyte antigen (HLA) class I and II genes (HLA-A, B, C, DRB1 and DQB1), luminex-based sequence-specific oligonucleotide (SSO) genotyping for a panel of 16 killer immunoglobulin-like receptor (KIR) genes and automated fragment size analysis for a panel of 15 short tandem repeat (STR) loci and amelogenin gene. For post-allogeneic hematopoietic cell transplantation (HCT) chimerism analysis, MDA-generated DNA appeared useful for enriching pre-transplant DNA but not for enriching post-transplant chimeric DNA. Overall, our results show that MDA-based WGA could generate DNA of high yield and fidelity that can be used for various clinical tests and research purposes.

X-chromosome inactivation pattern analysis for the assessment of cell clonality in cats.

X-chromosome inactivation pattern (XCIP) analysis has been widely used to assess cell clonality in various types of neoplasms in humans. In the present study, a polymerase chain reaction-based feline XCIP analysis using the feline androgen receptor gene was developed. To construct the system of the analysis, polymorphism in CAG tandem repeats within the feline androgen receptor gene was explored using somatic DNAs from 50 male and 103 female cats. CAG tandem repeats in exon 1 of the feline androgen receptor gene were found to be polymorphic, containing 15 to 22 CAG repeats. Of the 103 female cats, 70 (68%) were heterozygous for the number of CAG repeats, indicating the possible usefulness of XCIP analysis in cats. Application of the feline XCIP analysis to 3 feline mammary gland adenocarcinoma cell lines revealed distinctly skewed XCIPs in these cell lines, indicating their clonal origins. Twelve (80%) of the 15 primary tissue/cell samples obtained from cats with various neoplastic diseases showed skewed XCIPs. Moreover, bone marrow samples from 3 cats with myelodysplastic syndrome were also found to have skewed XCIPs. The polymerase chain reaction-based XCIP analysis developed in this study can provide information on cell clonality in female cats, potentially facilitating the differential diagnosis of various disorders in cats.

High resolution melting PCR to differentiate Mycobacterium avium subsp. paratuberculosis"cattle type" and "sheep type".

In this study, we describe a novel HRM-PCR assay that clearly differentiates the two main types of Mycobacterium avium subsp. paratuberculosis (cattle and sheep) based on the polymorphic variation of a previously described tandem repeat. This modern genotyping technique has several advantages over alternative methods, including cost, ease of use and rapidity.

Markov models of amino acid substitution to study proteins with intrinsically disordered regions.

BACKGROUND: Intrinsically disordered proteins (IDPs) or proteins with disordered regions (IDRs) do not have a well-defined tertiary structure, but perform a multitude of functions, often relying on their native disorder to achieve the binding flexibility through changing to alternative conformations. Intrinsic disorder is frequently found in all three kingdoms of life, and may occur in short stretches or span whole proteins. To date most studies contrasting the differences between ordered and disordered proteins focused on simple summary statistics. Here, we propose an evolutionary approach to study IDPs, and contrast patterns specific to ordered protein regions and the corresponding IDRs. RESULTS: Two empirical Markov models of amino acid substitutions were estimated, based on a large set of multiple sequence alignments with experimentally verified annotations of disordered regions from the DisProt database of IDPs. We applied new methods to detect differences in Markovian evolution and evolutionary rates between IDRs and the corresponding ordered protein regions. Further, we investigated the distribution of IDPs among functional categories, biochemical pathways and their preponderance to contain tandem repeats. CONCLUSIONS: We find significant differences in the evolution between ordered and disordered regions of proteins. Most importantly we find that disorder promoting amino acids are more conserved in IDRs, indicating that in some cases not only amino acid composition but the specific sequence is important for function. This conjecture is also reinforced by the observation that for of our data set IDRs evolve more slowly than the ordered parts of the proteins, while we still support the common view that IDRs in general evolve more quickly. The improvement in model fit indicates a possible improvement for various types of analyses e.g. de novo disorder prediction using a phylogenetic Hidden Markov Model based on our matrices showed a performance similar to other disorder predictors.

An MCMC algorithm for detecting short adjacent repeats shared by multiple sequences.

MOTIVATION: Repeats detection problems are traditionally formulated as string matching or signal processing problems. They cannot readily handle gaps between repeat units and are incapable of detecting repeat patterns shared by multiple sequences. This study detects short adjacent repeats with interunit insertions from multiple sequences. For biological sequences, such studies can shed light on molecular structure, biological function and evolution. RESULTS: The task of detecting short adjacent repeats is formulated as a statistical inference problem by using a probabilistic generative model. An Markov chain Monte Carlo algorithm is proposed to infer the parameters in a de novo fashion. Its applications on synthetic and real biological data show that the new method not only has a competitive edge over existing methods, but also can provide a way to study the structure and the evolution of repeat-containing genes. AVAILABILITY: The related C++ source code and datasets are available at http://ihome.cuhk.edu.hk/%7Eb118998/share/BASARD.zip. CONTACT: xfan@sta.cuhk.edu.hk

Assessment of DNA degradation and the genotyping success of highly degraded samples.

DNA becomes progressively more fragmented as biological tissue degrades resulting in decreasing ability to gain a complete DNA profile. Successful identification of samples exhibiting very high levels of DNA degradation may be complicated by presenting in minute quantities. The industry standard method for human DNA identification utilising short tandem repeats (STR) may produce partial or no DNA profile with such samples. We report a comparative study of genotyping using STRs, mini-STRs and single nucleotide polymorphisms (SNPs) with template at different levels of degradation in varying amounts. Two methods of assessing quantity and quality of a DNA sample prior to genotyping were investigated. The QIAxcel capillary gel electrophoresis system provided a rapid, cost effective screening method for assessing sample quality. A real-time quantitative PCR (qPCR) assay was able to simultaneously quantify total human DNA, male DNA, DNA degradation and PCR inhibition. The extent of DNA degradation could be assessed with reasonable accuracy to 62.5 pg and genomic targets could be quantified to a lower limit of 15.6 pg. The qPCR assay was able to detect male DNA to a lower limit of 20 pg in a 1:1,000 background of female DNA. By considering the amount of DNA and the degradation ratio of a sample, a general prediction of genotyping success using AmpFlSTR® Profiler Plus®, MiniFiler™ kits and SNP analysis can be made. The results indicate mini-STRs and SNP markers are usually more successful in typing degraded samples and in cases of extreme DNA degradation (≤200 bp) and template amounts below 250 pg, mini-STR and SNP analysis yielded significantly more complete profiles and lower match probabilities than corresponding STR profiles.