单细胞测序技术研究进展

单细胞测序技术是能够在单个细胞的水平上,对基因组进行高通量测序分析的一项新技术。与传统高通量测序相比,单细胞测序不仅能够分析相同表型细胞的遗传异质性,还能获取难以培养微生物的遗传信息,具有广阔的应用前景。单细胞测序技术的流程主要包括单细胞分离、细胞溶解与基因组DNA获取、全基因组扩增、测序与数据分析4个方面,以该技术流程为主线,分析了现有技术的优缺点,并对最新的改进办法进行了综述。     

Fine Breakpoint Mapping by Genome Sequencing Reveals the First Large X Inversion Disrupting the NHS Gene in a Patient with Syndromic Cataracts

Inversions are structural variants that are generally balanced. However, they could lead to gene disruptions or have positional effects leading to diseases. Mutations in the NHS gene cause Nance-Horan syndrome, an X-linked disorder characterised by congenital cataracts and dental anomalies. Here, we aimed to characterise a balanced pericentric inversion X(p22q27), maternally inherited, in a child with syndromic bilateral cataracts by breakpoint mapping using whole-genome sequencing (WGS). 30× Illumina paired-end WGS was performed in the proband, and breakpoints were confirmed by Sanger sequencing. EdU assays and FISH analysis were used to assess skewed X-inactivation patterns. RNA expression of involved genes in the breakpoint boundaries was evaluated by droplet-digital PCR. We defined the breakpoint position of the inversion at Xp22.13, with a 15 bp deletion, disrupting the unusually large intron 1 of the canonical NHS isoform, and also perturbing topologically-associated domains (TADs). Moreover, a microhomology region of 5 bp was found on both sides. RNA analysis confirmed null and reduced NHS expression in the proband and his unaffected mother, respectively. In conclusion, we report the first chromosomal inversion disrupting NHS, fine-mapped by WGS. Our data expand the clinical spectrum and the pathogenic mechanisms underlying the NHS defects.     

Data Interoperability of Whole Exome Sequencing (WES) Based Mutational Burden Estimates from Different Laboratories

Immune checkpoint inhibitors, which unleash a patient’s own T cells to kill tumors, are revolutionizing cancer treatment. Several independent studies suggest that higher non-synonymous mutational burden assessed by whole exome sequencing (WES) in tumors is associated with improved objective response, durable clinical benefit, and progression-free survival in immune checkpoint inhibitors treatment. Next-generation sequencing (NGS) is a promising technology being used in the clinic to direct patient treatment. Cancer genome WES poses a unique challenge due to tumor heterogeneity and sequencing artifacts introduced by formalin-fixed, paraffin-embedded (FFPE) tissue. In order to evaluate the data interoperability of WES data from different sources to survey tumor mutational landscape, we compared WES data of several tumor/normal matched samples from five commercial vendors. A large data discrepancy was observed from vendors’ self-reported data. Independent data analysis from vendors’ raw NGS data shows that whole exome sequencing data from qualified vendors can be combined and analyzed uniformly to derive comparable quantitative estimates of tumor mutational burden.     

Pre-Analytical Considerations for Successful Next-Generation Sequencing (NGS): Challenges and Opportunities for Formalin-Fixed and Paraffin-Embedded Tumor Tissue (FFPE) Samples

In cancer drug discovery, it is important to investigate the genetic determinants of response or resistance to cancer therapy as well as factors that contribute to adverse events in the course of clinical trials. Despite the emergence of new technologies and the ability to measure more diverse analytes (e.g., circulating tumor cell (CTC), circulating tumor DNA (ctDNA), etc.), tumor tissue is still the most common and reliable source for biomarker investigation. Because of its worldwide use and ability to preserve samples for many decades at ambient temperature, formalin-fixed, paraffin-embedded tumor tissue (FFPE) is likely to be the preferred choice for tissue preservation in clinical practice for the foreseeable future. Multiple analyses are routinely performed on the same FFPE samples (such as Immunohistochemistry (IHC), in situ hybridization, RNAseq, DNAseq, TILseq, Methyl-Seq, etc.). Thus, specimen prioritization and optimization of the isolation of analytes is critical to ensure successful completion of each assay. FFPE is notorious for producing suboptimal DNA quality and low DNA yield. However, commercial vendors tend to request higher DNA sample mass than what is actually required for downstream assays, which restricts the breadth of biomarker work that can be performed. We evaluated multiple genomics service laboratories to assess the current state of NGS pre-analytical processing of FFPE. Significant differences in pre-analytical capabilities were observed. Key aspects are highlighted and recommendations are made to improve the current practice in translational research.     

Comparative Mitogenomics of the Genus Odontobutis (Perciformes: Gobioidei: Odontobutidae) Revealed Conserved Gene Rearrangement and High Sequence Variations

To understand the molecular evolution of mitochondrial genomes (mitogenomes) in the genus Odontobutis, the mitogenome of Odontobutis yaluensis was sequenced and compared with those of another four Odontobutis species. Our results displayed similar mitogenome features among species in genome organization, base composition, codon usage, and gene rearrangement. The identical gene rearrangement of trnS-trnL-trnH tRNA cluster observed in mitogenomes of these five closely related freshwater sleepers suggests that this unique gene order is conserved within Odontobutis. Additionally, the present gene order and the positions of associated intergenic spacers of these Odontobutis mitogenomes indicate that this unusual gene rearrangement results from tandem duplication and random loss of large-scale gene regions. Moreover, these mitogenomes exhibit a high level of sequence variation, mainly due to the differences of corresponding intergenic sequences in gene rearrangement regions and the heterogeneity of tandem repeats in the control regions. Phylogenetic analyses support Odontobutis species with shared gene rearrangement forming a monophyletic group, and the interspecific phylogenetic relationships are associated with structural differences among their mitogenomes. The present study contributes to understanding the evolutionary patterns of Odontobutidae species.     

Structural and Thermal Characterization of A Novel High Nitrogen Energetic Material: (NH4)2DNMT

1,4‐Dihydro‐5H‐(dinitromethylene)‐tetrazole ammonium salt ((NH₄)₂DNMT), a high nitrogen energetic compound, was synthesized and structurally characterized by single‐crystal X‐ray diffraction. The thermal behavior of (NH₄)₂DNMT was studied with DSC and TG‐DTG methods. The kinetic equation of the thermal decomposition reaction is: dα/dT=10^13.17/3β(1−α)⁻² exp(−1.388×10⁵/RT). The critical temperature of thermal explosion is 182.7 °C. The specific heat capacity of (NH₄)₂DNMT was determined and the molar heat capacity is 301 J mol⁻¹ K⁻¹ at 298.15 K. The adiabatic time‐to‐explosion of (NH₄)₂DNMT was calculated to be 277 s. The detonation velocity and detonation pressure were also estimated. All results showed that (NH₄)₂DNMT presents good performance.     

A Novel Mitochondrial Genome Fragmentation Pattern in the Buffalo Louse Haematopinus tuberculatus (Psocodea: Haematopinidae)

Sucking lice are obligate ectoparasites of mammalian hosts, causing serious public health problems and economic losses worldwide. It is well known that sucking lice have fragmented mitochondrial (mt) genomes, but many remain undetermined. To better understand patterns of mt genome fragmentation in the sucking lice, we sequenced the mt genome of the buffalo louse Haematopinus tuberculatus using next-generation sequencing (NGS). The mt genome of H. tuberculatus has ten circular minichromosomes containing a total of 37 genes. Each minichromosome is 2.9–5.0 kb long and carries one to eight genes plus one large non-coding region. The number of mt minichromosomes of H. tuberculatus (ten) is different from those of congeneric species (horse louse H. asini, domestic pig louse H. suis and wild pig louse H. apri) and other sucking lice. Two events (gene translocation and merger of mt minichromosome) are observed in Haematopinus. Compared to other studies, our phylogeny generated from mt genome datasets showed a different topology, suggesting that inclusion of data other than mt genomes would be required to resolve phylogeny of sucking lice. To our knowledge, this is the first report of a ten mt minichromosomes genome in sucking lice, which opens a new outlook into unexplored mt genome fragmentation patterns in sucking lice.     

The Presentation of Two Unrelated Clinical Cases from the Republic of North Ossetia-Alania with the Same Previously Undescribed Variant in the COL6A2 Gene

Here, we described three affected boys from two unrelated families of Ossetian-Digor origin from the Republic of North Ossetia-Alania who were admitted to the Research Centre for Medical Genetics with unspecified muscular dystrophy. High-throughput sequencing was performed and revealed two novel frameshift variants in the COL6A2 gene (NM_001849.3) in a heterozygous state each in both cases: c.508_535delinsCTGTGG and c.1659_1660del (case 1) and c.1689del and c.1659_1660del (case 2). In two cases, the same nucleotide variant in the COL6A2 gene (c.1659_1660del) was observed. We have suggested that the variant c.1659_1660del may be common in the Ossetian-Digor population because two analyzed families have the same ancestry from the same subethnic group of Ossetians). The screening for an asymptomatic carriage of the nucleotide variant c.1659_1660del in 54 healthy donors from Ossetian-Digor population revealed that the estimated carrier frequency is 0.0093 (CI: 0.0002–0.0505), which is high for healthy carriers of the pathogenic variant. Molecular genetic, anamnestic data and clinical examination results allowed us to diagnose Ullrich muscular dystrophy in those affected boys. Genetic heterogeneity and phenotypic diversity of muscular dystrophies complicate diagnosis. It is important to make a differential diagnosis of such conditions and use HTS methods to determine the most accurate diagnosis.     

Single-Cell RNA Sequencing for the Detection of Clonotypic V(D)J Rearrangements in Multiple Myeloma

Multiple myeloma (MM) has a highly heterogeneous genetic background, which complicates its molecular tracking over time. Nevertheless, each MM patient’s malignant plasma cells (PCs) share unique V(D)J rearranged sequences at immunoglobulin loci, which represent ideal disease biomarkers. Because the tumor-specific V(D)J sequence is highly expressed in bulk RNA in MM patients, we wondered whether it can be identified by single-cell RNA sequencing (scRNA-seq). To this end we analyzed CD138⁺ cells purified from bone marrow aspirates of 19 samples with PC dyscrasias by both a standard method based on bulk DNA and by an implementation of the standard 10x Genomics protocol to detect expressed V(D)J sequences. A dominant clonotype was easily identified in each sample, accounting on average for 83.65% of V(D)J-rearranged cells. Compared with standard methods, scRNA-seq analysis proved highly concordant and even more effective in identifying clonal productive rearrangements, by-passing limitations related to the misannealing of consensus primers in hypermutated regions. We next validated its accuracy to track 5 clonal cells with absolute sensitivity in a virtual sample containing 3180 polyclonal cells. This shows that single-cell V(D)J analysis may be used to find rare clonal cells, laying the foundations for functional single-cell dissection of minimal residual disease.     

Development and Characterization of Novel Microsatellite Markers for the Peach Fruit Moth Carposina sasakii (Lepidoptera: Carposinidae) Using Next-Generation Sequencing

The peach fruit moth Carposina sasakii is an economically important pest on dozens of fruits from Rosaceae and Rhamnaceae in Northeast Asia. We developed novel microsatellite markers for C. sasakii from randomly sequenced regions of the genome using next-generation sequencing. In total, 95,153 microsatellite markers were isolated from 4.70 GB genomic sequences. Thirty-five polymorphic markers were developed by assessing in 63 individuals from two geographical populations. The allele numbers ranged from 2 to 9 with an average value of 4.60 per locus, while the polymorphism information content ranged from 0.075 to 0.696 with an average value of 0.407. Furthermore, the observed and expected heterozygosity varied from 0.000 to 0.677 and 0.062 to 0.771, respectively. The microsatellites developed provide abundant molecular markers for investigating genetic structure, genetic diversity, and existence of host-plant associated biotypes of C. sasakii.     

Fatty Imidazolines: A Novel Extractant for the Recovery of Palladium(II) from Hydrochloric Acid Solutions

The extraction of palladium(II) from hydrochloric acid solutions with a novel highly basic extractant, a mixture of homologous 1-[2-(alkanoylamino)ethyl]-2-alkyl-2-imidazolines (AAI) in toluene with 15% (v/v) of n-octanol was studied. Palladium(II) is rapidly and most effectively extracted with AAI hydrochloride at the low hydrochloric acid (chloride ions) concentration (up to 1 M) and can be completely separated from Fe(III), Cu(II), and Co(II). The palladium(II) extraction at the optimum acidity occurs via an anion-exchange mechanism with the formation of ionic associates (LH)₂PdCl₄ (K _ex = (1.5 ± 0.2) · 10⁴ at 0.5 M HCl) and is accompanied by the dimerization of palladium(II) in the organic phase with the formation of ionic associates (LH)₂Pd₂Cl₆ (K _dim = (3.9 ± 0.4) · 10^?4 at 0.5 M HCl). The anion-exchange extraction of palladium(II) at the acidity of 0.5 M HCl is temperature independent in the range 20–49°C. Complete stripping of palladium(II) can be performed using a 5% solution of thiourea in 0.1 M HCl.     

Genome Editing of Golden SNP-Carrying Lycopene Epsilon-Cyclase (LcyE) Gene Using the CRSPR-Cas9/HDR and Geminiviral Replicon System in Rice

Lycopene epsilon-cyclase (LcyE) is a key enzyme in the carotenoid biosynthetic pathway of higher plants. Using the CRSPR/Cas9 and the geminiviral replicon, we optimized a method for targeted mutagenesis and golden SNP replacement of the LcyE gene in rice. We have exploited the geminiviral replicon amplification as a means to provide a large amount of donor template for the repair of a CRISPR-Cas-induced DNA double-strand break (DSB) in the target gene via homology-directed repair (HDR). Mutagenesis experiments performed on the Donggin variety achieved precise modification of the LcyE loci with an efficiency of up to 90%. In HDR experiments, our target was the LcyE allele (LcyE-H523L) derived from anther culture containing a golden SNP replacement. The phenotype of the homologous recombination (HR) mutant obtained through the geminiviral replicon-based template delivery system was tangerine color, and the frequency was 1.32% of the transformed calli. In addition, the total carotenoid content of the LcyEsg2-HDR1 and LcyEsg2-HDR2 lines was 6.8–9.6 times higher than that of the wild-type (WT) calli, respectively. The reactive oxygen species content was lower in the LcyEsg2-HDR1 and LcyEsg2-HDR2 lines. These results indicate that efficient HDR can be achieved in the golden SNP replacement using a single and modular configuration applicable to different rice targets and other crops. This work demonstrates the potential to replace all genes with elite alleles within one generation and greatly expands our ability to improve agriculturally important traits.     

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.     

Mystery(n) Phenotypic Presentation in Europeans: Report of Three Further Novel Missense RNF213 Variants Leading to Severe Syndromic Forms of Moyamoya Angiopathy and Literature Review

Moyamoya angiopathy (MMA) is a rare cerebral vasculopathy in some cases occurring in children. Incidence is higher in East Asia, where the heterozygous p.Arg4810Lys variant in RNF213 (Mysterin) represents the major susceptibility factor. Rare variants in RNF213 have also been found in European MMA patients with incomplete penetrance and are today a recognized susceptibility factor for other cardiovascular disorders, from extracerebral artery stenosis to hypertension. By whole exome sequencing, we identified three rare and previously unreported missense variants of RNF213 in three children with early onset of bilateral MMA, and subsequently extended clinical and radiological investigations to their carrier relatives. Substitutions all involved highly conserved residues clustered in the C-terminal region of RNF213, mainly in the E3 ligase domain. Probands showed a de novo occurring variant, p.Phe4120Leu (family A), a maternally inherited heterozygous variant, p.Ser4118Cys (family B), and a novel heterozygous variant, p.Glu4867Lys, inherited from the mother, in whom it occurred de novo (family C). Patients from families A and C experienced transient hypertransaminasemia and stenosis of extracerebral arteries. Bilateral MMA was present in the proband’s carrier grandfather from family B. The proband from family C and her carrier mother both exhibited annular figurate erythema. Our data confirm that rare heterozygous variants in RNF213 cause MMA in Europeans as well as in East Asian populations, suggesting that substitutions close to positions 4118–4122 and 4867 of RNF213 could lead to a syndromic form of MMA showing elevated aminotransferases and extracerebral vascular involvement, with the possible association of peculiar skin manifestations.     

ABCBA Pentablock Copolymers Consisting of Poly(l‐lactide) (PLLA: A), Poly(d‐lactide) (PDLA: B), and Poly(butylene succinate) (PBS: C): Effects of Semicrystalline PBS Segments on the Stereo‐Crystallinity and Properties

A new stereo pentablock copolymer consisting of poly(l ‐lactide) (PLLA: A), poly‐d ‐lactide (PDLA: B), and poly(butylene succinate) (PBS: C) is synthesized by two‐step ring‐opening polymerization of d ‐ and l ‐lactides in the presence of bis‐hydroxyl‐terminated PBS prepolymer that has been prepared by the ordinary polycondensation. The pentablock copolymers (PLLA‐PDLA‐PBS‐PDLA‐PLLA) as well as the triblock copolymers (PLLA‐PBS‐PLLA) obtained as the intermediates show different properties depending on the polymer compositions. In the pentablock copolymers, the direct connection of the PLLA and PDLA blocks allows easy formation of the stereocomplex crystals, while the introduction of the semicrystalline PBS block is effective not only for changing the crystallization kinetics but also for imparting an elastomeric property.

     

Identification of the Biosynthetic Gene Cluster for Himeic Acid A: A Ubiquitin‐Activating Enzyme (E1) Inhibitor in Aspergillus japonicus MF275

Himeic acid A, which is produced by the marine fungus Aspergillus japonicus MF275, is a specific inhibitor of the ubiquitin‐activating enzyme E1 in the ubiquitin–proteasome system. To elucidate the mechanism of himeic acid biosynthesis, feeding experiments with labeled precursors have been performed. The long fatty acyl side chain attached to the pyrone ring is of polyketide origin, whereas the amide substituent is derived from leucine. These results suggest that a polyketide synthase–nonribosomal peptide synthase (PKS‐NRPS) is involved in himeic acid biosynthesis. A candidate gene cluster was selected from the results of genome sequencing analysis. Disruption of the PKS‐NRPS gene by Agrobacterium‐mediated transformation confirms that HimA PKS‐NRPS is involved in himeic acid biosynthesis. Thus, the him biosynthetic gene cluster for himeic acid in A. japonicus MF275 has been identified.     

Gene Expression of the D-Series Resolvin Pathway Predicts Activation of Anti-Tumor Immunity and Clinical Outcomes in Head and Neck Cancer

In human medicine, the progression from early neoplasia development to either complete resolution of tumorigenesis and associated inflammation or chronicity and fatal outcomes remain difficult to predict. Resolution of inflammation is an active process that stimulates the termination of the inflammatory response and promotes return to homeostasis, while failure in resolution contributes to the development of a number of diseases. To understand how resolution pathways contribute to tumorigenesis, we defined and employed a cumulative score based on the expression level of genes involved in synthesis, signaling, and metabolism of the D-series resolvin (RvD). This score was used for comparative analyses of clinical, cellular, and molecular features of tumors, based on RNA-sequencing (RNA-seq) datasets collected within The Cancer Genome Atlas (TCGA) program. Our results indicate that higher RvD scores are associated with better clinical outcome of patients with head and neck squamous cell carcinoma (HNSC), and with molecular and cellular signatures indicative of enhanced anti-tumor immunity and better response to immune-checkpoint inhibitors (ICI), also in human papilloma virus (HPV) negative HNSC subtypes. Thus, higher activity of the RvD pathway identifies patients with improved resolution and a more efficient immune reaction against cancer.     

The Immunomodulatory and Neuroprotective Effects of Mesenchymal Stem Cells (MSCs) in Experimental Autoimmune Encephalomyelitis (EAE): A Model of Multiple Sclerosis (MS)

Mesenchymal stem cells (MSCs) are multipotent cells that differentiate into the mesenchymal lineages of adipocytes, osteocytes and chondrocytes. MSCs can also transdifferentiate and thereby cross lineage barriers, differentiating for example into neurons under certain experimental conditions. MSCs have anti-proliferative, anti-inflammatory and anti-apoptotic effects on neurons. Therefore, MSCs were tested in experimental autoimmune encephalomyelitis (EAE), an animal model of multiple sclerosis (MS), for their effectiveness in modulating the pathogenic process in EAE to develop effective therapies for MS. The data in the literature have shown that MSCs can inhibit the functions of autoreactive T cells in EAE and that this immunomodulation can be neuroprotective. In addition, MSCs can rescue neural cells via a mechanism that is mediated by soluble factors, which provide a suitable environment for neuron regeneration, remyelination and cerebral blood flow improvement. In this review, we discuss the effectiveness of MSCs in modulating the immunopathogenic process and in providing neuroprotection in EAE.