The hereditary N363K POLE exonuclease mutant extends PPAP tumor spectrum to glioblastomas by causing DNA damage and aneuploidy in addition to increased mismatch mutagenicity.

The exonuclease domain of DNA polymerases epsilon's catalytic subunit (POLE) removes misincorporated nucleotides, called proofreading. POLE-exonuclease mutations cause colorectal- and endometrial cancers with an extreme burden of single nucleotide substitutions. We recently reported that particularly the hereditary POLE exonuclease mutation N363K predisposes in addition to aggressive giant cell glioblastomas. We knocked-in this mutation homozygously into human cell lines and compared its properties to knock-ins of the likewise hereditary POLE L424V mutation and to a complete proofreading-inactivating mutation (exo-null). We found that N363K cells have higher mutation rates as both L424V- or exo-null mutant cells. In contrast to L424V cells, N363K cells expose a growth defect, replication stress and DNA damage. In non-transformed cells, these burdens lead to aneuploidy but macroscopically normal nuclei. In contrast, transformed N363K cells phenocopy the enlarged and disorganized nuclei of giant cell glioblastomas. Taken together, our data characterize a POLE exonuclease domain mutant that not only causes single nucleotide hypermutation, but in addition DNA damage and chromosome instability, leading to an extended tumor spectrum. Our results expand the understanding of the polymerase exonuclease domain and suggest that an assessment of both the mutational potential and the genetic instability might refine classification and treatment of POLE-mutated tumors.

Evaluation of somatic and/or germline mosaicism in congenital malformation of the eye.

Microphthalmia, Anophthalmia and Coloboma (MAC) form a spectrum of congenital eye malformations responsible for severe visual impairment. Despite the exploration of hundreds of genes by High-Throughput Sequencing (HTS), most of the patients remain without genetic diagnosis. One explanation could be the not yet demonstrated involvement of somatic mosaicism (undetected by conventional analysis pipelines) in those patients. Furthermore, the proportion of parental germline mosaicism in presumed de novo variations is still unknown in ocular malformations. Thus, using dedicated bioinformatics pipeline designed to detect mosaic variants, we reanalysed the sequencing data obtained from a 119 ocular development genes panel performed on blood samples of 78 probands with sporadic MAC without genetic diagnosis. Using the same HTS strategy, we sequenced 80 asymptomatic parents of 41 probands carrying a disease-causing variant in an ocular development gene considered de novo after Sanger sequencing of both parents. Reanalysis of the previously sequencing data did not find any mosaic variant in probands without genetic diagnosis. However, HTS of parents revealed undetected SOX2 and PAX6 mosaic variants in two parents. Finally, this work, performed on two large cohorts of patients with MAC spectrum, provides for the first time an overview of the interest of looking for mosaicism in ocular development disorders. Somatic mosaicism does not appear to be frequent in MAC spectrum and might explain only few diagnoses. Thus, other approaches such as whole genome sequencing should be considered in those patients. Parental mosaicism is however not that rare (around 5%) and challenging for genetic counselling.

Resistance of B-Cell Lymphomas to CAR T-Cell Therapy Is Associated With Genomic Tumor Changes Which Can Result in Transdifferentiation.

Despite the impressive efficacy of chimeric antigen receptor (CAR) T-cell therapy (CART) in B-cell non-Hodgkin lymphomas, durable responses are uncommon. The histopathologic and molecular features associated with treatment failure are still largely unknown. Therefore, we have analyzed 19 sequential tumor samples from 9 patients, prior anti-CD19 CART (pre-CART) and at relapse (post-CART), using immunohistochemistry, fluorescence in situ hybridization, array comparative genomic hybridization, next-generation DNA and RNA sequencing, and genome-scale DNA methylation. The initial diagnosis was diffuse large B-cell lymphoma (n=6), double-hit high-grade B-cell lymphoma (n=1), and Burkitt lymphoma (n=2). Histopathologic features were mostly retained at relapse in 7/9 patients, except the frequent loss of 1 or several B-cell markers. The remaining 2 cases (1 diffuse large B-cell lymphoma and 1 Burkitt lymphoma) displayed a dramatic phenotypic shift in post-CART tumors, with the drastic downfall of B-cell markers and emergence of T-cell or histiocytic markers, despite the persistence of identical clonal immunoglobulin gene rearrangements. The post-CART tumor with aberrant T-cell phenotype showed reduced mRNA expression of most B-cell genes with increased methylation of their promoter. Fluorescence in situ hybridization and comparative genomic hybridization showed global stability of chromosomal alterations in all paired samples, including 17p/TP53 deletions. New pathogenic variants acquired in post-CART samples included mutations triggering the PI3K pathway (PIK3R1, PIK3R2, PIK3C2G) or associated with tumor aggressiveness (KRAS, INPP4B, SF3B1, SYNE1, TBL1XR1). These results indicate that CART-resistant B-cell non-Hodgkin lymphomas display genetic remodeling, which may result in profound dysregulation of B-cell differentiation. Acquired mutations in the PI3K and KRAS pathways suggest that some targeted therapies could be useful to overcome CART resistance.

Molecular diagnosis of T-cell lymphoma: a correlative study of PCR-based T-cell clonality assessment and targeted NGS.

Immunomorphological diagnosis of T-cell lymphoma (TCL) may be challenging, especially on needle biopsies. Multiplex polymerase chain reaction (PCR) assays to assess T-cell receptor (TCR) gene rearrangements are now widely used to detect T-cell clones and provide diagnostic support. However, PCR assays detect only 80% of TCL, and clonal lymphocyte populations may also appear in nonneoplastic conditions. More recently, targeted next-generation sequencing (t-NGS) technologies have been deployed to improve lymphoma classification. To the best of our knowledge, the comparison of these techniques' performance in TCL diagnosis has not been reported yet. In this study, 82 TCL samples and 25 nonneoplastic T-cell infiltrates were divided into 2 cohorts (test and validation) and analyzed with both multiplex PCR and t-NGS to investigate TCR gene rearrangements and somatic mutations, respectively. The detection of mutations appeared to be more specific (100.0%) than T-cell clonality assessment (41.7%-45.5%), whereas no differences were observed in terms of sensitivity (95.1%-97.4%). Furthermore, t-NGS provided a reliable basis for TCL diagnosis in samples with partially degraded DNA that was impossible to assess with PCR. Finally, although multiplex PCR assays appeared to be less specific than t-NGS, both techniques remain complementary, as PCR recovered some t-NGS negative cases.

High-grade Follicular Lymphomas Exhibit Clinicopathologic, Cytogenetic, and Molecular Diversity Extending Beyond Grades 3A and 3B.

Although follicular lymphoma (FL) is usually graded as FL1-2, FL3A, and FL3B, some borderline cases can be observed and led us to investigate the clinicopathologic diversity of grade 3 FL (FL3). Among 2449 FL patients enrolled in Lymphoma Study Association (LYSA) trials, 1921 cases with sufficient material underwent a central pathologic review. The resulting diagnoses comprised 89.6% FL1-2 (n=1723), 7.2% FL3A (n=138), and 0.5% purely follicular FL3B (n=9). The remaining 51 unclassifiable cases (2.7%) exhibited high-grade features but did not meet WHO criteria for either FL3A or FL3B; and were considered as "unconventional" high-grade FL (FL3U). FL3U morphological pattern consisted of nodular proliferation of large cleaved cells or small-sized to medium-sized blast cells. Compared with FL3A, FL3U exhibited higher MUM1 and Ki67 expression, less BCL2 breaks and more BCL6 rearrangements, together with a higher number of cases without any BCL2, BCL6 or MYC rearrangement. FL3U harbored less frequent mutations in BCL2, KMT2D, KMT2B, and CREBBP than FL3A. MYC and BCL2 were less frequently mutated in FL3U than FL3B. Rituximab cyclophosphamide, doxorubicin, vincristine, and prednisone treated FL3U patients had a worse survival than FL1-2 patients with similar follicular lymphoma international prognostic index and treatment. These results suggest that high-grade FLs encompass a heterogeneous spectrum of tumors with variable morphology and genomic alterations, including FL3U cases that do not strictly fit WHO criteria for either FL3A or FL3B, and display a worse outcome than FL1-2. The distinction of FL3U may be useful to allow a better comprehension of high-grade FLs and to design clinical trials.

Epstein-Barr Virus-Associated T- and NK-Cell Lymphoproliferative Diseases: A Review of Clinical and Pathological Features.

Epstein-Barr virus (EBV) is a ubiquitous virus detected in up to 95% of the general population. Most people are asymptomatic, while some may develop a wide range of EBV-associated lymphoproliferative disorders (LPD). Among them, EBV-positive T/NK LPD are uncommon diseases defined by the proliferation of T- or NK-cells infected by EBV. The 2017 World Health Organization (WHO) classification recognizes the following entities characterized by different outcomes: chronic active EBV infection of T- or NK-cell types (cutaneous and systemic forms), systemic EBV-positive T-cell lymphoma of childhood, EBV-positive aggressive NK-cell leukemia, extra nodal NK/T-cell lymphoma nasal type, and the new provisional entity known as primary EBV-positive nodal T/NK-cell lymphoma. In addition, EBV associated-hemophagocytic lymphohistiocytosis is part of EBV-positive T/NK LPD, but has not been included in the WHO classification due to its reactive nature. Despite novel insights from high-throughput molecular studies, EBV-positive NK/T-cell LPD diagnoses remain challenging, especially because of their rarity and overlap. Until now, an accurate EBV-positive NK/T LPD diagnosis has been based on its clinical presentation and course correlated with its histological features. This review aims to summarize clinical, pathological and molecular features of EBV-positive T/NK LPD subtypes and to provide an overview of new understandings regarding these rare disorders.

Management of T-Cell Lymphoma: In Quest of the Holy Grail.

T-cell lymphomas (TCL) represent a very heterogeneous group of lymphoid tumors which are clearly distinct from B-cell neoplasms [...].

Gene alterations in epigenetic modifiers and JAK-STAT signaling are frequent in breast implant-associated ALCL.

The oncogenic events involved in breast implant-associated anaplastic large cell lymphoma (BI-ALCL) remain elusive. To clarify this point, we have characterized the genomic landscape of 34 BI-ALCLs (15 tumor and 19 in situ subtypes) collected from 54 BI-ALCL patients diagnosed through the French Lymphopath network. Whole-exome sequencing (n = 22, with paired tumor/germline DNA) and/or targeted deep sequencing (n = 24) showed recurrent mutations of epigenetic modifiers in 74% of cases, involving notably KMT2C (26%), KMT2D (9%), CHD2 (15%), and CREBBP (15%). KMT2D and KMT2C mutations correlated with a loss of H3K4 mono- and trimethylation by immunohistochemistry. Twenty cases (59%) showed mutations in ≥1 member of the JAK/STAT pathway, including STAT3 (38%), JAK1 (18%), and STAT5B (3%), and in negative regulators, including SOCS3 (6%), SOCS1 (3%), and PTPN1 (3%). These mutations were more frequent in tumor-type samples than in situ samples (P = .038). All BI-ALCLs expressed pSTAT3, regardless of the mutational status of genes in the JAK/STAT pathway. Mutations in the EOMES gene (12%) involved in lymphocyte development, PI3K-AKT/mTOR (6%), and loss-of-function mutations in TP53 (12%) were also identified. Copy-number aberration (CNA) analysis identified recurrent alterations, including gains on chromosomes 2, 9p, 12p, and 21 and losses on 4q, 8p, 15, 16, and 20. Regions of CNA encompassed genes involved in the JAK/STAT pathway and epigenetic regulators. Our results show that the BI-ALCL genomic landscape is characterized by not only JAK/STAT activating mutations but also loss-of-function alterations of epigenetic modifiers.

MIAmS: microsatellite instability detection on NGS amplicons data.

MOTIVATION: Microsatellite instability (MSI) is a molecular marker of DNA mismatch repair deficiency frequently at play in oncogenesis. MSI testing is used for diagnostic, prognostic and therapeutic purposes in several cancers. The current gold standard analysis for microsatellite status is based on length distribution analysis of multiplex-PCR generated DNA fragments from tumor samples which is a laborious and time consuming method. Next generation sequencing (NGS) using amplicon panels is an easy, accurate and scalable technique to determine both the microsatellite status and tumor genome mutations at the same time. Here, we describe MIAmS, an application designed to tag microsatellite status from amplicon NGS of tumor samples. Interestingly, this tool does not require paired normal tissue for comparison. In addition, this scalable application provides a user-friendly report for the interpretation of the results by biologists and exhibits a strong accuracy and robustness for determination of the MSI status. AVAILABILITY: https://github.com/bialimed/miams. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online, evaluation data are available at http://www.ebi.ac.uk/ena/data/view/PRJEB31725.

EWSR1-PATZ1 gene fusion may define a new glioneuronal tumor entity.

We investigated the challenging diagnostic case of a ventricular cystic glioneuronal tumor with papillary features, by RNA sequencing using the Illumina TruSight RNA Fusion panel. We did not retrieve the SLC44A1-PRKCA fusion gene specific for papillary glioneuronal tumor, but an EWSR1-PATZ1 fusion transcript. RT-PCR followed by Sanger sequencing confirmed the EWSR1-PATZ1 fusion. It matched with canonic EWSR1 fusion oncogene, juxtaposing the entire N-terminal transcriptional activation domain of EWSR1 gene and the C-terminal DNA binding domain of a transcription factor gene, PATZ1. PATZ1 protein belongs to the BTB-ZF (broad-complex, tramtrack and bric-à-brac -zinc finger) family. It directly regulates Pou5f1 and Nanog and is essential to maintaining stemness by inhibiting neural differentiation. EWSR1-PATZ1 fusion is a rare event in tumors: it was only reported in six round cell sarcomas and in three gliomas of three exclusively molecular studies. The first reported glioma was a BRAFV600E negative ganglioglioma, the second a BRAFV600E negative glioneuronal tumor, not otherwise specified and the third, very recently reported, a high grade glioma, not otherwise specified. In our study, forty BRAFV600E negative gangliogliomas were screened by FISH using EWSR1 break-apart probes. We performed methylation profiling for the index case and for seven out of the ten FISH positive cases. The index case clustered apart from other pediatric low grade glioneuronal entities, and specifically from the well-defined ganglioglioma methylation group. An additional pediatric intraventricular ganglioglioma clustered slightly more closely with ganglioglioma, but showed differences from the main ganglioglioma group and similarities with the index case. Both cases harbored copy number variations at the PATZ1 locus. EWSR1-PATZ1 gene fusion might define a new type of glioneuronal tumors, distinct from gangliogliomas.

RREB1-MKL2 fusion in biphenotypic "oropharyngeal" sarcoma: New entity or part of the spectrum of biphenotypic sinonasal sarcomas?

An increasing number of sarcomas displaying a primitive, monomorphic spindle cell phenotype have been shown to harbor recurrent gene fusions, including biphenotypic sinonasal sarcoma (SNS). Occurring in the sinonasal area of middle-aged patients, SNS is a locally aggressive tumor harboring in 90% of cases recurrent gene fusions involving the PAX3 gene, in which the chimeric transcription factor induces an aberrant dual myogenic and neural phenotype. Here, we report an unusual oropharyngeal monomorphic spindle cell sarcoma in a 53-year-old man that revealed a novel RREB1-MKL2 gene fusion by RNA sequencing with the Illumina TruSight RNA Fusion Panel. The gene fusion was validated by RT-PCR. Although the tumor location is unusual (but head and neck seated), most of the other clinical, morphologic, immunophenotypic (focal combined expression of S100 protein, SMA, desmin, and myogenin) and oncogenic data suggest that this biphenotypic "oropharyngeal" sarcoma is closely related to the biphenotypic SNS spectrum. Notably, the RREB1-MKL2 chimeric transcription factor encoded by this fusion gene produced an increase in MKL2 expression, which regulates both neural and myogenic differentiation, mimicking the crucial role of PAX3 reported in SNS oncogenesis. NGS and especially RNA sequencing may be used to identify new candidate fusion oncogenes in soft tissue tumors, which would help in updating the existing classification. In turn, this would lead to better therapeutic management of patients.

FROGS: Find, Rapidly, OTUs with Galaxy Solution.

Motivation: Metagenomics leads to major advances in microbial ecology and biologists need user friendly tools to analyze their data on their own. Results: This Galaxy-supported pipeline, called FROGS, is designed to analyze large sets of amplicon sequences and produce abundance tables of Operational Taxonomic Units (OTUs) and their taxonomic affiliation. The clustering uses Swarm. The chimera removal uses VSEARCH, combined with original cross-sample validation. The taxonomic affiliation returns an innovative multi-affiliation output to highlight databases conflicts and uncertainties. Statistical results and numerous graphical illustrations are produced along the way to monitor the pipeline. FROGS was tested for the detection and quantification of OTUs on real and in silico datasets and proved to be rapid, robust and highly sensitive. It compares favorably with the widespread mothur, UPARSE and QIIME. Availability and implementation: Source code and instructions for installation: https://github.com/geraldinepascal/FROGS.git. A companion website: http://frogs.toulouse.inra.fr. Contact: geraldine.pascal@inra.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Root transcriptomic responses of grafted grapevines to heterogeneous nitrogen availability depend on rootstock genotype.

In many fruit species, including grapevine, grafting is used to improve scion productivity and quality and to adapt the plant to environmental conditions. However, the mechanisms underlying the rootstock control of scion development are still poorly understood. The ability of rootstocks to regulate nitrogen uptake and assimilation may contribute to this control. A split-root system was used to grow heterografted grapevines and to investigate the molecular responses to changes in nitrate availability of two rootstocks known to affect scion growth differently. Transcriptome profiling by RNA sequencing was performed on root samples collected 3 and 24 h after nitrogen supply. The results demonstrated a common response involving nitrogen-related genes, as well as a more pronounced transcriptomic reprogramming in the genotype conferring the lower scion growth. A weighted gene co-expression network analysis allowed the identification of co-regulated gene modules, suggesting a role for nitrate transporter 2 family genes and some transcription factors as main actors controlling this genotype-dependent response to heterogeneous nitrogen supply. The relationship between nitrate, ethylene, and strigolactone hormonal pathways was found to differ between the two genotypes. These findings indicated that the genotypes responded differently to heterogeneous nitrogen availability, and this may contribute to their contrasting effect on scion growth.

Compacting and correcting Trinity and Oases RNA-Seq de novo assemblies.

BACKGROUND: De novo transcriptome assembly of short reads is now a common step in expression analysis of organisms lacking a reference genome sequence. Several software packages are available to perform this task. Even if their results are of good quality it is still possible to improve them in several ways including redundancy reduction or error correction. Trinity and Oases are two commonly used de novo transcriptome assemblers. The contig sets they produce are of good quality. Still, their compaction (number of contigs needed to represent the transcriptome) and their quality (chimera and nucleotide error rates) can be improved. RESULTS: We built a de novo RNA-Seq Assembly Pipeline (DRAP) which wraps these two assemblers (Trinity and Oases) in order to improve their results regarding the above-mentioned criteria. DRAP reduces from 1.3 to 15 fold the number of resulting contigs of the assemblies depending on the read set and the assembler used. This article presents seven assembly comparisons showing in some cases drastic improvements when using DRAP. DRAP does not significantly impair assembly quality metrics such are read realignment rate or protein reconstruction counts. CONCLUSION: Transcriptome assembly is a challenging computational task even if good solutions are already available to end-users, these solutions can still be improved while conserving the overall representation and quality of the assembly. The de novo RNA-Seq Assembly Pipeline (DRAP) is an easy to use software package to produce compact and corrected transcript set. DRAP is free, open-source and available under GPL V3 license at http://www.sigenae.org/drap.

Jflow: a workflow management system for web applications.

SUMMARY: Biologists produce large data sets and are in demand of rich and simple web portals in which they can upload and analyze their files. Providing such tools requires to mask the complexity induced by the needed High Performance Computing (HPC) environment. The connection between interface and computing infrastructure is usually specific to each portal. With Jflow, we introduce a Workflow Management System (WMS), composed of jQuery plug-ins which can easily be embedded in any web application and a Python library providing all requested features to setup, run and monitor workflows. AVAILABILITY AND IMPLEMENTATION: Jflow is available under the GNU General Public License (GPL) at http://bioinfo.genotoul.fr/jflow. The package is coming with full documentation, quick start and a running test portal. CONTACT: Jerome.Mariette@toulouse.inra.fr.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds.

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.

Contrasted effects of natural complex mixtures of PAHs and metals on oxygen cycle in a microbial mat.

The contamination of polluted environments is often due to a complex mixture of pollutants sometimes at trace levels which nevertheless may have significant effects on the diversity and functioning of organisms. The aim of this study was to assess the functional responses of a microbial mat exposed to a natural complex mixture of PAHs and metals as a function of the maturation stage of the biofilm. Microbial mats sampled in a slightly polluted environment were exposed to contaminated water of a retention basin of an oil refinery. The responses of the microbial mats differed according to season. In spring 2012, strong inhibition of both oxygen production and respiration was observed relative to the control, with rates representing less than 5% of the control after 72 h of incubation. A decrease of microbial activities was followed by a decrease of the coupling between autotrophs and heterotrophs. In contrast, in autumn 2012, no significant changes for oxygen production and respiration were observed and the coupling between autotrophs and heterotrophs was not altered. The differences observed between the spring and autumn mats might be explained by the maturity of the microbial mat with dominance of heterotrophic bacteria in spring, and diatoms and cyanobacteria in autumn, as well as by the differences in the chemical composition of the complex mixture of PAHs and metals.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities.

BACKGROUND: Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform. RESULTS: The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations. CONCLUSIONS: Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.