A Kinesin-14 Motor Activates Neocentromeres to Promote Meiotic Drive in Maize.

Maize abnormal chromosome 10 (Ab10) encodes a classic example of true meiotic drive that converts heterochromatic regions called knobs into motile neocentromeres that are preferentially transmitted to egg cells. Here, we identify a cluster of eight genes on Ab10, called the Kinesin driver (Kindr) complex, that are required for both neocentromere motility and preferential transmission. Two meiotic drive mutants that lack neocentromere activity proved to be kindr epimutants with increased DNA methylation across the entire gene cluster. RNAi of Kindr induced a third epimutant and corresponding loss of meiotic drive. Kinesin gliding assays and immunolocalization revealed that KINDR is a functional minus-end-directed kinesin that localizes specifically to knobs containing 180 bp repeats. Sequence comparisons suggest that Kindr diverged from a Kinesin-14A ancestor ∼12 mya and has driven the accumulation of > 500 Mb of knob repeats and affected the segregation of thousands of genes linked to knobs on all 10 chromosomes.

Early transcriptome changes induced by the Geminivirus C4 oncoprotein: setting the stage for oncogenesis.

BACKGROUND: The Beet curly top virus C4 oncoprotein is a pathogenic determinant capable of inducing extensive developmental abnormalities. No studies to date have investigated how the transcriptional profiles differ between plants expressing or not expressing the C4 oncoprotein. RESULTS: We investigated early transcriptional changes in Arabidopsis associated with expression of the Beet curly top virus C4 protein that represent initial events in pathogenesis via a comparative transcriptional analysis of mRNAs and small RNAs. We identified 48 and 94 differentially expressed genes at 6- and 12-h post-induction versus control plants. These early time points were selected to focus on direct regulatory effects of C4 expression. Since previous evidence suggested that the C4 protein regulated the brassinosteroid (BR)-signaling pathway, differentially expressed genes could be divided into two groups: those responsive to alterations in the BR-signaling pathway and those uniquely responsive to C4. Early transcriptional changes that disrupted hormone homeostasis, 18 and 19 differentially expressed genes at both 6- and 12-hpi, respectively, were responsive to C4-induced regulation of the BR-signaling pathway. Other C4-induced differentially expressed genes appeared independent of the BR-signaling pathway at 12-hpi, including changes that could alter cell development (4 genes), cell wall homeostasis (5 genes), redox homeostasis (11 genes) and lipid transport (4 genes). Minimal effects were observed on expression of small RNAs. CONCLUSION: This work identifies initial events in genetic regulation induced by a geminivirus C4 oncoprotein. We provide evidence suggesting the C4 protein regulates multiple regulatory pathways and provides valuable insights into the role of the C4 protein in regulating initial events in pathogenesis.

Histone H3 Variant Regulates RNA Polymerase II Transcription Termination and Dual Strand Transcription of siRNA Loci in Trypanosoma brucei.

Base J, ß-D-glucosyl-hydroxymethyluracil, is a chromatin modification of thymine in the nuclear DNA of flagellated protozoa of the order Kinetoplastida. In Trypanosoma brucei, J is enriched, along with histone H3 variant (H3.V), at sites involved in RNA Polymerase (RNAP) II termination and telomeric sites involved in regulating variant surface glycoprotein gene (VSG) transcription by RNAP I. Reduction of J in T. brucei indicated a role of J in the regulation of RNAP II termination, where the loss of J at specific sites within polycistronic gene clusters led to read-through transcription and increased expression of downstream genes. We now demonstrate that the loss of H3.V leads to similar defects in RNAP II termination within gene clusters and increased expression of downstream genes. Gene derepression is intensified upon the subsequent loss of J in the H3.V knockout. mRNA-seq indicates gene derepression includes VSG genes within the silent RNAP I transcribed telomeric gene clusters, suggesting an important role for H3.V in telomeric gene repression and antigenic variation. Furthermore, the loss of H3.V at regions of overlapping transcription at the end of convergent gene clusters leads to increased nascent RNA and siRNA production. Our results suggest base J and H3.V can act independently as well as synergistically to regulate transcription termination and expression of coding and non-coding RNAs in T. brucei, depending on chromatin context (and transcribing polymerase). As such these studies provide the first direct evidence for histone H3.V negatively influencing transcription elongation to promote termination.

Abundance and Multilocus Sequence Analysis of Vibrio Bacteria Associated with Diseased Elkhorn Coral (Acropora palmata) of the Florida Keys.

The critically endangered elkhorn coral (Acropora palmata) is affected by white pox disease (WPX) throughout the Florida Reef Tract and wider Caribbean. The bacterium Serratia marcescens was previously identified as one etiologic agent of WPX but is no longer consistently detected in contemporary outbreaks. It is now believed that multiple etiologic agents cause WPX; however, to date, no other potential pathogens have been thoroughly investigated. This study examined the association of Vibrio bacteria with WPX occurrence from August 2012 to 2014 at Looe Key Reef in the Florida Keys, USA. The concentration of cultivable Vibrio was consistently greater in WPX samples than in healthy samples. The abundance of Vibrio bacteria relative to total bacteria was four times higher in samples from WPX lesions than in adjacent apparently healthy regions of diseased corals based on quantitative PCR (qPCR). Multilocus sequence analysis (MLSA) was used to assess the diversity of 69 Vibrio isolates collected from diseased and apparently healthy A. palmata colonies and the surrounding seawater. Vibrio species with known pathogenicity to corals were detected in both apparently healthy and diseased samples. While the causative agent(s) of contemporary WPX outbreaks remains elusive, our results suggest that Vibrio spp. may be part of a nonspecific heterotrophic bacterial bloom rather than acting as primary pathogens. This study highlights the need for highly resolved temporal sampling in situ to further elucidate the role of Vibrio during WPX onset and progression.IMPORTANCE Coral diseases are increasing worldwide and are now considered a major contributor to coral reef decline. In particular, the Caribbean has been noted as a coral disease hot spot, owing to the dramatic loss of framework-building acroporid corals due to tissue loss diseases. The pathogenesis of contemporary white pox disease (WPX) outbreaks in Acropora palmata remains poorly understood. This study investigates the association of Vibrio bacteria with WPX.

Whole-genome sequence and annotation of Penstemon davidsonii.

Penstemon is the most speciose flowering plant genus endemic to North America. Penstemon species' diverse morphology and adaptation to various environments have made them a valuable model system for studying evolution. Here, we report the first full reference genome assembly and annotation for Penstemon davidsonii. Using PacBio long-read sequencing and Hi-C scaffolding technology, we constructed a de novo reference genome of 437,568,744 bases, with a contig N50 of 40 Mb and L50 of 5. The annotation includes 18,199 gene models, and both the genome and transcriptome assembly contain over 95% complete eudicot BUSCOs. This genome assembly will serve as a valuable reference for studying the evolutionary history and genetic diversity of the Penstemon genus.

A detailed gene expression study of the Miscanthus genus reveals changes in the transcriptome associated with the rejuvenation of spring rhizomes.

BACKGROUND: The Miscanthus genus of perennial C4 grasses contains promising biofuel crops for temperate climates. However, few genomic resources exist for Miscanthus, which limits understanding of its interesting biology and future genetic improvement. A comprehensive catalog of expressed sequences were generated from a variety of Miscanthus species and tissue types, with an emphasis on characterizing gene expression changes in spring compared to fall rhizomes. RESULTS: Illumina short read sequencing technology was used to produce transcriptome sequences from different tissues and organs during distinct developmental stages for multiple Miscanthus species, including Miscanthus sinensis, Miscanthus sacchariflorus, and their interspecific hybrid Miscanthus × giganteus. More than fifty billion base-pairs of Miscanthus transcript sequence were produced. Overall, 26,230 Sorghum gene models (i.e., ~ 96% of predicted Sorghum genes) had at least five Miscanthus reads mapped to them, suggesting that a large portion of the Miscanthus transcriptome is represented in this dataset. The Miscanthus × giganteus data was used to identify genes preferentially expressed in a single tissue, such as the spring rhizome, using Sorghum bicolor as a reference. Quantitative real-time PCR was used to verify examples of preferential expression predicted via RNA-Seq. Contiguous consensus transcript sequences were assembled for each species and annotated using InterProScan. Sequences from the assembled transcriptome were used to amplify genomic segments from a doubled haploid Miscanthus sinensis and from Miscanthus × giganteus to further disentangle the allelic and paralogous variations in genes. CONCLUSIONS: This large expressed sequence tag collection creates a valuable resource for the study of Miscanthus biology by providing detailed gene sequence information and tissue preferred expression patterns. We have successfully generated a database of transcriptome assemblies and demonstrated its use in the study of genes of interest. Analysis of gene expression profiles revealed biological pathways that exhibit altered regulation in spring compared to fall rhizomes, which are consistent with their different physiological functions. The expression profiles of the subterranean rhizome provides a better understanding of the biological activities of the underground stem structures that are essentials for perenniality and the storage or remobilization of carbon and nutrient resources.

Effects of the probiotic Lactobacillus animalis in murine Mycobacterium avium subspecies paratuberculosis infection.

BACKGROUND: MAP is a suspected zoonotic pathogen and the causative agent of Johne's Disease in cattle and other ruminant animals. With over $1 billion dollars in loss to the dairy industry due to Johne's Disease, efforts to eliminate or reduce MAP from cattle are of importance. The purpose of this study was to determine if daily intake of probiotics could eliminate or reduce Johne's Disease associated symptoms and pathogenesis by MAP. Post infection, animals are often asymptomatic carriers with limited shedding of the pathogen, proving early detection to be difficult. Disease and symptoms often appear 3-4 years after infection with antibiotic treatment proving ineffective. Symptoms include chronic gastrointestinal inflammation leading to severe weight-loss from poor feed and water intake cause a wasting disease. These symptoms are similar to those found in individuals with Crohn's Disease (CD); MAP has been implicated by not proven to be the causative agent of CD. Probiotics administered to livestock animals, including dairy and beef cattle have demonstrated improvements in cattle performance and health. Our objectives included determining the benefits of Lactobacillus animalis (strain name: NP-51) in MAP infected BALB/c mice by evaluating systemic and gastrointestinal response by the host and gut microbiota. Male and female animals were fed 1×106 CFU/g probiotics in sterile, powdered mouse chow daily and infected with 1 × 107 CFU/ml MAP and compared to controls. Animals were evaluated for 180 days to assess acute and chronic stages of disease, with sample collection from animals every 45 days. MAP concentrations from liver and intestinal tissues were examined using real time-PCR methods and the expression of key inflammatory markers were measured during MAP infection (interferon-gamma [IFN-Υ], Interleukin-1α, IL-12, IL-10, IL-6, and Tumor necrosis factor alpha [TNF-α]). RESULTS: Our results demonstrate administration of probiotics reduces production of IFN-Υ and IL-6 while increasing TNF-α and IL-17 in chronic disease; healthful immune responses that reduce chronic inflammation associated to MAP infection. CONCLUSIONS: We observed that the immune system's response in the presence of probiotics to MAP contributes towards host health by influencing the activity of the immune system and gut microbial populations.

The species coalescent indicates possible bat and pangolin origins of the COVID-19 pandemic.

A consensus species tree is reconstructed from 11 gene trees for human, bat, and pangolin beta coronaviruses from samples taken early in the pandemic (prior to April 1, 2020). Using coalescent theory, the shallow (short branches relative to the hosts) consensus species tree provides evidence of recent gene flow events between bat and pangolin beta coronaviruses predating the zoonotic transfer to humans. The consensus species tree was also used to reconstruct the ancestral sequence of human SARS-CoV-2, which was 2 nucleotides different from the Wuhan sequence. The time to most recent common ancestor was estimated to be Dec 8, 2019 with a bat origin. Some human, bat, and pangolin coronavirus lineages found in China are phylogenetically distinct, a rare example of a class II phylogeography pattern (Avise et al. in Ann Rev Eco Syst 18:489-422, 1987). The consensus species tree is a product of evolutionary factors, providing evidence of repeated zoonotic transfers between bat and pangolin as a reservoir for future zoonotic transfers to humans.

Serum albumin alters the expression of iron-controlled genes in Pseudomonas aeruginosa.

Pseudomonas aeruginosa, which causes serious infections in immunocompromised patients, produces numerous virulence factors, including exotoxin A and the siderophore pyoverdine. As production of these virulence factors is influenced by the host environment, we examined the effect serum has on global transcription within P. aeruginosa strain PAO1 at different phases of growth in an iron-deficient medium. At early exponential phase, serum significantly enhanced expression of 138 genes, most of which are repressed by iron, including pvdS, regA and the pyoverdine synthesis genes. However, serum did not interfere with the repression of these genes by iron. Serum enhanced regA expression in a fur mutant of PAO1 but not in a pvdS mutant. The serum iron-binding protein apotransferrin, but not ferritin, enhanced regA and pvdS expression. However, in PAO1 grown in a chemically defined medium that contains no iron, serum but not apotransferrin enhanced pvdS and regA expression. While complement inactivation failed to eliminate this effect, albumin absorption reduced the effect of serum on pvdS and regA expression in the iron-deficient medium chelexed tryptic soy broth dialysate. Additionally, albumin absorption eliminated the effect of serum on pvdS and regA expression in the chemically defined medium. These results suggest that serum enhances the expression of P. aeruginosa iron-controlled genes by two mechanisms: one through apotransferrin and another one through albumin.

Diet composition and diversity does not explain fewer, smaller urban nestlings.

The reproductive success of animals breeding in cities is often lower compared to counterparts that inhabit rural, suburban, and peri-urban areas. Urban dwelling may be especially costly for offspring development and survival. Diet composition and diversity may underlie factors that lead to lower fitness, particularly if prey abundance and quality decline in modified environments. Moreover, breeding success may change over the course of a season, an effect that may be augmented in urban areas. In this study, we tested the hypothesis that habitat and date affected nestling house wren (Troglodytes aedon) body condition and survival, and examined whether diet explained differences in nestling success. We monitored urban and rural populations of house wrens breeding in nest boxes, and tested whether clutch size, nestling survivorship, and nestling body condition varied by habitat or by date, and then characterized the diet of a subset of nestlings with DNA metabarcoding of fecal samples. Urbanization had clear impacts on house wren nestling fitness: urban broods contained fewer, smaller nestlings. Early nestling survival decreased as the breeding season progressed, and this effect was more pronounced in the urban population. However, the diets of urban and rural nestlings were similar and did not explain differences in body condition. Instead, across populations, diet changed with date, becoming less diverse, with fewer Lepidoptera and more Orthoptera. Regardless of habitat, adult house wrens provide nestlings with similar types of foods, but other factors, such as quantity or quality of prey delivered, may lead to fitness disparities between urban and rural nestlings.

Analysis of genes associated with proinflammatory and profibrotic pathways upregulated in ischemia-induced chronic kidney disease in cats.

OBJECTIVE: To use RNA sequencing (RNAseq) to characterize renal transcriptional activities of genes associated with proinflammatory and profibrotic pathways in ischemia-induced chronic kidney disease (CKD) in cats. SAMPLES: Banked renal tissues from 6 cats with experimentally induced CKD (renal ischemia [RI] group) and 9 healthy cats (control group). PROCEDURES: Transcriptome analysis with RNAseq, followed by gene ontology and cluster analyses, were performed on banked tissue samples of the right kidneys (control kidneys) from cats in the control group and of both kidneys from cats in the RI group, in which unilateral (right) RI had been induced 6 months before the cats were euthanized and the ischemic kidneys (IKs) and contralateral nonischemic kidneys (CNIKs) were harvested. Results for the IKs, CNIKs, and control kidneys were compared to identify potential differentially expressed genes and overrepresented proinflammatory and profibrotic pathways. RESULTS: Genes from the gene ontology pathways of collagen binding (eg, transforming growth factor-ß1), metalloendopeptidase activity (eg, metalloproteinase [MMP]-7, MMP-9, MMP-11, MMP-13, MMP-16, MMP-23B, and MMP-28), chemokine activity, and T-cell migration were overrepresented as upregulated in tissue samples of the IKs versus control kidneys. Genes associated with the extracellular matrix (eg, TIMP-1, fibulin-1, secreted phosphoprotein-1, matrix Gla protein, and connective tissue growth factor) were upregulated in tissue samples from both the IKs and CNIKs, compared with tissues from the control kidneys. CONCLUSIONS AND CLINICAL RELEVANCE: Unilateral ischemic injury differentially altered gene expression in both kidneys, compared with control kidneys. Fibulin-1, secreted phosphoprotein-1, and matrix Gla protein may be candidate biomarkers of active kidney injury in cats.

Characterizing changes in soil microbiome abundance and diversity due to different cover crop techniques.

Soil-based microorganisms assume a direct and crucial role in the promotion of soil health, quality and fertility, all factors known to contribute heavily to the quality and yield of agricultural products. Cover cropping, used in both traditional and organic farming, is a particularly efficient and environmentally favorable tool for manipulating microbiome composition in agricultural soils and has had clear benefits for soil quality and crop output. Several long-term investigations have evaluated the influence of multi-mix (multiple species) cover crop treatments on soil health and microbial diversity. The present study investigated the short-term effects of a seven species multi-mix cover crop treatment on soil nutrient content and microbial diversity, compared to a single-mix cover crop treatment and control. Analysis of 16S sequencing data of isolated soil DNA revealed that the single-mix cover crop treatment decreased overall microbial abundance and diversity, whereas the control and multi-mix treatments altered the overall microbial composition in similar fluctuating trends. Furthermore, we observed significant changes in specific bacteria belonging to the phyla Acidobacteria, Actinobacteria, Planctomycetes, Proteobacteria and Verrucombicrobia for all treatments, but only the single-mix significantly decreased in abundance of the selected bacteria over time. Our findings indicate that the control and multi-mix treatments are better at maintaining overall microbial composition and diversity compared to the single-mix. Further study is required to elucidate the specific difference between the treatment effect of the multi-mix treatment and the control, given that their microbial composition changes over time were similar but they diverge into two populations of unique bacterial types by the end of this short-term study.

Gapless assembly of maize chromosomes using long-read technologies.

Creating gapless telomere-to-telomere assemblies of complex genomes is one of the ultimate challenges in genomics. We use two independent assemblies and an optical map-based merging pipeline to produce a maize genome (B73-Ab10) composed of 63 contigs and a contig N50 of 162 Mb. This genome includes gapless assemblies of chromosome 3 (236 Mb) and chromosome 9 (162 Mb), and 53 Mb of the Ab10 meiotic drive haplotype. The data also reveal the internal structure of seven centromeres and five heterochromatic knobs, showing that the major tandem repeat arrays (CentC, knob180, and TR-1) are discontinuous and frequently interspersed with retroelements.

Gene-rich islands for fiber development in the cotton genome.

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext).

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.

Double feature selection and cluster analyses in mining of microarray data from cotton.

BACKGROUND: Cotton fiber is a single-celled seed trichome of major biological and economic importance. In recent years, genomic approaches such as microarray-based expression profiling were used to study fiber growth and development to understand the developmental mechanisms of fiber at the molecular level. The vast volume of microarray expression data generated requires a sophisticated means of data mining in order to extract novel information that addresses fundamental questions of biological interest. One of the ways to approach microarray data mining is to increase the number of dimensions/levels to the analysis, such as comparing independent studies from different genotypes. However, adding dimensions also creates a challenge in finding novel ways for analyzing multi-dimensional microarray data. RESULTS: Mining of independent microarray studies from Pima and Upland (TM1) cotton using double feature selection and cluster analyses identified species-specific and stage-specific gene transcripts that argue in favor of discrete genetic mechanisms that govern developmental programming of cotton fiber morphogenesis in these two cultivated species. Double feature selection analysis identified the highest number of differentially expressed genes that distinguish the fiber transcriptomes of developing Pima and TM1 fibers. These results were based on the finding that differences in fibers harvested between 17 and 24 day post-anthesis (dpa) represent the greatest expressional distance between the two species. This powerful selection method identified a subset of genes expressed during primary (PCW) and secondary (SCW) cell wall biogenesis in Pima fibers that exhibits an expression pattern that is generally reversed in TM1 at the same developmental stage. Cluster and functional analyses revealed that this subset of genes are primarily regulated during the transition stage that overlaps the termination of PCW and onset of SCW biogenesis, suggesting that these particular genes play a major role in the genetic mechanism that underlies the phenotypic differences in fiber traits between Pima and TM1. CONCLUSION: The novel application of double feature selection analysis led to the discovery of species- and stage-specific genetic expression patterns, which are biologically relevant to the genetic programs that underlie the differences in the fiber phenotypes in Pima and TM1. These results promise to have profound impacts on the ongoing efforts to improve cotton fiber traits.

A carnivorous plant genetic map: pitcher/insect-capture QTL on a genetic linkage map of Sarracenia.

The study of carnivorous plants can afford insight into their unique evolutionary adaptations and their interactions with prokaryotic and eukaryotic species. For Sarracenia (pitcher plants), we identified 64 quantitative trait loci (QTL) for insect-capture traits of the pitchers, providing the genetic basis for differences between the pitfall and lobster-trap strategies of insect capture. The linkage map developed here is based upon the F2 of a cross between Sarracenia rosea and Sarracenia psittacina; we mapped 437 single nucleotide polymorphism and simple sequence repeat markers. We measured pitcher traits which differ between S. rosea and S. psittacina, mapping 64 QTL for 17 pitcher traits; there are hot-spot locations where multiple QTL map near each other. There are epistatic interactions in many cases where there are multiple loci for a trait. The QTL map uncovered the genetic basis for the differences between pitfall- and lobster-traps, and the changes that occurred during the divergence of these species. The longevity and clonability of Sarracenia plants make the F2 mapping population a resource for mapping more traits and for phenotype-to-genotype studies.

Widespread natural variation of DNA methylation within angiosperms.

BACKGROUND: DNA methylation is an important feature of plant epigenomes, involved in the formation of heterochromatin and affecting gene expression. Extensive variation of DNA methylation patterns within a species has been uncovered from studies of natural variation. However, the extent to which DNA methylation varies between flowering plant species is still unclear. To understand the variation in genomic patterning of DNA methylation across flowering plant species, we compared single base resolution DNA methylomes of 34 diverse angiosperm species. RESULTS: By analyzing whole-genome bisulfite sequencing data in a phylogenetic context, it becomes clear that there is extensive variation throughout angiosperms in gene body DNA methylation, euchromatic silencing of transposons and repeats, as well as silencing of heterochromatic transposons. The Brassicaceae have reduced CHG methylation levels and also reduced or loss of CG gene body methylation. The Poaceae are characterized by a lack or reduction of heterochromatic CHH methylation and enrichment of CHH methylation in genic regions. Furthermore, low levels of CHH methylation are observed in a number of species, especially in clonally propagated species. CONCLUSIONS: These results reveal the extent of variation in DNA methylation in angiosperms and show that DNA methylation patterns are broadly a reflection of the evolutionary and life histories of plant species.

Development of Transcriptomic Markers for Population Analysis Using Restriction Site Associated RNA Sequencing (RARseq).

We describe restriction site associated RNA sequencing (RARseq), an RNAseq-based genotype by sequencing (GBS) method. It includes the construction of RNAseq libraries from double stranded cDNA digested with selected restriction enzymes. To test this, we constructed six single- and six-dual-digested RARseq libraries from six F2 pitcher plant individuals and sequenced them on a half of a Miseq run. On average, the de novo approach of population genome analysis detected 544 and 570 RNA SNPs, whereas the reference transcriptome-based approach revealed an average of 1907 and 1876 RNA SNPs per individual, from single- and dual-digested RARseq data, respectively. The average numbers of RNA SNPs and alleles per loci are 1.89 and 2.17, respectively. Our results suggest that the RARseq protocol allows good depth of coverage per loci for detecting RNA SNPs and polymorphic loci for population genomics and mapping analyses. In non-model systems where complete genomes sequences are not always available, RARseq data can be analyzed in reference to the transcriptome. In addition to enriching for functional markers, this method may prove particularly useful in organisms where the genomes are not favorable for DNA GBS.

Genomic and small RNA sequencing of Miscanthus x giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses.

BACKGROUND: Miscanthus x giganteus (Mxg) is a perennial grass that produces superior biomass yields in temperate environments. The essentially uncharacterized triploid genome (3n = 57, x = 19) of Mxg is likely critical for the rapid growth of this vegetatively propagated interspecific hybrid. RESULTS: A survey of the complex Mxg genome was conducted using 454 pyrosequencing of genomic DNA and Illumina sequencing-by-synthesis of small RNA. We found that the coding fraction of the Mxg genome has a high level of sequence identity to that of other grasses. Highly repetitive sequences representing the great majority of the Mxg genome were predicted using non-cognate assembly for de novo repeat detection. Twelve abundant families of repeat were observed, with those related to either transposons or centromeric repeats likely to comprise over 95% of the genome. Comparisons of abundant repeat sequences to a small RNA survey of three Mxg organs (leaf, rhizome, inflorescence) revealed that the majority of observed 24-nucleotide small RNAs are derived from these repetitive sequences. We show that high-copy-number repeats match more of the small RNA, even when the amount of the repeat sequence in the genome is accounted for. CONCLUSIONS: We show that major repeats are present within the triploid Mxg genome and are actively producing small RNAs. We also confirm the hypothesized origins of Mxg, and suggest that while the repeat content of Mxg differs from sorghum, the sorghum genome is likely to be of utility in the assembly of a gene-space sequence of Mxg.