Genome size evolution is known to be related with transposable elements, yet such relation in incipient species remains poorly understood. For decades, the willistoni subgroup of Drosophila has been a model for evolutionary studies because of the different evolutionary stages and degrees of reproductive isolation its species present. Our main question here was how speciation influences genome size evolution and the fraction of repetitive elements, with a focus on transposable elements. We quantitatively compared the mobilome of four species and two subspecies belonging to this subgroup with their genome size, and performed comparative phylogenetic analyses. Our results showed that genome size and the fraction of repetitive elements evolved according to the evolutionary history of these species, but the content of transposable elements showed some discrepancies. Signals of recent transposition events were detected for different superfamilies. Their low genomic GC content suggests that in these species transposable element mobilization might be facilitated by relaxed natural selection. Additionally, a possible role of the superfamily DNA/TcMar-Tigger in the expansion of these genomes was also detected. We hypothesize that the undergoing process of speciation could be promoting the observed increase in the fraction of repetitive elements and, consequently, genome size.

Natural wild populations of C. rupestris and C. salonitana were studied to determine possible relationships between the volatile oil (VO) composition and ploidy level. The chemical composition of the volatile oil was investigated using the GC/MS technique. The predominant components of the VO of diploid and tetraploid C. salonitana were hexadecanoic acid and α-linoleic acids, while in C. rupestris they were germacrene D and β-caryophyllene in one population and heptacosane and germacrene D, in another. The nuclear DNA amounts (2 C DNA), determined by flow cytometry, were 3.54 pg for C. rupestris, 3.39 pg for the diploid and 6.79 pg for the tetraploid population of C. salonitana. Evidence that the degree of ploidy solely influences the chemical composition of the essential oil of C. salonitana was not found. The results presented are the first data to be reported on the DNA content of the studied Centaurea populations from Croatia, as well as on the chemical composition of C. salonitana volatile oil.

Polyploidy is a widespread phenomenon across angiosperms, and one of the main drivers of diversification. Whilst it frequently involves hybridisation, autopolyploidy is also an important feature of plant evolution. Minority cytotypes are frequently overlooked due to their lower frequency in populations, but the development of techniques such as flow cytometry, which enable the rapid screening of cytotype diversity across large numbers of individuals, is now providing a more comprehensive understanding of cytotype diversity within species. Senecio doronicum is a relatively common daisy found throughout European mountain grasslands from subalpine to almost nival elevations. We have carried out a population-level cytotype screening of 500 individuals from Tête Grosse (Alpes-de-Haute-Provence, France), confirming the coexistence of tetraploid (28.2%) and octoploid cytotypes (71.2%), but also uncovering a small number of hexaploid individuals (0.6%). The analysis of repetitive elements from short-read genome-skimming data combined with nuclear (ITS) and whole plastid DNA sequences support an autopolyploid origin of the polyploid S. doronicum individuals and provide molecular evidence regarding the sole contribution of tetraploids in the formation of hexaploid individuals. The evolutionary impact and resilience of the new cytotype have yet to be determined, although the coexistence of different cytotypes may indicate nascent speciation.

Genomes provide a platform for storage of chemical information that must be stable under the context in which an organism thrives. The 2\'-deoxyguanosine (G) nucleotide has the potential to provide additional chemical information beyond its Watson-Crick base-pairing capacity. Sequences with four or more runs of three G nucleotides each are potential G-quadruplex forming sequences (PQSs) that can adopt G-quadruplex folds. Herein, we analyzed sequenced genomes from the NCBI database to determine the PQS densities of the genome sequences. First, we found organisms with large genomes, including humans, alligators, and maize, have similar densities of PQSs (~300 PQSs/Mbp), and the genomes are significantly enriched in PQSs with more than four G tracks. Analysis of microorganism genomes found a greater diversity of PQS densities. In general, PQS densities positively tracked with the GC% of the genome. Exceptions to this observation were the genomes from thermophiles that had many more PQSs than expected by random chance. Analysis of the location of these PQSs in annotated genomes from the order Thermales showed these G-rich sequences to be randomly distributed; in contrast, in the order Deinococcales the PQSs were enriched and biased around transcription start sites of genes. Four representative PQSs, two each from the Thermales and Deinococcales, were studied by biophysical methods to establish the ability of them to fold to G-quadruplexes. The experiments found the two PQSs in the Thermales did not adopt G-quadruplex folds, while the two most common in the Deinococcales adopted stable parallel-stranded G-quadruplexes. The findings lead to a hypothesis that thermophilic organisms are enriched with PQSs as an unavoidable consequence to stabilize thermally their genomes to live at high temperature; in contrast, the genomes from stress-resistant bacteria found in the Deinococcales may utilize PQSs for gene regulatory purposes.

Parallel phylogenies between aphid and its obligate symbiont Buchnera are hot topics which always focused on aphid lower taxonomic levels. Symbionts in the subfamily Lachninae are special. Buchnera in many lachnine species has undergone functional and genome size reduction that was replaced by other co-obligate symbionts. In this study, we constructed the phylogenetic relationships of Lachninae with a combined dataset of five genes sequenced from Buchnera to estimate the effects of a dual symbiotic system in the aphid-Buchnera cospeciation association. The phylogeny of Buchnera in Lachninae was well-resolved in the combined dataset. Each of the genera formed strongly supported monophyletic groups, with the exception of the genus Cinara. The phylogeny based on sequences from Buchnera was divided into five tribes according to the clades of the Lachninae hosts tree, with the phylogenies of Buchnera and Lachninae being generally congruent. These results first provided evidence of parallel evolution at the aphid subfamily level comprehensively and supported the view that topological congruence between the phylogenies of Buchnera and Lachninae would not be interfered with the other co-obligate symbionts, such as Sarretia, in aphid-entosymbiont association. These results also provided new insight in understanding host-plant coevolution in lachnine lineages.

In this study we showed that constitutive heterochromatin, GC-rich DNA and rDNA are implicated in chromosomal rearrangements during the basic chromosome number changing (dysploidy) in Reichardia genus. This small Mediterranean genus comprises 8-10 species and presents three basic chromosome numbers (x = 9, 8 and 7). To assess genome evolution and differentiation processes, studies were conducted in a dysploid series of six species: R. dichotoma, R. macrophylla and R. albanica (2n = 18), R. tingitana and R. gaditana (2n = 16), and R. picroides (2n = 14). The molecular phylogeny reconstruction comprised three additional species (R. crystallina and R. ligulata, 2n = 16 and R. intermedia, 2n = 14). Our results indicate that the way of dysploidy is descending. During this process, a positive correlation was observed between chromosome number and genome size, rDNA loci number and pollen size, although only the correlation between chromosome number and genome size is still recovered significant once considering the phylogenetic effect. Fluorescent in situ hybridisation also evidenced changes in number, position and organisation of two rDNA families (35S and 5S), including the reduction of loci number and, consequently, reduction in the number of secondary constrictions and nuclear organising regions from three to one per diploid genome. The potential mechanisms of chromosomal and genome evolution, strongly implicating heterochromatin, are proposed and discussed, with particular consideration for Reichardia genus.

Transposable elements (TEs) are ubiquitous in eukaryotic genomes and their mobility impacts genome structure and function in myriad ways. Because of their abundance, activity, and repetitive nature, the characterization and analysis of TEs remain challenging, particularly from short-read sequencing projects. To overcome this difficulty, we have developed a method that estimates TE copy number from short-read sequences. To test the accuracy of our method, we first performed an in silico analysis of the reference Sorghum bicolor genome, using both reference-based and de novo approaches. The resulting TE copy number estimates were strikingly similar to the annotated numbers. We then tested our method on real short-read data by estimating TE copy numbers in several accessions of S. bicolor and its close relative S. propinquum. Both methods effectively identify and rank similar TE families from highest to lowest abundance. We found that de novo characterization was effective at capturing qualitative variation, but underestimated the abundance of some TE families, specifically families of more ancient origin. Also, interspecific reference-based mapping of S. propinquum reads to the S. bicolor database failed to fully describe TE content in S. propinquum, indicative of recent TE activity leading to changes in the respective repetitive landscapes over very short evolutionary timescales. We conclude that reference-based analyses are best suited for within-species comparisons, while de novo approaches are more reliable for evolutionarily distant comparisons.

BACKGROUND: Recent findings indicated that a correlation between genomic % AT and genome size within strains of microbial species was predominantly associated with the uptake of foreign DNA. One species however, Chlamydia trachomatis, defied any explanation. In the present study 79 fully sequenced C. trachomatis genomes, representing ocular- (nine strains), urogenital- (36 strains) and lymphogranuloma venereum strains (LGV, 22 strains), in three pathogroups, in addition to 12 laboratory isolates, were scrutinized with the intent of elucidating the positive correlation between genomic AT content and genome size.
RESULTS: The average size difference between the strains of each pathogroup was largely explained by the incorporation of genetic fragments. These fragments were slightly more AT rich than their corresponding host genomes, but not enough to justify the difference in AT content between the strains of the smaller genomes lacking the fragments. In addition, a genetic region predominantly found in the ocular strains, which had the largest genomes, was on average more GC rich than the host genomes of the urogenital strains (58.64% AT vs. 58.69% AT), which had the second largest genomes, implying that the foreign genetic regions cannot alone explain the association between genome size and AT content in C. trachomatis. 23,492 SNPs were identified for all 79 genomes, and although the SNPs were on average slightly GC rich (~47% AT), a significant association was found between genome-wide SNP AT content, for each pathogroup, and genome size (p < 0.001, R (2) = 0.86) in the C. trachomatis strains.
CONCLUSIONS: The correlation between genome size and AT content, with respect to the C. trachomatis pathogroups, was explained by the incorporation of genetic fragments unique to the ocular and/or urogenital strains into the LGV- and urogential strains in addition to the genome-wide SNP AT content differences between the three pathogroups.

Flow cytometry (FCM) is a commonly used method for estimating genome size in many organisms. The use of FCM in plants is influenced by endogenous fluorescence inhibitors and may cause an inaccurate estimation of genome size; thus, falsifying the relationship between genome size and phenotypic traits/ecological performance. Quantitative optimization of FCM methodology minimizes such errors, yet there are few studies detailing this methodology. We selected the genus Primulina, one of the most representative and diverse genera of the Old World Gesneriaceae, to evaluate the methodology effect on determining genome size. Our results showed that buffer choice significantly affected genome size estimation in six out of the eight species examined and altered the 2C-value (DNA content) by as much as 21.4%. The staining duration and propidium iodide (PI) concentration slightly affected the 2C-value. Our experiments showed better histogram quality when the samples were stained for 40 min at a PI concentration of 100 μg ml(-1). The quality of the estimates was not improved by 1-day incubation in the dark at 4°C or by centrifugation. Thus, our study determined an optimum protocol for genome size measurement in Primulina: LB01 buffer supplemented with 100 μg ml(-1) PI and stained for 40 min. This protocol also demonstrated a high universality in other Gesneriaceae genera. We report the genome size of nine Gesneriaceae species for the first time. The results showed substantial genome size variation both within and among the species, with the 2C-value ranging between 1.62 and 2.71 pg. Our study highlights the necessity of optimizing the FCM methodology prior to obtaining reliable genome size estimates in a given taxon.

The question on the patterns and limits of reduction of plastid genomes in nonphotosynthetic plants and the reasons of their conservation is one of the intriguing topics in plant genome evolution. Here, we report sequencing and analysis of plastid genome in nonphotosynthetic orchids Epipogium aphyllum and Epipogium roseum, which, with sizes of 31 and 19 kbp, respectively, represent the smallest plastid genomes characterized by now. Besides drastic reduction, which is expected, we found several unusual features of these \"minimal\" plastomes: Multiple rearrangements, highly biased nucleotide composition, and unprecedentedly high substitution rate. Only 27 and 29 genes remained intact in the plastomes of E. aphyllum and E. roseum-those encoding ribosomal components, transfer RNAs, and three additional housekeeping genes (infA, clpP, and accD). We found no signs of relaxed selection acting on these genes. We hypothesize that the main reason for retention of plastid genomes in Epipogium is the necessity to translate messenger RNAs (mRNAs) of accD and/or clpP proteins which are essential for cell metabolism. However, these genes are absent in plastomes of several plant species; their absence is compensated by the presence of a functional copy arisen by gene transfer from plastid to the nuclear genome. This suggests that there is no single set of plastid-encoded essential genes, but rather different sets for different species and that the retention of a gene in the plastome depends on the interaction between the nucleus and plastids.