Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation.

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

Review of evidence for bosentan therapy for treatment of Eisenmenger syndrome.

BACKGROUND AND PURPOSE: Eisenmenger syndrome (ES) is a rare condition caused by a right-to-left cyanotic shunt. To date, only heart-lung transplant has been shown to be curative. Bosentan is the only medication studied with a double-blind placebo-controlled randomized trial for management of this condition. The intent of this article is to explore the literature surrounding bosentan in ES and assess its efficacy. METHODS: A literature review was conducted with no limitation on date. Titles were scanned for applicability, and abstracts of those articles found to be pertinent were reviewed. Those articles considered relevant based on the abstract were read in entirety. CONCLUSIONS: Eisenmenger syndrome remains incurable except through heart-lung transplant. Although no specific medical treatment or algorithm exists, three pharmacological classes show promise in disease management: endothelin receptor antagonists, phosphodiesterase inhibitors, and prostacyclins. Combined therapy with these agents may improve cardiopulmonary function. Bosentan has not been proven as a monotherapy for ES and is not appropriate in all patients as side effects are commonly reported. IMPLICATIONS FOR PRACTICE: Further study is required to assess efficacy of combination therapy and utilization as a bridge to transplant or surgical correction of the underlying defect.

Ctf4-related protein recruits LHP1-PRC2 to maintain H3K27me3 levels in dividing cells in Arabidopsis thaliana.

Polycomb Repressive Complex (PRC) 2 catalyzes the H3K27me3 modification that warrants inheritance of a repressive chromatin structure during cell division, thereby assuring stable target gene repression in differentiated cells. It is still under investigation how H3K27me3 is passed on from maternal to filial strands during DNA replication; however, cell division can reinforce H3K27me3 coverage at target regions. To identify novel factors involved in the Polycomb pathway in plants, we performed a forward genetic screen for enhancers of the like heterochromatin protein 1 (lhp1) mutant, which shows relatively mild phenotypic alterations compared with other plant PRC mutants. We mapped enhancer of lhp1 (eol) 1 to a gene related to yeast Chromosome transmission fidelity 4 (Ctf4) based on phylogenetic analysis, structural similarities, physical interaction with the CMG helicase component SLD5, and an expression pattern confined to actively dividing cells. A combination of eol1 with the curly leaf (clf) allele, carrying a mutation in the catalytic core of PRC2, strongly enhanced the clf phenotype; furthermore, H3K27me3 coverage at target genes was strongly reduced in eol1 clf double mutants compared with clf single mutants. EOL1 physically interacted with CLF, its partially redundant paralog SWINGER (SWN), and LHP1. We propose that EOL1 interacts with LHP1-PRC2 complexes during replication and thereby participates in maintaining the H3K27me3 mark at target genes.

Improving and correcting the contiguity of long-read genome assemblies of three plant species using optical mapping and chromosome conformation capture data.

Long-read sequencing can overcome the weaknesses of short reads in the assembly of eukaryotic genomes; however, at present additional scaffolding is needed to achieve chromosome-level assemblies. We generated Pacific Biosciences (PacBio) long-read data of the genomes of three relatives of the model plant Arabidopsis thaliana and assembled all three genomes into only a few hundred contigs. To improve the contiguities of these assemblies, we generated BioNano Genomics optical mapping and Dovetail Genomics chromosome conformation capture data for genome scaffolding. Despite their technical differences, optical mapping and chromosome conformation capture performed similarly and doubled N50 values. After improving both integration methods, assembly contiguity reached chromosome-arm-levels. We rigorously assessed the quality of contigs and scaffolds using Illumina mate-pair libraries and genetic map information. This showed that PacBio assemblies have high sequence accuracy but can contain several misassemblies, which join unlinked regions of the genome. Most, but not all, of these misjoints were removed during the integration of the optical mapping and chromosome conformation capture data. Even though none of the centromeres were fully assembled, the scaffolds revealed large parts of some centromeric regions, even including some of the heterochromatic regions, which are not present in gold standard reference sequences.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Chromosome-level assembly of Arabidopsis thaliana Ler reveals the extent of translocation and inversion polymorphisms.

Resequencing or reference-based assemblies reveal large parts of the small-scale sequence variation. However, they typically fail to separate such local variation into colinear and rearranged variation, because they usually do not recover the complement of large-scale rearrangements, including transpositions and inversions. Besides the availability of hundreds of genomes of diverse Arabidopsis thaliana accessions, there is so far only one full-length assembled genome: the reference sequence. We have assembled 117 Mb of the A. thaliana Landsberg erecta (Ler) genome into five chromosome-equivalent sequences using a combination of short Illumina reads, long PacBio reads, and linkage information. Whole-genome comparison against the reference sequence revealed 564 transpositions and 47 inversions comprising ∼3.6 Mb, in addition to 4.1 Mb of nonreference sequence, mostly originating from duplications. Although rearranged regions are not different in local divergence from colinear regions, they are drastically depleted for meiotic recombination in heterozygotes. Using a 1.2-Mb inversion as an example, we show that such rearrangement-mediated reduction of meiotic recombination can lead to genetically isolated haplotypes in the worldwide population of A. thaliana Moreover, we found 105 single-copy genes, which were only present in the reference sequence or the Ler assembly, and 334 single-copy orthologs, which showed an additional copy in only one of the genomes. To our knowledge, this work gives first insights into the degree and type of variation, which will be revealed once complete assemblies will replace resequencing or other reference-dependent methods.

cDNA Library Enrichment of Full Length Transcripts for SMRT Long Read Sequencing.

The utility of genome assemblies does not only rely on the quality of the assembled genome sequence, but also on the quality of the gene annotations. The Pacific Biosciences Iso-Seq technology is a powerful support for accurate eukaryotic gene model annotation as it allows for direct readout of full-length cDNA sequences without the need for noisy short read-based transcript assembly. We propose the implementation of the TeloPrime Full Length cDNA Amplification kit to the Pacific Biosciences Iso-Seq technology in order to enrich for genuine full-length transcripts in the cDNA libraries. We provide evidence that TeloPrime outperforms the commonly used SMARTer PCR cDNA Synthesis Kit in identifying transcription start and end sites in Arabidopsis thaliana. Furthermore, we show that TeloPrime-based Pacific Biosciences Iso-Seq can be successfully applied to the polyploid genome of bread wheat (Triticum aestivum) not only to efficiently annotate gene models, but also to identify novel transcription sites, gene homeologs, splicing isoforms and previously unidentified gene loci.

Complementary Activities of TELOMERE REPEAT BINDING Proteins and Polycomb Group Complexes in Transcriptional Regulation of Target Genes.

In multicellular organisms, Polycomb Repressive Complex 1 (PRC1) and PRC2 repress target genes through histone modification and chromatin compaction. Arabidopsis thaliana mutants strongly compromised in the pathway cannot develop differentiated organs. LIKE HETEROCHROMATIN PROTEIN1 (LHP1) is so far the only known plant PRC1 component that directly binds to H3K27me3, the histone modification set by PRC2, and also associates genome-wide with trimethylation of lysine 27 of histone H3 (H3K27me3). Surprisingly, lhp1 mutants show relatively mild phenotypic alterations. To explain this paradox, we screened for genetic enhancers of lhp1 mutants to identify novel components repressing target genes together with, or in parallel to, LHP1. Two enhancing mutations were mapped to TELOMERE REPEAT BINDING PROTEIN1 (TRB1) and its paralog TRB3. We show that TRB1 binds to thousands of genomic sites containing telobox or related cis-elements with a significant increase of sites and strength of binding in the lhp1 background. Furthermore, in combination with lhp1, but not alone, trb1 mutants show increased transcription of LHP1 targets, such as floral meristem identity genes, which are more likely to be bound by TRB1 in the lhp1 background. By contrast, expression of a subset of LHP1-independent TRB1 target genes, many involved in primary metabolism, is decreased in the absence of TRB1 alone. Thus, TRB1 is a bivalent transcriptional modulator that maintains downregulation of Polycomb Group (PcG) target genes in lhp1 mutants, while it sustains high expression of targets that are regulated independently of PcG.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Fast isogenic mapping-by-sequencing of ethyl methanesulfonate-induced mutant bulks.

Mapping-by-sequencing (or SHOREmapping) has revitalized the powerful concept of forward genetic screens in plants. However, as in conventional genetic mapping approaches, mapping-by-sequencing requires phenotyping of mapping populations established from crosses between two diverged accessions. In addition to the segregation of the focal phenotype, this introduces natural phenotypic variation, which can interfere with the recognition of quantitative phenotypes. Here, we demonstrate how mapping-by-sequencing and candidate gene identification can be performed within the same genetic background using only mutagen-induced changes as segregating markers. Using a previously unknown suppressor of mutants of like heterochromatin protein1 (lhp1), which in its functional form is involved in chromatin-mediated gene repression, we identified three closely linked ethyl methanesulfonate-induced changes as putative candidates. In order to assess allele frequency differences between such closely linked mutations, we introduced deep candidate resequencing using the new Ion Torrent Personal Genome Machine sequencing platform to our mutant identification pipeline and thereby reduced the number of causal candidate mutations to only one. Genetic analysis of two independent additional alleles confirmed that this mutation was causal for the suppression of lhp1.