Duplicated Sequences Research Articles

Sequali: Efficient and Comprehensive Quality Control of Short- and Long-Read Sequencing Data

Abstract Motivation Quality control of sequencing data is the first step in many sequencing workflows. Short- and long-read sequencing technologies have many commonalities with regards to quality control. Several quality control programs exist, however none possess a feature set that is adequate for both technologies. Quality Control programs aimed at Oxford Nanopore Technologies sequencing lack vital features such as adapter searching, overrepresented sequence analysis and duplication analysis. Results Sequali was developed to provide sequencing quality control for both short- and long-read sequencing technologies. It features adapter search, overrepresented sequence analysis and duplication analysis and supports FASTQ and uBAM inputs. It is significantly faster than comparable sequencing quality control programs for both short- and long-read sequencing technologies. Availability and Implementation Sequali is an open source Python application using C extensions and is freely available under the AGPL-3.0 license at https://github.com/rhpvorderman/sequali. The source code for each release is archived at zenodo: https://zenodo.org/doi/10.5281/zenodo.10822485.

Order in Chaos: Lesser-Conserved and Repeat Structures in Dehydrins.

Dehydrins (Dhns) are a group of intrinsically disordered land plant proteins that are closely associated with tolerance of dehydrative stress. Dhns are recognized and classified by the presence and sequence of five different conserved segments, varying in length from 8 to 15 residues, separated by highly variable disordered regions. In addition to one or more copies of the diagnostic, fifteen-residue K segment, most Dhns can be classified into one of three major groups based on the mutually exclusive presence of three other conserved segments (H, Y, or F), with all three groups typically incorporating multi-serine S segments. Many Dhns also include repeat structures. From an input library of 8675 non-redundant candidate sequences, a specialized R script identified and classified 2658 complete and 236 partial Dhn sequences in all major green plant (Viridiplantae) lineages, including a few green algal genera. An examination of the connecting segments bridging the conserved segments identified additional conserved patterns, suggesting that multi-Y, S-K, and K-S domains may act as functional units. Dhn Decoder identified 857 Dhns with repeat structures, ranging from 3 short, simple repeats to elaborate variations with up to 45 repeats or repeats of up to 85 residues comprising 1 or more of the conserved segments, suggesting that internal sequence duplication is an important mode of evolution in Dhns.

Structural polymorphism and diversity of human segmental duplications.

Segmental duplications (SDs) contribute significantly to human disease, evolution and diversity but have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies (from 85 samples representing 38 Africans and 47 non-Africans) in which the majority of autosomal SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms and sex chromosomes, we identify 173.2 Mb of duplicated sequence (47.4 Mb not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable, with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy numbers than non-African samples. Comparison to a resource of 563 million full-length isoform sequencing reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.

Tailed molecular beacon probes: an approach for the detection of structured DNA and RNA analytes.

Molecular beacon (MB) probes have been extensively used for nucleic acid analysis. However, MB probes fail to hybridize with folded DNA or RNA. Here, we demonstrate that MB probes equipped with extra sequences complementary to the analyte, named 'tail', can increase the signal-to-background ratio by ∼40-fold and hybridization rates by ∼800-fold compared to conventional MB probes. Tailed MB probes can be used as mismatched-tolerant alternatives to traditional hairpin probes for fast assays.

Insights into Adaption and Growth Evolution: Genome-Wide Copy Number Variation Analysis in Chinese Hainan Yellow Cattle Using Whole-Genome Re-Sequencing Data.

Copy number variation (CNV) serves as a crucial source of genomic variation and significantly aids in the mining of genomic information in cattle. This study aims to analyze re-sequencing data from Chinese Hainan yellow cattle, to uncover breed CNV information, and to elucidate the resources of population genetic variation. We conducted whole-genome sequencing on 30 Chinese Hainan yellow cattle, thus generating 814.50 Gb of raw data. CNVs were called using CNVnator software, and subsequent filtering with Plink and HandyCNV yielded 197,434 high-quality CNVs and 5852 CNV regions (CNVRs). Notably, the proportion of deleted sequences (81.98%) exceeded that of duplicated sequences (18.02%), with the lengths of CNVs predominantly ranging between 20 and 500 Kb This distribution demonstrated a decrease in CNVR count with increasing fragment length. Furthermore, an analysis of the population genetic structure using CNVR databases from Chinese, Indian, and European commercial cattle breeds revealed differences between Chinese Bos indicus and Indian Bos indicus. Significant differences were also observed between Hainan yellow cattle and European commercial breeds. We conducted gene annotation for both Hainan yellow cattle and European commercial cattle, as well as for Chinese Bos indicus and Indian Bos indicus, identifying 206 genes that are expressed in both Chinese and Indian Bos indicus. These findings may provide valuable references for future research on Bos indicus. Additionally, selection signatures analysis based on Hainan yellow cattle and three European commercial cattle breeds identified putative pathways related to heat tolerance, disease resistance, fat metabolism, environmental adaptation, candidate genes associated with reproduction and the development of sperm and oocytes (CABS1, DLD, FSHR, HSD17B2, KDM2A), environmental adaptation (CNGB3, FAM161A, DIAPH3, EYA4, AAK1, ERBB4, ERC2), oxidative stress anti-inflammatory response (COMMD1, OXR1), disease resistance (CNTN5, HRH4, NAALADL2), and meat quality (EHHADH, RHOD, GFPT1, SULT1B1). This study provides a comprehensive exploration of CNVs at the molecular level in Chinese Hainan yellow cattle, offering theoretical support for future breeding and selection programs aimed at enhancing qualities of this breed.

The origin and maintenance of supergenes contributing to ecological adaptation in Atlantic herring

Chromosomal inversions are associated with local adaptation in many species. However, questions regarding how they are formed, maintained and impact various other evolutionary processes remain elusive. Here, using a large genomic dataset of long-read and short-read sequencing, we ask these questions in one of the most abundant vertebrates on Earth, the Atlantic herring. This species has four megabase-sized inversions associated with ecological adaptation that correlate with water temperature. The S and N inversion alleles at these four loci dominate in the southern and northern parts, respectively, of the species distribution in the North Atlantic Ocean. By determining breakpoint coordinates of the four inversions and the structural variations surrounding them, we hypothesize that these inversions are formed by ectopic recombination between duplicated sequences immediately outside of the inversions. We show that these are old inversions (>1 MY), albeit formed after the split between the Atlantic herring and its sister species, the Pacific herring. There is evidence for extensive gene flux between inversion alleles at all four loci. The large Ne of herring combined with the common occurrence of opposite homozygotes across the species distribution has allowed effective purifying selection to prevent the accumulation of genetic load and repeats within the inversions.

Comprehensive characterization of gibberellin oxidase gene family in Brassica napus reveals BnGA2ox15 involved in hormone signaling and response to drought stress

Brassica napus is a well-known allopolyploid oil crop with high commercial potential. Gibberellin oxidase (GAox) is an essential enzyme that activates gibberellins, which regulate plant growth, and development, and have a significant impact on plant responses to abiotic stress. However, the comprehensive understanding of GAox genes and their evolution in Brassica plants remains elusive. Using advanced bioinformatics tools, this study identified 125 candidate GAox genes from the whole genomes of three key Brassica species. This study also investigated sequence characteristics, conserved motifs, exon/intron structures, cis-acting elements, syntenic analysis, duplication events and expression patterns. Subcellular localization analysis showed that the BnGA2ox14 and BnGA2ox15 proteins are located in the nucleus, whereas BnGA2ox26 is specifically localized to the chloroplast. Yeast one-hybrid and dual-luciferase assays demonstrated that MYELOCYTOMATOSIS 4 (BnMYC4) and ABA-INDUCIBLE BHLH-TYPE TRANSCRIPTION FACTOR (BnAIB) bind to the BnGA2ox15 promoter and activate its transcription. Molecular docking analysis further elucidated their interaction structures and identified potential binding sites. Roots transformations show that overexpression of BnGA2ox15 increased sensitivity to PEG-6000 treatment in rapeseed. In brief, this study reveals that BnGA2ox15 is a downstream target in JA and ABA signaling pathways, functioning as a negative regulator in response to drought stress.

An innovative CRISPR/Cas9 mouse model of human isolated microtia indicates the potential contribution of CNVs near HMX1 gene

BackgroundMicrotia is a prevalent congenital malformation, the precise etiology and pathogenesis of which remain elusive. Mutations in the non-coding region of the HMX1 gene have been implicated in isolated cases of microtia, emerging as a significant focus of contemporary research. Several pathogenic copy number variations (CNVs) proximal to the HMX1 gene have been documented in wild animal populations, whereas only a single large segmental duplication in this region has been identified in humans. However, the absence of a gene-edited animal model has impeded the investigation of the unclear gene function associated with HMX1 mutations in human isolated microtia. In this study, we sought to precisely identify the pathogenic mutation by analyzing three pedigrees alongside population controls. Subsequently, our objective was to develop a CRISPR/Cas9 gene-edited mouse model to elucidate the functional implications of the identified mutation. MethodsGenomic DNA was collected from 32 affected individuals across three pedigrees, as well as from 2000 control subjects. Comprehensive genomic analyses, including genome-wide linkage analysis, targeted capture, second-generation sequencing, and copy number analysis, were conducted to identify potential mutations associated with congenital auricle malformation. CRISPR/Cas9 gene-edited murine models were generated in response to the identified mutation. The auricular phenotypes of these gene-edited mice were systematically monitored. Small-animal Micro-CT scanning was employed to identify potential craniofacial or skeletal abnormalities. Furthermore, the expression of the HMX1 gene in the PA2 region of mouse embryos was quantified using RT-qPCR. ResultsA co-segregated 600 base pair duplication located on chromosome 4 (chr4:8701900-8702500, hg19) was identified in affected individuals across three pedigrees, but was absent in healthy controls. Two types of CRISPR/Cas9 gene-edited mice were subsequently generated. The knock-in (KI) mouse model was engineered by inserting one copy of the duplicated sequence directly adjacent to the mutated site, whereas the knockout (KO) mouse model was created by excising the mutation sequence. The phenotypes of different group of CRISPR/Cas9 gene-edited mice demonstrated distinct auricular deformities. Furthermore, an increase in the copy number of the mutated sequence was associated with elevated expression levels of HMX1 in the gene-edited mouse model. ConclusionsIn this study, we further narrowed down and identified a 600 base pair copy number variation (CNV) located at chr4:8701900-8702500 (hg19), which is implicated in human bilateral, isolated microtia. Utilizing CRISPR/Cas9 technology, we developed novel mouse models harboring the identified mutation. These models serve as a robust platform for the comprehensive investigation of the underlying mechanisms of the disease.

Genome-wide investigation of caffeic acid O-methyltransferases and expression analysis under different abiotic stresses and anthracnose infection in pepper (Capsicum annuum L.)

ABSTRACT The caffeic acid O-methyl transferases (COMTs) facilitate the biosynthesis of lignin and play important roles in many processes, including growth, development and response to stresses. Pepper (Capsicum annuum L.) is a major horticultural crop with significant nutritional and medicinal values. Nevertheless, a comprehensive assessment of the COMT genes in pepper and their involvement in drought, salt tolerance and anthracnose resistance remains unestablished. In the present study, 80 pepper COMT genes were identified and classified into four major groups based on domain structure and phylogenetics. Additionally, the sequence analysis, gene structures, and duplication events were analysed. The genes were unevenly distributed across 10 chromosomes, and 10 pairs of tandem and 13 pairs of segmental duplicated genes were detected in the COMT gene family. Examination of cis-acting regulatory elements detected the presence of phytohormone-specific, stress-receptive and light responsive elements in the pepper COMT promoters. Moreover, qRT-PCR analysis revealed significant upregulated expression of CanCOMT48 to drought stress, CanCOMT14 to salinity stress and CanCOMT18 to anthracnose infection by Colletotrichum truncatum, suggesting their important functions in different stresses. Our results provide significant insights into the expression and probable function of selective COMTs which could be used as potential candidate genes for pepper tolerance breeding.

Endogenous viral elements are targeted by RNA silencing pathways in banana.

Endogenous banana streak virus (eBSV) integrants derived from three distinct species, present in Musa balbisiana (B) but not Musa acuminata (A) banana genomes are able to reconstitute functional episomal viruses causing banana streak disease in interspecific triploid AAB banana hybrids but not in the diploid (BB) parent line, which harbours identical eBSV loci. Here, we investigated the regulation of these eBSV. In-depth characterization of siRNAs, transcripts and methylation derived from eBSV using Illumina and bisulfite sequencing were carried out on eBSV-free Musa acuminata AAA plants and BB or AAB banana plants with eBSV. eBSV loci produce low-abundance transcripts covering most of the viral sequence and generate predominantly 24-nt siRNAs. siRNA accumulation is restricted to duplicated and inverted viral sequences present in eBSV. Both siRNA-accumulating and nonaccumulating sequences of eBSV in BB plants are heavily methylated in all three CG, CHG and CHH contexts. Our data suggest that eBSVs are controlled at the epigenetic level in BB diploids. This regulation not only prevents their awakening and systemic infection of the plant but is also probably involved in the inherent resistance of the BB plants to mealybug-transmitted viral infection. These findings are thus of relevance to other plant resources hosting integrated viruses.

Recurrent Independent Pseudogenization Events of the Sperm Fertilization Gene ZP3r in Apes and Monkeys.

Many reproductive proteins show signatures of rapid evolution through sequence divergence and duplication. These features of reproductive genes may complicate the detection of orthologs across taxa, making it difficult to connect studies in model systems to human biology. In mice, ZP3r/sp56 is a binding partner to the egg coat protein ZP3 and may mediate induction of the acrosome reaction, a crucial step in fertilization. In rodents, ZP3r, as a member of the Regulators of Complement Activation cluster, is surrounded by paralogs, some of which have been shown to be evolving under positive selection. Although primate egg coats also contain ZP3, sequence divergence paired with paralogous relationships with neighboring genes has complicated the accurate identification of the human ZP3r ortholog. Here, we phylogenetically and syntenically resolve that the human ortholog of ZP3r is the pseudogene C4BPAP1. We investigate the evolution of this gene within primates. We observe independent pseudogenization events of ZP3r in all Apes with the exception of Orangutans, and independent pseudogenization events in many monkey species. ZP3r in both primates that retain ZP3r and in rodents contains positively selected sites. We hypothesize that redundant mechanisms mediate ZP3 recognition in mammals and ZP3r's relative importance to ZP recognition varies across species.

An intronic copy number variation in Syntaxin 17 determines speed of greying and melanoma incidence in Grey horses

The Greying with age phenotype in horses involves loss of hair pigmentation whereas skin pigmentation is not reduced, and a predisposition to melanoma. The causal mutation was initially reported as a duplication of a 4.6 kb intronic sequence in Syntaxin 17. The speed of greying varies considerably among Grey horses. Here we demonstrate the presence of two different Grey alleles, G2 carrying two tandem copies of the duplicated sequence and G3 carrying three. The latter is by far the most common allele, probably due to strong selection for the striking white phenotype. Our results reveal a remarkable dosage effect where the G3 allele is associated with fast greying and high incidence of melanoma whereas G2 is associated with slow greying and low incidence of melanoma. The copy number expansion transforms a weak enhancer to a strong melanocyte-specific enhancer that underlies hair greying (G2 and G3) and a drastically elevated risk of melanoma (G3 only). Our direct pedigree-based observation of the origin of a G2 allele from a G3 allele by copy number contraction demonstrates the dynamic evolution of this locus and provides the ultimate evidence for causality of the copy number variation of the 4.6 kb intronic sequence.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References.

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

The Internal Extra Sequence Regions in Satellite RNA TA-Tb Are Important for Suppressing RNA Accumulations of Cucumber Mosaic Virus to Attenuate the Virulence of the Helper Virus

Cucumber mosaic virus (CMV) infection is often associated with satellite RNA (satRNA), which can sometimes interfere with the replication and symptom expression of CMV. However, the mechanism underlying symptom attenuation has remained unclear. We previously discovered a larger type (than the usual type) of satellite RNA (satRNA TA-Tb) of CMV that reduced the symptom severity of CMV. Herein, we show that satRNA TA-Tb is associated with a reduction in CMV RNA accumulation, and particularly, a strong reduction of RNA4 accumulation at later stages of infection. Deletion analysis showed that the deletion of ten nucleotides of 5′ and 3′ termini, but not the internal sequence regions proximal to the 5′- and 3′-terminal regions, abolished satRNA TA-Tb replication. The alignment of satRNA TA-Tb with usual satRNA isolates showed four internal extra sequence regions (exR1–4) in satRNA TA-Tb. A satRNA TA-Tb mutant with deletion in the exR1 region retained the ability to attenuate CMV symptoms, whereas deletion of the exR2–4 regions abolished the attenuating effect of satRNA TA-Tb, but did not affect its replication. Overall, these results suggest that some short, internal extra sequence regions are dispensable for satRNA TA-Tb replication, but important for symptom attenuation function, supporting the possibility that the RNA structure of satRNA TA-Tb is important for its function in symptom attenuation.

Duplications and Retrogenes Are Numerous and Widespread in Modern Canine Genomic Assemblies.

Recent years have seen a dramatic increase in the number of canine genome assemblies available. Duplications are an important source of evolutionary novelty and are also prone to misassembly. We explored the duplication content of nine canine genome assemblies using both genome self-alignment and read-depth approaches. We find that 8.58% of the genome is duplicated in the canFam4 assembly, derived from the German Shepherd Dog Mischka, including 90.15% of unplaced contigs. Highlighting the continued difficulty in properly assembling duplications, less than half of read-depth and assembly alignment duplications overlap, but the mCanLor1.2 Greenland wolf assembly shows greater concordance. Further study shows the presence of multiple segments that have alignments to four or more duplicate copies. These high-recurrence duplications correspond to gene retrocopies. We identified 3,892 candidate retrocopies from 1,316 parental genes in the canFam4 assembly and find that ∼8.82% of duplicated base pairs involve a retrocopy, confirming this mechanism as a major driver of gene duplication in canines. Similar patterns are found across eight other recent canine genome assemblies, with metrics supporting a greater quality of the PacBio HiFi mCanLor1.2 assembly. Comparison between the wolf and other canine assemblies found that 92% of retrocopy insertions are shared between assemblies. By calculating the number of generations since genome divergence, we estimate that new retrocopy insertions appear, on average, in 1 out of 3,514 births. Our analyses illustrate the impact of retrogene formation on canine genomes and highlight the variable representation of duplicated sequences among recently completed canine assemblies.

Search for genes gained by horizontal gene transfer in an entomopoxvirus, with special reference to the analysis of the transfer of an ABC transporter gene

Although it is generally believed that large DNA viruses capture genes by horizontal gene transfer (HGT), the detailed manner of such transfer has not been fully elucidated. Here, we searched for genes in the coleopteran entomopoxvirus (EV) Anomala cuprea entomopoxvirus (ACEV) that might have been gained by ACEV by HGT. We classified the potential source organisms for HGT into three categories: the host A. cuprea; other organisms, including viruses unrelated to EVs; and organisms with uncertain host attribution. Of the open reading frames (ORFs) of the ACEV genome, 2.1 % were suggested to have been gained from the host by ACEV or its recent ancestor via HGT; 8.7 % were possibly from organisms other than the host, and 3.7 % were possibly from the third category of organisms via HGT. The analysis showed that ACEV contains some interesting ORFs obtained by HGT, including a large ATP-binding cassette protein (ABC transporter) ORF and a tenascin ORF (IDs ACV025 and ACV123, respectively). We then performed a detailed analysis of the HGT of the ACEV large ABC transporter ORF—the largest of the ACEV ORFs. mRNA sequences obtained by RNA-seq from fat bodies—sites of ACEV replication—and midgut tissues—sites of initial infection—of the virus's host A. cuprea larvae were subjected to BLAST analysis. One type of ABC transporter ORF from the fat bodies and two types from the midgut tissues, one of which was identical to that in the fat bodies, had the greatest identity to the ABC transporter ORF of ACEV. The two types from the host had high levels of identity to each other (approximately 95 % nucleotide sequence identity), strongly suggesting that the host ABC transporter group consisting of the two types was the origin of ACV025. We then determined the sequence (12,381 bp) containing a full-length gene of the A. cuprea ABC transporter. It turned out to be a transcription template for the abovementioned mRNA found in both tissues. In addition, we determined a large part (ca. 6.9 kb) of the template sequence for the mRNA found only in the midgut tissues. The results showed that the ACEV ABC transporter ORF is missing parts corresponding to introns of the host ABC transporter genes, indicating that the ORF was likely acquired by HGT in the form of mRNA. The presence of definite duplicated sequences adjacent to the ACEV ABC transporter genes—a sign of LINE-1 retrotransposon-mediated HGT—was not observed. An approximately 2-month ACV025 transcription experiment suggested that the transporter sequence is presumed to be continuously functional. The amino acid sequence of ACV025 suggests that its product might function in the regulation of phosphatide in the host-cell membranes.

Rice Big Grain1 improves grain yield in ectopically expressing rice and heterologously expressing tobacco plants.

Functional genomics through transgenesis has provided faster and more reliable methods for identifying, characterizing, and utilizing genes or quantitative trait loci linked to agronomic traits to target yield. The present study explored the role of Big Grain1 (BG1) gene of rice (Oryza sativa L.) in yield improvement of crop plants. We aimed to identify the genetic variation of OsBG1 in various indica rice cultivars by studying the allelic polymorphism of the gene, while also investigating the gene's potential to increase crop yield through the transgenic approach. Our study reports the presence of an extra 393bp sequence having two 6bp enhancer elements in the 3' regulatory sequence of OsBG1 in the large-grain cultivar IR64 but not in the small-grain cultivar Badshahbhog. A single copy of the OsBG1 gene in both the cultivars and a 4.1-fold higher expression of OsBG1 in IR64 than in Badshahbhog imply that the grain size is positively correlated with the level of OsBG1 expression in rice. The ectopic expression of OsBG1 under the endosperm-specific glutelin C promoter in Badshahbhog enhanced the flag leaf length, panicle weight, and panicle length by an average of 33.2%, 33.7%, and 30.5%, respectively. The length of anthers, spikelet fertility, and grain yield per plant increased in transgenic rice lines by an average of 27.5%, 8.3%, and 54.4%, respectively. Heterologous expression of OsBG1 under the constitutive 2xCaMV35S promoter improved the number of seed pods per plant and seed yield per plant in transgenic tobacco lines by an average of 2.2-fold and 2.6-fold, respectively. Improving crop yield is crucial to ensure food security and socio-economic stability, and identifying suitable genetic factor is the essential step towards this endeavor. Our findings suggest that the OsBG1 gene is a promising candidate for improving the grain yield of monocot and dicot plant systems by molecular breeding and genetic engineering.

Emergence of a Novel G4P[6] Porcine Rotavirus with Unique Sequence Duplication in NSP5 Gene in China.

Rotavirus is a major causative agent of diarrhoea in children, infants, and young animals around the world. The associated zoonotic risk necessitates the serious consideration of the complete genetic information of rotavirus. A segmented genome makes rotavirus prone to rearrangement and the formation of a new viral strain. Monitoring the molecular epidemiology of rotavirus is essential for its prevention and control. The quantitative RT-PCR targeting the NSP5 gene was used to detect rotavirus group A (RVA) in pig faecal samples, and two pairs of universal primers and protocols were used for amplifying the G and P genotype. The genotyping and phylogenetic analysis of 11 genes were performed by RT-PCR and a basic bioinformatics method. A unique G4P[6] rotavirus strain, designated S2CF (RVA/Pig-tc/CHN/S2CF/2023/G4P[6]), was identified in one faecal sample from a piglet with severe diarrhoea in Guangdong, China. Whole genome sequencing and analysis suggested that the 11 segments of the S2CF strain showed a unique Wa-like genotype constellation and a typical porcine RVA genomic configuration of G4-P[6]-I1-R1-C1-M1-A8-N1-T1-E1-H1. Notably, 4 of the 11 gene segments (VP4, VP6, VP2, and NSP5) clustered consistently with human-like RVAs, suggesting independent human-to-porcine interspecies transmission. Moreover, a unique 344-nt duplicated sequence was identified for the first time in the untranslated region of NSP5. This study further reveals the genetic diversity and potential inter-species transmission of porcine rotavirus.

Identification of novel driver risk genes in CNV loci associated with neurodevelopmental disorders

Copy number variants (CNVs) are genome-wide structural variations involving the duplication or deletion of large nucleotide sequences. While these types of variations can be commonly found in humans, large and rare CNVs are known to contribute to the development of various neurodevelopmental disorders (NDDs), including autism spectrum disorder (ASD). Nevertheless, given that these NDD-risk CNVs cover broad regions of the genome, it is particularly challenging to pinpoint the critical gene(s) responsible for the manifestation of the phenotype. In this study, we performed a meta-analysis of CNV data from 11,614 patients with NDDs and 4,031 controls from SFARI database to identify 41 NDD-risk CNV loci, including 24 novel regions. We also found evidence for dosage-sensitive genes within these regions being significantly enriched for known NDD-risk genes and pathways. In addition, a significant proportion of these genes was found to i) converge in protein-protein interaction networks; ii) be among most expressed genes in the brain across all developmental stages; and iii) be hit by deletions that are significantly over-transmitted to individuals with ASD within multiplex ASD families from the iHART cohort. Finally, we conducted a burden analysis using 4,281 NDD cases from Decipher and iHART cohorts, and 2,504 neurotypical controls from 1,000 Genomes and iHART, that resulted in the validation of the association of 162 dosage sensitive genes driving risk for NDDs, including 22 novel NDD-risk genes. Importantly, most NDD-risk CNV loci entail multiple NDD-risk genes in agreement with a polygenic model associated with the majority of NDD cases.

Structural polymorphism and diversity of human segmental duplications.

Segmental duplications (SDs) contribute significantly to human disease, evolution, and diversity yet have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies where the majority of SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms, we identify 173.2 Mbp of duplicated sequence (47.4 Mbp not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy number when compared to non-African samples. A comparison to a resource of 563 million full-length Iso-Seq reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.