Search

Supplementary Figure from IFNα Potentiates Anti–PD-1 Efficacy by Remodeling Glucose Metabolism in the Hepatocellular Carcinoma Microenvironment

Supplementary Table 1 from Nicotine Withdrawal Sensitivity, Linkage to chr6q26, and Association of OPRM1 SNPs in the SMOking in FAMilies (SMOFAM) Sample

Supplementary Table 2 from Nicotine Withdrawal Sensitivity, Linkage to chr6q26, and Association of OPRM1 SNPs in the SMOking in FAMilies (SMOFAM) Sample

ObjectivesIn a dangerous future pandemic without effective vaccines and medicines, a reliable screening-and-isolation strategy can be the last opportunity to keep critical facilities and communities running and avoid a complete shutdown.MethodsIn this study, we introduced an epidemiological model that included essential parameters of infection transmission and screening. With varying parameters, we studied the dynamics of viral infection in the semi-isolated communities.ResultsIn the scenario with a periodic infection screening once per 3 days and a viral basic reproduction number 3.0, more than 85% of the infection waves have a duration less than 7 days and the infection size in each of the waves is generally less than 4 individuals when the efficiency of infection discovery is 0.9 in the screening. When the period of screening was elongated to once per 7 days, the cases of infection dramatically increased to 5 folds of that mentioned previously. Further, with a weak discovery efficiency of 0.7 and the aforementioned low screening frequency, the spread of infection would be out of control.ConclusionsOur study suggests that frequent periodic screening is capable of controlling a future epidemic in a semi-isolated community without vaccines and medicines.

DNA-hydrolyzing DNAs represent an attractive type of DNA-processing catalysts distinctive from the protein-based restriction enzymes. The innate DNA property has enabled them to readily join DNA-based manipulations to promote the development of DNA biotechnology. A major in vitro selection strategy to identify these DNA catalysts relies tightly on the isolation of linear DNAs processed from a circular single-stranded (ss) DNA sequence library by self-hydrolysis. Herein, we report that by programming a terminal hybridization stem in the library, other than the previously reported classes (I & II) of deoxyribozymes, two new classes (III & IV) were identified with the old selection strategy to site-specifically hydrolyze DNA in the presence of Zn2+. Their representatives own a catalytic core consisting of ∼20 conserved nucleotides and a half-life of ∼15 min at neutral pH. In a bimolecular construct, class III exhibits unique broad generality on the enzyme strand, which can be potentially harnessed to engineer DNA-responsive DNA hydrolyzers for detection of any target ssDNA sequence. Besides the new findings, this work should also provide an improved approach to select for DNA-hydrolyzing deoxyribozymes that use various molecules and ions as cofactors.

Background Jaw bones are the most common organs to be invaded by oral malignancies, such as oral squamous cell carcinoma (OSCC), because of their special anatomical relationship. Various serious complications, such as pathological fractures and bone pain can significantly decrease the quality of life or even survival outcomes for a patient. Although chemotherapy is a promising strategy for bone invasion treatment, its clinical applications are limited by the lack of tumor-specific targeting and poor permeability in bone tissue. Therefore, it is necessary to develop a smart bone and cancer dual targeting drug delivery platform. Results We designed a dual targeting nano-biomimetic drug delivery vehicle Asp8[H40-TPZ/IR780@(RBC-H)] that has excellent bone and cancer targeting as well as immune escape abilities to treat malignancies in jaw bones. These nanoparticles were camouflaged with a head and neck squamous cell carcinoma WSU-HN6 cell (H) and red blood cell (RBC) hybrid membrane, which were modified by an oligopeptide of eight aspartate acid (Asp8). The spherical morphology and typical core-shell structure of biomimetic nanoparticles were observed by transmission electron microscopy. These nanoparticles exhibited the same surface proteins as those of WSU-HN6 and RBC. Flow cytometry and confocal microscopy showed a greater uptake of the biomimetic nanoparticles when compared to bare H40-PEG nanoparticles. Biodistribution of the nanoparticles in vivo revealed that they were mainly localized in the area of bone invasion by WSU-HN6 cells. Moreover, the Asp8[H40-TPZ/IR780@(RBC-H)] nanoparticles exhibited effective cancer growth inhibition properties when compared to other TPZ or IR780 formulations. Conclusions Asp8[H40-TPZ/IR780@(RBC-H)] has bone targeting, tumor-homing and immune escape abilities, therefore, it is an efficient multi-targeting drug delivery platform for achieving precise anti-cancer therapy during bone invasion. Graphical Abstract

The overall response rate for anti–PD-1 therapy remains modest in hepatocellular carcinoma (HCC). We found that a combination of IFNα and anti–PD-1–based immunotherapy resulted in enhanced antitumor activity in patients with unresectable HCC. In both immunocompetent orthotopic and spontaneous HCC models, IFNα therapy synergized with anti–PD-1 and the combination treatment led to significant enrichment of cytotoxic CD27+CD8+ T cells. Mechanistically, IFNα suppressed HIF1α signaling by inhibiting FosB transcription in HCC cells, resulting in reduced glucose consumption capacity and consequentially establishing a high-glucose microenvironment that fostered transcription of the T-cell costimulatory molecule Cd27 via mTOR–FOXM1 signaling in infiltrating CD8+ T cells. Together, these data reveal that IFNα reprograms glucose metabolism within the HCC tumor microenvironment, thereby liberating T-cell cytotoxic capacities and potentiating the PD-1 blockade–induced immune response. Our findings suggest that IFNα and anti–PD-1 cotreatment is an effective novel combination strategy for patients with HCC. Significance: Our study supports a role of tumor glucose metabolism in IFNα-mediated antitumor immunity in HCC, and tumor-infiltrating CD27+CD8+ T cells may be a promising biomarker for stratifying patients for anti–PD-1 therapy. See related commentary by Kao et al., p. 1615. This article is highlighted in the In This Issue feature, p. 1599

Hepatitis B virus surface antigen (HBsAg) encoded by the S gene is highly expressed during the replication cycle of hepatitis B virus (HBV). However, the frequent usage of tryptophan in HBsAg, which leads to a high cost of biosynthesis, is inconsistent with the high expression level of this protein. Tryptophan-truncated mutation of HBsAg, that is, a tryptophan to stop codon mutation resulting in truncated HBsAg, might help to maintain its high expression with lower biosynthetic cost. We aimed to investigate the prevalence of tryptophan-truncated S quasispecies in treatment-naïve patients with chronic hepatitis B (CHB) by applying CirSeq as well as a site-by-site algorithm developed by us to identify variants at extremely low frequencies in the carboxyl terminus of HBsAg. A total of 730 mutations were identified in 27 patients with CHB, varying from seven to 56 mutations per sample. The number of synonymous mutations was much higher than that of nonsynonymous mutations in the reverse transcriptase (RT) coding region and vice versa in the S coding region, implying that the evolutionary constraints on the RT and S genes might be different. We showed that 25 (92.6 %) of 27 patients had at least one S-truncated mutation, most of which were derived from tryptophan, indicating a high prevalence of tryptophan-truncated S mutations in treatment-naïve patients with CHB. In terms of the RT gene, 21 (77.8 %) patients had pre-existing drug-resistant mutations, while no truncated mutations were detected. Our findings that tryptophan-truncated S quasispecies and drug-resistant RT mutants were highly prevalent in treatment-naïve patients with CHB provide new insights into the composition of the HBV population, which might help optimize the treatment and management of patients with CHB.

Background: The cancer molecular targeted therapy has achieved unprecedented progress in the past decade and is thought to be the most promising direction for cancer treatment in future. As the fast growing of the clinical trials of targeted anticancer agents for different cancer types, it is critical to collect and integrate such information to guide clinical practice. Methods: We constructed the Cancer Molecular Targeted Therapy database (CMTTdb) to store and retrieve molecular targeted therapy data about randomized clinical trials (RCTs) of targeted agents and also accompanied targets, biomarkers, targeted cancer subtypes, etc. Results: Different with some existing resources, CMTTdb particularly focuses on clinical application of the trails. Design of the trails, such as treatment modalities (monotherapy or combination with other therapies), as well as results on clinical efficacy parameters, adverse events are also collected. In this current version, CMTTdb contains data for 1,088 clinical trials which cover 165 agents, 80 targets, 15 cancer types (95 molecular subtypes and 56 histological or cytological subtypes) from public literatures. This database is freely available at http://www.biosino.org/CMTTdb. A user-friendly web interface was designed so that these data can be easily retrieved. Conclusions: CMTTdb will be a valuable source for providing access to information of clinical trials on the rapidly growing number of novel targeted agent and be useful in guiding oncologists for the optimization of the therapy strategy for cancer treatment.

Oral squamous cell carcinoma (OSCC) is a common pernicious tumor in the head and neck regions. However, the function of tumor extracellular matrix (ECM) has not been elucidated. A tissue engineering method was applied for remodeling ECM through decellularization. The cellular components were removed, and the biological composition was mostly preserved. Proteomics was performed to analyze the characterization between normal and tumor ECM. According to LC-MS/MS results, 26 proteins just showed in tumor ECM, and 14 proteins only showed in late-stage tumor ECM. KEGG pathway analysis showed that most variant proteins were linked to metabolic regulation and tumor immunity (such as SCC-Ag1, LOX). To affirm the influence of tumor ECM on the progression of OSCC, tumor cells and macrophages were co-cultured with ECM scaffold. Marked differences in proliferation, apoptosis, and migration of OSCC cells were observed between tumor and normal ECM. Tumor ECM polarized macrophages towards an anti-inflammatory phenotype (higher IL-10 and CD68, and relatively lower CD86 and IL1-β). Collectively, these findings suggest that tumor ECM served as a permissive role in OSCC progression. SIGNIFICANCE: The variation between OSCC ECM and normal ECM confirm tumor ECM plays a significant role in OSCC deterioration, which is conducive to exploring the occurrence and progression mechanisms of OSCC, and further improving the curative effect of this disease.

Sample preparation is an important step in the untargeted cell metabolomics workflow. In this study, a gas chromatography-mass spectrometry (GC-MS)-based metabolome preparation protocol was optimized with the aim of achieving global, reliable and reproducible intracellular metabolite profiling of leukemia cells. We have evaluated and compared 8 sample preparation methods: 4 extraction solvents (i.e., pure methanol, aqueous acetonitrile, methanol/methanol/water and aqueous methanol) coupled with 2 common trimethylsilylation reagents (i.e., N,O-bis(trimethylsilyl)trifluoroacetamide and N-methyl-N-(trimethylsilyl)trifluoroacetamide). The extraction solvents and derivatization agents were evaluated for their suitability for leukemia cells using criteria based on the number of total detected metabolites, extraction efficiency and repeatability. Based on the overall performance, methanol/methanol/water was chosen as the optimal extraction solvent for covering the highest number of intracellular metabolites with the best reproducibility. The derivatization agent N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) proved to be more suitable for profiling of the leukemia intracellular metabolome with better repeatability than N-methyl-N-(trimethylsilyl)trifluoroacetamide (MSTFA). This optimized method provides a robust, efficient and reproducible way to gather information on the metabolome of a leukemia cell. It also offers a solid basis to further understand the physiology of leukemia cells and unravel the mechanism of leukemia’s multidrug resistance from the perspective of metabolomics.

Background Nucleoside/nucleotide analogues are widely used to treat chronic HBV infection, but drug resistance is common. The role of HBV surface gene variants in drug resistance to nucleoside/nucleotide analogues is unknown. We are trying to investigate the dynamics of S gene mutations and how they relate to a patient's virological response in this study. Methods Thirty patients with chronic hepatitis B were enrolled and serum samples were collected at multiple time points during treatment with telbivudine (LdT). The coding regions of the small surface antigen (S-HBsAg) were amplified and sequenced using the 454 GS FLX+ System. Results Sequencing results revealed different dynamics of non-synonymous mutations, such as sL9P, sN40S, sG44E, sW172*, sW182* and sS187F, between patients with a complete virological response and those with a partial virological response. The viral population heterogeneity decreased at week 12 of LdT treatment in patients with a complete virological response, with a concomitant decline in non-synonymous mutations (from an average of 14 to 9.9 per sample) and an increase in the frequencies of major variants (from 14.3% to 40.4%). Conclusions Our findings suggest that the decrease in viral population heterogeneity at an early stage of LdT treatment was associated with the subsequent optimal virological response, and the early appearance of some specific mutations, such as sG44E, sW172* and sW182*, is a potential indicator of a partial virological response in continuing therapy.

Background Kluyveromyces marxianus, the known fastest-growing eukaryote on the earth, has remarkable thermotolerance and capacity to utilize various agricultural residues to produce low-cost bioethanol, and hence is industrially important to resolve the imminent energy shortage crisis. Currently, the poor ethanol tolerance hinders its operable application in the industry, and it is necessary to improve K. marxianus’ ethanol resistance and unravel the underlying systematical mechanisms. However, this has been seldom reported to date. Results We carried out a wild-type haploid K. marxianus FIM1 in adaptive evolution in 6% (v/v) ethanol. After 100-day evolution, the KM-100d population was obtained; its ethanol tolerance increased up to 10% (v/v). Interestingly, DNA analysis and RNA-seq analysis showed that KM-100d yeasts’ ethanol tolerance improvement was not due to ploidy change or meaningful mutations, but founded on transcriptional reprogramming in a genome-wide range. Even growth in an ethanol-free medium, many genes in KM-100d maintained their up-regulation. Especially, pathways of ethanol consumption, membrane lipid biosynthesis, anti-osmotic pressure, anti-oxidative stress, and protein folding were generally up-regulated in KM-100d to resist ethanol. Notably, enhancement of the secretory pathway may be the new strategy KM-100d developed to anti-osmotic pressure, instead of the traditional glycerol production way in S. cerevisiae. Inferred from the transcriptome data, besides ethanol tolerance, KM-100d may also develop the ability to resist osmotic, oxidative, and thermic stresses, and this was further confirmed by the cell viability test. Furthermore, under such environmental stresses, KM-100d greatly improved ethanol production than the original strain. In addition, we found that K. marxianus may adopt distinct routes to resist different ethanol concentrations. Trehalose biosynthesis was required for low ethanol, while sterol biosynthesis and the whole secretory pathway were activated for high ethanol. Conclusions This study reveals that ethanol-driven laboratory evolution could improve K. marxianus’ ethanol tolerance via significant up-regulation of multiple pathways including anti-osmotic, anti-oxidative, and anti-thermic processes, and indeed consequently raised ethanol yield in industrial high-temperature and high-ethanol circumstance. Our findings give genetic clues for further rational optimization of K. marxianus’ ethanol production, and also partly confirm the positively correlated relationship between yeast’s ethanol tolerance and production. Electronic supplementary material The online version of this article (10.1186/s13068-019-1393-z) contains supplementary material, which is available to authorized users.

Additional file 1: Fig. S1. The growth rate of KM and KM-100d at different ethanol concentrations. Fig. S2. Flow cytometer analysis for DNA content variation during K. marxianus evolution. Fig. S3. Heat map for log2ratio values with clustered genes and groups. Fig. S4. GO enrichment analysis for KM and KM-100d exposed in low and high ethanol. Fig. S5. Secretory pathway and cell wall biogenesis in K. marxianus pre- and post-evolution. Table S1. Primers for RT-qPCR analysis.

Human pigmentation is a highly diverse and complex trait among populations and has drawn particular attention from both academic and non-academic investigators for thousands of years. Previous studies detected selection signals in several human pigmentation genes, but few studies have integrated contribution from multiple genes to the evolution of human pigmentation. Moreover, none has quantified selective pressures on human pigmentation over epochs and between populations. Here, we dissect dynamics and differences of selective pressures during different periods and between distinct populations with new approaches. We use genotype data of 19 genes associated with human pigmentation from 17 publicly available datasets and obtain data for 2346 individuals of six representative population groups from across the world. Our results quantify the strength of natural selection on light pigmentation not only in modern Europeans (0.0259/generation) but also in proto-Eurasians (0.00650/generation). Our results also suggest that several derived alleles associated with human dark pigmentation may be under positive directional selection in some African populations. Our study provides the first attempt to quantitatively investigate the dynamics of selective pressures during different time periods in the evolution of human pigmentation. This article has an associated First Person interview with the first author of the article., Summary: This study developed a novel approach to estimate selection strength over time and applied it on the evolution of human pigmentation.

Human populations around the world encounter various environmental challenges and, consequently, develop genetic adaptations to different selection forces. Identifying the differences in natural selection between populations is critical for understanding the roles of specific genetic variants in evolutionary adaptation. Although numerous methods have been developed to detect genetic loci under recent directional selection, a probabilistic solution for testing and quantifying selection differences between populations is lacking. Here we report the development of a probabilistic method for testing and estimating selection differences between populations. By use of a probabilistic model of genetic drift and selection, we showed that logarithm odds ratios of allele frequencies provide estimates of the differences in selection coefficients between populations. The estimates approximate a normal distribution, and variance can be estimated using genome-wide variants. This allows us to quantify differences in selection coefficients and to determine the confidence intervals of the estimate. Our work also revealed the link between genetic association testing and hypothesis testing of selection differences. It therefore supplies a solution for hypothesis testing of selection differences. This method was applied to a genome-wide data analysis of Han and Tibetan populations. The results confirmed that both the EPAS1 and EGLN1 genes are under statistically different selection in Han and Tibetan populations. We further estimated differences in the selection coefficients for genetic variants involved in melanin formation and determined their confidence intervals between continental population groups. Application of the method to empirical data demonstrated the outstanding capability of this novel approach for testing and quantifying differences in natural selection.

The Y chromosome is one of the best genetic materials to explore the evolutionary history of human populations. Global analyses of Y chromosomal short tandem repeats (STRs) data can reveal very interesting world population structures and histories. However, previous Y-STR works tended to focus on small geographical ranges or only included limited sample sizes. In this study, we have investigated population structure and demographic history using 17 Y chromosomal STRs data of 979 males from 44 worldwide populations. The largest genetic distances have been observed between pairs of African and non-African populations. American populations with the lowest genetic diversities also showed large genetic distances and coancestry coefficients with other populations, whereas Eurasian populations displayed close genetic affinities. African populations tend to have the oldest time to the most recent common ancestors (TMRCAs), the largest effective population sizes and the earliest expansion times, whereas the American, Siberian, Melanesian, and isolated Atayal populations have the most recent TMRCAs and expansion times, and the smallest effective population sizes. This clear geographic pattern is well consistent with serial founder model for the origin of populations outside Africa. The Y-STR dataset presented here provides the most detailed view of worldwide population structure and human male demographic history, and additionally will be of great benefit to future forensic applications and population genetic studies.

Studies of natural selection, followed by functional validation, are shedding light on understanding of genetic mechanisms underlying human evolution and adaptation. Classic methods for detecting selection, such as the integrated haplotype score (iHS) and Fay and Wu's H statistic, are useful for candidate gene searching underlying positive selection. These methods, however, have limited capability to localize causal variants in selection target regions. In this study, we developed a novel method based on conditional coalescent tree to detect recent positive selection by counting unbalanced mutations on coalescent gene genealogies. Extensive simulation studies revealed that our method is more robust than many other approaches against biases due to various demographic effects, including population bottleneck, expansion, or stratification, while not sacrificing its power. Furthermore, our method demonstrated its superiority in localizing causal variants from massive linked genetic variants. The rate of successful localization was about 20-40% higher than that of other state-of-the-art methods on simulated data sets. On empirical data, validated functional causal variants of four well-known positive selected genes were all successfully localized by our method, such as ADH1B, MCM6, APOL1, and HBB. Finally, the computational efficiency of this new method was much higher than that of iHS implementations, that is, 24-66 times faster than the REHH package, and more than 10,000 times faster than the original iHS implementation. These magnitudes make our method suitable for applying on large sequencing data sets. Software can be downloaded from https://github.com/wavefancy/scct.

Autism is a highly heritable neurodevelopmental disorder, and known genetic variants, mostly rare, account only for a small proportion of cases. Here we report a genome-wide association study on autism using two Chinese cohorts as gene discovery (n=2150) and three data sets of European ancestry populations for replication analysis of top association signals. Meta-analysis identified three single-nucleotide polymorphisms, rs936938 (P=4.49 × 10(-8)), non-synonymous rs6537835 (P=3.26 × 10(-8)) and rs1877455 (P=8.70 × 10(-8)), and related haplotypes, AMPD1-NRAS-CSDE1, TRIM33 and TRIM33-BCAS2, associated with autism; all were mapped to a previously reported linkage region (1p13.2) with autism. These genetic associations were further supported by a cis-acting regulatory effect on the gene expressions of CSDE1, NRAS and TRIM33 and by differential expression of CSDE1 and TRIM33 in the human prefrontal cortex of post-mortem brains between subjects with and those without autism. Our study suggests TRIM33 and NRAS-CSDE1 as candidate genes for autism, and may provide a novel insight into the etiology of autism.

Sherpa population is an ethnic group living in south mountainside of Himalayas for hundreds of years. They are famous as extraordinary mountaineers and guides, considered as a good example for successful adaptation to low oxygen environment in Tibetan highlands. Mitochondrial DNA (mtDNA) variations might be important in the highland adaption given its role in coding core subunits of oxidative phosphorylation in mitochondria. In this study, we sequenced the complete mtDNA genomes of 76 unrelated Sherpa individuals. Generally, Sherpa mtDNA haplogroup constitution was close to Tibetan populations. However, we found three lineage expansions in Sherpas, two of which (C4a3b1 and A4e3a) were Sherpa-specific. Both lineage expansions might begin within the past hundreds of years. Especially, nine individuals carry identical Haplogroup C4a3b1. According to the history of Sherpas and Bayesian skyline plot, we constructed various demographic models and found out that it is unlikely for these lineage expansions to occur in neutral models especially for C4a3b1. Nonsynonymous mutations harbored in C4a3b1 (G3745A) and A4e3a (T4216C) are both ND1 mutants (A147T and Y304H, respectively). Secondary structure predictions showed that G3745A were structurally closing to other pathogenic mutants, whereas T4216C itself was reported as the primary mutation for Leber's hereditary optic neuropathy. Thus, we propose that these mutations had certain effect on Complex I function and might be important in the high altitude adaptation for Sherpa people.

Mutations of the influenza virus lead to antigenic changes that cause recurrent epidemics and vaccine resistance. Preventive measures would benefit greatly from the ability to predict the potential distribution of new antigenic sites in future strains. By leveraging the extensive historical records of HA sequences for 90 years, we designed a computational model to simulate the dynamic evolution of antigenic sites in A/H1N1. With templates of antigenic sequences, the model can effectively predict the potential distribution of future antigenic mutants. Validation on 10932 HA sequences from the last 16 years showing that the mutated antigenic sites of over 94% of reported strains fell in our predicted profile. Meanwhile, our model can successfully capture 96% of antigenic sites in those dominant epitopes. Similar results are observed on the complete set of H3N2 historical data, supporting the general applicability of our model to multiple sub-types of influenza. Our results suggest that the mutational profile of future antigenic sites can be predicted based on historical evolutionary traces despite the widespread, random mutations in influenza. Coupled with closely monitored sequence data from influenza surveillance networks, our method can help to forecast changes in viral antigenicity for seasonal flu and inform public health interventions.

A few genetic polymorphisms of TP53 are known to have a significant effect on cancer susceptibility. Intron 3 16 bp duplication polymorphism of TP53 has been reported to be associated with breast cancer, colorectal cancer, lung cancer and other cancers, but the reported results remain inconclusive. The present study, a meta-analysis including a total of 9801 cases and 10,391 controls from 26 studies, revealed that the 16 bp insertion (Ins) allele is significantly associated with an increased cancer risk in overall analysis [Ins/Ins + deletion (Del)/Ins versus Del/Del: odds ratio (OR) = 1.14, 95% confidence interval (CI) = 1.02-1.27, P = 0.02; Ins/Ins versus Del/Del: OR = 1.35, 95% CI = 1.11-1.63, P = 0.002; Del/Ins versus Del/Del: OR = 1.10, 95% CI = 0.98-1.23, P = 0.11.), particularly in breast cancer subgroup (Ins/Ins + Del/Ins versus Del/Del: OR = 1.16, 95% CI = 1.03-1.31, P = 0.02; Ins/Ins versus Del/Del: OR = 1.81, 95% CI = 1.30-2.52, P < 0.001; Del/Ins versus Del/Del: OR = 1.10, 95% CI = 0.97-1.25, P = 0.13). The relative risks to the colorectal and lung cancers increased but their association power was relatively weak, which may result from a limited number of studies of these two cancer types. These results suggest that intron 3 16 bp duplication polymorphism of TP53 is potentially an important and clinically relevant genetic marker contributing to cancer susceptibility.

The outbreak of a novel influenza A (H1N1) virus across the globe poses a threat to human health. It is of paramount importance to develop a rapid, reliable and inexpensive diagnostic procedure. Based on the bioinformatic information from public database, primers specific for influenza A virus surface protein haemagglutinin (HA) of several subtypes (including H1, H2, H3, H5, H7 and H9) were designed. Primer-specific PCR products were subjected to sequencing for accurately distinguishing H1 and H3 subtypes from others. This sequencing-based detection method will not only be applied to rapid detection and simultaneous subtype identification of new influenza A virus H1N1, but also provide the strategies to monitor other new types of influenza virus with explosive potential.

We utilized a cohort of 828 treatment-seeking self-identified white cigarette smokers (50% female) to rank candidate gene single nucleotide polymorphisms (SNPs) associated with the Fagerström Test for Nicotine Dependence (FTND), a measure of nicotine dependence which assesses quantity of cigarettes smoked and time- and place-dependent characteristics of the respondent's smoking behavior. A total of 1123 SNPs at 55 autosomal candidate genes, nicotinic acetylcholine receptors and genes involved in dopaminergic function, were tested for association to baseline FTND scores adjusted for age, depression, education, sex, and study site. SNP P-values were adjusted for the number of transmission models, the number of SNPs tested per candidate gene, and their intragenic correlation. DRD2, SLC6A3, and NR4A2 SNPs with adjusted P-values

Summary The power of statistical tests to measure effect sizes in the presence of population stratification is an important issue for the design and analysis of population-based association studies. Comparisons of statistical tests have shown that the power of different statistical approaches varies in different genetic scenarios. However, the impact of stratified population parameters on statistical power is not yet understood in a general statistical framework, particularly the impact of correlated population parameters. To investigate such impact in detail, we implemented a genetic model for population-based association studies with stratified samples and evaluated the impact on power with different genetic scenarios. The investigation shows that correlation between disease prevalence and risk allele frequency among subpopulations impacts statistical power. In a model with five subpopulations and moderate population divergence (Fst= 0.01), the correlation accounts for more than 85% of power difference. Our results also show that the estimation of genetic effect for candidate loci is biased by population divergence. Beneficial alleles could be wrongly characterized as risk alleles when prevalence differences and divergences of risk loci are large among subpopulations.

Summary The gene families in the AZFc region of the Y chromosome have been shown to be functionally important in human spermatogenesis. The gr/gr deletion, a partial AZFc deletion that reduces the copy numbers of all the AZFc gene families, was identified as a significant risk factor for spermatogenic impairment in Dutch, Spanish and Italians. However, the presence of this deletion in healthy French and Germans questioned its importance in male infertility. In this study, we have shown that the gr/gr deletion does not render an increased risk in Han Chinese. In fact, the gr/gr deletion is frequent (about 8%) in our survey of 886 East Asians from 8 ethnic groups. Furthermore, the DAZ1/DAZ2 deletion has been detected as the primary subtype of the gr/gr deletion in East Asians, though this doublet has been considered as crucial for normal spermatogenesis in Europeans. The different spermatogenic effects of various types of the partial AZFc deletion suggest that the functional difference between AZFc gene copies is a likely cause of inconsistent associations of the gr/gr deletion with spermatogenic impairment across populations.

Substantial increases of linkage disequilibrium (LD) both in magnitude and in range have been observed in recently admixed populations such as African-American (AfA). On the other hand, it has also been shown that LD in AfAs was very similar to that of African. In this study, we attempted to resolve these contradicting observations by conducting a systematic examination of the LD structure in AfAs by genotyping a sample of AfA individuals at 24,341 single nucleotide polymorphisms (SNPs) spanning almost the entire chromosome 21, with an average density of 1.5 kb/SNP. The overall LD in AfAs is similar to that in African populations and much less than that in European populations. Even when the ancestry-informative markers (AIMs) were used, extended LD in AfA was found to be limited to certain magnitude range (0.2 < or = r(2) < or = 0.8) and certain distance range, that is, between-marker distance more than 200 kb. Furthermore, the inclusion of AfA individuals with predominant African ancestry was found to reduce the overall magnitude of LD. Elevation of LD in the AfA population, compared with its parental populations, can only be observed at the markers with large allele frequency differences between 2 parental populations at limited scenario. AfA individuals of wholly African ancestry contribute little to the extended LD in the AfA population, and further genotyping or association analysis conducted using only admixed individuals may lead to higher statistical power and possibly reduced cost.

Reverse transcriptase (RT) mutations contribute to hepatitis B virus resistance during antiviral therapy with nucleos(t)ide analogs. However, the composition of the RT quasispecies and their interactions during antiviral treatment have not yet been thoroughly defined. In this report, 10 patients from each of 3 different virological response groups, i.e., complete virological response, partial virological response and virological breakthrough, were selected from a multicenter trial of Telbivudine treatment. Variations in the drug resistance-related critical RT regions in 107 serial serum samples from the 30 patients were examined by ultra-deep sequencing. A total of 496,577 sequence reads were obtained, with an average sequencing coverage of 4,641X per sample. The phylogenies of the quasispecies revealed the independent origins of two critical quasispecies, i.e., the rtA181T and rtM204I mutants. Data analyses and theoretical modeling showed a cooperative-competitive interplay among the quasispecies. In particular, rtM204I mutants compete against other quasispecies, which eventually leads to virological breakthrough. However, in the absence of rtM204I mutants, synergistic growth of the drug-resistant rtA181T mutants with the wild-type quasispecies could drive the composition of the viral population into a state of partial virological response. Furthermore, we demonstrated that the frequency of drug-resistant mutations in the early phase of treatment is important for predicting the virological response to antiviral therapy.

In this study, we present an analysis of Neanderthal introgression at the dipeptidase 1 gene, DPEP1. A Neanderthal origin for the putative introgressive haplotypes was demonstrated using an established three-step approach. This introgression was under positive natural selection, reached a frequency of >50%, and introduced a homocysteine level- and pigmentation-associated allele (rs460879-T) into East Asians. However, the same allele was also found in non-East Asians, but not from Neanderthal introgression. It is likely that rs460879-T was lost in East Asians and was reintroduced subsequently through Neanderthal introgression. Our findings suggest that Neanderthal introgression could reintroduce an important previously existing allele into populations where the allele had been lost. This study sheds new light on understanding the contribution of Neanderthal introgression to the adaptation of non-Africans.

Previous research reported that Papua New Guineans (PNG) and Australians contain introgressions from Denisovans. Here we present a genome-wide analysis of Denisovan introgressions in PNG and Australians. We firstly developed a two-phase method to detect Denisovan introgressions from whole-genome sequencing data. This method has relatively high detection power (79.74%) and low false positive rate (2.44%) based on simulations. Using this method, we identified 1.34 Gb of Denisovan introgressions from sixteen PNG and four Australian genomes, in which we identified 38,877 Denisovan introgressive alleles (DIAs). We found that 78 Denisovan introgressions were under positive selection. Genes located in the 78 introgressions are related to evolutionarily important functions, such as spermatogenesis, fertilization, cold acclimation, circadian rhythm, development of brain, neural tube, face, and olfactory pit, immunity, etc. We also found that 121 DIAs are missense. Genes harboring the 121 missense DIAs are also related to evolutionarily important functions, such as female pregnancy, development of face, lung, heart, skin, nervous system, and male gonad, visual and smell perception, response to heat, pain, hypoxia, and UV, lipid transport, metabolism, blood coagulation, wound healing, aging, etc. Taken together, this study suggests that Denisovan introgressions in PNG and Australians are evolutionarily important, and may help PNG and Australians in local adaptation. In this study, we also proposed a method that could efficiently identify archaic hominin introgressions in modern non-African genomes.