Loss of phospholipase PLAAT3 causes a mixed lipodystrophic and neurological syndrome due to impaired PPARγ signaling

Phospholipase A/acyltransferase 3 (PLAAT3) is a phospholipid-modifying enzyme predominantly expressed in neural and white adipose tissue (WAT). It is a potential drug target for metabolic syndrome, as Plaat3 deficiency in mice protects against diet-induced obesity. We identified seven patients from four unrelated consanguineous families, with homozygous loss-of-function variants in PLAAT3, who presented with a lipodystrophy syndrome with loss of fat varying from partial to generalized and associated with metabolic complications, as well as variable neurological features including demyelinating neuropathy and intellectual disability. Multi-omics analysis of mouse Plaat3-/- and patient-derived WAT showed enrichment of arachidonic acid-containing membrane phospholipids and a strong decrease in the signaling of peroxisome proliferator-activated receptor gamma (PPARγ), the master regulator of adipocyte differentiation. Accordingly, CRISPR-Cas9-mediated PLAAT3 inactivation in human adipose stem cells induced insulin resistance, altered adipocyte differentiation with decreased lipid droplet formation and reduced the expression of adipogenic and mature adipocyte markers, including PPARγ. These findings establish PLAAT3 deficiency as a hereditary lipodystrophy syndrome with neurological manifestations, caused by a PPARγ-dependent defect in WAT differentiation and function.

Syntaxin 18 Defects in Human and Zebrafish Unravel Key Roles in Early Cartilage and Bone Development

SNARE proteins comprise a conserved protein family responsible for catalyzing membrane fusion during vesicle traffic. Syntaxin18 (STX18) is a poorly characterized endoplasmic reticulum (ER)-resident t-SNARE. Recently, together with TANGO1 and SLY1, its involvement was shown in ER to Golgi transport of collagen II during chondrogenesis. We report a fetus with a severe osteochondrodysplasia in whom we identified a homozygous substitution of the highly conserved p.Arg10 to Pro of STX18. CRISPR/Cas9-mediated Stx18 deficiency in zebrafish reveals a crucial role for Stx18 in cartilage and bone development. Furthermore, increased expression of multiple components of the Stx18 SNARE complex and of COPI and COPII proteins suggests that Stx18 deficiency impairs antero- and retrograde vesicular transport in the crispant stx18 zebrafish. Taken together, our studies highlight a new candidate gene for a recessive form of osteochondrodysplasia, thereby possibly broadening the SNAREopathy phenotypic spectrum and opening new doors toward future research avenues. © 2023 American Society for Bone and Mineral Research (ASBMR).

ENIGMA CHEK2gether Project: A Comprehensive Study Identifies Functionally Impaired CHEK2 Germline Missense Variants Associated with Increased Breast Cancer Risk

PURPOSE

Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT).

EXPERIMENTAL DESIGN

We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1-CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case-control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls.

RESULTS

A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)-like (N = 226). We then examined their association with breast cancer risk in the case-control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35-3.41), 1.57 (95% CI, 1.41-1.75), and 1.19 (95% CI, 1.08-1.31), respectively. The meta-analysis of population-specific datasets showed similar results.

CONCLUSIONS

We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

A novel duplication involving PRDM13 in a Turkish family supports its role in North Carolina macular dystrophy (NCMD/MCDR1)

Purpose

To clinically and molecularly investigate a new family with North Carolina macular dystrophy (NCMD) from Turkey, a previously unreported geographic origin for this phenotype.

Methods

Clinical ophthalmic examinations, including fundus imaging and spectral domain-optical coherence tomography (SD-OCT), were performed on eight members of a two-generation non-consanguineous family from southern Turkey. Whole genome sequencing (WGS) was performed on two affected subjects, followed by variant filtering and copy number variant (CNV) analysis. Junction PCR and Sanger sequencing were used to confirm and characterize the duplication involving PRDM13 at the nucleotide level. The underlying mechanism was assessed with in silico analyses.

Results

The proband presented with lifelong bilateral vision impairment and displayed large grade 3 coloboma-like central macular lesions. Five of her six children showed similar macular malformations, consistent with autosomal dominant NCMD. The severity grades in the six affected individuals from two generations are not evenly distributed. CNV analysis of WGS data of the two affected family members, followed by junction PCR and Sanger sequencing, revealed a novel 56.2 kb tandem duplication involving PRDM13 (chr6:99560265-99616492dup, hg38) at the MCDR1 locus. This duplication cosegregates with the NCMD phenotype in the five affected children. No other (likely) pathogenic variants in known inherited retinal disease genes were found in the WGS data. Bioinformatics analyses of the breakpoints suggest a replicative-based repair mechanism underlying the duplication.

Conclusions

We report a novel tandem duplication involving the PRDM13 gene in a family with NCMD from a previously unreported geographic region. The duplication size is the smallest that has been reported thus far and may correlate with the particular phenotype.

Dealing with Pseudogenes in Molecular Diagnostics in the Next Generation Sequencing Era

Presence of pseudogenes is a dreadful issue in next generation sequencing (NGS), because their contamination can interfere with the detection of variants in the genuine gene and generate false positive and false negative variants.In this chapter we focus on issues related to the application of NGS strategies for analysis of genes with pseudogenes in a clinical setting. The degree to which a pseudogene impacts the ability to accurately detect and map variants in its parent gene depends on the degree of similarity (homology) with the parent gene itself. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences. For target enrichment, we describe advantages and disadvantages of PCR- and capture-based strategies. For mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some examples of genes associated with Mendelian disorders, for which interesting NGS approaches are described to avoid interference with pseudogene sequences.

Long-Read Sequencing to Unravel Complex Structural Variants of CEP78 Leading to Cone-Rod Dystrophy and Hearing Loss

Inactivating variants as well as a missense variant in the centrosomal CEP78 gene have been identified in autosomal recessive cone-rod dystrophy with hearing loss (CRDHL), a rare syndromic inherited retinal disease distinct from Usher syndrome. Apart from this, a complex structural variant (SV) implicating CEP78 has been reported in CRDHL. Here we aimed to expand the genetic architecture of typical CRDHL by the identification of complex SVs of the CEP78 region and characterization of their underlying mechanisms. Approaches used for the identification of the SVs are shallow whole-genome sequencing (sWGS) combined with quantitative polymerase chain reaction (PCR) and long-range PCR, or ExomeDepth analysis on whole-exome sequencing (WES) data. Targeted or whole-genome nanopore long-read sequencing (LRS) was used to delineate breakpoint junctions at the nucleotide level. For all SVs cases, the effect of the SVs on CEP78 expression was assessed using quantitative PCR on patient-derived RNA. Apart from two novel canonical CEP78 splice variants and a frameshifting single-nucleotide variant (SNV), two SVs affecting CEP78 were identified in three unrelated individuals with CRDHL: a heterozygous total gene deletion of 235 kb and a partial gene deletion of 15 kb in a heterozygous and homozygous state, respectively. Assessment of the molecular consequences of the SVs on patient's materials displayed a loss-of-function effect. Delineation and characterization of the 15-kb deletion using targeted LRS revealed the previously described complex CEP78 SV, suggestive of a recurrent genomic rearrangement. A founder haplotype was demonstrated for the latter SV in cases of Belgian and British origin, respectively. The novel 235-kb deletion was delineated using whole-genome LRS. Breakpoint analysis showed microhomology and pointed to a replication-based underlying mechanism. Moreover, data mining of bulk and single-cell human and mouse transcriptional datasets, together with CEP78 immunostaining on human retina, linked the CEP78 expression domain with its phenotypic manifestations. Overall, this study supports that the CEP78 locus is prone to distinct SVs and that SV analysis should be considered in a genetic workup of CRDHL. Finally, it demonstrated the power of sWGS and both targeted and whole-genome LRS in identifying and characterizing complex SVs in patients with ocular diseases.

Phenocopy of a heterozygous carrier of X-linked retinitis pigmentosa due to mosaicism for a RHO variant

We describe both phenotype and pathogenesis in two male siblings with typical retinitis pigmentosa (RP) and the potentially X-linked RP (XLRP) carrier phenotype in their mother. Two affected sons, two unaffected daughters, and their mother underwent detailed ophthalmological assessments including Goldmann perimetry, color vision testing, multimodal imaging and ISCEV-standard electroretinography. Genetic testing consisted of targeted next-generation sequencing (NGS) of known XLRP genes and whole exome sequencing (WES) of known inherited retinal disease genes (RetNet-WES). Variant validation and segregation analysis were performed by Sanger sequencing. The mutational load of the RHO variant in the mother was assessed in DNA from leucocytes, buccal cells and hair follicles using targeted NGS. Both affected sons showed signs of classical RP, while the mother displayed patches of hyperautofluorescence on blue light autofluorescence imaging and regional, intraretinal, spicular pigmentation, reminiscent of a carrier phenotype of XLRP. XLRP testing was negative. RetNet-WES testing revealed RHO variant c.404G > C p.(Arg135Pro) in a mosaic state (21% of the reads) in the mother and in a heterozygous state in both sons. Targeted NGQSS of the RHO variant in different maternal tissues showed a mutation load between 25.06% and 41.72%. We report for the first time that somatic mosaicism of RHO variant c.404G > C p.(Arg135Pro) mimics the phenotype of a female carrier of XLRP, in combination with heterozygosity for the variant in the two affected sons.

Functional characterization of the first missense variant in CEP78, a founder allele associated with cone-rod dystrophy, hearing loss, and reduced male fertility

Inactivating variants in the centrosomal CEP78 gene have been found in cone-rod dystrophy with hearing loss (CRDHL), a particular phenotype distinct from Usher syndrome. Here, we identified and functionally characterized the first CEP78 missense variant c.449T>C, p.(Leu150Ser) in three CRDHL families. The variant was found in a biallelic state in two Belgian families and in a compound heterozygous state-in trans with c.1462-1G>T-in a third German family. Haplotype reconstruction showed a founder effect. Homology modeling revealed a detrimental effect of p.(Leu150Ser) on protein stability, which was corroborated in patients' fibroblasts. Elongated primary cilia without clear ultrastructural abnormalities in sperm or nasal brushes suggest impaired cilia assembly. Two affected males from different families displayed sperm abnormalities causing infertility. One of these is a heterozygous carrier of a complex allele in SPAG17, a ciliary gene previously associated with autosomal recessive male infertility. Taken together, our data indicate that a missense founder allele in CEP78 underlies the same sensorineural CRDHL phenotype previously associated with inactivating variants. Interestingly, the CEP78 phenotype has been possibly expanded with male infertility. Finally, CEP78 loss-of-function variants may have an underestimated role in misdiagnosed Usher syndrome, with or without sperm abnormalities.

ABCA4-associated disease as a model for missing heritability in autosomal recessive disorders: novel noncoding splice, cis-regulatory, structural, and recurrent hypomorphic variants

Purpose

ABCA4-associated disease, a recessive retinal dystrophy, is hallmarked by a large proportion of patients with only one pathogenic ABCA4 variant, suggestive for missing heritability.

Methods

By locus-specific analysis of ABCA4, combined with extensive functional studies, we aimed to unravel the missing alleles in a cohort of 67 patients (p), with one (p = 64) or no (p = 3) identified coding pathogenic variants of ABCA4.

Results

We identified eight pathogenic (deep-)intronic ABCA4 splice variants, of which five are novel and six structural variants, four of which are novel, including two duplications. Together, these variants account for the missing alleles in 40.3% of patients. Furthermore, two novel variants with a putative cis-regulatory effect were identified. The common hypomorphic variant c.5603A>T p.(Asn1868Ile) was found as a candidate second allele in 43.3% of patients. Overall, we have elucidated the missing heritability in 83.6% of our cohort. In addition, we successfully rescued three deep-intronic variants using antisense oligonucleotide (AON)-mediated treatment in HEK 293-T cells and in patient-derived fibroblast cells.

Conclusion

Noncoding pathogenic variants, novel structural variants, and a common hypomorphic allele of the ABCA4 gene explain the majority of unsolved cases with ABCA4-associated disease, rendering this retinopathy a model for missing heritability in autosomal recessive disorders.

Accurate detection and quantification of epigenetic and genetic second hits in BRCA1 and BRCA2-associated hereditary breast and ovarian cancer reveals multiple co-acting second hits

BACKGROUND

This study characterizes the second hit spectrum in BRCA1 and BRCA2-associated breast and ovarian cancers at both gene loci to investigate if second hit mechanisms are mutually exclusive or able to coincide within the same tumor.

METHODS

Loss of heterozygosity, somatic point mutations and copy number alterations along with promoter methylation were studied in 56 breast and 15 ovarian cancers from BRCA1 and BRCA2 germline mutation carriers. A mathematical methodology was introduced to quantify the tumor cell population carrying a second hit.

RESULTS

Copy neutral LOH was the most prevalent LOH mechanism in this cohort (BC 69%, OC 67%). However, only 36% of BC and 47% of OC showed LOH in all cancerous cells. Somatic intragenic deletions and methylated subclones were also found in combination with (partial) loss of heterozygosity. Unequivocal deleterious somatic point mutations were not identified in this cohort.

CONCLUSION

Different mechanisms inactivating the wild type allele are present within the same tumor sample at various extents. Results indicate that BRCA1/2-linked breast and ovarian cancer cells are predominantly characterized by LOH, but harbor a complex combination of second hits at various frequencies.

Thorough in silico and in vitro cDNA analysis of 21 putative BRCA1 and BRCA2 splice variants and a complex tandem duplication in BRCA2 allowing the identification of activated cryptic splice donor sites in BRCA2 exon 11

For 21 putative BRCA1 and BRCA2 splice site variants, the concordance between mRNA analysis and predictions by in silico programs was evaluated. Aberrant splicing was confirmed for 12 alterations. In silico prediction tools were helpful to determine for which variants cDNA analysis is warranted, however, predictions for variants in the Cartegni consensus region but outside the canonical sites, were less reliable. Learning algorithms like Adaboost and Random Forest outperformed the classical tools. Further validations are warranted prior to implementation of these novel tools in clinical settings. Additionally, we report here for the first time activated cryptic donor sites in the large exon 11 of BRCA2 by evaluating the effect at the cDNA level of a novel tandem duplication (5' breakpoint in intron 4; 3' breakpoint in exon 11) and of a variant disrupting the splice donor site of exon 11 (c.6841+1G > C). Additional sites were predicted, but not activated. These sites warrant further research to increase our knowledge on cis and trans acting factors involved in the conservation of correct transcription of this large exon. This may contribute to adequate design of ASOs (antisense oligonucleotides), an emerging therapy to render cancer cells sensitive to PARP inhibitor and platinum therapies.

NR5A1 is a novel disease gene for 46,XX testicular and ovotesticular disorders of sex development

PURPOSE

We aimed to identify the genetic cause in a cohort of 11 unrelated cases and two sisters with 46,XX SRY-negative (ovo)testicular disorders of sex development (DSD).

METHODS

Whole-exome sequencing (n = 9), targeted resequencing (n = 4), and haplotyping were performed. Immunohistochemistry of sex-specific markers was performed on patients' gonads. The consequences of mutation were investigated using luciferase assays, localization studies, and RNA-seq.

RESULTS

We identified a novel heterozygous NR5A1 mutation, c.274C>T p.(Arg92Trp), in three unrelated patients. The Arg92 residue is highly conserved and located in the Ftz-F1 region, probably involved in DNA-binding specificity and stability. There were no consistent changes in transcriptional activation or subcellular localization. Transcriptomics in patient-derived lymphocytes showed upregulation of MAMLD1, a direct NR5A1 target previously associated with 46,XY DSD. In gonads of affected individuals, ovarian FOXL2 and testicular SRY-independent SOX9 expression observed.

CONCLUSIONS

We propose NR5A1, previously associated with 46,XY DSD and 46,XX primary ovarian insufficiency, as a novel gene for 46,XX (ovo)testicular DSD. We hypothesize that p.(Arg92Trp) results in decreased inhibition of the male developmental pathway through downregulation of female antitestis genes, thereby tipping the balance toward testicular differentiation in 46,XX individuals. In conclusion, our study supports a role for NR5A1 in testis differentiation in the XX gonad.Genet Med 19 4, 367-376.

Next-generation sequencing shows West Nile virus quasispecies diversification after a single passage in a carrion crow (Corvus corone) in vivo infection model

West Nile virus (WNV) occurs as a population of genetic variants (quasispecies) infecting a single animal. Previous low-resolution viral genetic diversity estimates in sampled wild birds and mosquitoes, and in multiple-passage adaptation studies in vivo or in cell culture, suggest that WNV genetic diversification is mostly limited to the mosquito vector. This study investigated genetic diversification of WNV in avian hosts during a single passage using next-generation sequencing. Wild-captured carrion crows were subcutaneously infected using a clonal Middle-East WNV. Blood samples were collected 2 and 4 days post-infection. A reverse-transcription (RT)-PCR approach was used to amplify the WNV genome directly from serum samples prior to next-generation sequencing resulting in an average depth of at least 700 ×  in each sample. Appropriate controls were sequenced to discriminate biologically relevant low-frequency variants from experimentally introduced errors. The WNV populations in the wild crows showed significant diversification away from the inoculum virus quasispecies structure. By contrast, WNV populations in intracerebrally infected day-old chickens did not diversify from that of the inoculum. Where previous studies concluded that WNV genetic diversification is only experimentally demonstrated in its permissive insect vector species, we have experimentally shown significant diversification of WNV populations in a wild bird reservoir species.

Evaluation of convenient pretreatment protocols for RNA virus metagenomics in serum and tissue samples

Viral metagenomic approaches are increasingly being used for viral discovery. Various strategies are applied to enrich viral sequences, but there is often a lack of knowledge about their effective influence on the viral discovery sensitivity. We evaluate some convenient and widely used approaches for RNA virus discovery in clinical samples in order to reveal their sensitivity and potential bias introduced by the enrichment or amplifications steps. An RNA virus was artificially spiked at a fixed titer in serum and lung tissue, respectively, low and high nucleic acid content matrices. For serum, a simple DNase treatment on the RNA extract gave the maximum gain in proportion of viral sequences (83×), and a subsequent ribosomal RNA removal nearly doubled once more the proportion of viral sequences. For lung tissue, a ribosomal RNA depletion step on the RNA extract had the biggest gain in proportion of viral sequences (32×). We show also that direct sequencing of cDNA is recommended above an extra random PCR amplification step, and a that the virion enrichment strategy (filtration and nuclease treatment) has a beneficial effect for sequencing-based virus discovery. Our findings provide sample-dependent guidelines for targeted virus discovery strategies.

Next-generation sequencing in veterinary medicine: how can the massive amount of information arising from high-throughput technologies improve diagnosis, control, and management of infectious diseases?

The development of high-throughput molecular technologies and associated bioinformatics has dramatically changed the capacities of scientists to produce, handle, and analyze large amounts of genomic, transcriptomic, and proteomic data. A clear example of this step-change is represented by the amount of DNA sequence data that can be now produced using next-generation sequencing (NGS) platforms. Similarly, recent improvements in protein and peptide separation efficiencies and highly accurate mass spectrometry have promoted the identification and quantification of proteins in a given sample. These advancements in biotechnology have increasingly been applied to the study of animal infectious diseases and are beginning to revolutionize the way that biological and evolutionary processes can be studied at the molecular level. Studies have demonstrated the value of NGS technologies for molecular characterization, ranging from metagenomic characterization of unknown pathogens or microbial communities to molecular epidemiology and evolution of viral quasispecies. Moreover, high-throughput technologies now allow detailed studies of host-pathogen interactions at the level of their genomes (genomics), transcriptomes (transcriptomics), or proteomes (proteomics). Ultimately, the interaction between pathogen and host biological networks can be questioned by analytically integrating these levels (integrative OMICS and systems biology). The application of high-throughput biotechnology platforms in these fields and their typical low-cost per information content has revolutionized the resolution with which these processes can now be studied. The aim of this chapter is to provide a current and prospective view on the opportunities and challenges associated with the application of massive parallel sequencing technologies to veterinary medicine, with particular focus on applications that have a potential impact on disease control and management.

False-positive results in metagenomic virus discovery: a strong case for follow-up diagnosis

A viral metagenomic approach using virion enrichment, random amplification and next-generation sequencing was used to investigate an undiagnosed cluster of dairy cattle presenting with high persistent fever, unresponsive to anti-microbial and anti-inflammatory treatment, diarrhoea and redness of nose and teat. Serum and whole blood samples were taken in the predicted hyperviraemic state of an animal that a few days later presented with these clinical signs. Bioinformatics analysis of the resulting data from the DNA virus identification workflow (a total of 32 757 sequences with average read length 335 bases) initially demonstrated the presence of parvovirus-like sequences in the tested blood sample. Thorough follow-up using specific real-time RT-PCR assays targeting the detected sequence fragments confirmed the presence of these sequences in the original sample as well as in a sample of an additional animal, but a contamination with an identical genetic signature in negative extraction controls was demonstrated. Further investigation using an alternative extraction method identified a contamination of the originally used Qiagen extraction columns with parvovirus-like nucleic acids or virus particles. Although we did not find any relevant virus that could be associated with the disease, these observations clearly illustrate the importance of using a proper control strategy and follow-up diagnostic tests in any viral metagenomic study.

Full-Genome Sequencing of Four Bluetongue Virus Serotype 11 Viruses

Recently, a contamination incident was described in which the challenge inoculum used in a bluetongue virus serotype 8 (BTV-8) vaccination trial was contaminated with a BTV-11 virus that was closely related to the Belgian BTV-11 virus from 2008. This study reports the first complete genome sequences of four BTV-11 viruses: the BTV-11 contaminant, BTV-11 reference strain, BTV-11 vaccine strain and a recently isolated BTV-11 field strain from Martinique. Full-genome analysis showed that these viruses belong to serotype 11/nucleotype A and cluster together with other western topotype bluetongue viruses. Detailed comparisons of the genomes further indicated that the contaminant was derived from the BTV-11 reference strain, as they were distinguished by a single synonymous nucleotide substitution. The previously reported partial sequence of genome segment 2 of the Belgian BTV-11 was found to be identical to that of the BTV-11 vaccine strain, indicating that it most likely was the BTV-11 vaccine strain. These findings also suggest that the BTV-11 contaminant and the Belgian BTV-11 are not the same viruses. Finally, comparison of the reference and vaccine strain did not allow determining the amino acid substitutions that contribute to the attenuated phenotype.

Phylogeographic analysis of avian influenza viruses isolated from Charadriiformes in Belgium confirms intercontinental reassortment in gulls

Nine influenza viruses isolated from gulls and shorebirds in Belgium (2008-2010), including H3N8, H5N2, H6N1, H11N9, H13N6, H13N8, and H16N3 subtypes, were targeted using random amplification and next-generation sequencing. The gene segments of these viruses segregated into three phylogeographic lineage types: (1) segments circulating in waterfowl in Eurasia with sporadic introduction in other species and in the Americas ("Eurasian avian"), (2) segments circulating in American waterfowl with sporadic introduction to other species and regions ("American avian"), and (3) segments circulating exclusively in gulls and shorebirds and having increased connectivity between the two hemispheres ("Charadriiformes specific"). Notably, an H6N1 and an H5N2 isolated from L. argentatus had mainly Eurasian avian genes but shared a matrix segment of American avian origin (first documentation in European gulls of transhemispheric reassortment). These data support the growing evidence of an important role of Charadriiformes birds in the dynamic nature of avian influenza ecology.

Identification and complete genome sequencing of paramyxoviruses in mallard ducks (Anas platyrhynchos) using random access amplification and next generation sequencing technologies

Background

During a wildlife screening program for avian influenza A viruses (AIV) and avian paramyxoviruses (APMV) in Belgium, we isolated two hemagglutinating agents from pools of cloacal swabs of wild mallards (Anas platyrhynchos) caught in a single sampling site at two different times. AIV and APMV1 were excluded using hemagglutination inhibition (HI) testing and specific real-time RT-PCR tests.

Methods

To refine the virological identification of APMV2-10 realized by HI subtyping tests and in lack of validated molecular tests for APMV2-10, random access amplification was used in combination with next generation sequencing for the sequence independent identification of the viruses and the determination of their genomes.

Results

Three different APMVs were identified. From one pooled sample, the complete genome sequence (15054 nucleotides) of an APMV4 was assembled from the random sequences. From the second pooled sample, the nearly complete genome sequence of an APMV6 (genome size of 16236 nucleotides) was determined, as well as a partial sequence for an APMV4. This APMV4 was closely related but not identical to the APMV4 isolated from the first sample. Although a cross-reactivity with other APMV subtypes did not allow formal identification, the HI subtyping revealed APMV4 and APMV6 in the respective pooled samples but failed to identify the co-infecting APMV4 in the APMV6 infected pool.

Conclusions

These data further contribute to the knowledge about the genetic diversity within the serotypes APMV4 and 6, and confirm the limited sensitivity of the HI subtyping test. Moreover, this study demonstrates the value of a random access nucleic acid amplification method in combination with massive parallel sequencing. Using only a moderate and economical sequencing effort, the characterization and full genome sequencing of APMVs can be obtained, including the identification of viruses in mixed infections.