Diploid Genome Research Articles

Convergence and divergence of diploid and tetraploid cotton genomes.

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

A tip of the iceberg: genome survey indicated a complex evolutionary history of Garuga Roxb. species

Background Garuga Roxb. is a genus endemic to southwest China and other tropical regions in Southeast Asia facing risk of extinction due to the loss of tropical forests and changes in land use. Conducting a genome survey of G. forrestii contribute to a deeper understanding and conservation of the genus.ResultsThis study utilized genome survey of G. forrestii generated approximately 54.56 GB of sequence data, with approximately 112 × coverage. K-mer analysis indicated a genome size of approximately 0.48 GB, smaller than 0.52GB estimated by flow cytometry. The heterozygosity is of about 0.54%, and a repeat rate of around 51.54%. All the shotgun data were assembled into 339,729 scaffolds, with an N50 of 17,344 bp. The average content of guanine and cytosine was approximately 35.16%. A total of 330,999 SSRs were detected, with mononucleotide repeats being the most abundant at 70.16%, followed by dinucleotide repeats at 20.40%. We conducted a preliminary ploidy assessment using Smudgeplot and observed a clear bimodal distribution in G. forrestii at 1/2 relative coverage depth and total coverage depth (2n), suggesting a potential diploid genome structure. A pseudo chromosome of G. forrestii and a gemone of Boswellia sacra were used as reference genome to perform a primer population resequencing analysis within three Garuga species. Principal component analysis (PCA) indicated three distinct groups, but genome wide phylogenetics represented conflicting both between the dataset of different reference genomes and between maternal and nuclear genome.ConclusionIn summary, the genome of G. forrestii is small, and the phylogenetic relationships within the Garuga genus are complex. The genetic data presented in this study holds significant value for comprehensive whole-genome analyses, the evaluation of population genetic diversity, investigations into adaptive evolution, the advancement of artificial breeding efforts, and the support of species conservation and restoration initiatives. Ultimately, this research contributes to reinforcing the conservation and management of natural ecosystems, promoting biodiversity conservation, and advancing sustainable development.

Telomere-to-telomere assembly by preserving contained reads.

Automated telomere-to-telomere (T2T) de novo assembly of diploid and polyploid genomes remains a formidable task. A string graph is a commonly used assembly graph representation in the assembly algorithms. The string graph formulation employs graph simplification heuristics, which drastically reduce the count of vertices and edges. One of these heuristics involves removing the reads contained in longer reads. In practice, this heuristic occasionally introduces gaps in the assembly by removing all reads that cover one or more genome intervals. The factors contributing to such gaps remain poorly understood. In this work, we mathematically derived the frequency of observing a gap near a germline and a somatic heterozygous variant locus. Our analysis shows that (i) an assembly gap due to contained read deletion is an order of magnitude more frequent in Oxford Nanopore reads than PacBio HiFi reads due to differences in their read-length distributions, and (ii) this frequency decreases with an increase in the sequencing depth. Drawing cues from these observations, we addressed the weakness of the string graph formulation by developing the RAFT assembly algorithm. RAFT addresses the issue of contained reads by fragmenting reads and producing a more uniform read-length distribution. The algorithm retains spanned repeats in the reads during the fragmentation. We empirically demonstrate that RAFT significantly reduces the number of gaps using simulated datasets. Using real Oxford Nanopore and PacBio HiFi datasets of the HG002 human genome, we achieved a twofold increase in the contig NG50 and the number of haplotype-resolved T2T contigs compared to Hifiasm.

Graphasing: phasing diploid genome assembly graphs with single-cell strand sequencing

Haplotype information is crucial for biomedical and population genetics research. However, current strategies to produce de novo haplotype-resolved assemblies often require either difficult-to-acquire parental data or an intermediate haplotype-collapsed assembly. Here, we present Graphasing, a workflow which synthesizes the global phase signal of Strand-seq with assembly graph topology to produce chromosome-scale de novo haplotypes for diploid genomes. Graphasing readily integrates with any assembly workflow that both outputs an assembly graph and has a haplotype assembly mode. Graphasing performs comparably to trio phasing in contiguity, phasing accuracy, and assembly quality, outperforms Hi-C in phasing accuracy, and generates human assemblies with over 18 chromosome-spanning haplotypes.

A Dominance Hypothesis Argument for Historical Genetic Gains and the Fixation of Heterosis in Octoploid Strawberry.

Heterosis was the catalyst for the domestication of cultivated strawberry (Fragaria × ananassa), an interspecific hybrid species that originated in the 1700s. The hybrid origin was discovered because the phenotypes of spontaneous hybrids transgressed those of their parent species. The transgressions included fruit yield increases and other genetic gains in the twentieth century that sparked the global expansion of strawberry production. The importance of heterosis to the agricultural success of the hybrid species, however, has remained a mystery. Here we show that heterosis has disappeared (become fixed) among improved hybrids within a population (the California population) that has been under long-term selection for increased fruit yield, weight, and firmness. We found that the highest yielding hybrids are among the most highly inbred (59-79%), which seems counterintuitive for a highly heterozygous, outbreeder carrying heavy genetic loads. Although faint remnants of heterosis were discovered, the between-parent allele frequency differences and dispersed favorable dominant alleles necessary for heterosis have decreased nearly genome-wide within the California population. Conversely, heterosis was prevalent and significant among wide hybrids, especially for fruit count, a significant driver of genetic gains for fruit yield. We attributed the disappearance (fixation) of heterosis within the California population to increased homozygosity of favorable dominant alleles and inbreeding associated with selection, random genetic drift, and selective sweeps. Despite historical inbreeding, the highest yielding hybrids reported to-date are estimated to be heterozygous for 20,370-44,280 of 97,000-108,000 genes in the octoploid genome, the equivalent of an entire diploid genome or more.

A chromosome level reference genome of Diviner’s sage (Salvia divinorum) provides insight into salvinorin A biosynthesis

BackgroundDiviner’s sage (Salvia divinorum; Lamiaceae) is the source of the powerful hallucinogen salvinorin A (SalA). This neoclerodane diterpenoid is an agonist of the human Κ-opioid receptor with potential medical applications in the treatment of chronic pain, addiction, and post-traumatic stress disorder. Only two steps of the approximately twelve step biosynthetic sequence leading to SalA have been resolved to date.ResultsTo facilitate pathway elucidation in this ethnomedicinal plant species, here we report a chromosome level genome assembly. A high-quality genome sequence was assembled with an N50 value of 41.4 Mb and a BUSCO completeness score of 98.4%. The diploid (2n = 22) genome of ~ 541 Mb is comparable in size and ploidy to most other members of this genus. Two diterpene biosynthetic gene clusters were identified and are highly enriched in previously unidentified cytochrome P450s as well as crotonolide G synthase, which forms the dihydrofuran ring early in the SalA pathway. Coding sequences for other enzyme classes with likely involvement in downstream steps of the SalA pathway (BAHD acyl transferases, alcohol dehydrogenases, and O-methyl transferases) were scattered throughout the genome with no clear indication of clustering. Differential gene expression analysis suggests that most of these genes are not inducible by methyl jasmonate treatment.ConclusionsThis genome sequence and associated gene annotation are among the highest resolution in Salvia, a genus well known for the medicinal properties of its members. Here we have identified the cohort of genes responsible for the remaining steps in the SalA pathway. This genome sequence and associated candidate genes will facilitate the elucidation of SalA biosynthesis and enable an exploration of its full clinical potential.

Genealogical asymmetry under the IM model and a two-taxon test for gene flow.

Methods for detecting gene flow between populations often rely on asymmetry in the average length of particular genealogical branches, with the ABBA-BABA test being a well known example. Currently, asymmetry-based methods cannot be applied to a pair of populations and such analyses are instead performed using model-based methods. Here we investigate genealogical asymmetry under a two-population Isolation with Migration model. We focus on genealogies where the first coalescence event is between lineages sampled from different populations, as the external branches of these genealogies have equal expected length as long as there is no post-divergence gene flow. We show that unidirectional gene flow breaks this symmetry and results in the recipient population having longer external branches. We derive expectations for the probability of this genealogical asymmetry and propose a simple statistic (Am) to detect it from genome sequence data. Am provides a two-taxon test for gene flow that only requires a single unphased diploid genome from each population, with no outgroup information. We use analytic expectations and simulations to explore how recombination, unequal effective population sizes, bidirectional gene flow and background selection influence Am and find that the statistic provides unambiguous evidence for gene flow under a continent-island history. We estimate Am for genome sequence data from Heliconius butterflies and Odocoileus deer, generating results consistent with previous model-based analyses. Our work highlights a signal of gene flow overlooked to date and provides a method that complements existing approaches for investigating the demographic history of recently diverged populations.

A chromosome phased diploid genome assembly of African hunting dog (Lycaon pictus).

The African hunting dog (Lycaon pictus, 2n=78) once ranged over most sub-Saharan ecosystems except its deserts and rainforests. However as a result of (still ongoing) population declines, today they remain only as small fragmented populations. Furthermore, the future of the species remains unclear, due to both anthropogenic pressure as well as interactions with domestic dogs, thus their preservation is a conservation priority. On the tree of life, the hunting dog is basal to Canis and Cuon and forms a crown group with them, making it a useful species for comparative genomic studies. Here, we present a diploid chromosome level assembly of an African hunting dog. Assembled according to VGP guidelines from a combination of PacBio HiFi reads and HiC data, it is phased at the level of individual chromosomes. The maternal (pseudo)haplotype (mat) of our assembly has a length of 2.38 Gbp, and 99.36 % of the sequence is encompassed by 39 chromosomal scaffolds. The rest is included in only 36 unplaced short scaffolds. At the contig level, mat consists of only 166 contigs with an N50 of 39 Mbp. BUSCO analysis showed 95.4 % completeness based on Сarnivora conservative genes (carnivora_odb10). When compared to other available genomes from subtribe Canina, the quality of the assembly is excellent, typically between the 1st and 3rd depending on the parameter used, and a significant improvement on previously published genomes for the species. We hope this assembly will play an important role in future conservation efforts and comparative studies of canid genomes.

Replacing sperm with genotyped haploid androgenetic blastomeres to generate cattle with predetermined paternal genomes

Abstract Although meiosis plays an essential role for the survival of species in natural selection, the genetic diversity resulting from sexual reproduction impedes human-driven strategies to transmit the most suitable genomes for genetic improvement, forcing breeders to select diploid genomes generated after fertilization, that is, after the encounter of sperm and oocytes carrying unknown genomes. To determine whether genomic assessment could be used before fertilization, some androgenetic haploid morula-stage bovine embryos derived from individual sperm were biopsied for genomic evaluation and others used to reconstruct “semi-cloned” (SC) diploid zygotes by the intracytoplasmic injection into parthenogenetically activated oocytes, and the resulting embryos were transferred to surrogate females to obtain gestations. Compared to controls, in vitro development to the blastocyst stage was lower and fewer surrogates became pregnant from the transfer of SC embryos. However, fetometric measurements of organs and placental membranes of all SC conceptuses were similar to controls, suggesting a normal post-implantation development. Moreover, transcript amounts of imprinted genes IGF2, IGF2R, PHLDA2, SNRPN and KCNQ1OT1 and methylation pattern of the KCNQ1 DMR were unaltered in SC conceptuses. Overall, this study shows that sperm can be replaced by genotyped haploid embryonic-derived cells to produce bovine embryos carrying a predetermined paternal genome and viable first trimester fetuses after transfer to female recipients.

A genome assembly for the Chryxus Arctic (Oeneis chryxus), the highest butterfly in North America.

We describe a highly contiguous and complete diploid genome assembly for the Chryxus Arctic, Oeneis chryxus (E. Doubleday, [1849]), a butterfly species complex spanning much of northern and western North America. One subspecies, the Ivallda Arctic (O. c. ivallda), is endemic to California's Sierra Nevada and of particular biogeographic interest and conservation concern. Extreme alpine habitats occupied by this subspecies include the summit of Mt. Whitney, California, representing the highest elevation butterfly population in North America. The assembly presented here consists of two haplotypes, 738.92 and 770.85Mb in length, with contig N50 values of 10.49 and 10.13Mb, scaffold N50 values of 25.35 and 25.69Mb, scaffold L50 values of 13 and 14, and BUSCO completeness scores of 96.5 and 98.3%, respectively. More than 97% of the assembly is organized into 29 scaffolds, which likely represent whole chromosomes. This assembly is the first major genomic resource for Oeneis, providing a foundational reference for future genomic studies on the taxonomy, evolutionary history, and conservation of the genus. As part of the California Conservation Genomics Project, we will use this assembly in conjunction with short-read resequencing to resolve patterns of evolutionary differentiation, adaptive genomic variation, and gene flow among remaining O. c. ivallda populations. These data can and will be used to inform the subspecies' conservation as warming climatic conditions continue to lead to the loss and fragmentation of alpine habitats. We also provide genome assemblies for the O. chryxus mitochondrion and a Wolbachia endosymbiont.

Analysis of CYP701A1 genes in Gossypium species and functional characterization through gene silencing

Gibberellins (GA) are known to play crucial roles in various aspects of plant growth and development. The cytochrome P450 enzyme family is recognized for its significance in plant metabolic processes. Specifically, CYP701s, a subgroup of CYP71, encode ent-kaurene oxidase in the gibberellin synthesis pathway. In this study, we analyzed genomic data from 30 Gossypium species, including nine allotetraploid genomes (AD1-AD7, with two each for AD1 and AD2), 21 diploid genomes (A-G, K, with two A-genomes and 12 D-genomes), and Gossypioides kirkii genome as an outgroup for evolutionary analysis, totaling 31 genomes. Subsequently, 40 CYP701A1 genes were identified from various genomes and conducted a comprehensive analysis of their structure and evolution. Virus-induced gene silencing (VIGS) technology was utilized to knock out the GhCYP701A1 gene in Gossypium hirsutum ac TM-1. Subsequent analysis revealed changes in hormone content, with decreased gibberellin levels and notable increases in auxin, cytokinin, and jasmonic acid contents. Conversely, salicylic acid content decreased, while the precursor for ethylene synthesis, 1-aminocyclopropane-1-carboxylic acid (ACC), remained relatively stable. Transcriptome analysis of the gene silencing plants identified 15,962 differentially expressed genes, including 8,376 upregulated and 7,586 downregulated genes. Enrichment analysis through KEGG pathway analysis highlighted 'Plant hormone signal transduction' as a prominent pathway with 234 differentially expressed genes. The study provided insights into the function and regulatory network of the gene.

Consistent Delivery of Adeno-Associated Virus via Lateral Tail-Vein Injection in Adult Mice.

Many disorders affect multiple organs or involve different regions of the body, so it is critical to deliver therapeutics systemically to target the affected cells located in different sites. Intravenous injection is a widely used systemic delivery route in preclinical studies that assess treatments intended for body-wide administration. In adult mice, it involves the intravenous administration of the therapeutic agent into the mouse's lateral tail veins. When mastered, tail-vein injections are safe and fast, and only require simple and commonly available tools. However, tail-vein injections are technically challenging and require extensive training and continuous practice to ensure the accurate delivery of the intended dose. Here we describe a detailed, optimized, lateral tail-vein injection protocol that we have developed based on our experience and on recommendations that had been previously reported by other groups. Other than the mouse restrainers and insulin syringes, this protocol requires only reagents and equipment that are readily available in most labs. We found that following this protocol results in consistently successful intravenous delivery of adeno-associated virus (AAV) into the tail veins of unsedated 7-9 week-old mice. Additionally, we describe the optimized protocols for the histological detection of fluorescent reporter proteins and vector genome per diploid genome (vg/dg) quantification used to assess AAV transduction and biodistribution. The goal of this protocol is to aid experimenters in easily performing tail-vein injections successfully and consistently, which can reduce the practice time needed to master the technique.

VolcanoSV enables accurate and robust structural variant calling in diploid genomes from single-molecule long read sequencing

Structural variants (SVs) significantly contribute to human genome diversity and play a crucial role in precision medicine. Although advancements in single-molecule long-read sequencing offer a groundbreaking resource for SV detection, identifying SV breakpoints and sequences accurately and robustly remains challenging. We introduce VolcanoSV, an innovative hybrid SV detection pipeline that utilizes both a reference genome and local de novo assembly to generate a phased diploid assembly. VolcanoSV uses phased SNPs and unique k-mer similarity analysis, enabling precise haplotype-resolved SV discovery. VolcanoSV is adept at constructing comprehensive genetic maps encompassing SNPs, small indels, and all types of SVs, making it well-suited for human genomics studies. Our extensive experiments demonstrate that VolcanoSV surpasses state-of-the-art assembly-based tools in the detection of insertion and deletion SVs, exhibiting superior recall, precision, F1 scores, and genotype accuracy across a diverse range of datasets, including low-coverage (10x) datasets. VolcanoSV outperforms assembly-based tools in the identification of complex SVs, including translocations, duplications, and inversions, in both simulated and real cancer data. Moreover, VolcanoSV is robust to various evaluation parameters and accurately identifies breakpoints and SV sequences.

A familial, telomere-to-telomere reference for human de novo mutation and recombination from a four-generation pedigree.

Using five complementary short- and long-read sequencing technologies, we phased and assembled >95% of each diploid human genome in a four-generation, 28-member family (CEPH 1463) allowing us to systematically assess de novo mutations (DNMs) and recombination. From this family, we estimate an average of 192 DNMs per generation, including 75.5 de novo single-nucleotide variants (SNVs), 7.4 non-tandem repeat indels, 79.6 de novo indels or structural variants (SVs) originating from tandem repeats, 7.7 centromeric de novo SVs and SNVs, and 12.4 de novo Y chromosome events per generation. STRs and VNTRs are the most mutable with 32 loci exhibiting recurrent mutation through the generations. We accurately assemble 288 centromeres and six Y chromosomes across the generations, documenting de novo SVs, and demonstrate that the DNM rate varies by an order of magnitude depending on repeat content, length, and sequence identity. We show a strong paternal bias (75-81%) for all forms of germline DNM, yet we estimate that 17% of de novo SNVs are postzygotic in origin with no paternal bias. We place all this variation in the context of a high-resolution recombination map (~3.5 kbp breakpoint resolution). We observe a strong maternal recombination bias (1.36 maternal:paternal ratio) with a consistent reduction in the number of crossovers with increasing paternal (r=0.85) and maternal (r=0.65) age. However, we observe no correlation between meiotic crossover locations and de novo SVs, arguing against non-allelic homologous recombination as a predominant mechanism. The use of multiple orthogonal technologies, near-telomere-to-telomere phased genome assemblies, and a multi-generation family to assess transmission has created the most comprehensive, publicly available "truth set" of all classes of genomic variants. The resource can be used to test and benchmark new algorithms and technologies to understand the most fundamental processes underlying human genetic variation.

Simultaneous de novo calling and phasing of genetic variants at chromosome-scale using NanoStrand-seq

The successful accomplishment of the first telomere-to-telomere human genome assembly, T2T-CHM13, marked a milestone in achieving completeness of the human reference genome. The upcoming era of genome study will focus on fully phased diploid genome assembly, with an emphasis on genetic differences between individual haplotypes. Most existing sequencing approaches only achieved localized haplotype phasing and relied on additional pedigree information for further whole-chromosome scale phasing. The short-read-based Strand-seq method is able to directly phase single nucleotide polymorphisms (SNPs) at whole-chromosome scale but falls short when it comes to phasing structural variations (SVs). To shed light on this issue, we developed a Nanopore sequencing platform-based Strand-seq approach, which we named NanoStrand-seq. This method allowed for de novo SNP calling with high precision (99.52%) and acheived a superior phasing accuracy (0.02% Hamming error rate) at whole-chromosome scale, a level of performance comparable to Strand-seq for haplotype phasing of the GM12878 genome. Importantly, we demonstrated that NanoStrand-seq can efficiently resolve the MHC locus, a highly polymorphic genomic region. Moreover, NanoStrand-seq enabled independent direct calling and phasing of deletions and insertions at whole-chromosome level; when applied to long genomic regions of SNP homozygosity, it outperformed the strategy that combined Strand-seq with bulk long-read sequencing. Finally, we showed that, like Strand-seq, NanoStrand-seq was also applicable to primary cultured cells. Together, here we provided a novel methodology that enabled interrogation of a full spectrum of haplotype-resolved SNPs and SVs at whole-chromosome scale, with broad applications for species with diploid or even potentially polypoid genomes.

The haplotype-resolved Prymnesium parvum (type B) microalga genome reveals the genetic basis of its fish-killing toxins

The catastrophic loss of aquatic life in the Central European Oder River in 2022, caused by a toxic bloom of the haptophyte microalga Prymnesium parvum (in a wide sense, s.l.), underscores the need to improve our understanding of the genomic basis of the toxin. Previous morphological, phylogenetic, and genomic studies have revealed cryptic diversity within P.parvum s.l. and uncovered three clade-specific (types A, B, and C) prymnesin toxins. Here, we used state-of-the-art long-read sequencing and assembled the first haplotype-resolved diploid genome of a P.parvum type B from the strain responsible for the Oder disaster. Comparative analyses with type A genomes uncovered a genome-size expansion driven by repetitive elements in type B. We also found conserved synteny but divergent evolution in several polyketide synthase (PKS) genes, which are known to underlie toxin production in combination with environmental cues. We identified an approximately 20-kbp deletion in the largest PKS gene of type B that we link to differences in the chemical structure of types A and B prymnesins. Flow cytometry and electron microscopy analyses confirmed diploidy in the Oder River strain and revealed differences to closely related strains in both ploidy and morphology. Our results provide unprecedented resolution of strain diversity in P.parvum s.l. and a better understanding of the genomic basis of toxin variability in haptophytes. The reference-quality genome will enable us to better understand changes in microbial diversity in the face of increasing environmental pressures and provides a basis for strain-level monitoring of invasive Prymnesium in the future.

Estimating Haplotype Structure and Frequencies: A Bayesian Approach to Unknown Design in Pooled Genomic Data.

The estimation of haplotype structure and frequencies provides crucial information about the composition of genomes. Techniques, such as single-individual haplotyping, aim to reconstruct individual haplotypes from diploid genome sequencing data. However, our focus is distinct. We address the challenge of reconstructing haplotype structure and frequencies from pooled sequencing samples where multiple individuals are sequenced simultaneously. A frequentist method to address this issue has recently been proposed. In contrast to this and other methods that compute point estimates, our proposed Bayesian hierarchical model delivers a posterior that permits us to also quantify uncertainty. Since matching permutations in both haplotype structure and corresponding frequency matrix lead to the same reconstruction of their product, we introduce an order-preserving shrinkage prior that ensures identifiability with respect to permutations. For inference, we introduce a blocked Gibbs sampler that enforces the required constraints. In a simulation study, we assessed the performance of our method. Furthermore, by using our approach on two distinct sets of real data, we demonstrate that our Bayesian approach can reconstruct the dominant haplotypes in a challenging, high-dimensional set-up.

O-066 Parental inheritance of monopronuclear blastocysts from ICSI cycles and its association with specific morphokinetic traits: a retrospective observational study

Abstract Study question Are monopronuclear-originated blastocysts from ICSI prone to have uniparental genome and does it correlate with specific morphokinetic parameters? Summary answer Extended compaction was notably observed in parthenogenetic monopronuclear blastocysts, suggesting a potential correlation between morphokinetic traits and parental origin. What is known already Previous studies indicated monopronuclear embryos with larger pronucleus exhibit a propensity to develop into blastocysts as normal fertilized embryos, which was considered as evidence supporting the presence of a diploid genome enclosed within a singular pronucleus.However, diploid blastocysts may also arise from haploid spontaneous diploidization following parthenogenetic or androgenetic activation, particularly in ICSI-derived blastocysts. Although certain characteristics were discerned in prior morphokinetic studies on monopronuclear blastocysts, these embryos were not proceeded to comprehensive chromosomal analysis to authenticate the association with parental inheritance, which requires further clarification and investigation. Study design, size, duration This retrospective cohort study, conducted in a university-affiliated reproductive medical center, focused on ICSI-derived monopronuclear blastocysts (n = 151) obtained from 145 PGT cycles spanning January 2013 to December 2022. Among these, nineteen were subject to morphokinetic analysis within a time-lapse system, while 35 bipronuclear-originated sibling blastocysts from the same cycles comprised the control group. Notably, the morphokinetic analysis incorporated 10 uniparental monopronuclear blastocysts, with 7 originating from parthenogenesis and the remaining from androgenesis. Participants/materials, setting, methods The data was collected from patients who had at least one biopsied ICSI-derived monopronuclear blastocyst between 2013 and 2022 at the reproductive center. The genetic analysis of blastocysts were performed using SNP array or NGS. Morphokinetic recordings were obtained from blastocysts cultured in time-lapse system along with euploid bipronuclear sibling blastocysts in the same cycles. All annotations underwent manual scrutiny and standardization by the same embryologist to ensure consistency. Main results and the role of chance This study explored links between genetic origin and morphokinetics in ICSI-derived monopronuclear blastocyst through combine analysis. Among 151 blastocysts, 53 (35.1%) were biparental euploidies, while 20% demonstrated uniparental patterns, with 18 confirmed as parthenogenetic and 7 as androgenetic. The ratio of parthenogenesis to androgenesis was 2.6:1, aligning with reported figures. Pronuclear metrics revealed that pronuclei area of monopronuclear blastocysts were notably larger (91∼171 μm2) than that of maternal-origin pronuclei but equivalent to paternal-origin pronuclei from bipronuclear blastocysts.Kinetic analysis indicated that pronuclei of monopronuclear blastocysts vanished 1.9 hours earlier than those of bipronuclear blastocysts, and furtherly biparental monopronuclear blastocysts had a 2.74-hour shorter visible duration compared to bipronuclear blastocysts, whereas uniparental monopronulcei showed no significant difference. Following pronucleus disappearance, the compaction period of monopronuclear blastocysts lasted 2.85 hours longer than bipronuclear blastocysts. Notably, in parthenogenetic uniparental monopronuclear blastocysts, this gap expanded to 10.92 hours. Additionally, there was a significantly higher proportion of asymmetrical halo and dysmorphism in monopronuclear blastocysts (78.95%, 63.16%) compared to bipronuclear blastocysts (40.00%, 34.29%), but neither differences was related to parental origins. The study corroborated a high incidence of uniparental diploidy and identified the discrepancies in pronuclear parameters and compaction before blastulation of monopronuclear blastocysts with different parental origins. Limitations, reasons for caution The scarcity of ICSI-derived monopronuclear blastocysts undergoing PGT and time-lapse culture posed challenges in discerning additional morphokinetic distinctions related to parental inheritance and conducting in-depth analyses for androgenetic blastocysts. It is noted that parental ploidy results from NGS were deduced from targeted chromosome, which may not represent the whole genome. Wider implications of the findings Our results confirmed the high risk of uniparental diploidy in ICSI-derived monopronuclear blastocysts, suggesting the necessity of comprehensive screening. We pinpointed morphokinetic features in monopronuclear embryos associated with parental origin. While it remains unreliable to gauge parental compositions through morphokinetics, these distinctions could serve as supplementary indicators evaluating monopronuclear embryos. Trial registration number not applicable

MethPhaser: methylation-based long-read haplotype phasing of human genomes

The assignment of variants across haplotypes, phasing, is crucial for predicting the consequences, interaction, and inheritance of mutations and is a key step in improving our understanding of phenotype and disease. However, phasing is limited by read length and stretches of homozygosity along the genome. To overcome this limitation, we designed MethPhaser, a method that utilizes methylation signals from Oxford Nanopore Technologies to extend Single Nucleotide Variation (SNV)-based phasing. We demonstrate that haplotype-specific methylations extensively exist in Human genomes and the advent of long-read technologies enabled direct report of methylation signals. For ONT R9 and R10 cell line data, we increase the phase length N50 by 78%-151% at a phasing accuracy of 83.4-98.7% To assess the impact of tissue purity and random methylation signals due to inactivation, we also applied MethPhaser on blood samples from 4 patients, still showing improvements over SNV-only phasing. MethPhaser further improves phasing across HLA and multiple other medically relevant genes, improving our understanding of how mutations interact across multiple phenotypes. The concept of MethPhaser can also be extended to non-human diploid genomes. MethPhaser is available at https://github.com/treangenlab/methphaser.

Phasing Diploid Genome Assembly Graphs with Single-Cell Strand Sequencing.

Haplotype information is crucial for biomedical and population genetics research. However, current strategies to produce de-novo haplotype-resolved assemblies often require either difficult-to-acquire parental data or an intermediate haplotype-collapsed assembly. Here, we present Graphasing, a workflow which synthesizes the global phase signal of Strand-seq with assembly graph topology to produce chromosome-scale de-novo haplotypes for diploid genomes. Graphasing readily integrates with any assembly workflow that both outputs an assembly graph and has a haplotype assembly mode. Graphasing performs comparably to trio-phasing in contiguity, phasing accuracy, and assembly quality, outperforms Hi-C in phasing accuracy, and generates human assemblies with over 18 chromosome-spanning haplotypes.