A comparison between whole transcript and 3' RNA sequencing methods using Kapa and Lexogen library preparation methods.

BACKGROUND: 3' RNA sequencing provides an alternative to whole transcript analysis. However, we do not know a priori the relative advantage of each method. Thus, a comprehensive comparison between the whole transcript and the 3' method is needed to determine their relative merits. To this end, we used two commercially available library preparation kits, the KAPA Stranded mRNA-Seq kit (traditional method) and the Lexogen QuantSeq 3' mRNA-Seq kit (3' method), to prepare libraries from mouse liver RNA. We then sequenced and analyzed the libraries to determine the advantages and disadvantages of these two approaches. RESULTS: We found that the traditional whole transcript method and the 3' RNA-Seq method had similar levels of reproducibility. As expected, the whole transcript method assigned more reads to longer transcripts, while the 3' method assigned roughly equal numbers of reads to transcripts regardless of their lengths. We found that the 3' RNA-Seq method detected more short transcripts than the whole transcript method. With regard to differential expression analysis, we found that the whole transcript method detected more differentially expressed genes, regardless of the level of sequencing depth. CONCLUSIONS: The 3' RNA-Seq method was better able to detect short transcripts, while the whole transcript RNA-Seq was able to detect more differentially expressed genes. Thus, both approaches have relative advantages and should be selected based on the goals of the experiment.

Systematic inference of copy-number genotypes from personal genome sequencing data reveals extensive olfactory receptor gene content diversity.

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95-99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting ~15% and ~20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies involving high-throughput sequencing.

High-resolution copy-number variation map reflects human olfactory receptor diversity and evolution.

Olfactory receptors (ORs), which are involved in odorant recognition, form the largest mammalian protein superfamily. The genomic content of OR genes is considerably reduced in humans, as reflected by the relatively small repertoire size and the high fraction ( approximately 55%) of human pseudogenes. Since several recent low-resolution surveys suggested that OR genomic loci are frequently affected by copy-number variants (CNVs), we hypothesized that CNVs may play an important role in the evolution of the human olfactory repertoire. We used high-resolution oligonucleotide tiling microarrays to detect CNVs across 851 OR gene and pseudogene loci. Examining genomic DNA from 25 individuals with ancestry from three populations, we identified 93 OR gene loci and 151 pseudogene loci affected by CNVs, generating a mosaic of OR dosages across persons. Our data suggest that approximately 50% of the CNVs involve more than one OR, with the largest CNV spanning 11 loci. In contrast to earlier reports, we observe that CNVs are more frequent among OR pseudogenes than among intact genes, presumably due to both selective constraints and CNV formation biases. Furthermore, our results show an enrichment of CNVs among ORs with a close human paralog or lacking a one-to-one ortholog in chimpanzee. Interestingly, among the latter we observed an enrichment in CNV losses over gains, a finding potentially related to the known diminution of the human OR repertoire. Quantitative PCR experiments performed for 122 sampled ORs agreed well with the microarray results and uncovered 23 additional CNVs. Importantly, these experiments allowed us to uncover nine common deletion alleles that affect 15 OR genes and five pseudogenes. Comparison to the chimpanzee reference genome revealed that all of the deletion alleles are human derived, therefore indicating a profound effect of human-specific deletions on the individual OR gene content. Furthermore, these deletion alleles may be used in future genetic association studies of olfactory inter-individual differences.

Genetic elucidation of human hyperosmia to isovaleric acid.

The genetic basis of odorant-specific variations in human olfactory thresholds, and in particular of enhanced odorant sensitivity (hyperosmia), remains largely unknown. Olfactory receptor (OR) segregating pseudogenes, displaying both functional and nonfunctional alleles in humans, are excellent candidates to underlie these differences in olfactory sensitivity. To explore this hypothesis, we examined the association between olfactory detection threshold phenotypes of four odorants and segregating pseudogene genotypes of 43 ORs genome-wide. A strong association signal was observed between the single nucleotide polymorphism variants in OR11H7P and sensitivity to the odorant isovaleric acid. This association was largely due to the low frequency of homozygous pseudogenized genotype in individuals with specific hyperosmia to this odorant, implying a possible functional role of OR11H7P in isovaleric acid detection. This predicted receptor-ligand functional relationship was further verified using the Xenopus oocyte expression system, whereby the intact allele of OR11H7P exhibited a response to isovaleric acid. Notably, we also uncovered another mechanism affecting general olfactory acuity that manifested as a significant inter-odorant threshold concordance, resulting in an overrepresentation of individuals who were hyperosmic to several odorants. An involvement of polymorphisms in other downstream transduction genes is one possible explanation for this observation. Thus, human hyperosmia to isovaleric acid is a complex trait, contributed to by both receptor and other mechanisms in the olfactory signaling pathway.

Tissue-specific pathways and networks underlying sexual dimorphism in non-alcoholic fatty liver disease.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) encompasses benign steatosis and more severe conditions such as non-alcoholic steatohepatitis (NASH), cirrhosis, and liver cancer. This chronic liver disease has a poorly understood etiology and demonstrates sexual dimorphisms. We aim to examine the molecular mechanisms underlying sexual dimorphisms in NAFLD pathogenesis through a comprehensive multi-omics study. We integrated genomics (DNA variations), transcriptomics of liver and adipose tissue, and phenotypic data of NAFLD derived from female mice of ~ 100 strains included in the hybrid mouse diversity panel (HMDP) and compared the NAFLD molecular pathways and gene networks between sexes. RESULTS: We identified both shared and sex-specific biological processes for NAFLD. Adaptive immunity, branched chain amino acid metabolism, oxidative phosphorylation, and cell cycle/apoptosis were shared between sexes. Among the sex-specific pathways were vitamins and cofactors metabolism and ion channel transport for females, and phospholipid, lysophospholipid, and phosphatidylinositol metabolism and insulin signaling for males. Additionally, numerous lipid and insulin-related pathways and inflammatory processes in the adipose and liver tissue appeared to show more prominent association with NAFLD in male HMDP. Using data-driven network modeling, we identified plausible sex-specific and tissue-specific regulatory genes as well as those that are shared between sexes. These key regulators orchestrate the NAFLD pathways in a sex- and tissue-specific manner. Gonadectomy experiments support that sex hormones may partially underlie the sexually dimorphic genes and pathways involved in NAFLD. CONCLUSIONS: Our multi-omics integrative study reveals sex- and tissue-specific genes, processes, and networks underlying sexual dimorphism in NAFLD and may facilitate sex-specific precision medicine.

Multi-omics approaches to disease.

High-throughput technologies have revolutionized medical research. The advent of genotyping arrays enabled large-scale genome-wide association studies and methods for examining global transcript levels, which gave rise to the field of "integrative genetics". Other omics technologies, such as proteomics and metabolomics, are now often incorporated into the everyday methodology of biological researchers. In this review, we provide an overview of such omics technologies and focus on methods for their integration across multiple omics layers. As compared to studies of a single omics type, multi-omics offers the opportunity to understand the flow of information that underlies disease.

Analysis of genetic polymorphisms in acetylcholinesterase as reflected in different populations.

Acetylcholinesterase (AChE) plays a crucial physiological role in termination of impulse transmission at cholinergic synapses through rapid hydrolysis of acetylcholine. In addition, it was implicated in amyloid plaque formation, a hallmark of Alzheimer's disease (AD), and most of the drugs used in AD treatment are AChE inhibitors. Thus ACHE is an obvious candidate gene for pharmacogenetic study of AD treatment. However, AChE is a highly conserved molecule, and only a few naturally occurring genetic polymorphisms have been reported in the human gene. The goals of this study were to make a systematic effort to identify natural single nucleotide polymorphisms (SNPs) in the human ACHE gene, and to reveal their population specific architecture. To this end, the genomic coding sequences for AChE of 96 unrelated control individuals from three distinct ethnic groups, African Americans, Ashkenazi Jews and Israeli Arabs, were analyzed. Thirteen ACHE SNPs were identified, ten of which are newly described, and five of which should produce amino-acid substitutions (Arg34Gln, Gly57Arg, Glu344Gly, His353Asn and Pro592Arg). Population frequencies of 11 of the 13 SNPs were established in four different populations, African Americans, Ashkenazi Jews, Sephardic Jews and Israeli Arabs; 17 haplotypes and 5 ethno-specific alleles were identified, and a cladogram of ACHE haplotypes was constructed. Among the SNPs resulting in an amino-acid substitution, three are within the mature protein, mapping on its external surface; they are thus unlikely to affect its catalytic properties, yet could have antigenic consequences or affect putative protein-protein interactions. Furthermore, the newly identified SNPs open the door to a study of the possible association of AChE with deleterious phenotypes - such as adverse drug responses to AChE inhibitors employed in treatment of AD patients and hypersensitivity to pesticides.

A paradigm for single nucleotide polymorphism analysis: the case of the acetylcholinesterase gene.

Acetylcholinesterase (AChE) plays a crucial physiological role in termination of impulse transmission at cholinergic synapses through rapid hydrolysis of acetylcholine. It is a highly conserved molecule, and only a few naturally occurring genetic polymorphisms have been reported in the human gene. The goal of the present study was to make a systematic effort to identify natural single nucleotide polymorphisms (SNPs) in the human ACHE gene. To this end, the genomic coding sequences for acetylcholinesterase of 96 unrelated control individuals from three distinct ethnic groups were analyzed. A total of 13 ACHE SNPs were identified, 10 of which are newly described, and five that should produce amino acid substitutions [c.101G>A (p.Arg34Gln), c.169G>A (p.Gly57Arg), c.1031A>G (p.Glu344Gly), c.1057C>A (p.His353Asn), and c.1775C>G (p.Pro592Arg)]. Population frequencies of 11 of the 13 SNPs were established in four different populations: African Americans, Ashkenazi Jews, Sephardic Jews, and Israeli Arabs; 15 haplotypes and five ethnospecific alleles were identified. The low number of SNPs identified until now in the ACHE gene is ascribed to technical hurdles arising from the high GC content and the presence of numerous repeat sequences, and does not reflect its intrinsic heterozygosity. Among the SNPs resulting in an amino acid substitution, three are within the mature protein, mapping on its external surface: they are thus unlikely to affect its catalytic properties, yet could have antigenic consequences or affect putative protein-protein interactions. Furthermore, the newly identified SNPs open the door to a study of the possible association of AChE with deleterious phenotypes-such as adverse drug responses to AChE inhibitors employed in treatment of Alzheimer patients and hypersensitivity to pesticides.

Limited RNA editing in exons of mouse liver and adipose.

Several studies have investigated RNA-DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.