Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross‐species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next‐generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross‐genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter‐ and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.     

Inference of the phylogenetic position of the phylum Deferribacteres from gene order comparison

The phylogenetic placement of the phylum Deferribacteres was investigated on the basis of gene order comparisons of completely sequenced bacterial genomes. Two completely sequenced Deferribacteres species share five sets of gene arrangements with a group of phyla, Proteobacteria, Aquificae, Planctomycetes, Spirochaetes, Bacteroidetes, Chlorobi, Acidobacteria, Verrucomicrobia, Elusimicrobia and Nitrospirae, while the other group of phyla, Synergistetes, Firmicutes, Actinobacteria, Thermotogae, Chloroflexi and Deinococcus-Thermus, Fusobacteria, shares alternative sets of gene arrangements, suggesting that the Deferribacteres is classified in the former group of phyla. Gene transfers that are thought to have occurred in a common ancestor of the Deferribacteres, Deltaproteobacteria and Nitrospirae exclusive of virtually all other phyla were identified, which suggests that the Deferribacteres is phylogenetically proximal to the Proteobacteria and Nitrospirae.     

Inter‐varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter‐varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high‐throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter‐varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.     

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non‐symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi‐parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5–8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7–8 5S loci) by the Rubiginosa ribotype. The co‐localized 5S/18S loci occurred on either bivalent‐forming (R. canina) or univalent‐forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here‐observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.     

Comprehensive analyses of hox gene expression in Xenopus laevis embryos and adult tissues

From whole genome sequencing of an allotetraploid frog, Xenopus laevis, two homeologous sets (L and S) of four Hox clusters A through D (HoxA.L/S, HoxB.L/S, HoxC.L/S, and HoxD.L/S) and 13 paralogous groups (PGs) with 76 genes were identified, allowing us to carry out the first comprehensive analyses of hox gene expression in vertebrates. Expression of all hox genes during development and in adult tissues was analyzed by RNA‐sequencing. The expression levels of most hox genes were similar between homeologs, but in some pairs, large differences were observed and several of these were confirmed by RT‐PCR and whole mount in situ hybridization experiments. These results indicate that subfunctionalization of hox genes may have occurred since allotetraploidization. Furthermore, comprehensive analysis of hox gene expression during early development did not agree with the hypothesis of temporal collinearity especially in genes belonging to PG2 to PG10 .     

Nucleosomal arrays self‐assemble into supramolecular globular structures lacking 30‐nm fibers

The existence of a 30‐nm fiber as a basic folding unit for DNA packaging has remained a topic of active discussion. Here, we characterize the supramolecular structures formed by reversible Mg²⁺‐dependent self‐association of linear 12‐mer nucleosomal arrays using microscopy and physicochemical approaches. These reconstituted chromatin structures, which we call “oligomers”, are globular throughout all stages of cooperative assembly and range in size from ~50 nm to a maximum diameter of ~1,000 nm. The nucleosomal arrays were packaged within the oligomers as interdigitated 10‐nm fibers, rather than folded 30‐nm structures. Linker DNA was freely accessible to micrococcal nuclease, although the oligomers remained partially intact after linker DNA digestion. The organization of chromosomal fibers in human nuclei in situ was stabilized by 1 mM MgCl₂, but became disrupted in the absence of MgCl₂, conditions that also dissociated the oligomers in vitro. These results indicate that a 10‐nm array of nucleosomes has the intrinsic ability to self‐assemble into large chromatin globules stabilized by nucleosome–nucleosome interactions, and suggest that the oligomers are a good in vitro model for investigating the structure and organization of interphase chromosomes.     

De novo assembly of a wild pear (Pyrus betuleafolia) genome

China is the origin and evolutionary centre of Oriental pears. Pyrus betuleafolia is a wild species native to China and distributed in the northern region, and it is widely used as rootstock. Here, we report the de novo assembly of the genome of P. betuleafolia‐Shanxi Duli using an integrated strategy that combines PacBio sequencing, BioNano mapping and chromosome conformation capture (Hi‐C) sequencing. The genome assembly size was 532.7 Mb, with a contig N50 of 1.57 Mb. A total of 59 552 protein‐coding genes and 247.4 Mb of repetitive sequences were annotated for this genome. The expansion genes in P. betuleafolia were significantly enriched in secondary metabolism, which may account for the organism's considerable environmental adaptability. An alignment analysis of orthologous genes showed that fruit size, sugar metabolism and transport, and photosynthetic efficiency were positively selected in Oriental pear during domestication. A total of 573 nucleotide‐binding site (NBS)‐type resistance gene analogues (RGAs) were identified in the P. betuleafolia genome, 150 of which are TIR‐NBS‐LRR (TNL)‐type genes, which represented the greatest number of TNL‐type genes among the published Rosaceae genomes and explained the strong disease resistance of this wild species. The study of flavour metabolism‐related genes showed that the anthocyanidin reductase (ANR) metabolic pathway affected the astringency of pear fruit and that sorbitol transporter (SOT) transmembrane transport may be the main factor affecting the accumulation of soluble organic matter. This high‐quality P. betuleafolia genome provides a valuable resource for the utilization of wild pear in fundamental pear studies and breeding.     

Chromosomal‐level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white‐backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome‐level genome resources has hindered in‐depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome‐level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re‐sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex‐linked genes, and all chromosomes showed high synteny. The chromosome‐level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.     

Drought,post‐dispersal seed predation,and the establishment of epiphytic bromeliads (Tillandsia spp.)

Post‐dispersal, epiphyte seed predation is poorly documented. Our study on the in situ germination of two Tillandsia species in a deciduous forest in Yucatan, Mexico showed an average post‐dispersal seed predation rate by the Yucatan deer mouse (Peromyscus yucatanicus) of 90%. Post‐dispersal predation was thus more limiting than drought.     

Syncephalastrum contaminatum,a new species in the Mucorales from Australia

A new species is described in the Mucorales family Syncephalastraceae: Syncephalastrum contaminatum, isolated as an in vitro culture from a laboratory contaminant. The species has variable copies of the internal transcribed spacer (ITS) regions, requiring cloning of these regions prior to Sanger sequencing before subsequent use in phylogenetic comparisons with other fungi. The genome of the strain was sequenced using short paired-reads to yield a draft genome of 28.6 Mb. Syncephalastrum contaminatum is distinguished by diverse DNA sequences at several loci from the other species of Syncephalastrum, including only 81% sequence identity with its ITS regions to that of S. racemosum. Its merosporangium produces four or more asexual spores and the genome sequencing information suggests that the species is heterothallic. The identification of this species highlights the limited knowledge about the early lineages of fungi both in Australia and globally.     

Nuclear rDNA instability in in vitro‐generated plants is amplified after sexual reproduction with conspecific wild individuals

Using micropropagation through tissue culture has become the most used approach worldwide for mass production for the conservation of endangered species. However, the screening of somaclonal variations generated using in vitro culture is usually restricted to the first generation of micropropagated plants, when they have not yet been released in the field. Accordingly, the fate of genetically modified regenerants after sexual reproduction is usually not assessed and changes in the genetic structures of species are unknown. In this work, we assess the cytogenetic stability of two rDNA gene families in the offspring of experimental crosses between accessions generated after in vitro culture and wild individuals of Cistus heterophyllus (Cistaceae). The cytogenetic rDNA profiles (45S rDNA, 5S rDNA) of 118 accessions including wild and in vitro micropropagated individuals and bi‐directional artificial crosses between wild and in vitro‐generated plants were assessed by fluorescence in situ hybridization (FISH) and Ag‐NOR staining. Plants regenerated by micropropagation showed a lower size of the FISH signals in a 45S rDNA site, but this condition was not present in the wild accessions. Three new cytogenetic and cytological variants were present in 36% of the experimental progeny, involving the amplification of one additional 45S rDNA site and the presence of heteromorphic nucleoli. rDNA‐based genomic instability was present after sexual reproduction between wild and in vitro‐generated plants. The results of this study discourage the use of micropropagated materials for plant conservation unless comprehensive surveys of the genetic integrity and stability of the regenerants are performed after crossing between wild and micropropagated plants.     

Compartmentalisation of [FeFe]‐hydrogenase maturation in Chlamydomonas reinhardtii

Molecular hydrogen (H₂) can be produced in green microalgae by [FeFe]‐hydrogenases as a direct product of photosynthesis. The Chlamydomonas reinhardtii hydrogenase HYDA1 contains a catalytic site comprising a classic [4Fe4S] cluster linked to a unique 2Fe sub‐cluster. From in vitro studies it appears that the [4Fe4S] cluster is incorporated first by the housekeeping FeS cluster assembly machinery, followed by the 2Fe sub‐cluster, whose biosynthesis requires the specific maturases HYDEF and HYDG. To investigate the maturation process in vivo, we expressed HYDA1 from the C. reinhardtii chloroplast and nuclear genomes (with and without a chloroplast transit peptide) in a hydrogenase‐deficient mutant strain, and examined the cellular enzymatic hydrogenase activity, as well as in vivo H₂ production. The transformants expressing HYDA1 from the chloroplast genome displayed levels of H₂ production comparable to the wild type, as did the transformants expressing full‐length HYDA1 from the nuclear genome. In contrast, cells equipped with cytoplasm‐targeted HYDA1 produced inactive enzyme, which could only be activated in vitro after reconstitution of the [4Fe4S] cluster. This indicates that the HYDA1 FeS cluster can only be built by the chloroplastic FeS cluster assembly machinery. Further, the expression of a bacterial hydrogenase gene, CPI, from the C. reinhardtii chloroplast genome resulted in H₂‐producing strains, demonstrating that a hydrogenase with a very different structure can fulfil the role of HYDA1 in vivo and that overexpression of foreign hydrogenases in C. reinhardtii is possible. All chloroplast transformants were stable and no toxic effects were seen from HYDA1 or CPI expression.     

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high‐throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole‐genome sequencing project, while the others were obtained by means of long‐range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference‐based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes (‘Brindisino’ artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty‐nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short ‘variable’ cp regions. The evaluation of the molecular resources obtained from this study led us to support the ‘super‐barcode’ theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.     

Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.