Linked‐read sequencing enables haplotype‐resolved resequencing at population scale

The feasibility to sequence entire genomes of virtually any organism provides unprecedented insights into the evolutionary history of populations and species. Nevertheless, many population genomic inferences – including the quantification and dating of admixture, introgression and demographic events, and inference of selective sweeps – are still limited by the lack of high‐quality haplotype information. The newest generation of sequencing technology now promises significant progress. To establish the feasibility of haplotype‐resolved genome resequencing at population scale, we investigated properties of linked‐read sequencing data of songbirds of the genus Oenanthe across a range of sequencing depths. Our results based on the comparison of downsampled (25×, 20×, 15×, 10×, 7×, and 5×) with high‐coverage data (46–68×) of seven bird genomes mapped to a reference suggest that phasing contiguities and accuracies adequate for most population genomic analyses can be reached already with moderate sequencing effort. At 15× coverage, phased haplotypes span about 90% of the genome assembly, with 50% and 90% of phased sequences located in phase blocks longer than 1.25–4.6 Mb (N50) and 0.27–0.72 Mb (N90). Phasing accuracy reaches beyond 99% starting from 15× coverage. Higher coverages yielded higher contiguities (up to about 7 Mb/1 Mb [N50/N90] at 25× coverage), but only marginally improved phasing accuracy. Phase block contiguity improved with input DNA molecule length; thus, higher‐quality DNA may help keeping sequencing costs at bay. In conclusion, even for organisms with gigabase‐sized genomes like birds, linked‐read sequencing at moderate depth opens an affordable avenue towards haplotype‐resolved genome resequencing at population scale.     

Evolutionary processes of melanomas from giant congenital melanocytic nevi

Melanoma can develop in a congenital melanocytic nevus (CMN). In fact, a large CMN is associated with a high risk of developing melanoma. Although melanomas arising from CMNs are thought to have a pathogenesis distinct from conventional melanomas, no studies have been conducted on the evolution or tumor heterogeneity of CMN melanomas. We applied multi‐region whole‐exome sequencing to investigate the clonal nature of driver events and evolutionary processes in CMNs and melanomas arising from CMNs. In two patients, we observed an independent subclonal evolution in cancerized fields of CMNs and chromosome 8q amplification in both melanomas arising from CMNs. The amplification of MYC, located in chromosome 8q, was correlated with the percentage of tumor cells expressing high levels of MYC protein detected in melanoma cells by immunohistochemistry. Our analysis suggests that each CMN cell may evolve sporadically and that amplification of MYC might be a key event for melanoma development in CMNs.     

Life cycle assessment of forest biomass energy feedstock in the Northeast United States

Life cycle assessment (LCA) was combined with primary data from nine forest harvesting operations in New York, Maine, Massachusetts, and Vermont, from 2013 to 2019 where forest biomass (FB) for bioenergy was one of several products. The objective was to conduct a data‐driven study of greenhouse gas emissions associated with FB feedstock harvesting operations in the Northeast United States. Deterministic and stochastic LCA models were built to simulate the current FB bioenergy feedstock supply chain in the Northeast US with a cradle‐to‐gate scope (forest harvest through roadside loading) and a functional unit of 1.0 Mg of green FB feedstock at a 50% moisture content. Baseline LCA, sensitivity analysis, and uncertainty analyses were conducted for three different FB feedstock types—dirty chips, clean chips, and grindings—enabling an empirically driven investigation of differences between feedstock types, individual harvesting process contributions, and literature comparisons. The baseline LCA average impacts were lower for grindings (4.57 kg CO_2eq/Mg) and dirty chips (7.16 kg CO_2eq/Mg) than for clean chips (23.99 kg CO_2eq/Mg) under economic allocation, but impacts were of similar magnitude under mass allocation, ranging from 24.42 to 27.89 kg CO_2eq/Mg. Uncertainty analysis showed a wider range of probable results under mass allocation compared to economic allocation. Sensitivity analysis revealed the impact of variations in the production masses and total economic values of primary products of forest harvests on the LCA results due to allocation of supply chain emissions. The high variability in fuel use between logging contractors also had a distinct influence on LCA results. The results of this study can aid decision‐makers in energy policy and guide emissions reductions efforts while informing future LCAs that expand the system boundary to regional FB energy pathways, including electricity generation, transportation fuels, pellets for heat, and combined heat and power.     

Predicting soil carbon changes in switchgrass grown on marginal lands under climate change and adaptation strategies

The United States Great Lakes Region (USGLR) is a critical geographic area for future bioenergy production. Switchgrass (Panicum virgatum) is widely considered a carbon (C)‐neutral or C‐negative bioenergy production system, but projected increases in air temperature and precipitation due to climate change might substantially alter soil organic C (SOC) dynamics and storage in soils. This study examined long‐term SOC changes in switchgrass grown on marginal land in the USGLR under current and projected climate, predicted using a process‐based model (Systems Approach to Land‐Use Sustainability) extensively calibrated with a wealth of plant and soil measurements at nine experimental sites. Simulations indicate that these soils are likely a net C sink under switchgrass (average gain 0.87 Mg C ha^?1 year^?1), although substantial variation in the rate of SOC accumulation was predicted (range: 0.2–1.3 Mg C ha^?1 year^?1). Principal component analysis revealed that the predicted intersite variability in SOC sequestration was related in part to differences in climatic characteristics, and to a lesser extent, to heterogeneous soils. Although climate change impacts on switchgrass plant growth were predicted to be small (4%–6% decrease on average), the increased soil respiration was predicted to partially negate SOC accumulations down to 70% below historical rates in the most extreme scenarios. Increasing N fertilizer rate and decreasing harvest intensity both had modest SOC sequestration benefits under projected climate, whereas introducing genotypes better adapted to the longer growing seasons was a much more effective strategy. Best‐performing adaptation scenarios were able to offset >60% of the climate change impacts, leading to SOC sequestration 0.7 Mg C ha^?1 year^?1 under projected climate. On average, this was 0.3 Mg C ha^?1 year^?1 more C sequestered than the no adaptation baseline. These findings provide crucial knowledge needed to guide policy and operational management for maximizing SOC sequestration of future bioenergy production on marginal lands in the USGLR.     

Melanocortin‐1 receptor (MC1R) genotypes do not correlate with size in two cohorts of medium‐to‐giant congenital melanocytic nevi

Congenital melanocytic nevi (CMN) are cutaneous malformations whose prevalence is inversely correlated with projected adult size. CMN are caused by somatic mutations, but epidemiological studies suggest that germline genetic factors may influence CMN development. In CMN patients from the U.K., genetic variants in MC1R, such as p.V92M and loss‐of‐function variants, have been previously associated with larger CMN. We analyzed the association of MC1R variants with CMN characteristics in two distinct cohorts of medium‐to‐giant CMN patients from Spain (N = 113) and from France, Norway, Canada, and the United States (N = 53), similar at the clinical and phenotypical level except for the number of nevi per patient. We found that the p.V92M or loss‐of‐function MC1R variants either alone or in combination did not correlate with CMN size, in contrast to the U.K. CMN patients. An additional case–control analysis with 259 unaffected Spanish individuals showed a higher frequency of MC1R compound heterozygous or homozygous variant genotypes in Spanish CMN patients compared to the control population (15.9% vs. 9.3%; p = .075). Altogether, this study suggests that MC1R variants are not associated with CMN size in these non‐UK cohorts. Additional studies are required to define the potential role of MC1R as a risk factor in CMN development.     

Induced pluripotent stem cell‐derived myeloid cells expressing OX40 ligand amplify antigen‐specific T cells in advanced melanoma

Immune checkpoint inhibitors improved the survival rate of patients with unresectable melanoma. However, some patients do not respond, and variable immune‐related adverse events have been reported. Therefore, more effective and antigen‐specific immune therapies are urgently needed. We previously reported the efficacy of an immune cell therapy with immortalized myeloid cells derived from induced pluripotent stem cells (iPS‐ML). In this study, we generated OX40L‐overexpressing iPS‐ML (iPS‐ML‐Zsgreen‐OX40L) and investigated their characteristics and in vivo efficacy against mouse melanoma. We found that iPS‐ML‐Zsgreen‐OX40L suppressed the progression of B16‐BL6 melanoma, and prolonged survival of mice with ovalbumin (OVA)‐expressing B16 melanoma (MO4). The number of antigen‐specific CD8⁺ T cells was higher in spleen cells treated with OVA peptide‐pulsed iPS‐ML‐Zsgreen‐OX40L than in those without OX40L. The OVA peptide‐pulsed iPS‐ML‐Zsgreen‐OX40L significantly increased the number of tumor‐infiltrating T lymphocytes (TILs) in MO4 tumor. Flow cytometry showed decreased regulatory T cells but increased effector and effector memory T cells among the TILs. Although we plan to use allogeneic iPS‐ML in the clinical applications, iPS‐ML showed the tumorgenicity in the syngeneic mice model. Incorporating the suicide gene is necessary to ensure the safety in the future study. Collectively, these results indicate that iPS‐ML‐Zsgreen‐OX40L therapy might be a new method for antigen‐specific cancer immunotherapy.     

Distinct genomic traits of acral and mucosal melanomas revealed by targeted mutational profiling

The incidence of melanoma is rising globally including China. Comparing to Caucasians, the incidence of non‐cutaneous melanomas is significantly higher in Chinese. Herein, we performed genomic profiling of 89 Chinese surgically resected primary melanomas, including acral (n = 54), cutaneous (n = 22), and mucosal (n = 13), by hybrid capture‐based next‐generation sequencing. We show that mucosal melanomas tended to harbor more pathogenic mutations than other types of melanoma, though the biological significance of this finding remains uncertain. Chromosomal arm‐level alterations including 6q, 9p, and 10p/q loss were highly recurrent in all subtypes, but mucosal melanoma was significantly associated with increased genomic instability. Importantly, 7p gain significantly correlated with unfavorable clinical outcomes in non‐cutaneous melanomas, representing an intriguing prognostic biomarker of those subtypes. Furthermore, focal amplification of 4q12 (KIT, KDR, and PDGFRα) and RAD51 deletion were more abundant in mucosal melanoma, while NOTCH2 amplification was enriched in acral melanoma. Additionally, cutaneous melanomas had higher mutation load than acral melanomas, while mucosal melanomas did not differ from other subtypes in mutation burden. Together, our data revealed important features of acral and mucosal melanomas in Chinese including distinctive driver mutation pattern and increased genomic instability. These findings highlight the possibilities of combination therapies in the clinical management of melanoma.     

Do innate killing mechanisms activated by inflammasomes have a role in treating melanoma?

Melanoma, as for many other cancers, undergoes a selection process during progression that limits many innate and adaptive tumor control mechanisms. Immunotherapy with immune checkpoint blockade overcomes one of the escape mechanisms but if the tumor is not eliminated other escape mechanisms evolve that require new approaches for tumor control. Some of the innate mechanisms that have evolved against infections with microorganisms and viruses are proving to be active against cancer cells but require better understanding of how they are activated and what inhibitory mechanisms may need to be targeted. This is particularly so for inflammasomes which have evolved against many different organisms and which recruit a number of cytotoxic mechanisms that remain poorly understood. Equally important is understanding of where these mechanisms will fit into existing treatment strategies and whether existing strategies already involve the innate killing mechanisms.     

Chromosome‐level genome assembly of the greenhouse whitefly (Trialeurodes vaporariorum Westwood)

The greenhouse whitefly, Trialeurodes vaporariorum Westwood, is an agricultural pest of global importance. Here we report a 787‐Mb high‐quality draft genome sequence of T. vaporariorum assembled from PacBio long reads and Hi‐C chromatin interaction maps, which has scaffold and contig N50 lengths of 70 Mb and 500 kb, respectively, and contains 18,275 protein‐coding genes. About 98.8% of the assembled contigs were placed onto the 11 T. vaporariorum chromosomes. Comparative genomic analysis reveals significantly expanded gene families such as aspartyl proteases in T. vaporariorum compared to Bemisia tabaci Mediterranean (MED) and Middle East‐Asia Minor 1 (MEAM1). Furthermore, the cytochrome CYP6 subfamily shows significant expansion in T. vaporariorum and several genes in this subfamily display developmental stage‐specific expression patterns. The high‐quality T. vaporariorum genome provides a valuable resource for research in a broad range of areas such as fundamental molecular ecology, insect–plant/insect–microorganism or virus interactions and pest resistance management.