Analysis of genomic distributions of SARS-CoV-2 reveals a dominant strain type with strong allelic associations.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of COVID 19, continues to evolve since its first emergence in December 2019. Using the complete sequences of 1,932 SARS-CoV-2 genomes, various clustering analyses consistently identified six types of the strains. Independent of the dendrogram construction, 13 signature variations in the form of single nucleotide variations (SNVs) in protein coding regions and one SNV in the 5' untranslated region (UTR) were identified and provided a direct interpretation for the six types (types I to VI). The six types of the strains and their underlying signature SNVs were validated in two subsequent analyses of 6,228 and 38,248 SARS-CoV-2 genomes which became available later. To date, type VI, characterized by the four signature SNVs C241T (5'UTR), C3037T (nsp3 F924F), C14408T (nsp12 P4715L), and A23403G (Spike D614G), with strong allelic associations, has become the dominant type. Since C241T is in the 5' UTR with uncertain significance and the characteristics can be captured by the other three strongly associated SNVs, we focus on the other three. The increasing frequency of the type VI haplotype 3037T-14408T-23403G in the majority of the submitted samples in various countries suggests a possible fitness gain conferred by the type VI signature SNVs. The fact that strains missing one or two of these signature SNVs fail to persist implies possible interactions among these SNVs. Later SNVs such as G28881A, G28882A, and G28883C have emerged with strong allelic associations, forming new subtypes. This study suggests that SNVs may become an important consideration in SARS-CoV-2 classification and surveillance.

Covariate-adjusted heatmaps for visualizing biological data via correlation decomposition.

Motivation: Heatmap is a popular visualization technique in biology and related fields. In this study, we extend heatmaps within the framework of matrix visualization (MV) by incorporating a covariate adjustment process through the estimation of conditional correlations. MV can explore the embedded information structure of high-dimensional large-scale datasets effectively without dimension reduction. The benefit of the proposed covariate-adjusted heatmap is in the exploration of conditional association structures among the subjects or variables that cannot be done with conventional MV. Results: For adjustment of a discrete covariate, the conditional correlation is estimated by the within and between analysis. This procedure decomposes a correlation matrix into the within- and between-component matrices. The contribution of the covariate effects can then be assessed through the relative structure of the between-component to the original correlation matrix while the within-component acts as a residual. When a covariate is of continuous nature, the conditional correlation is equivalent to the partial correlation under the assumption of a joint normal distribution. A test is then employed to identify the variable pairs which possess the most significant differences at varying levels of correlation before and after a covariate adjustment. In addition, a z-score significance map is constructed to visualize these results. A simulation and three biological datasets are employed to illustrate the power and versatility of our proposed method. Availability and implementation: GAP is available to readers and is free to non-commercial applications. The installation instructions, the user's manual, and the detailed tutorials can be found at http://gap.stat.sinica.edu.tw/Software/GAP. Supplementary information: Supplementary Data are available at Bioinformatics online.

Symptom patterns and phenotypic subgrouping of women with polycystic ovary syndrome: association between endocrine characteristics and metabolic aberrations.

STUDY QUESTION: What are the potential endocrine characteristics related to risk and severity of metabolic disturbances in women with polycystic ovary syndrome (PCOS)? SUMMARY ANSWER: Women with PCOS could be subtyped into four subgroups according to heterogeneous endocrine characteristics and the major predictive endocrine factors for metabolic aberrations among different subgroups were free androgen index (FAI) and luteinizing hormone (LH) levels. WHAT IS KNOWN ALREADY: Women diagnosed with PCOS present with highly heterogeneous phenotypes, including endocrine and metabolic aberrations. Different strategies have been proposed to predict the metabolic outcomes but whether the endocrine factors can solely predict the metabolic aberrations is still inconclusive. STUDY DESIGN, SIZE, DURATION: A cross-sectional study including 460 patients recruited from a reproductive endocrinology outpatient clinic of a tertiary medical center. PARTICIPANTS/MATERIALS, SETTING, METHODS: Patients with PCOS diagnosed according to the 2003 Rotterdam criteria were studied. Clinical history recorded by questionnaires, anthropometric measurements, biochemistry tests after an overnight fast, and pelvic ultrasonography were collected from all patients. MAIN RESULTS AND THE ROLE OF CHANCE: Applying a matrix visualization and clustering approach (generalized association plots), the patients were divided into four distinct clusters according to the correlation with four endocrine parameters. Each cluster exhibited specific endocrine characteristics and the prevalence of metabolic syndrome (MS) was significantly different among the clusters (P < 0.0001). The high-risk subgroups for MS included one cluster with higher mean (SD) FAI (39.6 (14.7) in cluster 4), and another one with lower mean (SD) FAI (10 (6.4) in cluster 2). A common endocrine characteristic of these two metabolically unhealthy clusters was relatively lower LH level. Contrarily, higher LH level (≧15 mIU/ml) during early follicular phase was found to be the best indicator of the metabolically healthy cluster (cluster 1). While high FAI level did correlate with more severe metabolic aberrations, high LH level showed better predictive value than low FAI level to become a metabolically healthy cluster. LIMITATIONS, REASONS FOR CAUTION: The results should be applied to other populations with caution due to racial or environmental differences. Another limitation is a lack of normal non-PCOS control in our study. WIDER IMPLICATIONS OF THE FINDINGS: Stratifying women with PCOS into meaningful subtypes could provide a better understanding of related risk factors and potentially enable the design and delivery of more effective screening and treatment intervention. STUDY FUNDING/COMPETING INTERESTS: This study was supported by grant NSC 100-2314-B002-027-MY3 from the National Science Council of Taiwan. TRIAL REGISTRATION NUMBER: Nil.

Analysis of protein-protein interactions in cross-talk pathways reveals CRKL protein as a novel prognostic marker in hepatocellular carcinoma.

Deciphering the network of signaling pathways in cancer via protein-protein interactions (PPIs) at the cellular level is a promising approach but remains incomplete. We used an in situ proximity ligation assay to identify and quantify 67 endogenous PPIs among 21 interlinked pathways in two hepatocellular carcinoma (HCC) cells, Huh7 (minimally migratory cells) and Mahlavu (highly migratory cells). We then applied a differential network biology analysis and determined that the novel interaction, CRKL-FLT1, has a high centrality ranking, and the expression of this interaction is strongly correlated with the migratory ability of HCC and other cancer cell lines. Knockdown of CRKL and FLT1 in HCC cells leads to a decrease in cell migration via ERK signaling and the epithelial-mesenchymal transition process. Our immunohistochemical analysis shows high expression levels of the CRKL and CRKL-FLT1 pair that strongly correlate with reduced disease-free and overall survival in HCC patient samples, and a multivariate analysis further established CRKL and the CRKL-FLT1 as novel prognosis markers. This study demonstrated that functional exploration of a disease network with interlinked pathways via PPIs can be used to discover novel biomarkers.

Using an in situ proximity ligation assay to systematically profile endogenous protein-protein interactions in a pathway network.

Signal transduction pathways in the cell require protein-protein interactions (PPIs) to respond to environmental cues. Diverse experimental techniques for detecting PPIs have been developed. However, the huge amount of PPI data accumulated from various sources poses a challenge with respect to data reliability. Herein, we collected ∼ 700 primary antibodies and employed a highly sensitive and specific technique, an in situ proximity ligation assay, to investigate 1204 endogenous PPIs in HeLa cells, and 557 PPIs of them tested positive. To overview the tested PPIs, we mapped them into 13 PPI public databases, which showed 72% of them were annotated in the Human Protein Reference Database (HPRD) and 8 PPIs were new PPIs not in the PubMed database. Moreover, TP53, CTNNB1, AKT1, CDKN1A, and CASP3 were the top 5 proteins prioritized by topology analyses of the 557 PPI network. Integration of the PPI-pathway interaction revealed that 90 PPIs were cross-talk PPIs linking 17 signaling pathways based on Reactome annotations. The top 2 connected cross-talk PPIs are MAPK3-DAPK1 and FAS-PRKCA interactions, which link 9 and 8 pathways, respectively. In summary, we established an open resource for biological modules and signaling pathway profiles, providing a foundation for comprehensive analysis of the human interactome.

AI-enhanced integration of genetic and medical imaging data for risk assessment of Type 2 diabetes.

Type 2 diabetes (T2D) presents a formidable global health challenge, highlighted by its escalating prevalence, underscoring the critical need for precision health strategies and early detection initiatives. Leveraging artificial intelligence, particularly eXtreme Gradient Boosting (XGBoost), we devise robust risk assessment models for T2D. Drawing upon comprehensive genetic and medical imaging datasets from 68,911 individuals in the Taiwan Biobank, our models integrate Polygenic Risk Scores (PRS), Multi-image Risk Scores (MRS), and demographic variables, such as age, sex, and T2D family history. Here, we show that our model achieves an Area Under the Receiver Operating Curve (AUC) of 0.94, effectively identifying high-risk T2D subgroups. A streamlined model featuring eight key variables also maintains a high AUC of 0.939. This high accuracy for T2D risk assessment promises to catalyze early detection and preventive strategies. Moreover, we introduce an accessible online risk assessment tool for T2D, facilitating broader applicability and dissemination of our findings.

Fatty liver classification via risk controlled neural networks trained on grouped ultrasound image data.

Ultrasound imaging is a widely used technique for fatty liver diagnosis as it is practically affordable and can be quickly deployed by using suitable devices. When it is applied to a patient, multiple images of the targeted tissues are produced. We propose a machine learning model for fatty liver diagnosis from multiple ultrasound images. The machine learning model extracts features of the ultrasound images by using a pre-trained image encoder. It further produces a summary embedding on these features by using a graph neural network. The summary embedding is used as input for a classifier on fatty liver diagnosis. We train the machine learning model on a ultrasound image dataset collected by Taiwan Biobank. We also carry out risk control on the machine learning model using conformal prediction. Under the risk control procedure, the classifier can improve the results with high probabilistic guarantees.

dbDNV: a resource of duplicated gene nucleotide variants in human genome.

Gene duplications are scattered widely throughout the human genome. A single-base difference located in nearly identical duplicated segments may be misjudged as a single nucleotide polymorphism (SNP) from individuals. This imperfection is undistinguishable in current genotyping methods. As the next-generation sequencing technologies become more popular for sequence-based association studies, numerous ambiguous SNPs are rapidly accumulated. Thus, analyzing duplication variations in the reference genome to assist in preventing false positive SNPs is imperative. We have identified >10% of human genes associated with duplicated gene loci (DGL). Through meticulous sequence alignments of DGL, we systematically designated 1,236,956 variations as duplicated gene nucleotide variants (DNVs). The DNV database (dbDNV) (http://goods.ibms.sinica.edu.tw/DNVs/) has been established to promote more accurate variation annotation. Aside from the flat file download, users can explore the gene-related duplications and the associated DNVs by DGL and DNV searches, respectively. In addition, the dbDNV contains 304,110 DNV-coupled SNPs. From DNV-coupled SNP search, users observe which SNP records are also variants among duplicates. This is useful while ∼58% of exonic SNPs in DGL are DNV-coupled. Because of high accumulation of ambiguous SNPs, we suggest that annotating SNPs with DNVs possibilities should improve association studies of these variants with human diseases.

Morus alba and active compound oxyresveratrol exert anti-inflammatory activity via inhibition of leukocyte migration involving MEK/ERK signaling.

BACKGROUND: Morus alba has long been used in traditional Chinese medicine to treat inflammatory diseases; however, the scientific basis for such usage and the mechanism of action are not well understood. This study investigated the action of M. alba on leukocyte migration, one key step in inflammation. METHODS: Gas chromatography-mass spectrometry (GC-MS) and cluster analyses of supercritical CO2 extracts of three Morus species were performed for chemotaxonomy-aided plant authentication. Phytochemistry and CXCR4-mediated chemotaxis assays were used to characterize the chemical and biological properties of M. alba and its active compound, oxyresveratrol. fluorescence-activated cell sorting (FACS) and Western blot analyses were conducted to determine the mode of action of oxyresveratrol. RESULTS: Chemotaxonomy was used to help authenticate M. alba. Chemotaxis-based isolation identified oxyresveratrol as an active component in M. alba. Phytochemical and chemotaxis assays showed that the crude extract, ethyl acetate fraction and oxyresveratrol from M. alba suppressed cell migration of Jurkat T cells in response to SDF-1. Mechanistic study indicated that oxyresveratrol diminished CXCR4-mediated T-cell migration via inhibition of the MEK/ERK signaling cascade. CONCLUSIONS: A combination of GC-MS and cluster analysis techniques are applicable for authentication of the Morus species. Anti-inflammatory benefits of M. alba and its active compound, oxyresveratrol, may involve the inhibition of CXCR-4-mediated chemotaxis and MEK/ERK pathway in T and other immune cells.

A Deep Generative Model for Reordering Adjacency Matrices.

Depending on the node ordering, an adjacency matrix can highlight distinct characteristics of a graph. Deriving a "proper" node ordering is thus a critical step in visualizing a graph as an adjacency matrix. Users often try multiple matrix reorderings using different methods until they find one that meets the analysis goal. However, this trial-and-error approach is laborious and disorganized, which is especially challenging for novices. This paper presents a technique that enables users to effortlessly find a matrix reordering they want. Specifically, we design a generative model that learns a latent space of diverse matrix reorderings of the given graph. We also construct an intuitive user interface from the learned latent space by creating a map of various matrix reorderings. We demonstrate our approach through quantitative and qualitative evaluations of the generated reorderings and learned latent spaces. The results show that our model is capable of learning a latent space of diverse matrix reorderings. Most existing research in this area generally focused on developing algorithms that can compute "better" matrix reorderings for particular circumstances. This paper introduces a fundamentally new approach to matrix visualization of a graph, where a machine learning model learns to generate diverse matrix reorderings of a graph.

Exploring Factors Underlying Poorly-Controlled Asthma in Adults by Integrating Phenotypes and Genotypes Associated with Obesity and Asthma: A Case-Control Study.

Background: Uncontrolled asthma in adults leads to poor clinical outcome, while the clinical heterogeneity of phenotypes interferes the applicable genetic determinants. This study aimed to identify phenotypes and genetic impact on poorly-controlled asthma to optimize individualized treatment strategies. Methods: This propensity score-matched case-control study included 340 and 1020 asthmatics with poorly-controlled asthma and well-controlled asthma, respectively. Data were obtained from the 2008-2015 Taiwan Biobank Database and linked to the National Health Insurance Research Database. All asthmatics were aged ≥30 years, without cancer history, and each completed a questionnaire, physical examination, and genome-wide single nucleotide polymorphisms (SNPs). Multivariate adjusted odds ratios (ORs) for genetic risk scores were calculated using conditional logistic regression, stratified by age and sex. A model integrating obesity- and asthma-associated phenotypes and genotypes was applied for poorly-controlled asthma risk prediction. Results: General obesity with body mass index (BMI) ≥27 kg/m2 (OR:1.49, 95% confidence interval (CI) 1.09-2.03), central obesity with waist-to-height ratio (WHtR) ≥0.5 (OR:1.62, 95% CI 1.22-2.15), and parental history of asthma (OR:1.65, and 1.68; for BMI model and WHtR model, respectively) were significantly associated with poorly-controlled asthma in adults, and the combination effect of both obesity phenotypes was 1.66 (95% CI 1.17-2.35). A total of 16 obesity-associated SNPs and 9 asthma-associated SNPs were converted into genetic scores, and the aforementioned phenotypes were incorporated into the risk prediction model for poorly-controlled asthma, with an area under curve 0.72 in the receiver operating characteristic curve. The potential biological functions of genes are involved in immunity pathways. Conclusion: The prediction model integrating obesity-asthma phenotypes and genotypes for poorly-controlled asthma can facilitate the prediction of high-risk asthma and provide potential targets for novel treatment.

LOHAS: loss-of-heterozygosity analysis suite.

Detection of loss of heterozygosity (LOH) plays an important role in genetic, genomic and cancer research. We develop computational methods to estimate the proportion of homozygous SNP calls, identify samples with structural alterations and/or unusual genotypic patterns, cluster samples with close LOH structures and map the genomic segments bearing LOH by analyzing data of genome-wide SNP arrays or customized SNP arrays. In addition to cancer genetics/genomics, we also apply the methods to study long contiguous stretches of homozygosity (LCSH) in general populations. The LCSH analysis aids in the identification of samples with complex LCSH patterns indicative of nonrandom mating and/or meiotic recombination cold spots, separation of samples with different genetic backgrounds and sex, and mapping of regions of LCSH. Affymetrix Human Mapping 500K Set SNP data from an acute lymphoblastic leukemia study containing 304 cancer patients and 50 normal controls and from the HapMap Project containing 30 African trios, 30 Caucasian trios and 90 independent Asian samples were analyzed. We identified common gene regions of LOH, e.g., ETV6 and CDKN1B, and identified frequent regions of LCSH, e.g., the region that encompasses the centromeric gene desert region of chromosome 16. Unsupervised analysis separated cancer subtypes and ethnic subpopulations by patterns of LOH/LCSH. Simulation studies considering LOH width, effect size and heterozygous interference fraction were performed, and the results show that the proposed LOH association test has good test power and controls type 1 error well. The developed algorithms are packaged into LOHAS written in R and R GUI.

Integrative analysis of single nucleotide polymorphisms and gene expression efficiently distinguishes samples from closely related ethnic populations.

BACKGROUND: Ancestry informative markers (AIMs) are a type of genetic marker that is informative for tracing the ancestral ethnicity of individuals. Application of AIMs has gained substantial attention in population genetics, forensic sciences, and medical genetics. Single nucleotide polymorphisms (SNPs), the materials of AIMs, are useful for classifying individuals from distinct continental origins but cannot discriminate individuals with subtle genetic differences from closely related ancestral lineages. Proof-of-principle studies have shown that gene expression (GE) also is a heritable human variation that exhibits differential intensity distributions among ethnic groups. GE supplies ethnic information supplemental to SNPs; this motivated us to integrate SNP and GE markers to construct AIM panels with a reduced number of required markers and provide high accuracy in ancestry inference. Few studies in the literature have considered GE in this aspect, and none have integrated SNP and GE markers to aid classification of samples from closely related ethnic populations. RESULTS: We integrated a forward variable selection procedure into flexible discriminant analysis to identify key SNP and/or GE markers with the highest cross-validation prediction accuracy. By analyzing genome-wide SNP and/or GE markers in 210 independent samples from four ethnic groups in the HapMap II Project, we found that average testing accuracies for a majority of classification analyses were quite high, except for SNP-only analyses that were performed to discern study samples containing individuals from two close Asian populations. The average testing accuracies ranged from 0.53 to 0.79 for SNP-only analyses and increased to around 0.90 when GE markers were integrated together with SNP markers for the classification of samples from closely related Asian populations. Compared to GE-only analyses, integrative analyses of SNP and GE markers showed comparable testing accuracies and a reduced number of selected markers in AIM panels. CONCLUSIONS: Integrative analysis of SNP and GE markers provides high-accuracy and/or cost-effective classification results for assigning samples from closely related or distantly related ancestral lineages to their original ancestral populations. User-friendly BIASLESS (Biomarkers Identification and Samples Subdivision) software was developed as an efficient tool for selecting key SNP and/or GE markers and then building models for sample subdivision. BIASLESS was programmed in R and R-GUI and is available online at http://www.stat.sinica.edu.tw/hsinchou/genetics/prediction/BIASLESS.htm.

Medium-term course and outcome of schizophrenia depicted by the sixth-month subtype after an acute episode.

BACKGROUND/PURPOSE: The intermediate course of schizophrenia is a complex intertwined with the heterogeneity of the illness. This article attempts to simplify this complexity using a hypothetical tripartite based on the profile of symptoms at 6 months after acute treatment. METHODS: This is a prospective 5-year follow-up study including 163 schizophrenic inpatients in northern Taiwan comparing patients' demographic data at index admission, scores on the Positive and Negative Syndrome Scale (PANSS) for schizophrenia and social function scale measured at admission, 6-month follow-up, and annually, and scores on a neuropsychologic test battery measured approximately 5 years after recruitment. RESULTS: Patients were grouped into three subtypes based on their sixth-month symptomatology by Generalized Association Plots, designated as remitted (RM), persistent delusion/hallucination (PDH), and markedly blunting (MB) groups. These three subtypes presented with similar positive symptom profiles at recruitment, yet during follow-up, the PDH group tended to maintain the highest risk of having worse clinical symptomatology, social functioning, and neuropsychologic functioning, and the RM was the best outcome group. CONCLUSION: This three-subtype model provides a practical reference to predict medium-term outcomes by the subject's response to acute treatment and serves as a model to sort out part of the heterogeneous nature of schizophrenia that still should be examined by further psychopharmacological, neurobiological, and genetic studies.

Mixed sequence reader: a program for analyzing DNA sequences with heterozygous base calling.

The direct sequencing of PCR products generates heterozygous base-calling fluorescence chromatograms that are useful for identifying single-nucleotide polymorphisms (SNPs), insertion-deletions (indels), short tandem repeats (STRs), and paralogous genes. Indels and STRs can be easily detected using the currently available Indelligent or ShiftDetector programs, which do not search reference sequences. However, the detection of other genomic variants remains a challenge due to the lack of appropriate tools for heterozygous base-calling fluorescence chromatogram data analysis. In this study, we developed a free web-based program, Mixed Sequence Reader (MSR), which can directly analyze heterozygous base-calling fluorescence chromatogram data in .abi file format using comparisons with reference sequences. The heterozygous sequences are identified as two distinct sequences and aligned with reference sequences. Our results showed that MSR may be used to (i) physically locate indel and STR sequences and determine STR copy number by searching NCBI reference sequences; (ii) predict combinations of microsatellite patterns using the Federal Bureau of Investigation Combined DNA Index System (CODIS); (iii) determine human papilloma virus (HPV) genotypes by searching current viral databases in cases of double infections; (iv) estimate the copy number of paralogous genes, such as ß-defensin 4 (DEFB4) and its paralog HSPDP3.

Subtyping of major SARS-CoV-2 variants reveals different transmission dynamics based on 10 million genomes.

SARS-CoV-2 continues to evolve, causing waves of the pandemic. Up to May 2022, 10 million genome sequences have accumulated, which are classified into five major variants of concern. With the growing number of sequenced genomes, analysis of the big dataset has become increasingly challenging. Here we developed systematic approaches based on sets of correlated single nucleotide variations (SNVs) for comprehensive subtyping and pattern recognition of transmission dynamics. The approach outperformed single-SNV and spike-centric scans. Moreover, the derived subtypes elucidate the relationship of signature SNVs and transmission dynamics. We found that different subtypes of the same variant, including Delta and Omicron exhibited distinct temporal trajectories. For example, some Delta and Omicron subtypes did not spread rapidly, while others did. We identified sets of characteristic SNVs that appeared to enhance transmission or decrease efficacy of antibodies for some subtypes. We also identified a set of SNVs that appeared to suppress transmission or increase viral sensitivity to antibodies. For the Omicron variant, the dominant type in the world, we identified the subtypes with enhanced and suppressed transmission in an analysis of eight million genomes as of March 2022 and further confirmed the findings in a later analysis of ten million genomes as of May 2022. While the "enhancer" SNVs exhibited an enriched presence on the spike protein, the "suppressor" SNVs are mainly elsewhere. Disruption of the SNV correlation largely destroyed the enhancer-suppressor phenomena. These results suggest the importance of fine subtyping of variants, and point to potential complex interactions among SNVs.

Relationship among genetic variants, obesity traits and asthma in the Taiwan Biobank.

BACKGROUND AND OBJECTIVE: Obesity and asthma impose a heavy health and economic burden on millions of people around the world. The complex interaction between genetic traits and phenotypes caused the mechanism between obesity and asthma is still vague. This study investigates the relationship among obesity-related polygenic risk score (PRS), obesity phenotypes and the risk of having asthma. METHODS: This is a matched case-control study, with 4 controls (8288 non-asthmatic) for each case (2072 asthmatic). Data were obtained from the 2008-2015 Taiwan Biobank Database and linked to the 2000-2016 National Health Insurance Research Database. All participants were ≥30 years old with no history of cancer and had a complete questionnaire, as well as physical examination, genome-wide single nucleotide polymorphisms and clinical diagnosis data. Environmental exposure, PM2.5, was also considered. Multivariate adjusted ORs and 95% CIs were calculated using conditional logistic regression stratified by age and sex. Mediation analysis was also assessed, using a generalised linear model. RESULTS: We found that the obese phenotype was associated with significantly increased odds of asthma by approximately 26%. Four obesity-related PRS, including body mass index (OR=1.07 (1.01-1.13)), waist circumference (OR=1.10 (1.04-1.17)), central obesity as defined by waist-to-height ratio (OR=1.09 (1.03-1.15)) and general-central obesity (OR=1.06 (1.00-1.12)), were associated with increased odds of asthma. Additional independent risk factors for asthma included lower educational level, family history of asthma, certain chronic diseases and increased PM2.5 exposure. Obesity-related PRS is an indirect risk factor for asthma, the link being fully mediated by the trait of obesity. CONCLUSIONS: Obese phenotypes and obesity-related PRS are independent risk factors for having asthma in adults in the Taiwan Biobank. Overall, genetic risk for obesity increases the risk of asthma by affecting the obese phenotype.

SAQC: SNP array quality control.

BACKGROUND: Genome-wide single-nucleotide polymorphism (SNP) arrays containing hundreds of thousands of SNPs from the human genome have proven useful for studying important human genome questions. Data quality of SNP arrays plays a key role in the accuracy and precision of downstream data analyses. However, good indices for assessing data quality of SNP arrays have not yet been developed. RESULTS: We developed new quality indices to measure the quality of SNP arrays and/or DNA samples and investigated their statistical properties. The indices quantify a departure of estimated individual-level allele frequencies (AFs) from expected frequencies via standardized distances. The proposed quality indices followed lognormal distributions in several large genomic studies that we empirically evaluated. AF reference data and quality index reference data for different SNP array platforms were established based on samples from various reference populations. Furthermore, a confidence interval method based on the underlying empirical distributions of quality indices was developed to identify poor-quality SNP arrays and/or DNA samples. Analyses of authentic biological data and simulated data show that this new method is sensitive and specific for the detection of poor-quality SNP arrays and/or DNA samples. CONCLUSIONS: This study introduces new quality indices, establishes references for AFs and quality indices, and develops a detection method for poor-quality SNP arrays and/or DNA samples. We have developed a new computer program that utilizes these methods called SNP Array Quality Control (SAQC). SAQC software is written in R and R-GUI and was developed as a user-friendly tool for the visualization and evaluation of data quality of genome-wide SNP arrays. The program is available online (http://www.stat.sinica.edu.tw/hsinchou/genetics/quality/SAQC.htm).

Analysis of human meiotic recombination events with a parent-sibling tracing approach.

BACKGROUND: Meiotic recombination ensures that each child inherits distinct genetic materials from each parent, but the distribution of crossovers along meiotic chromosomes remains difficult to identify. In this study, we developed a parent-sibling tracing (PST) approach from previously reported methods to identify meiotic crossover sites of GEO GSE6754 data set. This approach requires only the single nucleotide polymorphism (SNP) data of the pedigrees of both parents and at least two of children. RESULTS: Compared to other SNP-based algorithms (identity by descent or pediSNP), fewer uninformative SNPs were derived with the use of PST. Analysis of a GEO GSE6754 data set containing 2,145 maternal and paternal meiotic events revealed that the pattern and distribution of paternal and maternal recombination sites vary along the chromosomes. Lower crossover rates near the centromeres were more prominent in males than in females. Based on analysis of repetitive sequences, we also showed that recombination hotspots are positively correlated with SINE/MIR repetitive elements and negatively correlated with LINE/L1 elements. The number of meiotic recombination events was positively correlated with the number of shorter tandem repeat sequences. CONCLUSIONS: The advantages of the PST approach include the ability to use only two-generation pedigrees with two siblings and the ability to perform gender-specific analyses of repetitive elements and tandem repeat sequences while including fewer uninformative SNP regions in the results.

Microarray labeling extension values: laboratory signatures for Affymetrix GeneChips.

Interlaboratory comparison of microarray data, even when using the same platform, imposes several challenges to scientists. RNA quality, RNA labeling efficiency, hybridization procedures and data-mining tools can all contribute variations in each laboratory. In Affymetrix GeneChips, about 11-20 different 25-mer oligonucleotides are used to measure the level of each transcript. Here, we report that 'labeling extension values (LEVs)', which are correlation coefficients between probe intensities and probe positions, are highly correlated with the gene expression levels (GEVs) on eukaryotic Affymetrix microarray data. By analyzing LEVs and GEVs in the publicly available 2414 cel files of 20 Affymetrix microarray types covering 13 species, we found that correlations between LEVs and GEVs only exist in eukaryotic RNAs, but not in prokaryotic ones. Surprisingly, Affymetrix results of the same specimens that were analyzed in different laboratories could be clearly differentiated only by LEVs, leading to the identification of 'laboratory signatures'. In the examined dataset, GSE10797, filtering out high-LEV genes did not compromise the discovery of biological processes that are constructed by differentially expressed genes. In conclusion, LEVs provide a new filtering parameter for microarray analysis of gene expression and it may improve the inter- and intralaboratory comparability of Affymetrix GeneChips data.