Segmental Duplication Regions Research Articles

Detection of an 8p23.1 Inversion Using High-Resolution Optical Genome Mapping

Abstract Objective To evaluate the performance of optical genome mapping (OGM) in identifying an inversion located in the short arm of chromosome 8 (8p, 8p23.1), flanked by regions of complex segmental duplication (SD), using the GRCh38 and telomere-to-telomere (T2T) genome references. Methods We investigated a couple suspected of carrying the 8p23.1 inversion due to a terminal deletion combined with an interstitial duplication of 8p found in their abortus. OGM was performed on both individuals. The data were mapped to the current GRCh38 and the updated T2T genome references, respectively. Results The 8p23.1 inversion was observed in the female when mapping OGM data to the T2T assembly. In contrast, under the GRCh38 reference, the orientation between the suspected breakpoints within the SD regions could not be distinguished. Additional variants of uncertain significance were also identified in both individuals. Conclusion Our findings highlight the superiority of the T2T reference in recognizing structural variations involving SD regions. The enhanced SV detection using the T2T reference may contribute to a better understanding of genome instability and human diseases.

Optical Genome Mapping for Chromosomal Aberrations Detection-False-Negative Results and Contributing Factors.

Optical genome mapping (OGM) has been known as an all-in-one technology for chromosomal aberration detection. However, there are also aberrations beyond the detection range of OGM. This study aimed to report the aberrations missed by OGM and analyze the contributing factors. OGM was performed by taking both GRCh37 and GRCh38 as reference genomes. The OGM results were analyzed in blinded fashion and compared to standard assays. Quality control (QC) metrics, sample types, reference genome, effective coverage and classes and locations of aberrations were then analyzed. In total, 154 clinically reported variations from 123 samples were investigated. OGM failed to detect 10 (6.5%, 10/154) aberrations with GRCh37 assembly, including five copy number variations (CNVs), two submicroscopic balanced translocations, two pericentric inversion and one isochromosome (mosaicism). All the samples passed pre-analytical and analytical QC. With GRCh38 assembly, the false-negative rate of OGM fell to 4.5% (7/154). The breakpoints of the CNVs, balanced translocations and inversions undetected by OGM were located in segmental duplication (SD) regions or regions with no DLE-1 label. In conclusion, besides variations with centromeric breakpoints, structural variations (SVs) with breakpoints located in large repetitive sequences may also be missed by OGM. GRCh38 is recommended as the reference genome when OGM is performed. Our results highlight the necessity of fully understanding the detection range and limitation of OGM in clinical practice.

A near-complete assembly of asparagus bean provides insights into anthocyanin accumulation in pods.

Asparagus bean (Vigna unguiculata ssp. sesquipedialis), a subspecies of V. unguiculata, is a vital legume crop widely cultivated in Asia for its tender pods consumed as vegetables. However, the existing asparagus bean assemblies still contain numerous gaps and unanchored sequences, which presents challenges to functional genomics research. Here, we present an improved reference genome sequence of an elite asparagus bean variety, Fengchan 6, achieved through the integration of nanopore ultra-long reads, PacBio high-fidelity reads, and Hi-C technology. The improved assembly is 521.3 Mb in length and demonstrates several enhancements, including a higher N50 length (46.4 Mb), an anchor ratio of 99.8%, and the presence of only one gap. Furthermore, we successfully assembled 14 telomeres and all 11 centromeres, including four telomere-to-telomere chromosomes. Remarkably, the centromeric regions cover a total length of 38.1 Mb, providing valuable insights into the complex architecture of centromeres. Among the 30 594 predicted protein-coding genes, we identified 2356 genes that are tandemly duplicated in segmental duplication regions. These findings have implications for defence responses and may contribute to evolutionary processes. By utilizing the reference genome, we were able to effectively identify the presence of the gene VuMYB114, which regulates the accumulation of anthocyanins, thereby controlling the purple coloration of the pods. This discovery holds significant implications for understanding the underlying mechanisms of color determination and the breeding process. Overall, the highly improved reference genome serves as crucial resource and lays a solid foundation for asparagus bean genomic studies and genetic improvement efforts.

Characterization of large-scale genomic differences in the first complete human genome

BackgroundThe first telomere-to-telomere (T2T) human genome assembly (T2T-CHM13) release is a milestone in human genomics. The T2T-CHM13 genome assembly extends our understanding of telomeres, centromeres, segmental duplication, and other complex regions. The current human genome reference (GRCh38) has been widely used in various human genomic studies. However, the large-scale genomic differences between these two important genome assemblies are not characterized in detail yet.ResultsHere, in addition to the previously reported “non-syntenic” regions, we find 67 additional large-scale discrepant regions and precisely categorize them into four structural types with a newly developed website tool called SynPlotter. The discrepant regions (~ 21.6 Mbp) excluding telomeric and centromeric regions are highly structurally polymorphic in humans, where the deletions or duplications are likely associated with various human diseases, such as immune and neurodevelopmental disorders. The analyses of a newly identified discrepant region—the KLRC gene cluster—show that the depletion of KLRC2 by a single-deletion event is associated with natural killer cell differentiation in ~ 20% of humans. Meanwhile, the rapid amino acid replacements observed within KLRC3 are probably a result of natural selection in primate evolution.ConclusionOur study provides a foundation for understanding the large-scale structural genomic differences between the two crucial human reference genomes, and is thereby important for future human genomics studies.

Variants detection in regions of segmental duplication facilitates molecular diagnosis for IEI such as chronic granulomatous disease caused by NCF1 mutations

With the expanding application of next-generation sequencing (NGS) in clinical practice, rapid progress has been made in molecular diagnoses for diseases such as inborn errors of immunity (IEI). However, due to characteristic short read length and fragment size for NGS, genomic regions with segmental duplication pose a big challenge for accurate mutation detection. Reads originating from segmental duplications are mapped with high ambiguity and are often discarded by most variant callers, leading to poor variant identification. Dozens of known IEI causal genes are located in regions of segmental duplications, such as NCF1, IGLL1, FCGR3A, etc. It is important to increase the sensitivity of mutation detections in such regions.In this study, we develop a strategy to significantly increase the sensitivity of identifying rare germline variants in regions of segmental duplication for IEI patients while controlling the number of false positives. By forcing multi-aligned reads from segmental duplications to map to functionally significant coding regions in the region concerned, all possible small variants are being identified. This is followed by the identification of false positives originating from intrinsic sequence differences among genomic duplications and the identification of potential common variants by matching with an in-house database.Using simulated data and real-world data from the 'Genome In a Bottle’ project, we demonstrate a significant improvement in variant detection sensitivity, ranging from 80% to 98% and 75% to 94%, respectively. This is accompanied by a limited increase in false positives that can be dealt with by downstream analysis. Applying our method to whole exome sequencing data led to the identification of the common ΔGT mutation in NCF1 in three patients with the chronic granulomatous disease, demonstrating the power of this novel approach.In summary, we developed SDrecall, a bioinformatic method that is by far the most robust open-source approach to detect rare, disease-causing variants in regions of segmental duplication. The method should apply to both IEI and other Mendelian diseases with causal mutations located in such regions and is suited to WGS, WES, and targeted sequencing data.

Promoter evolution of mammalian gene duplicates

BackgroundGene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.ResultsWe observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture—paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set.Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.ConclusionsIn this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

NPGREAT: assembly of human subtelomere regions with the use of ultralong nanopore reads and linked-reads

BackgroundHuman subtelomeric DNA regulates the length and stability of adjacent telomeres that are critical for cellular function, and contains many gene/pseudogene families. Large evolutionarily recent segmental duplications and associated structural variation in human subtelomeres has made complete sequencing and assembly of these regions difficult to impossible for many loci, complicating or precluding a wide range of genetic analyses to investigate their function.ResultsWe present a hybrid assembly method, NanoPore Guided REgional Assembly Tool (NPGREAT), which combines Linked-Read data with mapped ultralong nanopore reads spanning subtelomeric segmental duplications to potentially overcome these difficulties. Linked-Read sets of DNA sequences identified by matches with 1-copy subtelomere sequence adjacent to segmental duplications are assembled and extended into the segmental duplication regions using Regional Extension of Assemblies using Linked-Reads (REXTAL). Mapped telomere-containing ultralong nanopore reads are then used to provide contiguity and correct orientation for matching REXTAL sequence contigs as well as identification/correction of any misassemblies. Our method was tested for a subset of representative subtelomeres with ultralong nanopore read coverage in the haploid human cell line CHM13. A 10X Linked-Read dataset from CHM13 was combined with ultralong nanopore reads from the same genome to provide improved subtelomere assemblies. Comparison of Nanopore-only assemblies using SHASTA with our NPGREAT assemblies in the distal-most subtelomere regions showed that NPGREAT produced higher-quality and more complete assemblies than SHASTA alone when these regions had low ultralong nanopore coverage (such as cases where large segmental duplications were immediately adjacent to (TTAGGG) tracts).ConclusionIn genomic regions with large segmental duplications adjacent to telomeres, NPGREAT offers an alternative economical approach to improving assembly accuracy and coverage using linked-read datasets when more expensive HiFi datasets of 10–20 kb reads are unavailable.

VP.01 Array comparative genomic hybridisation and droplet digital PCR uncover recurrent copy number variation of the titin segmental duplication region

The sarcomeric giants nebulin and titin both contain intragenic segmental duplication regions. Segmental duplication regions are prevalent throughout the genome and are known potential hotspots for recurrent copy number variation and may harbour pathogenic aberrations. Using our custom comparative genomic hybridization array, we have shown that gains and losses of more than one copy of the repeated block of the nebulin triplicate region are a recessive pathogenic mutation contributing to the spectrum of nemaline myopathy-causing changes. Similar copy number variants have been detected in the segmental duplication region of titin, harboured within its PEVK domain. The exact copy number determination of the segmental duplication region in titin is challenging due to the limitations of the methodology and repetitiveness of the region. We developed three complementary custom droplet digital PCR assays for the titin segmental duplication region to confirm true variation. Two of the assays target the repeated region, and one targets a diploid sequence adjacent to the segmental duplication. By comparing the data acquired by the custom comparative genomic hybridization array and the custom droplet digital PCR assays, we assessed the copy number variation of the region in 55 samples from both healthy individuals and samples from patients with different congenital muscle disorders. Our combined methods show that the titin segmental duplication region is subject to recurrent copy number variation. Gains and losses of the region were detected in both samples from healthy individuals and in samples from patients. While the copy number variation observed in our cohort is mostly unlikely to be pathogenic, it is possible that this region of titin harbours pathogenic copy number variation in a fashion similar to that of nebulin. Putative pathogenic thresholds for gains and losses remain to be studied, but the region should no longer be neglected in genetic analyses.

Array Comparative Genomic Hybridisation and Droplet Digital PCR Uncover Recurrent Copy Number Variation of the TTN Segmental Duplication Region.

Intragenic segmental duplication regions are potential hotspots for recurrent copy number variation and possible pathogenic aberrations. Two large sarcomeric genes, nebulin and titin, both contain such segmental duplication regions. Using our custom Comparative Genomic Hybridisation array, we have previously shown that a gain or loss of more than one copy of the repeated block of the nebulin triplicate region constitutes a recessive pathogenic mutation. Using targeted array-CGH, similar copy number variants can be detected in the segmental duplication region of titin. Due to the limitations of the array-CGH methodology and the repetitiveness of the region, the exact copy numbers of the blocks could not be determined. Therefore, we developed complementary custom Droplet Digital PCR assays for the titin segmental duplication region to confirm true variation. Our combined methods show that the titin segmental duplication region is subject to recurrent copy number variation. Gains and losses were detected in samples from healthy individuals as well as in samples from patients with different muscle disorders. The copy number variation observed in our cohort is likely benign, but pathogenic copy number variants in the segmental duplication region of titin cannot be excluded. Further investigations are needed, however, this region should no longer be neglected in genetic analyses.

A custom ddPCR method for the detection of copy number variations in the nebulin triplicate region.

The human genome contains repetitive regions, such as segmental duplications, known to be prone to copy number variation. Segmental duplications are highly identical and homologous sequences, posing a specific challenge for most mutation detection methods. The giant nebulin gene is expressed in skeletal muscle. It harbors a large segmental duplication region composed of eight exons repeated three times, the so-called triplicate region. Mutations in nebulin are known to cause nemaline myopathy and other congenital myopathies. Using our custom targeted Comparative Genomic Hybridization arrays, we have previously shown that copy number variations in the nebulin triplicate region are pathogenic when the copy number of the segmental duplication block deviates two or more copies from the normal number, which is three per allele. To complement our Comparative Genomic Hybridization arrays, we have established a custom Droplet Digital PCR method for the detection of copy number variations within the nebulin triplicate region. The custom Droplet Digital PCR assays allow sensitive, rapid, high-throughput, and cost-effective detection of copy number variations within this region and is ready for implementation a screening method for disease-causing copy number variations of the nebulin triplicate region. We suggest that Droplet Digital PCR may also be used in the study and diagnostics of other segmental duplication regions of the genome.

Comprehensive identification of the VQ family genes in cucumber and their roles in response to abiotic and biotic stresses

Valine-glutamine (VQ) motif-containing proteins are plant-specific transcription factors with critical functions in plant growth, development and defense reactions. However, no systematic study has been performed for VQ family genes in cucumber to date. In this work, a total of 32 VQ family genes (CsVQ1–CsVQ32) were identified in cucumber, which were unevenly distributed on seven cucumber chromosomes (chromosomes 1–7). The CsVQs were divided into nine distinct groups (group I–IX) according to their phylogenetic relationships, with the members of the same group sharing highly conserved motif compositions. Gene structure analysis suggested that the CsVQ genes were generally lack of introns. Gene duplication analysis revealed that two and 11 CsVQ genes were located in the tandem and segmental duplication regions of cucumber chromosomes, respectively. Tissue expression profiles obtained from RNA-seq data showed that most CsVQ genes were differentially expressed in different tissues and developmental stages. Quantitative real-time PCR (qRT-PCR) of six selected CsVQ genes (CsVQ6, CsVQ9, CsVQ12, CsVQ13, CsVQ16 and CsVQ20) showed that all of them were regulated by cold, drought and salt stresses. In addition, RNA-seq analysis revealed changes in the expression of some CsVQ genes (such as CsVQ2, CsVQ9, CsVQ12, CsVQ13, CsVQ20 and CsVQ22) under downy mildew (DM) and powdery mildew (PM) treatments. Our findings provide a theoretical basis for studying the potential roles of CsVQ genes in the development and stress response of cucumber in the future.

Characterization of Germin-like Proteins (GLPs) and Their Expression in Response to Abiotic and Biotic Stresses in Cucumber

Germins and germin-like proteins (GLPs) are glycoproteins closely associated with plant development and stress response in the plant kingdom. Here, we carried out genome-wide identification and expression analysis of the GLP gene family in cucumber to study their possible functions. A total of 38 GLP genes were identified in cucumber, which could be mapped to six out of the seven cucumber chromosomes. A phylogenetic analysis of the GLP members from cucumber, Arabidopsis and rice showed that these GLPs could be divided into six groups, and cucumber GLPs in the same group had highly similar conserved motif distribution and gene structure. Gene duplication analysis revealed that six cucumber GLP genes were located in the segmental duplication regions of cucumber chromosomes, while 14 genes were associated with tandem duplications. Tissue expression profiles of cucumber GLP genes showed that many genes were preferentially expressed in specific tissues. In addition, some cucumber GLP genes were differentially expressed under salt, drought and ABA treatments, as well as under DM inoculation. Our results provide important information for the functional identification of GLP genes in the growth, development and stress response of cucumber.

Using de novo assembly to identify structural variation of eight complex immune system gene regions.

Driven by the necessity to survive environmental pathogens, the human immune system has evolved exceptional diversity and plasticity, to which several factors contribute including inheritable structural polymorphism of the underlying genes. Characterizing this variation is challenging due to the complexity of these loci, which contain extensive regions of paralogy, segmental duplication and high copy-number repeats, but recent progress in long-read sequencing and optical mapping techniques suggests this problem may now be tractable. Here we assess this by using long-read sequencing platforms from PacBio and Oxford Nanopore, supplemented with short-read sequencing and Bionano optical mapping, to sequence DNA extracted from CD14+ monocytes and peripheral blood mononuclear cells from a single European individual identified as HV31. We use this data to build a de novo assembly of eight genomic regions encoding four key components of the immune system, namely the human leukocyte antigen, immunoglobulins, T cell receptors, and killer-cell immunoglobulin-like receptors. Validation of our assembly using k-mer based and alignment approaches suggests that it has high accuracy, with estimated base-level error rates below 1 in 10 kb, although we identify a small number of remaining structural errors. We use the assembly to identify heterozygous and homozygous structural variation in comparison to GRCh38. Despite analyzing only a single individual, we find multiple large structural variants affecting core genes at all three immunoglobulin regions and at two of the three T cell receptor regions. Several of these variants are not accurately callable using current algorithms, implying that further methodological improvements are needed. Our results demonstrate that assessing haplotype variation in these regions is possible given sufficiently accurate long-read and associated data. Continued reductions in the cost of these technologies will enable application of these methods to larger samples and provide a broader catalogue of germline structural variation at these loci, an important step toward making these regions accessible to large-scale genetic association studies.

General practice after COVID-19: an introduction to our special series.

<h3>Abstract</h3> Intragenic segmental duplication regions are potential hotspots for recurrent copy number variation and possible pathogenic aberrations. Two large sarcomeric genes, nebulin and titin, both contain such segmental duplication regions. Using our custom Comparative Genomic Hybridization array, we have previously shown that a gain or loss of more than one copy of the repeated block of the nebulin triplicate region constitutes a recessive pathogenic mutation. Using targeted array-CGH, similar copy number variants can be detected in the segmental duplication region of titin. Due to the limitations of the array-CGH methodology and the repetitiveness of the region, the exact copy numbers of the blocks could not be determined. Therefore, we developed complementary custom Droplet Digital PCR assays for the titin segmental duplication region to confirm true variation. Our combined methods show that the titin segmental duplication region is subject to recurrent copy number variants, which is surprisingly common. Gains and losses were detected in samples from healthy individuals as well as in samples from patients with different muscle disorders. The copy number variation observed in our cohort is likely benign, but pathogenic copy number variants in the segmental duplication region of titin cannot be excluded. Further investigations are needed, however, this region should no longer be neglected in genetic analyses.

Focused Revision: ACMG practice resource: Genetic evaluation of short stature

Focused Revision: ACMG practice resource: Genetic evaluation of short stature

Generation of induced pluripotent stem cell GZHMCi005-A from amniotic fluid-derived cells with duplication of chromosome 8p.

Cases in which the duplication of chromosome 8p (dup 8p) is observed are characterized by facial dysmorphism, agenesis/hypoplasia of the corpus callosum, heart defects and severe mental retardation. The frequency of dup 8p cases is higher compared to other chromosomes because of the Non-allelic homologous recombination (NAHR) between two segmental duplication regions (SDs) containing olfactory receptor gene clusters, REPD (repeat-distal) and REPP (repeat-proximal), located in chromosome 8p23.1. Here we generated a human iPSC line from a patient's amniotic fluid cells with a 18Mb duplication in 8p23.3p22, which will serve as useful tools for studying dup 8p syndrome.

Segmental Duplication as Potential Biomarkers for Non-Invasive Prenatal Testing of Aneuploidies

Background: Segmental duplication (SD) regions are distinct targets for aneuploidy detection owing to the virtual elimination of amplification bias. The difficulty of searching SD sequences for assay design has hampered their applications. Methods:We developed a computational program, ChAPDes, which integrates SD searching, refinement, and design of specific PCR primer/probe sets in a pipeline to remove most of the manual work. The generated primer/probe sets were first tested in a multiplex multicolor melting curve analysis for the detection of five common aneuploidies. The primer/probe sets were then tested in a digital PCR assay for the detection of trisomy 21. Finally, a digital PCR protocol was established to quantify maternal plasma DNA sequences for the non-invasive prenatal detection of fetal trisomy 21. Findings: ChAPDes could output 21,772 candidate primer/probe sets for trisomy 13, 18, 21 and sex chromosome aneuploidies within 2 working days. Clinical evaluation of the multiplex multicolour melting curve analysis involving 463 fetal genomic DNA samples revealed a sensitivity of 100% and specificity of 99.64% in comparison with the reference methods. Using the established digital PCR protocol, we correctly identified two trisomy 21 fetuses and thirteen euploid foetuses from the maternal plasma samples. Interpretation: The combination of ChAPDes with digital PCR detection could facilitate the use of SD as potential biomarkers for the non-invasive prenatal testing of fetal chromosomal aneuploidies. Funding Information: This work was supported by the National Key Research & Developmental Program of China (2018YFC1003601), the National Natural Science Foundation of China (8167210 to Q.H.), the Science and Technology Project of Xiamen City (3502Z20183013), and the Guangxi Key Laboratory Thalassemia Research Project (20-065-20). Declaration of Interests: The authors declare no competing interests. Ethics Approval Statement: The collection of specimens and associated clinical data used in this study were approved by the Affiliated Suzhou Hospital of Nanjing Medical University and Suzhou Municipal Hospital. All patients who participated in this study provided written informed consent prior to enrollment and specimen collection. This study was performed in accordance with the principles of the Helsinki Declaration.

Analysis of Subtelomeric REXTAL Assemblies Using QUAST.

Genomic regions of high segmental duplication content and/or structural variation have led to gaps and misassemblies in the human reference sequence, and are refractory to assembly from whole-genome short-read datasets. Human subtelomere regions are highly enriched in both segmental duplication content and structural variations, and as a consequence are both impossible to assemble accurately and highly variable from individual to individual. Recently, we developed a pipeline for improved region-specific assembly called Regional Extension of Assemblies Using Linked-Reads (REXTAL). In this study, we evaluate REXTAL and genome-wide assembly (Supernova) approaches on 10X Genomics linked-reads data sets partitioned and barcoded using the Gel Bead in Emulsion (GEM) microfluidic method. Our results describe the accuracy and relative performance of these two approaches using the reference-based assessment module of QUAST. We show that REXTAL dramatically outperforms the Supernova whole genome assembler in subtelomeric segmental duplication regions, and results in highly accurate assemblies. Nearly all of the REXTAL "misassemblies" identified using default QUAST parameters simply pinpoint locations of tandem repeat arrays in the reference sequence where the repeat array length differs from that in the cognate REXTAL assembly by 1000 bp.

Genome-wide identification of the peptide transporter family in rice and analysis of the PTR expression modulation in two near-isogenic lines with different nitrogen use efficiency

BackgroundNitrogen (N) is a major nutrient element for crop growth. In plants, the members of the peptide transporter (PTR) gene family may involve in nitrate uptake and transport. Here, we identified PTR gene family in rice and analyzed their expression profile in near-isogenic lines.ResultsWe identified 96, 85 and 78 PTR genes in Nipponbare, R498 and Oryza glaberrima, and the phylogenetic trees were similar in Asian cultivated rice and African cultivated rice. The number of PTR genes was higher in peanut (125) and soybean (127). The 521 PTR genes in rice, maize, sorghum, peanut, soybean and Arabidopsis could be classified into 4 groups, and their distribution was different between monocots and dicots. In Nipponbare genome, the 25 PTR genes were distributed in 5 segmental duplication regions on chromosome 1, 2, 3, 4, 5, 7, 8, 9, and 10. The PTR genes in rice have 0–11 introns and 1–12 exons, and 16 of them have the NPF (NRT1/PTR family) domain. The results of RNA-seq showed that the number of differentially expressed genes (DEGs) between NIL15 and NIL19 at three stages were 928, 1467, and 1586, respectively. Under low N conditions, the number of differentially expressed PTR genes increased significantly. The RNA-seq data was analyzed using WGCNA to predict the potential interaction between genes. We classified the genes with similar expression pattern into one module, and obtained 25 target modules. Among these modules, three modules may be involved in rice N uptake and utilization, especially the brown module, in which hub genes were annotated as protein kinase that may regulate rice N metabolism.ConclusionsIn this study, we comprehensively analyzed the PTR gene family in rice. 96 PTR genes were identified in Nippobare genome and 25 of them were located on five large segmental duplication regions. The Ka/Ks ratio indicated that many PTR genes had undergone positive selection. The RNA-seq results showed that many PTR genes were involved in rice nitrogen use efficiency (NUE), and protein kinases might play an important role in this process. These results provide a fundamental basis to improve the rice NUE via molecular breeding.

Improved assembly and variant detection of a haploid human genome using single-molecule, high-fidelity long reads.

The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6kbp, CLR 191.5kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective standalone technology for de novo assembly of human genomes.