Deciphering genetic basis of developmental and agronomic traits by integrating high-throughput optical phenotyping and genome-wide association studies in wheat

Dissecting the genetic basis of complex traits such as dynamic growth and yield potential is a major challenge in crops. Monitoring the growth throughout growing season in a large wheat population to uncover the temporal genetic controls for plant growth and yield-related traits has so far not been explored. In this study, a diverse wheat panel composed of 288 lines was monitored by a non-invasive and high-throughput phenotyping platform to collect growth traits from seedling to grain filling stage and their relationship with yield-related traits was further explored. Whole genome re-sequencing of the panel provided 12.64 million markers for a high-resolution genome-wide association analysis using 190 image-based traits and 17 agronomic traits. A total of 8327 marker-trait associations were detected and clustered into 1605 quantitative trait loci (QTLs) including a number of known genes or QTLs. We identified 277 pleiotropic QTLs controlling multiple traits at different growth stages which revealed temporal dynamics of QTLs action on plant development and yield production in wheat. A candidate gene related to plant growth that was detected by image traits was further validated. Particularly, our study demonstrated that the yield-related traits are largely predictable using models developed based on i-traits and provide possibility for high-throughput early selection, thus to accelerate breeding process. Our study explored the genetic architecture of growth and yield-related traits by combining high-throughput phenotyping and genotyping, which further unravelled the complex and stage-specific contributions of genetic loci to optimize growth and yield in wheat.     

Density of common complex ocular traits in the aging eye: analysis of secondary traits in genome-wide association studies

Genetic association studies are identifying genetic risks for common complex ocular traits such as age-related macular degeneration (AMD). The subjects used for discovery of these loci have been largely from clinic-based, case-control studies. Typically, only the primary phenotype (e.g., AMD) being studied is systematically documented and other complex traits (e.g., affecting the eye) are largely ignored. The purpose of this study was to characterize these other or secondary complex ocular traits present in the cases and controls of clinic-based studies being used for genetic study of AMD. The records of 100 consecutive new patients (of any diagnosis) age 60 or older for which all traits affecting the eye had been recorded systematically were reviewed. The average patient had 3.5 distinct diagnoses. A subset of 10 complex traits was selected for further study because they were common and could be reliably diagnosed. The density of these 10 complex ocular traits increased by 0.017 log-traits/year (P = 0.03), ranging from a predicted 2.74 at age 60 to 4.45 at age 90. Trait-trait association was observed only between AMD and primary vitreomacular traction (P = 0.0009). Only 1% of subjects age 60 or older had no common complex traits affecting the eye. Extrapolations suggested that a study of 2000 similar subjects would have sufficient power to detect genetic association with an odds ratio of 2.0 or less for 4 of these 10 traits. In conclusion, the high prevalence of complex traits affecting the aging eye and the inherent biases in referral patterns leads to the potential for confounding by undocumented secondary traits within case-control studies. In addition to the primary trait, other common ocular phenotypes should be systematically documented in genetic association studies so that adjustments for potential trait-trait associations and other bias can be made and genetic risk variants identified in secondary analyses.     

羊重要性状全基因组关联分析研究进展

全基因组关联分析(genome-wide association study, GWAS)是一种复杂性状功能基因鉴定的分析策略,已成为挖掘畜禽重要经济性状候选基因的重要手段。随着绵羊和山羊基因组完成和公布,以及不同密度的SNP (single nucleotide polymorphism)芯片的推出并进行商业化推广,不仅大大丰富了羊标记辅助选择可利用的分子标记,而且还为开展重要性状的分子机理的探索提供了重要技术支撑。本文主要针对羊角、羊毛、羊奶、生长发育、肉质、繁殖和疾病等重要性状的GWAS研究所用的群体、主要研究方法和研究结果进行了综述,并对GWAS方法研究现状进行了归纳,以期为进一步利用GWAS进行羊的各种性状的遗传基础研究提供参考。     

SCOPE-Seq: a scalable technology for linking live cell imaging and single-cell RNA sequencing

Optically decodable beads link the identity of a sample to a measurement through an optical barcode, enabling libraries of biomolecules to be captured on beads in solution and decoded by fluorescence. This approach has been foundational to microarray, sequencing, and flow-based expression profiling technologies. We combine microfluidics with optically decodable beads and show that phenotypic analysis of living cells can be linked to single-cell sequencing. As a proof-of-concept, we demonstrate the accuracy and scalability of our tool called Single Cell Optical Phenotyping and Expression sequencing (SCOPE-Seq) to combine live cell imaging with single-cell RNA sequencing.     

MISC: missing imputation for single-cell RNA sequencing data

Background

Single-cell RNA sequencing (scRNA-seq) technology provides an effective way to study cell heterogeneity. However, due to the low capture efficiency and stochastic gene expression, scRNA-seq data often contains a high percentage of missing values. It has been showed that the missing rate can reach approximately 30% even after noise reduction. To accurately recover missing values in scRNA-seq data, we need to know where the missing data is; how much data is missing; and what are the values of these data.

Methods

To solve these three problems, we propose a novel model with a hybrid machine learning method, namely, missing imputation for single-cell RNA-seq (MISC). To solve the first problem, we transformed it to a binary classification problem on the RNA-seq expression matrix. Then, for the second problem, we searched for the intersection of the classification results, zero-inflated model and false negative model results. Finally, we used the regression model to recover the data in the missing elements.

Results

We compared the raw data without imputation, the mean-smooth neighbor cell trajectory, MISC on chronic myeloid leukemia data (CML), the primary somatosensory cortex and the hippocampal CA1 region of mouse brain cells. On the CML data, MISC discovered a trajectory branch from the CP-CML to the BC-CML, which provides direct evidence of evolution from CP to BC stem cells. On the mouse brain data, MISC clearly divides the pyramidal CA1 into different branches, and it is direct evidence of pyramidal CA1 in the subpopulations. In the meantime, with MISC, the oligodendrocyte cells became an independent group with an apparent boundary.

Conclusions

Our results showed that the MISC model improved the cell type classification and could be instrumental to study cellular heterogeneity. Overall, MISC is a robust missing data imputation model for single-cell RNA-seq data.

     

Genome-wide association studies for common diseases and complex traits

Genetic factors strongly affect susceptibility to common diseases and also influence disease-related quantitative traits. Identifying the relevant genes has been difficult, in part because each causal gene only makes a small contribution to overall heritability. Genetic association studies offer a potentially powerful approach for mapping causal genes with modest effects, but are limited because only a small number of genes can be studied at a time. Genome-wide association studies will soon become possible, and could open new frontiers in our understanding and treatment of disease. However, the execution and analysis of such studies will require great care.     

scAB detects multiresolution cell states with clinical significance by integrating single-cell genomics and bulk sequencing data

Although single-cell sequencing has provided a powerful tool to deconvolute cellular heterogeneity of diseases like cancer, extrapolating clinical significance or identifying clinically-relevant cells remains challenging. Here, we propose a novel computational method scAB, which integrates single-cell genomics data with clinically annotated bulk sequencing data via a knowledge- and graph-guided matrix factorization model. Once combined, scAB provides a coarse- and fine-grain multiresolution perspective of phenotype-associated cell states and prognostic signatures previously not visible by single-cell genomics. We use scAB to enhance live cancer single-cell RNA-seq data, identifying clinically-relevant previously unrecognized cancer and stromal cell subsets whose signatures show a stronger poor-survival association. The identified fine-grain cell subsets are associated with distinct cancer hallmarks and prognosis power. Furthermore, scAB demonstrates its utility as a biomarker identification tool, with the ability to predict immunotherapy, drug responses and survival when applied to melanoma single-cell RNA-seq datasets and glioma single-cell ATAC-seq datasets. Across multiple single-cell and bulk datasets from different cancer types, we also demonstrate the superior performance of scAB in generating prognosis signatures and survival predictions over existing models. Overall, scAB provides an efficient tool for prioritizing clinically-relevant cell subsets and predictive signatures, utilizing large publicly available databases to improve prognosis and treatments.     

GWAPower: a statistical power calculation software for genome-wide association studies with quantitative traits

Background

In designing genome-wide association (GWA) studies it is important to calculate statistical power. General statistical power calculation procedures for quantitative measures often require information concerning summary statistics of distributions such as mean and variance. However, with genetic studies, the effect size of quantitative traits is traditionally expressed as heritability, a quantity defined as the amount of phenotypic variation in the population that can be ascribed to the genetic variants among individuals. Heritability is hard to transform into summary statistics. Therefore, general power calculation procedures cannot be used directly in GWA studies. The development of appropriate statistical methods and a user-friendly software package to address this problem would be welcomed.

Results

This paper presents GWAPower, a statistical software package of power calculation designed for GWA studies with quantitative traits, where genetic effect is defined as heritability. Based on several popular one-degree-of-freedom genetic models, this method avoids the need to specify the non-centrality parameter of the F-distribution under the alternative hypothesis. Therefore, it can use heritability information directly without approximation. In GWAPower, the power calculation can be easily adjusted for adding covariates and linkage disequilibrium information. An example is provided to illustrate GWAPower, followed by discussions.

Conclusions

GWAPower is a user-friendly free software package for calculating statistical power based on heritability in GWA studies with quantitative traits. The software is freely available at: http://dl.dropbox.com/u/10502931/GWAPower.zip     

Multiple locus genome-wide association studies for important economic traits of oil palm

Palm oil has a balanced fatty acid composition and has no trans fat. As a result, its use in food has increased as food-labeling laws have changed to specify trans fat content. Increasing oil production is the main goal in oil palm breeding. Genetic mapping and genomic studies in palm trees are necessary to understand the genetic architecture of economic traits of importance for palm oil production. To help achieve this, we sampled 422 oil palms from MPOB (Malaysian Palm Oil Board)­Angola germplasm collection and measured 13 economic traits from these palms. Multi-locus genome-wide association studies (GWAS) were conducted using least absolute shrinkage and selection operator (LASSO) and genome-wide efficient mixed model analysis. We identified 19 quantitative trait loci (QTLs) for 8 traits. Of these, four Angola-specific QTLs associated with bunch components were detected on chromosomes 4, 8, and 11. These QTLs are potentially useful for introgression of desirable genes from the Angola palms to advanced breeding populations for improvement of bunch and oil yield traits. The majority of the QTLs were detected by LASSO-A, in which the p values of individual markers were calculated based on bootstrapped standard errors. Many of the detected QTLs are nearby known QTLs detected from linkage studies reported by other research groups. We also conducted genomic selection (GS) for the 13 traits and concluded that GS can be an effective tool for oil palm breeding. This is the first GWAS and GS study conducted on oil palm germplasm from Angola, and the results can be very useful in oil palm genetic studies and breeding.     

Progress and promise of genome-wide association studies for human complex trait genetics

Enormous progress in mapping complex traits in humans has been made in the last 5 yr. There has been early success for prevalent diseases with complex phenotypes. These studies have demonstrated clearly that, while complex traits differ in their underlying genetic architectures, for many common disorders the predominant pattern is that of many loci, individually with small effects on phenotype. For some traits, loci of large effect have been identified. For almost all complex traits studied in humans, the sum of the identified genetic effects comprises only a portion, generally less than half, of the estimated trait heritability. A variety of hypotheses have been proposed to explain why this might be the case, including untested rare variants, and gene-gene and gene-environment interaction. Effort is currently being directed toward implementation of novel analytic approaches and testing rare variants for association with complex traits using imputed variants from the publicly available 1000 Genomes Project resequencing data and from direct resequencing of clinical samples. Through integration with annotations and functional genomic data as well as by in vitro and in vivo experimentation, mapping studies continue to characterize functional variants associated with complex traits and address fundamental issues such as epistasis and pleiotropy. This review focuses primarily on the ways in which genome-wide association studies (GWASs) have revolutionized the field of human quantitative genetics.     

TedSim: temporal dynamics simulation of single-cell RNA sequencing data and cell division history

Recently, lineage tracing technology using CRISPR/Cas9 genome editing has enabled simultaneous readouts of gene expressions and lineage barcodes, which allows for the reconstruction of the cell division tree and makes it possible to reconstruct ancestral cell types and trace the origin of each cell type. Meanwhile, trajectory inference methods are widely used to infer cell trajectories and pseudotime in a dynamic process using gene expression data of present-day cells. Here, we present TedSim (single-cell temporal dynamics simulator), which simulates the cell division events from the root cell to present-day cells, simultaneously generating two data modalities for each single cell: the lineage barcode and gene expression data. TedSim is a framework that connects the two problems: lineage tracing and trajectory inference. Using TedSim, we conducted analysis to show that (i) TedSim generates realistic gene expression and barcode data, as well as realistic relationships between these two data modalities; (ii) trajectory inference methods can recover the underlying cell state transition mechanism with balanced cell type compositions; and (iii) integrating gene expression and barcode data can provide more insights into the temporal dynamics in cell differentiation compared to using only one type of data, but better integration methods need to be developed.     

Genetic association studies of complex traits: design and analysis issues

Most common diseases and many important quantitative traits are complex genetic traits, with multiple genetic and environmental variables contributing to the observed phenotype. Because of the multi-factorial nature of complex traits, each individual genetic variant generally has only a modest effect, and the interaction of genetic variants with each other or with environmental factors can potentially be quite important in determining the observed phenotype. It remains largely unknown what sort of genetic variants explain inherited variation in complex traits, but recent evidence suggests that common genetic variants will explain at least some of the inherited variation in susceptibility to common disease. Genetic association studies, in which the allele or genotype frequencies at markers are determined in affected individuals and compared with those of controls (either population- or family-based), may be an effective approach to detecting the effects of common variants with modest effects. With the explosion in single nucleotide polymorphism (SNP) discovery and genotyping technologies, large-scale association studies have become feasible, and small-scale association studies have become plentiful. We review the different types of association studies and discuss issues that are important to consider when performing and interpreting association studies of complex genetic traits. Heritable and accurately measured phenotypes, carefully matched large samples, well-chosen genetic markers, and adequate standards in genotyping, analysis, and interpretation are all integral parts of a high-quality association study.     

RS-SNP: a random-set method for genome-wide association studies

Background

Common carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.

Result

To develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.

Conclusion

BAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.     

SNPTransformer: a lightweight toolkit for genome-wide association studies

High-throughput genotyping chips have produced huge datasets for genome-wide association studies(GWAS)that have contributed greatly to discovering susceptibility genes for complex diseases.There are two strategies for performing data analysis for GWAS.One strategy is to use open-source or commercial packages that are designed for GWAS.The other is to take advantage of classic genetic programs with specific functions,such as linkage disequilibrium mapping,haplotype inference and transmission disequilibrium tests.However,most classic programs that are available are not suitable for analyzing chip data directly and require custom-made input,which results in the inconvenience of converting raw genotyping files into various data formats.We developed a powerful,user-friendly,lightweight program named SNPTransformer for GWAS that includes five major modules (Transformer,Operator,Previewer,Coder and Simulator).The toolkit not only works for transforming the genotyping files into ten input formats for use with classic genetics packages,but also carries out useful functions such as relational operations on IDs,previewing data files,recoding data formats and simulating marker files,among other functions.It bridges upstream raw genotyping data with downstream genetic programs,and can act as an in-hand toolkit for human geneticists,especially for non-programmers.SNPTransformer is freely available at http://snptransformer.sourceforge.net.     

Modelling human regulatory variation in mouse: finding the function in genome-wide association studies and whole-genome sequencing

An increasing body of literature from genome-wide association studies and human whole-genome sequencing highlights the identification of large numbers of candidate regulatory variants of potential therapeutic interest in numerous diseases. Our relatively poor understanding of the functions of non-coding genomic sequence, and the slow and laborious process of experimental validation of the functional significance of human regulatory variants, limits our ability to fully benefit from this information in our efforts to comprehend human disease. Humanized mouse models (HuMMs), in which human genes are introduced into the mouse, suggest an approach to this problem. In the past, HuMMs have been used successfully to study human disease variants; e.g., the complex genetic condition arising from Down syndrome, common monogenic disorders such as Huntington disease and β-thalassemia, and cancer susceptibility genes such as BRCA1. In this commentary, we highlight a novel method for high-throughput single-copy site-specific generation of HuMMs entitled High-throughput Human Genes on the X Chromosome (HuGX). This method can be applied to most human genes for which a bacterial artificial chromosome (BAC) construct can be derived and a mouse-null allele exists. This strategy comprises (1) the use of recombineering technology to create a human variant-harbouring BAC, (2) knock-in of this BAC into the mouse genome using Hprt docking technology, and (3) allele comparison by interspecies complementation. We demonstrate the throughput of the HuGX method by generating a series of seven different alleles for the human NR2E1 gene at Hprt. In future challenges, we consider the current limitations of experimental approaches and call for a concerted effort by the genetics community, for both human and mouse, to solve the challenge of the functional analysis of human regulatory variation.