Advantages and diagnostic effectiveness of the two most widely used resequencing approaches, whole exome (WES) and whole genome (WGS) sequencing, are often debated. WES dominated large-scale resequencing projects because of lower cost and easier data storage and processing. Rapid development of 3^rd generation sequencing methods and novel exome sequencing kits predicate the need for a robust statistical framework allowing informative and easy performance comparison of the emerging methods. In our study we developed a set of statistical tools to systematically assess coverage of coding regions provided by several modern WES platforms, as well as PCR-free WGS. We identified a substantial problem in most previously published comparisons which did not account for mappability limitations of short reads. Using regression analysis and simple machine learning, as well as several novel metrics of coverage evenness, we analyzed the contribution from the major determinants of CDS coverage. Contrary to a common view, most of the observed bias in modern WES stems from mappability limitations of short reads and exome probe design rather than sequence composition. We also identified the ~ 500 kb region of human exome that could not be effectively characterized using short read technology and should receive special attention during variant analysis. Using our novel metrics of sequencing coverage, we identified main determinants of WES and WGS performance. Overall, our study points out avenues for improvement of enrichment-based methods and development of novel approaches that would maximize variant discovery at optimal cost.

Cardiomyopathies represent a rare group of disorders often of genetic origin. While approximately 50% of genetic causes are known for other types of cardiomyopathies, the genetic spectrum of restrictive cardiomyopathy (RCM) is largely unknown. The aim of the present study was to identify the genetic background of idiopathic RCM and to compile the obtained genetic variants to the novel signalling pathways using in silico protein network analysis. We used Illumina MiSeq setup to screen for 108 cardiomyopathy and arrhythmia-associated genes in 24 patients with idiopathic RCM. Pathogenicity of genetic variants was classified according to American College of Medical Genetics and Genomics classification. Pathogenic and likely-pathogenic variants were detected in 13 of 24 patients resulting in an overall genotype-positive rate of 54%. Half of the genotype-positive patients carried a combination of pathogenic, likely-pathogenic variants and variants of unknown significance. The most frequent combination included mutations in sarcomeric and cytoskeletal genes (38%). A bioinformatics approach underlined the mechanotransducing protein networks important for RCM pathogenesis. Multiple gene mutations were detected in half of the RCM cases, with a combination of sarcomeric and cytoskeletal gene mutations being the most common. Mutations of genes encoding sarcomeric, cytoskeletal, and Z-line-associated proteins appear to have a predominant role in the development of RCM.

Malformations of the left ventricular outflow tract are one of the most common forms of congenital heart disorders. Recently, it has been shown that mutations in the NOTCH1 gene can lead to bicuspid aortic valve, aortic aneurysm, and hypoplastic left heart syndrome. The aim of our study was to estimate the frequency of NOTCH1 gene mutations/substitutions in patients with aortic coarctation, isolated or combined with bicuspid aortic valve. The study included 51 children with coarctation. Detailed family history was obtained for every study subject, and echocardiographic data were obtained for the relatives when available. We applied a strategy of targeted mutation screening for 10 out of 34 exons of the NOTCH1 gene by direct sequencing. Control DNA was obtained from 200 healthy donors. In more than half of the cases, coarctation was combined with bicuspid aortic valve, and in approximately half of the cases, it was combined with hypoplasia of the aortic arch or descending aorta. Familial history of congenital heart disease was observed in 34.3% of the cases. In total, 29 variants of the NOTCH1 gene were identified in the patient group and in the control subjects. Four of those variants led to amino acid exchange, of which only one, R1279H, was identified in both the patient group and in the controls. This variant was significantly overrepresented in the patients with aortic coarctation compared with those in the control group (P < .05). We conclude that the R1279H substitution in the NOTCH1 gene is significantly overrepresented in patients with aortic coarctation and, therefore, may represent a disease-susceptibility allele.

Ribosomal RNA (rRNA) genes, whose activity results in nucleolus formation, constitute an extremely important part of genome. Despite the extensive exploration into avian genomes, no complete description of avian rRNA gene primary structure has been offered so far. We publish a complete chicken rRNA gene cluster sequence here, including 5’ETS (1836 bp), 18S rRNA gene (1823 bp), ITS1 (2530 bp), 5.8S rRNA gene (157 bp), ITS2 (733 bp), 28S rRNA gene (4441 bp) and 3’ETS (343 bp). The rRNA gene cluster sequence of 11863 bp was assembled from raw reads and deposited to GenBank under KT445934 accession number. The assembly was validated through in situ fluorescent hybridization analysis on chicken metaphase chromosomes using computed and synthesized specific probes, as well as through the reference assembly against de novo assembled rRNA gene cluster sequence using sequenced fragments of BAC-clone containing chicken NOR (nucleolus organizer region). The results have confirmed the chicken rRNA gene cluster validity.

Background: Exosomes are involved in intercellular communication and can transfer regulatory molecules between cells. Consequently, they can participate in host immune response regulation. For the influenza A virus (IAV), there is very limited information on changes in exosome composition during cell infection shedding light on the potential role of these extracellular membrane vesicles. Thus, the aim of our work was to study changes in exosomal composition following IAV infection of cells, as well as to evaluate their effect on uninfected cells. Methods: To characterize changes in the composition of cellular miRNAs and mRNAs of exosomes during IAV infection of A549 cells, NGS was used, as well as PCR to identify viral genes. Naïve A549 cells were stimulated with infected-cell-secreted exosomes for studying their activity. Changes in the expression of genes associated with the cell’s immune response were shown using PCR. The effect of exosomes on IAV replication was shown in MDCK cells using In-Cell ELISA and PCR of the supernatants. Results: A change in the miRNA composition (miR-21-3p, miR-26a-5p, miR-23a-5p, miR-548c-5p) and mRNA composition (RPL13A, MKNK2, TRIB3) of exosomes under the influence of the IAV was shown. Many RNAs were involved in the regulation of the immune response of the cell, mainly by suppressing it. After exosome stimulation of naïve cells, a significant decrease in the expression of genes involved in the immune response was shown (RIG1, IFIT1, MDA5, COX2, NFκB, AnxA1, PKR, IL6, IL18). When infecting MDCK cells, a significant decrease in nucleoprotein levels was observed in the presence of exosomes secreted by mock-infected cells. Viral levels in supernatants also decreased. Conclusions: Exosomes secreted by IAV-infected cells could reduce the immune response of neighboring intact cells, leading to more effective IAV replication. This may be associated both with regulatory functions of cellular miRNAs and mRNAs carried by exosomes, or with the presence of viral mRNAs encoding proteins with an immunosuppressive function.