Identification and quantitation of multiple variants in RNA virus genomes.

The goal of the study was to identify and characterize RNA virus variants containing mutations spread over genomic distances >5 kb. As proof of concept, high-quality viral RNA of the Dengue 2 component of Takeda's tetravalent dengue vaccine candidate (TDV-2) was used to develop a reverse transcription-polymerase chain reaction protocol to amplify a ∼5.3 kb cDNA segment that contains the three genetic determinants of TDV-2 attenuation. Unique molecular identifiers were incorporated into each viral cDNA molecule for PacBio library preparation to improve the quantitative precision of the observed variants at the attenuation loci. Following assay optimization, PacBio long-read sequencing was validated with multiple clone-derived TDV-2 revertant variants and four complex revertant mixtures containing various compositions of TDV-2 and revertant viruses. PacBio sequencing analysis correctly identified and quantified variant composition in all tested samples, demonstrating that TDV-2 revertants could be identified and characterized and supporting the use of this method in the differentiation and quantification of complex variants of other RNA viruses. Long-read sequencing can identify complex RNA virus variants containing multiple mutations on a single-genome molecule, which is useful for in-depth genetic stability and revertant detection of live-attenuated viral vaccines, as well as research in virus evolution to reveal mechanisms of immune evasion and host cell adaption.

Annotated genome sequence of a fast-growing diploid clone of red alder (Alnus rubra Bong.).

Red alder (Alnus rubra Bong.) is an ecologically significant and important fast-growing commercial tree species native to western coastal and riparian regions of North America, having highly desirable wood, pigment, and medicinal properties. We have sequenced the genome of a rapidly growing clone. The assembly is nearly complete, containing the full complement of expected genes. This supports our objectives of identifying and studying genes and pathways involved in nitrogen-fixing symbiosis and those related to secondary metabolites that underlie red alder's many interesting defense, pigmentation, and wood quality traits. We established that this clone is most likely diploid and identified a set of SNPs that will have utility in future breeding and selection endeavors, as well as in ongoing population studies. We have added a well-characterized genome to others from the order Fagales. In particular, it improves significantly upon the only other published alder genome sequence, that of Alnus glutinosa. Our work initiated a detailed comparative analysis of members of the order Fagales and established some similarities with previous reports in this clade, suggesting a biased retention of certain gene functions in the vestiges of an ancient genome duplication when compared with more recent tandem duplications.

De novo sequencing and native mass spectrometry revealed hetero-association of dirigent protein homologs and potential interacting proteins in Forsythia × intermedia.

The discovery of dirigent proteins (DPs) and their functions in plant phenol biochemistry was made over two decades ago with Forsythia × intermedia. Stereo-selective, DP-guided, monolignol-derived radical coupling in vitro was then reported to afford the optically active lignan, (+)-pinoresinol from coniferyl alcohol, provided one-electron oxidase/oxidant capacity was present. It later became evident that DPs have several distinct sub-families, presumably with different functions. Some known DPs require other essential enzymes/proteins (e.g. oxidases) for their functions. However, the lack of a fully sequenced genome for Forsythia × intermedia made it difficult to profile other components co-purified with the (+)-pinoresinol forming DP. Herein, we used an integrated bottom-up, top-down, and native mass spectrometry (MS) approach to de novo sequence the extracted proteins via adaptation of our initial report of DP solubilization and purification. Using publicly available transcriptome and genomic data from closely related species, we identified 14 proteins that were putatively associated with either DP function or the cell wall. Although their co-occurrence after extraction and chromatographic separation is suggestive for potential protein-protein interactions, none were found to form stable protein complexes with DPs in native MS under the specific experimental conditions we have explored. Interestingly, two new DP homologs were found and they formed hetero-trimers. Molecular dynamics simulations suggested that similar hetero-trimers were possible between Arabidopsis DP homologs with comparable sequence similarities. Nevertheless, our integrated mass spectrometry method development helped prepare for future investigations directed to the discovery of novel proteins and protein-protein interactions. These advantages can be highly beneficial for plant and microbial research where fully sequenced genomes may not be readily available.

Contiguous Genome Sequence of Frankia sp. Strain ArI3, Isolated from Root Nodules of Alnus rubra Bong.

We report the genome sequence of Frankia sp. strain ArI3, recovered as a single contig from one run of the Oxford Nanopore Technologies (ONT) MinION instrument. The genome has a G+C content of 72%, is 7,541,222 bp long, and contains 5,427 predicted protein-coding genes.

New Insights Into Lignification via Network and Multi-Omics Analyses of Arogenate Dehydratase Knock-Out Mutants in Arabidopsis thaliana.

Multiple Arabidopsis arogenate dehydratase (ADT) knock-out (KO) mutants, with phenotypes having variable lignin levels (up to circa 70% reduction), were studied to investigate how differential reductions in ADTs perturb its overall plant systems biology. Integrated "omics" analyses (metabolome, transcriptome, and proteome) of wild type (WT), single and multiple ADT KO lines were conducted. Transcriptome and proteome data were collapsed into gene ortholog (GO) data, with this allowing for enzymatic reaction and metabolome cross-comparisons to uncover dominant or likely metabolic biosynthesis reactions affected. Network analysis of enzymes-highly correlated to stem lignin levels-deduced the involvement of novel putative lignin related proteins or processes. These included those associated with ribosomes, the spliceosome, mRNA transport, aminoacyl tRNA biosynthesis, and phosphorylation. While prior work helped explain lignin biosynthesis regulation at the transcriptional level, our data here provide support for a new hypothesis that there are additional post-transcriptional and translational level processes that need to be considered. These findings are anticipated to lead to development of more accurate depictions of lignin/phenylpropanoid biosynthesis models in situ, with new protein targets identified for further biochemical analysis and/or plant bioengineering. Additionally, using KEGG defined functional categorization of proteomics and transcriptomics analyses, we detected significant changes to glucosinolate, α-linolenic acid, nitrogen, carotenoid, aromatic amino acid, phenylpropanoid, and photosynthesis-related metabolic pathways in ADT KO mutants. Metabolomics results also revealed that putative carotenoid and galactolipid levels were generally increased in amount, whereas many glucosinolates and phenylpropanoids (including flavonoids and lignans) were decreased in the KO mutants.

Unique Molecular Identifiers reveal a novel sequencing artefact with implications for RNA-Seq based gene expression analysis.

Attaching Unique Molecular Identifiers (UMI) to RNA molecules in the first step of sequencing library preparation establishes a distinct identity for each input molecule. This makes it possible to eliminate the effects of PCR amplification bias, which is particularly important where many PCR cycles are required, for example, in single cell studies. After PCR, molecules sharing a UMI are assumed to be derived from the same input molecule. In our single cell RNA-Seq studies of Physcomitrella patens, we discovered that reads sharing a UMI, and therefore presumed to be derived from the same mRNA molecule, frequently map to different, but closely spaced locations. This behaviour occurs in all such libraries that we have produced, and in multiple other UMI-containing RNA-Seq data sets in the public domain. This apparent paradox, that reads of identical origin map to distinct genomic coordinates may be partially explained by PCR stutter, which is often seen in low-entropy templates and those containing simple tandem repeats. In the absence of UMI this artefact is undetectable. We show that the common assumption that sequence reads having different mapping coordinates are derived from different starting molecules does not hold. Unless taken into account, this artefact is likely to result in over-estimation of certain transcript abundances, depending on the counting method employed.

Draft Genome Sequence of a Gordonia sp. Isolated from the Soil of a Red Alder Plant.

A Gordonia species was cultured from soil of a red alder (Alnus rubra) plant. Here we present the assembled and annotated genome sequence to aid investigations into the potential of this organism as a symbiont and comparative studies of the genus Gordonia.

Gene-enriched draft genome of the cattle tick Rhipicephalus microplus: assembly by the hybrid Pacific Biosciences/Illumina approach enabled analysis of the highly repetitive genome.

The genome of the cattle tick Rhipicephalus microplus, an ectoparasite with global distribution, is estimated to be 7.1Gbp in length and consists of approximately 70% repetitive DNA. We report the draft assembly of a tick genome that utilized a hybrid sequencing and assembly approach to capture the repetitive fractions of the genome. Our hybrid approach produced an assembly consisting of 2.0Gbp represented in 195,170 scaffolds with a N50 of 60,284bp. The Rmi v2.0 assembly is 51.46% repetitive with a large fraction of unclassified repeats, short interspersed elements, long interspersed elements and long terminal repeats. We identified 38,827 putative R. microplus gene loci, of which 24,758 were protein coding genes (≥100 amino acids). OrthoMCL comparative analysis against 11 selected species including insects and vertebrates identified 10,835 and 3,423 protein coding gene loci that are unique to R. microplus or common to both R. microplus and Ixodes scapularis ticks, respectively. We identified 191 microRNA loci, of which 168 have similarity to known miRNAs and 23 represent novel miRNA families. We identified the genomic loci of several highly divergent R. microplus esterases with sequence similarity to acetylcholinesterase. Additionally we report the finding of a novel cytochrome P450 CYP41 homolog that shows similar protein folding structures to known CYP41 proteins known to be involved in acaricide resistance.

Probing the evolution, ecology and physiology of marine protists using transcriptomics.

Protists, which are single-celled eukaryotes, critically influence the ecology and chemistry of marine ecosystems, but genome-based studies of these organisms have lagged behind those of other microorganisms. However, recent transcriptomic studies of cultured species, complemented by meta-omics analyses of natural communities, have increased the amount of genetic information available for poorly represented branches on the tree of eukaryotic life. This information is providing insights into the adaptations and interactions between protists and other microorganisms and macroorganisms, but many of the genes sequenced show no similarity to sequences currently available in public databases. A better understanding of these newly discovered genes will lead to a deeper appreciation of the functional diversity and metabolic processes in the ocean. In this Review, we summarize recent developments in our understanding of the ecology, physiology and evolution of protists, derived from transcriptomic studies of cultured strains and natural communities, and discuss how these novel large-scale genetic datasets will be used in the future.

The mitochondrial genome of a Texas outbreak strain of the cattle tick, Rhipicephalus (Boophilus) microplus, derived from whole genome sequencing Pacific Biosciences and Illumina reads.

The cattle fever tick, Rhipicephalus (Boophilus) microplus is one of the most significant medical veterinary pests in the world, vectoring several serious livestock diseases negatively impacting agricultural economies of tropical and subtropical countries around the world. In our study, we assembled the complete R. microplus mitochondrial genome from Illumina and Pac Bio sequencing reads obtained from the ongoing R. microplus (Deutsch strain from Texas, USA) genome sequencing project. We compared the Deutsch strain mitogenome to the mitogenome from a Brazilian R. microplus and from an Australian cattle tick that has recently been taxonomically designated as Rhipicephalus australis after previously being considered R. microplus. The sequence divergence of the Texas and Australia ticks is much higher than the divergence between the Texas and Brazil ticks. This is consistent with the idea that the Australian ticks are distinct from the R. microplus of the Americas.

Composition, Diversity and Abundance of Gut Microbiome in Prediabetes and Type 2 Diabetes.

Association between type 2 diabetes (T2DM) and compositional changes in the gut micro biota is established, however little is known about the dysbiosis in early stages of Prediabetes (preDM). The purpose of this investigation is to elucidate the characteristics of the gut micro biome in preDM and T2DM, compared to Non-Diabetic (nonDM) subjects. Forty nine subjects were recruited for this study, 15 nonDM, 20 preDM and 14 T2DM. Bacterial community composition and diversity were investigated in fecal DNA samples using Illumina sequencing of the V4 region within the 16S rRNA gene. The five most abundant phyla identified were: Bacteroidetes, Firmicutes, Proteobacteria, Verrucomicrobia, and Actinobacteria. Class Chloracido bacteria was increased in preDM compared to T2DM (p = 0.04). An unknown genus from family Pseudonocardiaceae was significantly present in preDM group compared to the others (p = 0.04). Genus Collinsella, and an unknown genus belonging to family Enterobacteriaceae were both found to be significantly increased in T2DM compared to the other groups (Collinsella, and p = 0.03, Enterobacteriaceae genus p = 0.02). PERMANOVA and Mantel tests performed did not reveal a relationship between overall composition and diagnosis group or HbA1C level. This study identified dysbiosis associated with both preDM and T2DM, specifically at the class and genus levels suggesting that earlier treatment in preDM could possibly have an impact on the intestinal micro flora transitioning to T2DM.

Deep Panning: steps towards probing the IgOme.

BACKGROUND: Polyclonal serum consists of vast collections of antibodies, products of differentiated B-cells. The spectrum of antibody specificities is dynamic and varies with age, physiology, and exposure to pathological insults. The complete repertoire of antibody specificities in blood, the IgOme, is therefore an extraordinarily rich source of information-a molecular record of previous encounters as well as a status report of current immune activity. The ability to profile antibody specificities of polyclonal serum at exceptionally high resolution has been an important and serious challenge which can now be overcome. METHODOLOGY/PRINCIPAL FINDINGS: Here we illustrate the application of Deep Panning, a method that combines the flexibility of combinatorial phage display of random peptides with the power of high-throughput deep sequencing. Deep Panning is first applied to evaluate the quality and diversity of naïve random peptide libraries. The production of very large data sets, hundreds of thousands of peptides, has revealed unexpected properties of combinatorial random peptide libraries and indicates correctives to ensure the quality of the libraries generated. Next, Deep Panning is used to analyze a model monoclonal antibody in addition to allowing one to follow the dynamics of biopanning and peptide selection. Finally Deep Panning is applied to profile polyclonal sera derived from HIV infected individuals. CONCLUSIONS/SIGNIFICANCE: The ability to generate and characterize hundreds of thousands of affinity-selected peptides creates an effective means towards the interrogation of the IgOme and understanding of the humoral response to disease. Deep Panning should open the door to new possibilities for serological diagnostics, vaccine design and the discovery of the correlates of immunity to emerging infectious agents.

Carrier testing for severe childhood recessive diseases by next-generation sequencing.

Of 7028 disorders with suspected Mendelian inheritance, 1139 are recessive and have an established molecular basis. Although individually uncommon, Mendelian diseases collectively account for ~20% of infant mortality and ~10% of pediatric hospitalizations. Preconception screening, together with genetic counseling of carriers, has resulted in remarkable declines in the incidence of several severe recessive diseases including Tay-Sachs disease and cystic fibrosis. However, extension of preconception screening to most severe disease genes has hitherto been impractical. Here, we report a preconception carrier screen for 448 severe recessive childhood diseases. Rather than costly, complete sequencing of the human genome, 7717 regions from 437 target genes were enriched by hybrid capture or microdroplet polymerase chain reaction, sequenced by next-generation sequencing (NGS) to a depth of up to 2.7 gigabases, and assessed with stringent bioinformatic filters. At a resultant 160x average target coverage, 93% of nucleotides had at least 20x coverage, and mutation detection/genotyping had ~95% sensitivity and ~100% specificity for substitution, insertion/deletion, splicing, and gross deletion mutations and single-nucleotide polymorphisms. In 104 unrelated DNA samples, the average genomic carrier burden for severe pediatric recessive mutations was 2.8 and ranged from 0 to 7. The distribution of mutations among sequenced samples appeared random. Twenty-seven percent of mutations cited in the literature were found to be common polymorphisms or misannotated, underscoring the need for better mutation databases as part of a comprehensive carrier testing strategy. Given the magnitude of carrier burden and the lower cost of testing compared to treating these conditions, carrier screening by NGS made available to the general population may be an economical way to reduce the incidence of and ameliorate suffering associated with severe recessive childhood disorders.

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis.

Monozygotic or 'identical' twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4(+) lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

A highly annotated whole-genome sequence of a Korean individual.

Recent advances in sequencing technologies have initiated an era of personal genome sequences. To date, human genome sequences have been reported for individuals with ancestry in three distinct geographical regions: a Yoruba African, two individuals of northwest European origin, and a person from China. Here we provide a highly annotated, whole-genome sequence for a Korean individual, known as AK1. The genome of AK1 was determined by an exacting, combined approach that included whole-genome shotgun sequencing (27.8x coverage), targeted bacterial artificial chromosome sequencing, and high-resolution comparative genomic hybridization using custom microarrays featuring more than 24 million probes. Alignment to the NCBI reference, a composite of several ethnic clades, disclosed nearly 3.45 million single nucleotide polymorphisms (SNPs), including 10,162 non-synonymous SNPs, and 170,202 deletion or insertion polymorphisms (indels). SNP and indel densities were strongly correlated genome-wide. Applying very conservative criteria yielded highly reliable copy number variants for clinical considerations. Potential medical phenotypes were annotated for non-synonymous SNPs, coding domain indels, and structural variants. The integration of several human whole-genome sequences derived from several ethnic groups will assist in understanding genetic ancestry, migration patterns and population bottlenecks.

A genome-wide linkage scan of insulin level derived traits: the Amish Family Diabetes Study.

OBJECTIVE: Serum insulin levels are altered in insulin resistance and insulin deficiency, states that are associated with the development of type 2 diabetes. The goal of our study was to identify chromosomal regions that are likely to harbor genetic determinants of these traits. RESEARCH DESIGN AND METHODS: We conducted a series of genetic analyses, including genome-wide and fine-mapping linkage studies, based on insulin levels measured during an oral glucose tolerance test (OGTT) in 552 nondiabetic participants in the Amish Family Diabetes Study. Indices of insulin secretion included the insulinogenic index and insulin at 30 min postglucose load (insulin 30), while indices of insulin resistance included homeostasis model assessment of insulin resistance (HOMA-IR) and fasting insulin. Insulin area under the curve, a measure of both insulin secretion and insulin resistance, was also examined. RESULTS: All traits were modestly heritable, with heritability estimates ranging from 0.1 to 0.4 (all P < 0.05). There was significant genetic correlation between fasting insulin and HOMA-IR (rho(G) > 0.86, P < 0.05), as well as insulin 30 and insulinogenic index (rho(G) = 0.81, P < 0.0001), suggesting that common genes influence variation in these pairs of traits. Suggestive linkage signals in the genome scan were to insulin 30 on chromosome 15q23 (logarithm of odds [LOD] 2.53, P = 0.00032) and to insulinogenic index on chromosome 2p13 (LOD 2.51, P = 0.00034). Fine-mapping study further refined our signal for insulin 30 on chromosome 15 (LOD 2.38 at 68 cM). CONCLUSIONS: These results suggest that there may be different genes influencing variation in OGTT measures of insulin secretion and insulin resistance.

Genome-wide and fine-mapping linkage studies of type 2 diabetes and glucose traits in the Old Order Amish: evidence for a new diabetes locus on chromosome 14q11 and confirmation of a locus on chromosome 1q21-q24.

We conducted a genome scan using a 10-cM map to search for genes linked to type 2 diabetes in 691 individuals from a founder population, the Old Order Amish. We then saturated two regions on chromosomes 1 and 14 showing promising linkage signals with additional markers to produce a approximately 2-cM map for fine mapping. Analyses of both discrete traits (type 2 diabetes and the composite trait of type 2 diabetes and/or impaired glucose homeostasis [IGH]), and quantitative traits (glucose levels during a 75-g oral glucose challenge, designated glucose 0-180 and HbA(1c)) were performed. We obtained significant evidence for linkage to type 2 diabetes in a novel region on chromosome 14q11 (logarithm of odds [LOD] for diabetes = 3.48, P = 0.00005). Furthermore, we observed evidence for the existence of a diabetes-related locus on chromosome 1q21-q24 (LOD for type 2 diabetes/IGH = 2.35, P = 0.0008), a region shown to be linked to diabetes in several other studies. Suggestive evidence for linkage to glucose traits was observed on three other regions: 14q11-q13 (telomeric to that above with LOD = 1.82-1.85 for glucose 150 and 180), 1p31 (LOD = 1.28-2.30 for type 2 diabetes and glucose 120-180), and 18p (LOD = 3.07, P = 0.000085 for HbA(1c) and LOD = 1.50 for glucose 0). In conclusion, our findings provide evidence that type 2 diabetes susceptibility genes reside on chromosomes 1, 14, and 18.