Effects of Cryopreservation and Thawing on Single-Cell Transcriptomes of Human T Cells

Cryopreservation and thawing of PBMCs are inevitable processes in expanding the scale of experiments in human immunology. Here, we carried out a fundamental study to investigate the detailed effects of PBMC cryopreservation and thawing on transcriptomes. We sorted Tregs from fresh and cryopreserved/thawed PBMCs from an identical donor and performed single-cell RNA-sequencing (scRNA-seq). We found that the cryopreservation and thawing process minimally affects the key molecular features of Tregs, including FOXP3. However, the cryopreserved and thawed sample had a specific cluster with up-regulation of genes for heat shock proteins. Caution may be warranted in interpreting the character of any cluster of cells with heat shock-related properties when cryopreserved and thawed samples are used for scRNA-seq.

Gene mapping study for constitutive skin color in an isolated Mongolian population

To elucidate the genes responsible for constitutive human skin color, we measured the extent of skin pigmentation in the buttock, representative of lifelong non-sun-exposed skin, and conducted a gene mapping study on skin color in an isolated Mongolian population composed of 344 individuals from 59 families who lived in Dashbalbar, Mongolia. The heritability of constitutive skin color was 0.82, indicating significant genetic association on this trait. Through the linkage analysis using 1,039 short tandem repeat (STR) microsatellite markers, we identified a novel genomic region regulating constitutive skin color on 11q24.2 with an logarithm of odds (LOD) score of 3.39. In addition, we also found other candidate regions on 17q23.2, 6q25.1, and 13q33.2 (LOD > or = 2). Family-based association tests on these regions with suggestive linkage peaks revealed ten and two significant single nucleotide polymorphisms (SNPs) on the linkage regions of chromosome 11 and 17, respectively. We were able to discover four possible candidate genes that would be implicated to regulate human skin color: ETS1, UBASH3B, ASAM, and CLTC.

Detection of hydin Gene Duplication in Personal Genome Sequence Data

Human personal genome sequencing can be done with high efficiency by aligning a huge number of short reads derived from various next generation sequencing (NGS) technologies to the reference genome sequence. One of the major obstacles is the incompleteness of human reference genome. We tried to analyze the effect of hidden gene duplication on the NGS data using the known example of hydin gene. Hydin2 , a duplicated copy of hydin on chromosome 16q22, has been recently found to be localized to chromosome 1q21, and is not included in the current version of standard human genome reference. We found that all of eight personal genome data published so far do not contain hydin2, and there is large number of nsSNPs in hydin. The heterozygosity of those nsSNPs was significantly higher than expected. The sequence coverage depth in hydin gene was about two fold of average depth. We believe that these unique finding of hydin can be used as useful indicators to discover new hidden multiplication in human genome.

Higher mitochondrial DNA copy number is associated with lower prevalence of microalbuminuria

It has been suggested that mitochondrial dysfunction contributes to the initiation and development of atherosclerosis and cardiovascular disease. We examined the association between mitochondrial DNA (mtDNA) copy number and microalbuminuria in a cross-sectional community-based study. We measured peripheral blood mtDNA copy number in 694 adults without chronic kidney disease by a real-time PCR method. The overall prevalence of microalbuminuria (defined as an albumin creatinine ratio of 30 to 299 mg/g) was 4.5%. The prevalence of microalbuminuria decreased progressively from the lower to the upper quartiles of mtDNA copy number (6.9%, 5.7%, 2.9%, and 2.3% in quartiles 1, 2, 3, and 4, respectively, P = 0.017 for trend). Multiple logistic regression models showed that the quartile of mtDNA copy number was independently associated with the prevalence of microalbuminuria (P = 0.01 for trend). Compared with the lowest quartile, the highest quartile had an odds ratio of 0.22 for microalbuminuria (95% confidence interval, 0.05 to 0.87; P = 0.03). Higher mtDNA copy number was associated with the lower prevalence of microalbuminuria in a community-based population.

Genome-wide Linkage Study for Plasma HDL Cholesterol Level in an Isolated Population of Mongolia

High-density lipoprotein (HDL) whose primary role is to transport cholesterol from peripheral tissues to the liver, is associated with the incidence of coronary heart disease. We analyzed HDL cholesterol levels in a genetically isolated population of extended Mongolian families. A total of 1002 individuals (54.5% women) from 95 families were enrolled. After genotyping by use of 1000 microsatellite markers, we performed a genome-wide linkage search with variance component analysis. The estimated heritability of HDL cholesterol was 0.45, revealing that HDL cholesterol was under significant genetic influence. We found peak evidence of linkage (LOD score=1.88) for HDL cholesterol level on chromosome 6(nearest marker D6S1660) and potential evidences for linkage on chromosomes 1, 12 and 19 with the LOD scores of 1.32, 1.44 and 1.14, respectively. These results should pave the way for the discovery of the relevant genes by fine mapping and association analysis.

Genome-wide Linkage Study for Plasma HDL Cholesterol Level in an Isolated Population of Mongolia

High-density lipoprotein (HDL) whose primary role is to transport cholesterol from peripheral tissues to the liver, is associated with the incidence of coronary heart disease. We analyzed HDL cholesterol levels in a genetically isolated population of extended Mongolian families. A total of 1002 individuals (54.5% women) from 95 families were enrolled. After genotyping by use of 1000 microsatellite markers, we performed a genome-wide linkage search with variance component analysis. The estimated heritability of HDL cholesterol was 0.45, revealing that HDL cholesterol was under significant genetic influence. We found peak evidence of linkage (LOD score=1.88) for HDL cholesterol level on chromosome 6(nearest marker D6S1660) and potential evidences for linkage on chromosomes 1, 12 and 19 with the LOD scores of 1.32, 1.44 and 1.14, respectively. These results should pave the way for the discovery of the relevant genes by fine mapping and association analysis.

Heritability and linkage study on heart rates in a Mongolian population

Elevated heart rate has been proposed as an independent risk factor for cardiovascular diseases, but their interrelationships are not well understood. In this study, we performed a genome-wide linkage scan in 1,026 individuals (mean age 30.6 years, 54.5% women) from 73 extended families of Mongolia and determined quantitative trait loci that influence heart rate. The DNA samples were genotyped using deCODE 1,039 microsatellite markers for 3 cM density genome-wide linkage scan. Correlation analysis was carried out to evaluate the correlation of the covariates and the heart rate. T-tests of the heart rate were also performed on sex, smoking and alcohol intake. Consequently, this model was used in a nonparametric genome-wide linkage analysis using variance component model to create a multipoint logarithm of odds (LOD) score and a corresponding P value. In the adjusted model, the heritability of heart rate was estimated as 0.32 (P<.0001) and a maximum multipoint LOD score of 2.03 was observed in 77 cM region at chromosome 18. The second largest LOD score of 1.52 was seen on chromosome 5 at 216 cM. Genes located on the specified locations in chromosomes 5 and 18 may be involved in the regulation of heart rate.

Single-cell Transcriptome of Bronchoalveolar Lavage Fluid Reveals Dynamic Change of Macrophages During SARS-CoV-2 Infection in Ferrets

Although the profile of immune cells changes during the natural course of SARS-CoV-2 inflection in human patients, few studies have used a longitudinal approach to reveal their dynamic features. Here, we performed single-cell RNA sequencing of bronchoalveolar lavage fluid cells longitudinally obtained from SARS-CoV-2-infected ferrets. Landscape analysis of the lung immune microenvironment showed dynamic changes in cell proportions and characteristics in uninfected control, at 2 days post-infection (dpi) (early stage of SARS-CoV-2 infection with peak viral titer), and 5 dpi (resolution phase). NK cells and CD8+ T cells exhibited activated subclusters with interferon-stimulated features, which were peaked at 2 dpi. Intriguingly, macrophages were classified into 10 distinct subpopulations, and their relative proportions changed over the time. We observed prominent transcriptome changes among monocyte-derived infiltrating macrophages and differentiated M1/M2 macrophages, especially at 2 dpi. Moreover, trajectory analysis revealed gene expression changes from monocyte-derived infiltrating macrophages toward M1 or M2 macrophages and identified the distinct macrophage subpopulation that had rapidly undergone SARS-CoV-2-mediated activation of inflammatory responses. Finally, we found that different spectrums of M1 or M2 macrophages showed distinct patterns of gene modules downregulated by immune-modulatory drugs. Overall, these results elucidate fundamental aspects of the immune response dynamics provoked by SARS-CoV-2 infection.

Relative infectivity of the SARS-CoV-2 Omicron variant in human alveolar cells

With the emergence of multiple highly transmissible SARS-CoV-2 variants during the recent pandemic, the comparison of their infectivity has become a substantially critical issue for public health. However, a direct assessment of these viral characteristics has been challenging due to the lack of appropriate experimental models and efficient methods. Here, we integrated human alveolar organoids and single-cell transcriptome sequencing techniques to facilitate the evaluation. In a proof-of-concept study using the assay with four highly transmissible SARS-CoV-2 variants, including GR (B.1.1.119), Alpha (B.1.1.7), Delta (B.1.617.2), and Omicron (BA.1), a rapid evaluation of the relative infectivity was possible. Our results demonstrate that the Omicron (BA.1) variant is 3-5-fold more infectious to human alveolar cells than the other SARS-CoV-2 variants at the early phase of infection. To our knowledge, this study provides the first direct measurement of the infectivity of the Omicron variant and new experimental procedures that can be applied for monitoring newly emerging viral variants.

Robust three-dimensional expansion of human adult alveolar stem cells and SARS-CoV-2 infection

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), which is the cause of a present global pandemic, infects human lung alveolar cells (hACs). Characterising the pathogenesis is crucial for developing vaccines and therapeutics. However, the lack of models mirroring the cellular physiology and pathology of hACs limits the study. Here, we develop a feeder-free, long-term three-dimensional (3D) culture technique for human alveolar type 2 (hAT2) cells, and investigate infection response to SARS-CoV-2. By imaging-based analysis and single-cell transcriptome profiling, we reveal rapid viral replication and the increased expression of interferon-associated genes and pro-inflammatory genes in infected hAT2 cells, indicating robust endogenous innate immune response. Further tracing of viral mutations acquired during transmission identifies full infection of individual cells effectively from a single viral entry. Our study provides deep insights into the pathogenesis of SARS-CoV-2, and the application of long-term 3D hAT2 cultures as models for respiratory diseases.