Supplementary data are available at Bioinformatics online.

Our understanding of human gut microbiota in health and disease depends on accurate and reproducible microbial data acquisition. The critical step in this process is to apply an appropriate methodology to extract microbial DNA, since biases introduced during the DNA extraction process may result in inaccurate microbial representation. In this study, we attempted to find a DNA extraction protocol which could be effectively used to analyze both the bacterial and fungal community. We evaluated the effect of five DNA extraction methods (QIAamp DNA Stool Mini Kit, PureLink TM Microbiome DNA Purification Kit, ZR Fecal DNA MiniPrep TM Kit, NucleoSpin ® DNA Stool Kit, and IHMS protocol Q) on bacterial and fungal gut microbiome recovery using (i) a defined system of germ-free mice feces spiked with bacterial or fungal strains, and (ii) non-spiked human feces. In our experimental setup, we confirmed that the examined methods significantly differed in efficiency and quality, which affected the identified stool microbiome composition. In addition, our results indicated that fungal DNA extraction might be prone to be affected by reagent/kit contamination, and thus an appropriate blank control should be included in mycobiome research. Overall, standardized IHMS protocol Q, recommended by the International Human Microbiome Consortium, performed the best when considering all the parameters analyzed, and thus could be applied not only in bacterial, but also in fungal microbiome research.

Rumex acetosa is a dioecious plant with the XY1Y2 sex chromosome system. Both Y chromosomes are heterochromatic and are thought to be degenerated. We performed low-pass 454 sequencing and similarity-based clustering of male and female genomic 454 reads to identify and characterize major groups of R. acetosa repetitive DNA. We found that Copia and Gypsy retrotransposons dominated, followed by DNA transposons and nonlong terminal repeat retrotransposons. CRM and Tat/Ogre retrotransposons dominated the Gypsy superfamily, whereas Maximus/Sireviruses were most abundant among Copia retrotransposons. Only one Gypsy subfamily had accumulated on Y1 and Y2 chromosomes, whereas many retrotransposons were ubiquitous on autosomes and the X chromosome, but absent on Y1 and Y2 chromosomes, and others were depleted from the X chromosome. One group of CRM Gypsy was specifically localized to centromeres. We also found that majority of previously described satellites (RAYSI, RAYSII, RAYSIII, and RAE180) are accumulated on the Y chromosomes where we identified Y chromosome-specific variant of RAE180. We discovered two novel satellites—RA160 satellite dominating on the X chromosome and RA690 localized mostly on the Y1 chromosome. The expression pattern obtained from Illumina RNA sequencing showed that the expression of transposable elements is similar in leaves of both sexes and that satellites are also expressed. Contrasting patterns of transposable elements (TEs) and satellite localization on sex chromosomes in R. acetosa, where not only accumulation but also depletion of repetitive DNA was observed, suggest that a plethora of evolutionary processes can shape sex chromosomes.

Retrotransposons with long terminal repeats (LTR) form a significant proportion of eukaryotic genomes, especially in plants. They have gag and pol genes and several regulatory regions necessary for transcription and reverse transcription. We searched for potential quadruplex-forming sequences (PQSs) and potential triplex-forming sequences (PTSs) in 18 377 full-length LTR retrotransposons collected from 21 plant species. We found that PQSs were often located in LTRs, both upstream and downstream of promoters from which the whole retrotransposon is transcribed. Upstream-located guanine PQSs were dominant in the minus DNA strand, whereas downstream-located guanine PQSs prevailed in the plus strand, indicating their role both at transcriptional and post-transcriptional levels. Our circular dichroism spectroscopy measurements confirmed that these PQSs readily adopted guanine quadruplex structures—some of them were paralell-stranded, while others were anti-parallel-stranded. The PQS often formed doublets at a mutual distance of up to 400 bp. PTSs were most abundant in 3′UTR (but were also present in 5′UTR). We discuss the potential role of quadruplexes and triplexes as the regulators of various processes participating in LTR retrotransposon life cycle and as potential recombination sites during post-insertional retrotransposon-based genome rearrangements.