Universal germline genetic testing in patients with hematologic malignancies using DNA isolated from nail clippings.

Not available.

Cell-free DNA from nail clippings as source of normal control for genomic studies in hematologic malignancies.

Comprehensive genomic sequencing is becoming a critical component in the assessment of hematologic malignancies, with broad implications for patient management. In this context, unequivocally discriminating somatic from germline events is challenging but greatly facilitated by matched analysis of tumor:normal pairs. In contrast to solid tumors, conventional sources of normal control (peripheral blood, buccal swabs, saliva) could be highly involved by the neoplastic process, rendering them unsuitable. In this work we describe our real-world experience using cell free DNA (cfDNA) isolated from nail clippings as an alternate source of normal control, through the dedicated review of 2,610 tumor:nail pairs comprehensively sequenced by MSK-IMPACT-heme. Overall, we find nail cfDNA is a robust source of germline control for paired genomic studies. In a subset of patients, nail DNA may have tumor DNA contamination, reflecting unique attributes of the hematologic disease and transplant history. Contamination is generally low level, but significantly more common among patients with myeloid neoplasms (20.5%; 304/1482) compared to lymphoid diseases (5.4%; 61/1128) and particularly enriched in myeloproliferative neoplasms with marked myelofibrosis. When identified in patients with lymphoid and plasma-cell neoplasms, mutations commonly reflected a myeloid profile and correlated with a concurrent/evolving clonal myeloid neoplasm. For nails collected after allogeneic stem-cell transplantation, donor DNA was identified in 22% (11/50). In this cohort, an association with recent history of graft-vs-host disease was identified. These findings should be considered as a potential limitation for the use of nail as normal control but could also provide important diagnostic information regarding the disease process.

Enhanced clinical assessment of hematologic malignancies through routine paired tumor and normal sequencing.

Genomic profiling of hematologic malignancies has augmented our understanding of variants that contribute to disease pathogenesis and supported development of prognostic models that inform disease management in the clinic. Tumor only sequencing assays are limited in their ability to identify definitive somatic variants, which can lead to ambiguity in clinical reporting and patient management. Here, we describe the MSK-IMPACT Heme cohort, a comprehensive data set of somatic alterations from paired tumor and normal DNA using a hybridization capture-based next generation sequencing platform. We highlight patterns of mutations, copy number alterations, and mutation signatures in a broad set of myeloid and lymphoid neoplasms. We also demonstrate the power of appropriate matching to make definitive somatic calls, including in patients who have undergone allogeneic stem cell transplant. We expect that this resource will further spur research into the pathobiology and clinical utility of clinical sequencing for patients with hematologic neoplasms.

Leptomeningeal disease in histone-mutant gliomas.

Background: The 2016 WHO classification described a subtype of midline gliomas harboring histone 3 (H3) K27M alterations, and the 2021 edition added a new subtype of hemispheric diffuse gliomas with H3 G34R/V mutations. The incidence and clinical behavior of leptomeningeal disease (LMD) in these patients is not well defined. Methods: Retrospective study of patients with H3-altered gliomas diagnosed from 01/2012 to 08/2021; histone mutations were identified through next-generation sequencing (NGS) of tumor biopsy and/or cerebrospinal fluid (CSF). Results: We identified 42 patients harboring H3 mutations (K27M mutations in 33 patients, G34R/V in 8, and both in one). Median age was 21 (4-70); 27 were male. LMD was diagnosed in 21/42 (50%) patients, corresponding to a 3-year cumulative incidence of 44.7% (95% confidence interval (CI): 26.1%-63.4%) for the K27-mutant group and a 1-year cumulative incidence of 37.5% in the G34-mutant group (95% CI: 0.01%-74.4%; no events after 1 year). Median time from tumor diagnosis to LMD was 12.9 months for H3-K27 patients and 5.6 months for H3-G34 patients. H3 mutation was detected in CSF in all patients with LMD who had NGS (8 H3-K27-mutant patients). In the H3-K27-mutant group, modeled risk of death was increased in patients who developed LMD (hazard ratio: 7.37, 95% CI: 2.98-18.23, P < .0001). Conclusions: In our cohort, 50% of patients developed LMD. Although further studies are needed, CSF ctDNA characterization may aid in identifying molecular tumor profiles in glioma patients with LMD, and neuroaxis imaging and CSF NGS should be considered for early LMD detection.

Tumoral Intraductal Neoplasms of the Bile Ducts Comprise Morphologically and Genetically Distinct Entities.

CONTEXT.­: Tumoral (grossly visible) intraductal neoplasms of the bile ducts are still being characterized. OBJECTIVE.­: To investigate their morphologic, immunohistochemical, and molecular features. DESIGN.­: Forty-one cases were classified as gastric-, intestinal-, pancreatobiliary-type intraductal papillary neoplasm (IPN), intraductal oncocytic papillary neoplasm (IOPN), or intraductal tubulopapillary neoplasm (ITPN) on the basis of histology. All neoplasms were subjected to targeted next-generation sequencing. RESULTS.­: The mean age at diagnosis was 69 years (42-81 years); male to female ratio was 1.3. Most neoplasms (n = 23, 56%) were extrahepatic/large (mean size, 4.6 cm). The majority (n = 32, 78%) contained high-grade dysplasia, and 68% (n = 28) revealed invasion. All gastric-type IPNs (n = 9) and most ITPNs/IOPNs showed consistent colabeling for CK7/MUC6, which was less common among others (P = .004). Intestinal-type IPNs (n = 5) showed higher rates of CK20 expression than others (P < .001). Overall, the most commonly mutated genes included TP53 and APC, while copy number variants affected ELF3 and CDKN2A/B. All gastric-type IPNs contained an alteration affecting the Wnt signaling pathway; 7 of 9 (78%) showed aberrations in the MAPK pathway. Mutations in APC and KRAS were common in gastric-type IPNs as compared with others (P = .01 for both). SMAD4 was more frequently mutated in intestinal-type IPNs (P = .02). Pancreatobiliary-type IPNs (n = 14) exhibited frequent alterations in tumor suppressor genes including TP53, CDKN2A/B, and ARID2 (P = .04, P = .01 and P = .002, respectively). Of 6 IOPNs analyzed, 3 (50%) revealed ATP1B1-PRKACB fusion. ITPNs (n = 6) showed relatively few recurrent genetic aberrations. Follow-up information was available for 38 patients (median, 58.5 months). The ratio of disease-related deaths was higher for the cases with invasion (56% versus 10%). CONCLUSIONS.­: Tumoral intraductal neoplasms of the bile ducts, similar to their counterparts in the pancreas, are morphologically and genetically heterogeneous.

Loss of CDKN2A/B is a Molecular Marker of High-grade Histology and is Associated with Aggressive Behavior in Acinic Cell Carcinoma.

Acinic cell carcinoma (AciCC) is a rare salivary gland cancer with excellent prognosis in most cases. However, a subset of patients will develop distant metastasis and die of disease. Recently, a 2-tiered grading scheme in AciCC was proposed to recognize patients at risk of poor outcome. We performed a genetic analysis of AciCC to explore the underlying molecular correlates of the tumor grade and examined programmed death ligand 1 (PD-L1) expression to identify potential candidates for immunotherapy. A retrospective cohort of 55 patients included 34 high-grade (HG) and 21 low-grade AciCCs. Forty-three cases were subjected to targeted exome sequencing by Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets. PD-L1 immunohistochemistry was performed in 33 cases. Tumor mutation burden was low with a median of 1 and 2 mutations in low-grade and HG AciCCs, respectively. CDKN2A/B was the most frequently altered gene, and loss-of-function mutations were found only in HG but not in low-grade AciCCs (18/31 [58.1%] vs 0/12 [0%], P < .001). CDKN2A/B alterations were significantly associated with distant metastasis, which occurred in 16/18 (88.9%) CDKN2A/B mutants versus 11/25 (44%) wild-type cases (P = .004, Fisher exact test). Sequential profiling of multiple temporally distant samples from the same patient demonstrated intratumor heterogeneity, including the detection of CDKN2A/B deletion in the second, in HG metastasis only. ATM and PTEN mutations were detected in 6/31 (19.4%) and 5/31 (16.1%); ARID2, BIRC3, and FBXW7 mutations each in 4/31 (12.9%); and TP53, MTAP, and FAT1 each in 3/31 (9.7%) HG AciCC. PD-L1-positive labeling was more common in HG AciCC (9/17, 52.9% vs 3/16, 18.9%, P = .071). CDKN2A/B mutations in AciCC represent a molecular marker of HG histology and disease progression, providing a rationale for further studies to determine their prognostic and therapeutic significance in this salivary gland cancer. AciCC with ATM mutations may be amenable to targeted therapy. Immunotherapy can be considered to be a treatment option for a subset of patients with AciCC.

Clinical sequencing of soft tissue and bone sarcomas delineates diverse genomic landscapes and potential therapeutic targets.

The genetic, biologic, and clinical heterogeneity of sarcomas poses a challenge for the identification of therapeutic targets, clinical research, and advancing patient care. Because there are > 100 sarcoma subtypes, in-depth genetic studies have focused on one or a few subtypes. Herein, we report a comparative genetic analysis of 2,138 sarcomas representing 45 pathological entities. This cohort is prospectively analyzed using targeted sequencing to characterize subtype-specific somatic alterations in targetable pathways, rates of whole genome doubling, mutational signatures, and subtype-agnostic genomic clusters. The most common alterations are in cell cycle control and TP53, receptor tyrosine kinases/PI3K/RAS, and epigenetic regulators. Subtype-specific associations include TERT amplification in intimal sarcoma and SWI/SNF alterations in uterine adenosarcoma. Tumor mutational burden, while low compared to other cancers, varies between and within subtypes. This resource will improve sarcoma models, motivate studies of subtype-specific alterations, and inform investigations of genetic factors and their correlations with treatment response.

Primary cutaneous SMARCB1-deficient carcinoma.

BACKGROUND: SMARCB1-deficient malignancies can arise in various sites. We describe a novel primary SMARCB1-deficient carcinoma of skin (SDCS) and characterize SMARCB1 mutations in non-melanoma skin cancers (NMSC). METHODS: Cases underwent immunophenotyping and targeted exome sequencing (MSK-IMPACT) assay interrogating somatic mutations in 468 cancer-related genes. The MSK-IMPACT database from 2014 to 2020 encompassing 55, 000 cases was searched for NMSC with SMARCB1 mutations. RESULTS: SDCS arose on the scalp of an 18-year-old woman showing homozygous SMARCB1 deletion with a LATS2 G963E variant. Another case arose on the temple of a 76-year-old man harboring a SMARCB1 W206* mutation associated with loss of heterozygosity (LOH), 59 concurrent mutations, and a UV mutation signature (UV-MS). Both tumors exhibited INI1 loss, positive CK5/6, p40, p63, and claudin-4 with negative CD34. Of 378 NMSC cases, including 370 carcinomas, 7 SMARCB1-mutated tumors were identified: 3 squamous cell, 3 Merkel cell, and one basal cell carcinoma. Six showed UV-MS. Five INI1-interrogated cases retained protein expression suggesting they were SMARCB1-proficient. CONCLUSIONS: SDCS can be clinically aggressive, harbor SMARCB1 homozygous deletions or truncating SMARCB1 mutations associated with LOH, and can occur with or without UV-MS. Overall, SMARCB1 mutations in NMSC are rare with most being of undetermined significance and associated with retained INI1 and UV-MS.

The association between tumor mutational burden and prognosis is dependent on treatment context.

In multiple cancer types, high tumor mutational burden (TMB) is associated with longer survival after treatment with immune checkpoint inhibitors (ICIs). The association of TMB with survival outside of the immunotherapy context is poorly understood. We analyzed 10,233 patients (80% non-ICI-treated, 20% ICI-treated) with 17 cancer types before/without ICI treatment or after ICI treatment. In non-ICI-treated patients, higher TMB (higher percentile within cancer type) was not associated with better prognosis; in fact, in many cancer types, higher TMB was associated with poorer survival, in contrast to ICI-treated patients in whom higher TMB was associated with longer survival.

A Pan-Cancer Study of Somatic TERT Promoter Mutations and Amplification in 30,773 Tumors Profiled by Clinical Genomic Sequencing.

TERT gene promoter mutations are known in multiple cancer types. Other TERT alterations remain poorly characterized. Sequencing data from 30,773 tumors analyzed by a hybridization capture next-generation sequencing assay (Memorial Sloan Kettering Cancer Center Integrated Mutation Profiling of Actionable Cancer Targets) were analyzed for the presence of TERT alterations. Promoter rearrangements (500 bases upstream of the transcriptional start site), hypermethylation (n = 57), and gene expression (n = 155) were evaluated for a subset of cases. Mutually exclusive and recurrent promoter mutations were identified at three hot spots upstream of the transcriptional start site in 11.3% of cases (-124: 74%; -146: 24%; and -138: <2%). Mutually exclusive amplification events were identified in another 2.3% of cases, whereas mutually exclusive rearrangements proximal to the TERT gene were seen in 24 cases. The highest incidence of TERT promoter mutations was seen in cutaneous melanoma (82%), whereas amplification events significantly outnumbered promoter mutations in well-differentiated/dedifferentiated liposarcoma (14.1% versus 2.4%) and adrenocortical carcinoma (13.6% versus 4.5%). Gene expression analysis suggests that the highest levels of gene expression are seen in cases with amplifications and rearrangements. Hypermethylation events upstream of the TERT coding sequence were not mutually exclusive with known pathogenic alterations. Studies aimed at defining the prevalence and prognostic impact of TERT alterations should incorporate other pathogenic TERT alterations as these may impact telomerase function.

E-cadherin immunohistochemical expression in invasive lobular carcinoma of the breast: correlation with morphology and CDH1 somatic alterations.

E-cadherin (ECAD) immunohistochemical (IHC) expression is lost in ∼90% of invasive lobular carcinomas (ILCs) owing to genomic alterations of CDH1. We examined morphologic features and ECAD IHC expression in invasive breast carcinomas (BCs) with known CDH1 alterations. Between January 2014 and May 2018, 202 cases of BC with a CDH1 somatic alteration were identified. ECAD expression was lost in 77% (155/202) of cases and was retained in 23% (47/202) cases. Most (90%, 139/155) ECAD-negative cases were morphologically classified as ILC, while the remaining (10%, 16/155) were invasive mammary carcinoma with mixed ductal and lobular features (IMC). Of 47 cases with ECAD staining, 62% (29/47) were classified as ILC, 23% (11/47) were classified as IMC, and 15% (7/47) were classified as invasive ductal carcinoma (IDC). Of note, 51% (24/47) of ECAD-positive cases were initially diagnosed as IDC or IMC based on ECAD expression alone. For ECAD-negative BCs, 98% (152/155) of CDH1 alterations were truncating, and 2% (3/155) were variants of unknown significance (VUS). Truncating CDH1 alterations were identified in the majority of ECAD-positive BCs (72%, 34/47); however, VUS-type CDH1 alterations were more prevalent (28%, 13/47) in ECAD-positive BCs than in ECAD-negative BCs. Although 90% of ECAD-negative tumors were compatible with ILC in this study, 17% (29/168) of ILC cases were ECAD positive. In addition, CDH1 truncating alterations were seen in ECAD-positive ILC, supporting the notion of aberrant ECAD staining. Therefore, ECAD IHC expression must be interpreted in conjunction with morphology, and BC with classic histologic features of ILC should not be reclassified as IDC/IMC based solely on the status of ECAD IHC expression.

Retained mismatch repair protein expression occurs in approximately 6% of microsatellite instability-high cancers and is associated with missense mutations in mismatch repair genes.

Immunohistochemistry for mismatch repair protein expression is widely used as a surrogate for microsatellite instability status-an important signature for immunotherapy and germline testing. There are no systematic analyses examining the sensitivity of immunohistochemistry for microsatellite instability-high status. Mismatch repair immunohistochemistry and microsatellite instability testing were performed routinely as clinically validated assays. We classified germline/somatic mutation types as truncating (nonsense, frameshift, and in/del) versus missense and predicted pathogenicity of the latter. Discordant cases were compared with concordant groups: microsatellite instability-high/mismatch repair-deficient for mutation comparison and microsatellite stable/mismatch repair-proficient for immunohistochemical comparison. 32 of 443 (7%) microsatellite instability-high cases had immunohistochemistry. Four additional microsatellite instability-high research cases had discordant immunohistochemistry. Of 36 microsatellite instability-high cases with discordant immunohistochemistry, 30 were mismatch repair-proficient, while six (five MLH1 and one MSH2) retained expression of the defective mismatch repair protein and lost its partner. In microsatellite instability-high tumors with discordant immunohistochemistry, we observed an enrichment in deleterious missense mutations over truncating mutations, with 69% (25/36) of cases having pathogenic germline or somatic missense mutations, as opposed to only 19% (7/36) in a matched microsatellite instability-high group with concordant immunohistochemistry (p = 0.0007).  In microsatellite instability-high cases with discordant immunohistochemistry and MLH1 or PMS2 abnormalities, less cells showed expression (p = 0.015 and p = 0.00095, respectively) compared with microsatellite stable/mismatch repair-proficient cases. Tumor mutation burden, MSIsensor score, and truncating mismatch repair gene mutations were similar between microsatellite instability-high cases with concordant versus discordant immunohistochemical expression. Approximately 6% of microsatellite instability-high cases have retained mismatch repair protein expression and would be missed by immunohistochemistry-based testing, hindering patient access to immunotherapy. Another 1% of microsatellite instability-high cases show isolated loss of the defective gene's dimerization partner, which may lead to germline testing of the wrong gene. These cases are enriched for pathogenic mismatch repair missense mutations.

Cancer Susceptibility Mutations in Patients With Urothelial Malignancies.

PURPOSE: Urothelial cancers (UCs) have a substantial hereditary component, but, other than their association with Lynch syndrome, the contribution of genetic risk factors to UC pathogenesis has not been systematically defined. We sought to determine the prevalence of pathogenic/likely pathogenic (P/LP) germline variants in patients with UC and identify associated clinical factors. PATIENTS AND METHODS: Overall, 586 patients with UC underwent prospective, matched tumor-normal DNA sequencing. Seventy-seven genes associated with cancer predisposition were analyzed; allele frequencies were compared with publicly available database. RESULTS: P/LP germline variants were identified in 80 (14%) of 586 individuals with UC. The most common P/LP variants in high- or moderate-penetrance genes were BRCA2 (n = 9; 1.5%), MSH2 (n = 8; 1.4%), BRCA1 (n = 8; 1.4%), CHEK2 (n = 6; 1.0%), ERCC3 (n = 4; 0.7%), and NBN and RAD50 (n = 3; 0.5% each). Sixty-six patients (83%) had germline P/LP variants in DNA-damage repair (DDR) genes, of which 28 (42%) had biallelic inactivation. Patients with P/LP variants were more commonly diagnosed at an early age (22% v 6% in those without variants; P = .01). BRCA2 and MSH2 were significantly associated with an increased risk for UC (odds ratio, 3.7 [P = .004] and 4.6 [P = .001], respectively). Current clinical guidelines for referral for genetic testing failed to identify 6 (26%) patients with high-penetrance variants. CONCLUSION: Clinically significant P/LP germline variants in DDR genes frequently are present in patients with advanced UC. The presence of DDR germline variants could guide cancer screening for patients and their families and serve as predictive biomarkers of response to targeted or immunotherapies. Family history-based criteria to identify patients with hereditary UC susceptibility are insensitive. Broader germline testing in UC, particularly in those of young ages, should be considered.

Distinctive mechanisms underlie the loss of SMARCB1 protein expression in renal medullary carcinoma: morphologic and molecular analysis of 20 cases.

Renal medullary carcinoma is a rare but highly aggressive type of renal cancer occurring in patients with sickle cell trait or rarely with other hemoglobinopathies. Loss of SMARCB1 protein expression, a core subunit of the switch/sucrose nonfermentable (SWI/SNF) chromatin remodeling complex, has emerged as a key diagnostic feature of these tumors. However, the molecular mechanism underlying this loss remains unclear. We retrospectively identified 20 patients diagnosed with renal medullary carcinoma at two institutions from 1996 to 2017. All patients were confirmed to have sickle cell trait, and all tumors exhibited a loss of SMARCB1 protein expression by immunohistochemistry. The status of SMARCB1 locus was examined by fluorescence in situ hybridization (FISH) using 3-color probes, and somatic alterations were detected by targeted next-generation sequencing platforms. FISH analysis of all 20 cases revealed 11 (55%) with concurrent hemizygous loss and translocation of SMARCB1, 6 (30%) with homozygous loss of SMARCB1, and 3 (15%) without structural or copy number alterations of SMARCB1 despite protein loss. Targeted sequencing revealed a pathogenic somatic mutation of SMARCB1 in one of these 3 cases that were negative by FISH. Tumors in the 3 subsets with different FISH findings largely exhibited similar clinicopathologic features, however, homozygous SMARCB1 deletion was found to show a significant association with the solid growth pattern, whereas tumors dominated by reticular/cribriform growth were enriched for SMARCB1 translocation. Taken together, we demonstrate that different molecular mechanisms underlie the loss of SMARCB1 expression in renal medullary carcinoma. Biallelic inactivation of SMARCB1 occurs in a large majority of cases either via concurrent hemizygous loss and translocation disrupting SMARCB1 or by homozygous loss.

Majority of B2M-Mutant and -Deficient Colorectal Carcinomas Achieve Clinical Benefit From Immune Checkpoint Inhibitor Therapy and Are Microsatellite Instability-High.

PURPOSE: Microsatellite instability-high (MSI-H) colorectal carcinomas (CRCs) show high rates of response to immune checkpoint inhibitors (IOs). B2M mutations and protein loss have been proposed as causes of resistance to IOs, yet they are enriched in MSI-H CRC. We aimed to characterize B2M-mutant, IO-naive CRC. PATIENTS AND METHODS: All CRCs with results for Memorial Sloan Kettering Integrated Mutation Profiling of Actionable Cancer Targets, a next-generation sequencing assay that interrogates > 400 genes for mutations as well as MSI status, were surveyed for B2M mutations. All B2M-mutant CRCs were assessed for expression of B2M, major histocompatibility complex class I, and programmed death-1 ligand (PD-L1) via immunohistochemistry and average CD3+ and CD8+ tumor-infiltrating lymphocyte counts against a control group of MSI-H B2M wild-type CRCs. RESULTS: Fifty-nine (3.4%) of 1,751 patients with CRC harbored B2M mutations, with 84% (77 of 92) of the mutations predicted to be truncating. B2M mutations were significantly enriched in MSI-H CRCs, with 44 (24%) of 182 MSI-H CRCs harboring B2M mutations (P < .001). Thirty-two of 44 B2M-mutant CRCs with available material (73%) had complete loss of B2M expression, whereas all 26 CRCs with wild-type B2M retained expression (P < .001). B2M mutation status was not associated with major histocompatibility complex class I expression, KRAS or BRAF mutation, tumor-infiltrating lymphocyte level, or PD-L1 expression after adjustment for MSI status. Of 13 patients with B2M-mutant CRC who received programmed death-1 or PD-L1 IOs, 11 (85%) achieved clinical benefit, defined as stable disease or partial response using Response Evaluation Criteria in Solid Tumors criteria. CONCLUSION: B2M mutations occur in approximately 24% of MSI-H CRCs and are usually associated with loss of B2M expression. Most patients with B2M-mutant MSI-H CRC with loss of protein expression obtain clinical benefit from IOs.

Comparison of De Novo Transcriptome Assemblers and k-mer Strategies Using the Killifish, Fundulus heteroclitus.

BACKGROUND: De novo assembly of non-model organism's transcriptomes has recently been on the rise in concert with the number of de novo transcriptome assembly software programs. There is a knowledge gap as to what assembler software or k-mer strategy is best for construction of an optimal de novo assembly. Additionally, there is a lack of consensus on which evaluation metrics should be used to assess the quality of de novo transcriptome assemblies. RESULT: Six different assembly strategies were evaluated from four different assemblers. The Trinity assembly was used in its default 25 single k-mer value while Bridger, Oases, and SOAPdenovo-Trans were performed with multiple k-mer strategies. Bridger, Oases, and SOAPdenovo-Trans used a small multiple k-mer (SMK) strategy consisting of the k-mer lengths of 21, 25, 27, 29, 31, and 33. Additionally, Oases and SOAPdenovo-Trans were performed using a large multiple k-mer (LMK) strategy consisting of k-mer lengths of 25, 35, 45, 55, 65, 75, and 85. Eleven metrics were used to evaluate each assembly strategy including three genome related evaluation metrics (contig number, N50 length, Contigs >1 kb, reads) and eight transcriptome evaluation metrics (mapped back to transcripts (RMBT), number of full length transcripts, number of open reading frames, Detonate RSEM-EVAL score, and percent alignment to the southern platyfish, Amazon molly, BUSCO and CEGMA databases). The assembly strategy that performed the best, that is it was within the top three of each evaluation metric, was the Bridger assembly (10 of 11) followed by the Oases SMK assembly (8 of 11), the Oases LMK assembly (6 of 11), the Trinity assembly (4 of 11), the SOAP LMK assembly (4 of 11), and the SOAP SMK assembly (3 of 11). CONCLUSION: This study provides an in-depth multi k-mer strategy investigation concluding that the assembler itself had a greater impact than k-mer size regardless of the strategy employed. Additionally, the comprehensive performance transcriptome evaluation metrics utilized in this study identified the need for choosing metrics centered on user defined research goals. Based on the evaluation metrics performed, the Bridger assembly was able to construct the best assembly of the testis transcriptome in Fundulus heteroclitus.