Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next‐generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single‐nucleotide polymorphisms detected using RNA‐seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.     

Next‐generation sequencing,FISH mapping and synteny‐based modeling reveal mechanisms of decreasing dysploidy in Cucumis

In the large Cucurbitaceae genus Cucumis, cucumber (C. sativus) is the only species with 2n = 2x = 14 chromosomes. The majority of the remaining species, including melon (C. melo) and the sister species of cucumber, C. hystrix, have 2n = 2x = 24 chromosomes, implying a reduction from n = 12 to n = 7. To understand the underlying mechanisms, we investigated chromosome synteny among cucumber, C. hystrix and melon using integrated and complementary approaches. We identified 14 inversions and a C. hystrix lineage‐specific reciprocal inversion between C. hystrix and melon. The results reveal the location and orientation of 53 C. hystrix syntenic blocks on the seven cucumber chromosomes, and allow us to infer at least 59 chromosome rearrangement events that led to the seven cucumber chromosomes, including five fusions, four translocations, and 50 inversions. The 12 inferred chromosomes (AK1–AK12) of an ancestor similar to melon and C. hystrix had strikingly different evolutionary fates, with cucumber chromosome C1 apparently resulting from insertion of chromosome AK12 into the centromeric region of translocated AK2/AK8, cucumber chromosome C3 originating from a Robertsonian‐like translocation between AK4 and AK6, and cucumber chromosome C5 originating from fusion of AK9 and AK10. Chromosomes C2, C4 and C6 were the result of complex reshuffling of syntenic blocks from three (AK3, AK5 and AK11), three (AK5, AK7 and AK8) and five (AK2, AK3, AK5, AK8 and AK11) ancestral chromosomes, respectively, through 33 fusion, translocation and inversion events. Previous results (Huang, S., Li, R., Zhang, Z. et al., 2009 , Nat. Genet. 41, 1275–1281; Li, D., Cuevas, H.E., Yang, L., Li, Y., Garcia‐Mas, J., Zalapa, J., Staub, J.E., Luan, F., Reddy, U., He, X., Gong, Z., Weng, Y. 2011a, BMC Genomics, 12, 396) showing that cucumber C7 stayed largely intact during the entire evolution of Cucumis are supported. Results from this study allow a fine‐scale understanding of the mechanisms of dysploid chromosome reduction that has not been achieved previously.     

Precise marker excision system using an animal‐derived piggyBac transposon in plants

Accurate and effective positive marker excision is indispensable for the introduction of desired mutations into the plant genome via gene targeting (GT) using a positive/negative counter selection system. In mammals, the moth‐derived piggyBac transposon system has been exploited successfully to eliminate a selectable marker from a GT locus without leaving a footprint. Here, we present evidence that the piggyBac transposon also functions in plant cells. To demonstrate the use of the piggyBac transposon for effective marker excision in plants, we designed a transposition assay system that allows the piggyBac transposition to be visualized as emerald luciferase (Eluc) luminescence in rice cells. The Eluc signal derived from piggyBac excision was observed in hyperactive piggyBac transposase‐expressing rice calli. Polymerase chain reaction, Southern blot analyses and sequencing revealed the efficient and precise transposition of piggyBac in these calli. Furthermore, we have demonstrated the excision of a selection marker from a reporter locus in T₀ plants without concomitant re‐integration of the transposon and at a high frequency (44.0% of excision events), even in the absence of negative selection.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour‐joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within‐prairie genetic diversity (92%). Using Bayenv 2, 14 top‐ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScan F_ST outliers were in common with Bayenv 2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.     

Origins and functional diversification of salinity‐responsive Na+, K+ ATPase α1 paralogs in salmonids

The Salmoniform whole‐genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity‐responsive paralogs of the Na⁺, K⁺ ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue‐specific expression patterns. We found that these paralogs arose during a small‐scale duplication event prior to the Salmoniform, but after the teleost, whole‐genome duplication. The ‘freshwater paralog’ (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the ‘saltwater paralog’ (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na⁺, K⁺ ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity‐responsive Na⁺, K⁺ ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole‐genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity‐responsive Na⁺, K⁺ ATPase α1 paralogs in fishes.     

De novo MECP2 duplication derived from paternal germ line result in dysmorphism and developmental delay

Xq28 duplications encompassing the methyl CpG binding protein 2 (MECP2) in males exhibit a distinct phenotype, including developmental delay, facial dysmorphism, muscular hypotonia, intellectual disability, poor or absent speech, recurrent infections and early death. The vast majority of affected males inherit the MECP2 duplication from their usually asymptomatic carrier mothers. Only a few cases with Xq28 duplication originating from de novo unbalanced X/Y translocation have been reported and the paternal origin of the aberration has only been validated in three males in the related literature. Here we present a karyotypically normal male with features characteristic of the MECP2 duplication syndrome. The genome-wide SNP genotyping shows a de novo 2.26-Mb duplication from Xq28 to the terminus. The genotypes of the SNPs within the duplicated region indicated a paternal origin. Furthermore, the results of fluorescence in situ hybridization (FISH) indicated a novel Xq:Yp translocation, characterized as der(Y)t(Y;X)(p11.32;q28), which suggests an aberrant that occurred during spermatogenesis. The phenotype is compared to the previously reported cases with Xq28 duplication originated from an unbalanced X/Y translocation, and there was no specific part of the phenotype that could be contributed to the origin of parental imbalances. This report further highlights the capacity of high-molecular cytogenetic methods, such as SNP array and FISH, in the identification of submicroscopic rearrangement, structural configuration and parental origin of aberrant while in the evaluation of children with idiopathic developmental delay and intellectual disability.     

The analysis of microsatellites and compound microsatellites in 56 complete genomes of Herpesvirales

Simple sequence repeats (SSRs), or microsatellites, are special DNA/RNA sequences with repeated unit of 1–6 bp. The genomes of Herpesvirales have many repeating structures, which is an excellent system to study the evolution and roles of microsatellites and compound microsatellites in viruses. Therefore, 56 genomes of Herpesvirales were selected and the occurrence, composition and complexity of different repeats were investigated in the genomes. A total of 63,939 microsatellites and 5825 compound microsatellites were extracted from 56 genomes. It found that GC content has a significant strong correlation with both the counts of microsatellites (CM) and the counts of compound microsatellites (CCM). However, genome size has a moderate correlation only with CM and almost no correlation with CCM. The compound microsatellites occurring in genic regions are obviously more than that in intergenic regions. In general, the number of compound microsatellite decreases with the increase of complexity (C) (the count of individual microsatellites being part of a compound microsatellite) and the complexity hardly exceeds C = 4. The vast majority of compound microsatellites exist in intergenic regions, when C ≥ 10. The distributions of SSRs tend to be organism-specific rather than host-specific in herpesvirus genomes. The diversity of microsatellites and compound microsatellites may be helpful for a better understanding of the viral genetic diversity, genotyping, and evolutionary biology in herpesviruses genomes.     

Complete chloroplast genome sequences of Praxelis (Eupatorium catarium Veldkamp), an important invasive species

Praxelis (Eupatorium catarium Veldkamp) is a new hazardous invasive plant species that has caused serious economic losses and environmental damage in the Northern hemisphere tropical and subtropical regions. Although previous studies focused on detecting the biological characteristics of this plant to prevent its expansion, little effort has been made to understand the impact of Praxelis on the ecosystem in an evolutionary process. The genetic information of Praxelis is required for further phylogenetic identification and evolutionary studies. Here, we report the complete Praxelis chloroplast (cp) genome sequence. The Praxelis chloroplast genome is 151,410 bp in length including a small single-copy region (18,547 bp) and a large single-copy region (85,311 bp) separated by a pair of inverted repeats (IRs; 23,776 bp). The genome contains 85 unique and 18 duplicated genes in the IR region. The gene content and organization are similar to other Asteraceae tribe cp genomes. We also analyzed the whole cp genome sequence, repeat structure, codon usage, contraction of the IR and gene structure/organization features between native and invasive Asteraceae plants, in order to understand the evolution of organelle genomes between native and invasive Asteraceae. Comparative analysis identified the 14 markers containing greater than 2% parsimony-informative characters, indicating that they are potential informative markers for barcoding and phylogenetic analysis. Moreover, a sister relationship between Praxelis and seven other species in Asteraceae was found based on phylogenetic analysis of 28 protein-coding sequences. Complete cp genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.