A chromosome-length genome assembly and annotation of blackberry (Rubus argutus, cv. "Hillquist").

Blackberries (Rubus spp.) are the fourth most economically important berry crop worldwide. Genome assemblies and annotations have been developed for Rubus species in subgenus Idaeobatus, including black raspberry (R. occidentalis), red raspberry (R. idaeus), and R. chingii, but very few genomic resources exist for blackberries and their relatives in subgenus Rubus. Here we present a chromosome-length assembly and annotation of the diploid blackberry germplasm accession "Hillquist" (R. argutus). "Hillquist" is the only known source of primocane-fruiting (annual-fruiting) in tetraploid fresh-market blackberry breeding programs and is represented in the pedigree of many important cultivars worldwide. The "Hillquist" assembly, generated using Pacific Biosciences long reads scaffolded with high-throughput chromosome conformation capture sequencing, consisted of 298 Mb, of which 270 Mb (90%) was placed on 7 chromosome-length scaffolds with an average length of 38.6 Mb. Approximately 52.8% of the genome was composed of repetitive elements. The genome sequence was highly collinear with a novel maternal haplotype-resolved linkage map of the tetraploid blackberry selection A-2551TN and genome assemblies of R. chingii and red raspberry. A total of 38,503 protein-coding genes were predicted, of which 72% were functionally annotated. Eighteen flowering gene homologs within a previously mapped locus aligning to an 11.2 Mb region on chromosome Ra02 were identified as potential candidate genes for primocane-fruiting. The utility of the "Hillquist" genome has been demonstrated here by the development of the first genotyping-by-sequencing-based linkage map of tetraploid blackberry and the identification of possible candidate genes for primocane-fruiting. This chromosome-length assembly will facilitate future studies in Rubus biology, genetics, and genomics and strengthen applied breeding programs.

Divergent subgenome evolution in the allotetraploid frog Xenopus calcaratus.

Allopolyploid genomes are divided into compartments called subgenomes that are derived from lower ploidy ancestors. In African clawed frogs of the subgenus Xenopus (genus Xenopus), allotetraploid species have two subgenomes (L and S) with morphologically distinct homoeologous chromosomes. In allotetraploid species of the sister subgenus Silurana, independently evolved subgenomes also exist, but their cytogenetics has not been investigated in detail. We used a diverse suite of cytogenetic and molecular FISH techniques on an allotetraploid species in Silurana-Xenopus calcaratus-to explore evolutionary dynamics of chromosome morphology and rearrangements. We find that the subgenomes of X. calcaratus have distinctive characteristics, with a more conserved a-subgenome resembling the closely related genome of the diploid species X. tropicalis, and a more rapidly evolving b-subgenome having more pronounced changes in chromosome structure, including diverged heterochromatic blocks, repetitive sequences, and deletion of a nucleolar secondary constriction. Based on these cytogenetic differences, we propose a chromosome nomenclature for X. calcaratus that may apply to other allotetraploids in subgenus Silurana, depending on as yet unresolved details of their evolutionary origins. These findings highlight the potential for large-scale asymmetry in subgenome evolution following allopolyploidization.

The genome sequence of the scarce swallowtail, Iphiclides podalirius.

The scarce swallowtail, Iphiclides podalirius (Linnaeus, 1758), is a species of butterfly in the family Papilionidae. Here, we present a chromosome-level genome assembly for Iphiclides podalirius as well as gene and transposable element annotations. We investigate how the density of genomic features differs between the 30 Iphiclides podalirius chromosomes. We find that shorter chromosomes have higher heterozygosity at four-fold-degenerate sites and a greater density of transposable elements. While the first result is an expected consequence of differences in recombination rate, the second suggests a counter-intuitive relationship between recombination and transposable element evolution. This high-quality genome assembly, the first for any species in the tribe Leptocircini, will be a valuable resource for population genomics in the genus Iphiclides and comparative genomics more generally.

Chromosome length ratio as a biomarker of DNA damage in cells exposed to high dose ionizing radiation.

The premature chromosome condensation (PCC) assay is considered as complementary bio-dosimetry tool for chromosome aberration assay and the PCC assay can be used to estimate high dose exposure. Though the PCC ring is considered as prospective biomarker, chromosome length ratio (ratio of longest and shortest chromosome length in PCC spreads) of chemically induced PCC is shown to be very good indicator of ionizing radiation. In view of this, an in-vitro study has been performed using PCC assay to suggest chromosome length ratio (LR) as potential bio-dosimeter induced by high dose ionizing radiation. Blood samples were collected from healthy subjects (n = 3) after prior consent and irradiated to ten different doses ranging between 0 and 20 Gy using 6 MV LINAC X-rays with dose rate of 5.6 Gy/min. Irradiated lymphocytes were cultured and calyculin induced PCC spreads were prepared. PCC spreads were captured using image analysis system and chromosome lengths were measured using open-source ImageJ software. For each dose, LR for 50 chromosome spreads were computed and mean LR value was calculated. LR varies between 6.0 ± 0.08 and 23.6 ± 0.55 for the dose range between 2 and 20 Gy. The dose response curve for LR was observed to be linear with y = 1.02x + 3.36, R2 = 0.97. Linear dose response relationship obtained in the present study confirms the prospective use of LR measurement. This study is first of its kind to examine chromosome length ratio as a biomarker of DNA damage in cells exposed to high dose X-ray exposure.

Genome collinearity analysis illuminates the evolution of donkey chromosome 1 and horse chromosome 5 in perissodactyls: A comparative study.

BACKGROUND: It is important to resolve the evolutionary history of species genomes as it has affected both genome organization and chromosomal architecture. The rapid innovation in sequencing technologies and the improvement in assembly algorithms have enabled the creation of highly contiguous genomes. DNA Zoo, a global organization dedicated to animal conservation, offers more than 150 chromosome-length genome assemblies. This database has great potential in the comparative genomics field. RESULTS: Using the donkey (Equus asinus asinus, EAS) genome provided by DNA Zoo as an example, the scaffold N50 length and Benchmarking Universal Single-Copy Ortholog score reached 95.5 Mb and 91.6%, respectively. We identified the cytogenetic nomenclature, corrected the direction of the chromosome-length sequence of the donkey genome, analyzed the genome-wide chromosomal rearrangements between the donkey and horse, and illustrated the evolution of the donkey chromosome 1 and horse chromosome 5 in perissodactyls. CONCLUSIONS: The donkey genome provided by DNA Zoo has relatively good continuity and integrity. Sequence-based comparative genomic analyses are useful for chromosome evolution research. Several previously published chromosome painting results can be used to identify the cytogenetic nomenclature and correct the direction of the chromosome-length sequence of new assemblies. Compared with the horse genome, the donkey chromosomes 1, 4, 20, and X have several obvious inversions, consistent with the results of previous studies. A 4.8 Mb inverted structure was first discovered in the donkey chromosome 25 and plains zebra chromosome 11. We speculate that the inverted structure and the tandem fusion of horse chromosome 31 and 4 are common features of non-caballine equids, which supports the correctness of the existing Equus phylogeny to an extent.

Delineating the Tnt1 Insertion Landscape of the Model Legume Medicago truncatula cv. R108 at the Hi-C Resolution Using a Chromosome-Length Genome Assembly.

Legumes are of great interest for sustainable agricultural production as they fix atmospheric nitrogen to improve the soil. Medicago truncatula is a well-established model legume, and extensive studies in fundamental molecular, physiological, and developmental biology have been undertaken to translate into trait improvements in economically important legume crops worldwide. However, M. truncatula reference genome was generated in the accession Jemalong A17, which is highly recalcitrant to transformation. M. truncatula R108 is more attractive for genetic studies due to its high transformation efficiency and Tnt1-insertion population resource for functional genomics. The need to perform accurate synteny analysis and comprehensive genome-scale comparisons necessitates a chromosome-length genome assembly for M. truncatula cv. R108. Here, we performed in situ Hi-C (48×) to anchor, order, orient scaffolds, and correct misjoins of contigs in a previously published genome assembly (R108 v1.0), resulting in an improved genome assembly containing eight chromosome-length scaffolds that span 97.62% of the sequenced bases in the input assembly. The long-range physical information data generated using Hi-C allowed us to obtain a chromosome-length ordering of the genome assembly, better validate previous draft misjoins, and provide further insights accurately predicting synteny between A17 and R108 regions corresponding to the known chromosome 4/8 translocation. Furthermore, mapping the Tnt1 insertion landscape on this reference assembly presents an important resource for M. truncatula functional genomics by supporting efficient mutant gene identification in Tnt1 insertion lines. Our data provide a much-needed foundational resource that supports functional and molecular research into the Leguminosae for sustainable agriculture and feeding the future.

Inferring chromosome radial organization from Hi-C data.

BACKGROUND: The nonrandom radial organization of eukaryotic chromosome territories (CTs) inside the nucleus plays an important role in nuclear functional compartmentalization. Increasingly, chromosome conformation capture (Hi-C) based approaches are being used to characterize the genome structure of many cell types and conditions. Computational methods to extract 3D arrangements of CTs from this type of pairwise contact data will thus increase our ability to analyze CT organization in a wider variety of biological situations. RESULTS: A number of full-scale polymer models have successfully reconstructed the 3D structure of chromosome territories from Hi-C. To supplement such methods, we explore alternative, direct, and less computationally intensive approaches to capture radial CT organization from Hi-C data. We show that we can infer relative chromosome ordering using PCA on a thresholded inter-chromosomal contact matrix. We simulate an ensemble of possible CT arrangements using a force-directed network layout algorithm and propose an approach to integrate additional chromosome properties into our predictions. Our CT radial organization predictions have a high correlation with microscopy imaging data for various cell nucleus geometries (lymphoblastoid, skin fibroblast, and breast epithelial cells), and we can capture previously documented changes in senescent and progeria cells. CONCLUSIONS: Our analysis approaches provide rapid and modular approaches to screen for alterations in CT organization across widely available Hi-C data. We demonstrate which stages of the approach can extract meaningful information, and also describe limitations of pairwise contacts alone to predict absolute 3D positions.

Pacific Biosciences assembly with Hi-C mapping generates an improved, chromosome-level goose genome.

BACKGROUND: The domestic goose is an economically important and scientifically valuable waterfowl; however, a lack of high-quality genomic data has hindered research concerning its genome, genetics, and breeding. As domestic geese breeds derive from both the swan goose (Anser cygnoides) and the graylag goose (Anser anser), we selected a female Tianfu goose for genome sequencing. We generated a chromosome-level goose genome assembly by adopting a hybrid de novo assembly approach that combined Pacific Biosciences single-molecule real-time sequencing, high-throughput chromatin conformation capture mapping, and Illumina short-read sequencing. FINDINGS: We generated a 1.11-Gb goose genome with contig and scaffold N50 values of 1.85 and 33.12 Mb, respectively. The assembly contains 39 pseudo-chromosomes (2n = 78) accounting for ∼88.36% of the goose genome. Compared with previous goose assemblies, our assembly has more continuity, completeness, and accuracy; the annotation of core eukaryotic genes and universal single-copy orthologs has also been improved. We have identified 17,568 protein-coding genes and a repeat content of 8.67% (96.57 Mb) in this genome assembly. We also explored the spatial organization of chromatin and gene expression in the goose liver tissues, in terms of inter-pseudo-chromosomal interaction patterns, compartments, topologically associating domains, and promoter-enhancer interactions. CONCLUSIONS: We present the first chromosome-level assembly of the goose genome. This will be a valuable resource for future genetic and genomic studies on geese.

Inference of Chromosome-Length Haplotypes Using Genomic Data of Three or a Few More Single Gametes.

Compared with genomic data of individual markers, haplotype data provide higher resolution for DNA variants, advancing our knowledge in genetics and evolution. Although many computational and experimental phasing methods have been developed for analyzing diploid genomes, it remains challenging to reconstruct chromosome-scale haplotypes at low cost, which constrains the utility of this valuable genetic resource. Gamete cells, the natural packaging of haploid complements, are ideal materials for phasing entire chromosomes because the majority of the haplotypic allele combinations has been preserved. Therefore, compared with the current diploid-based phasing methods, using haploid genomic data of single gametes may substantially reduce the complexity in inferring the donor's chromosomal haplotypes. In this study, we developed the first easy-to-use R package, Hapi, for inferring chromosome-length haplotypes of individual diploid genomes with only a few gametes. Hapi outperformed other phasing methods when analyzing both simulated and real single gamete cell sequencing data sets. The results also suggested that chromosome-scale haplotypes may be inferred by using as few as three gametes, which has pushed the boundary to its possible limit. The single gamete cell sequencing technology allied with the cost-effective Hapi method will make large-scale haplotype-based genetic studies feasible and affordable, promoting the use of haplotype data in a wide range of research.

The Incidence of Mosaicism for Individual Chromosome in Human Blastocysts Is Correlated With Chromosome Length.

Mosaicism, known as partial aneuploidies, mostly originates from mitotic errors during the post-zygotic stage; it consists of different cell lineages within a human embryo. The incidence of mosaicism has not been shown to correlate with maternal age, and its correlation with individual chromosome characteristics has not been well investigated. In this study, the results of preimplantation genetic testing for aneuploidy (PGT-A) derived from 4,036 blastocysts (930 IVF couples) were collected from 2015 to 2017. Via next-generation sequencing for comprehensive chromosome screening, embryo ploidy was identified as aneuploid, mosaic, and euploid. Total mosaicism was classified into two categories: "mosaic euploid/aneuploidy" (with mosaic aneuploidy between 20 and 80%) and "mosaic and aneuploidy" (a uniformly abnormal embryo superimposed with mosaic aneuploidies). Frequency of mosaicism was analyzed according to the function of chromosomal lengths, which divides involved chromosomes into three groups: group A (156-249 Mb), group B (102-145 Mb), and group C (51-90 Mb). The results show that the aneuploidy was more frequent in group C than in group A and group B (A: 23.7%, B: 35.1, 41.2%, p < 0.0001), while the mosaicism was more frequent in group A and group B than in group C [(Mosaic euploid/aneuploid) A: 14.6%, B: 12.4%, C: 9.9%, p < 0.0001; (mosaic and aneuploid) A: 21.3%, B: 22.9%, C: 18.9%, p < 0.0001; (Total mosaicism) A: 35.9%, B: 35.3%, C: 28.8%, p < 0.0001]. The significantly higher frequency of aneuploidy was on the shorter chromosome (< 90 Mb), and that of mosaicism was on the longer chromosomes (> 100 Mb). The length association did not reach significance in the patients with advanced age (≥ 36 years), and of the chromosome-specific mosaicism rate, the highest prevalence was on chromosome 14 (5.8%), 1 (5.7%), and 9 (5.6%). Although the length association was observed via group comparison, there may be affecting mechanisms other than chromosomes length. Eventually, twenty patients with mosaic embryo cryotransfers resulted in six live births. No significant correlation was observed between the transfer outcomes and chromosome length; however, the analysis was limited by small sample size.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics.

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

Chromosome-Range Whole-Genome High-Throughput Experimental Haplotyping by Single-Chromosome Microdissection.

Haplotype is fundamental genetic information; it provides essential information for deciphering the functional and etiological roles of genetic variants. As haplotype information is closely related to the functional and etiological impact of genetic variants, it is widely anticipated that haplotype information will be extremely valuable in a wide spectra of applications, including academic research, clinical diagnosis of genetic disease and in the pharmaceutical industry. Haplotyping is essential for LD (linkage disequilibrium) mapping, functional studies on cis-interactions, big data imputation, association studies, population studies, and evolutionary studies. Unfortunately, current sequencing technologies and genotyping arrays do not routinely deliver this information for each individual, but yield only unphased genotypes. Here, we describe a high-throughput and cost-effective experimental protocol to obtain high-resolution chromosomal haplotypes of each individual diploid (including human) genome by the single-chromosome microdissection and sequencing approach.

Electrophoretic karyotyping of Hypsizygus marmoreus and evaluation of variation among its basidiospores.

The molecular karyotype of Hypsizygus marmoreus was explored by contour-clamped homogeneous electric field gel electrophoresis. Eleven chromosomal bands were separated from the dikaryotic mycelia of H. marmoreus (strain Hm 3-10), and the chromosomes ranged in size from 1.9 to 5.8 Mb. The total genome size of the strain was estimated to be 36.3 Mb. The chromosome numbers were also confirmed by telomere fingerprinting, and 22 telomeric bands were identified. This result suggests that 11 chromosomes exist in Hm 3-10. The marker sequences for each chromosome were determined and were applied to identify each chromosome. Karyotyping and Southern blot analysis revealed that the size of chromosomes in the basidiospores were greatly different from those of parental dikaryon Hm 3-10 cells.