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.

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.

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.

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.

Eating behavior in the Old Order Amish: heritability analysis and a genome-wide linkage analysis.

BACKGROUND: Eating behavior and thus dietary intake affect the development of obesity-related diseases such as diabetes, hypertension, and hyperlipidemia. OBJECTIVE: We investigated the genetic underpinnings of eating behavior. DESIGN: We administered a standardized eating behavior inventory to 624 adults from 28 families participating in the Amish Family Diabetes Study. Three quantifiable components of eating behavior were measured: restraint, disinhibition, and hunger. Associations between eating behavior scores and physical characteristics were evaluated. Heritability analysis and a genome-wide multipoint linkage analysis were performed. RESULTS: Eating behavior scores were associated with obesity and obesity-related phenotypes. Heritability estimates were 0.28 +/- 0.09 for restraint, 0.40 +/- 0.10 for disinhibition, and 0.23 +/- 0.09 for hunger (P < 0.001). The linkage analysis showed 4 regions of suggestive linkage. We observed suggestive evidence for linkage of restraint scores to 2 chromosomal regions, near markers D3S1304 [LOD (log of odds) = 2.5, P = 0.0003] and D6S276 (LOD = 2.3, P = 0.0006). We previously reported that D3S1304 is linked to a locus influencing percentage body fat in this same population (LOD = 1.6), suggesting that this behavioral phenotype may be secondary to obesity. The maximum LOD scores for disinhibition were 1.6 (P = 0.003) near marker D7S657 and 1.4 (P = 0.005) near marker D16S752. The maximum LOD score for hunger was 1.4 (P = 0.005) near marker D3S1278. CONCLUSION: Significant familial effects on eating behavior and suggestive genetic linkage were found in Amish adults.