The monoclonal antibody combination REGEN-COV protects against SARS-CoV-2 mutational escape in preclinical and human studies.

Monoclonal antibodies against SARS-CoV-2 are a clinically validated therapeutic option against COVID-19. Because rapidly emerging virus mutants are becoming the next major concern in the fight against the global pandemic, it is imperative that these therapeutic treatments provide coverage against circulating variants and do not contribute to development of treatment-induced emergent resistance. To this end, we investigated the sequence diversity of the spike protein and monitored emergence of virus variants in SARS-COV-2 isolates found in COVID-19 patients treated with the two-antibody combination REGEN-COV, as well as in preclinical in vitro studies using single, dual, or triple antibody combinations, and in hamster in vivo studies using REGEN-COV or single monoclonal antibody treatments. Our study demonstrates that the combination of non-competing antibodies in REGEN-COV provides protection against all current SARS-CoV-2 variants of concern/interest and also protects against emergence of new variants and their potential seeding into the population in a clinical setting.

Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca).

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

Use of Whole-Genome Sequencing to Estimate the Contribution of Immune Evasion and Waning Immunity on Decreasing COVID-19 Vaccine Effectiveness.

BACKGROUND: The impact variant-specific immune evasion and waning protection have on declining coronavirus disease 2019 (COVID-19) vaccine effectiveness (VE) remains unclear. Using whole-genome sequencing (WGS), we examined the contribution these factors had on the decline that followed the introduction of the Delta variant. Furthermore, we evaluated calendar-period-based classification as a WGS alternative. METHODS: We conducted a test-negative case-control study among people tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between 1 April and 24 August 2021. Variants were classified using WGS and calendar period. RESULTS: We included 2029 cases (positive, sequenced samples) and 343 727 controls (negative tests). VE 14-89 days after second dose was significantly higher against Alpha (84.4%; 95% confidence interval [CI], 75.6%-90.0%) than Delta infection (68.9%; 95% CI, 58.0%-77.1%). The odds of Delta infection were significantly higher 90-149 than 14-89 days after second dose (P value = .003). Calendar-period-classified VE estimates approximated WGS-classified estimates; however, calendar-period-based classification was subject to misclassification (35% Alpha, 4% Delta). CONCLUSIONS: Both waning protection and variant-specific immune evasion contributed to the lower effectiveness. While calendar-period-classified VE estimates mirrored WGS-classified estimates, our analysis highlights the need for WGS when variants are cocirculating and misclassification is likely.

Common genetic risk factors in ASD and ADHD co-occurring families.

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

The lineage-specific transcription factor CDX2 navigates dynamic chromatin to control distinct stages of intestine development.

Lineage-restricted transcription factors, such as the intestine-specifying factor CDX2, often have dual requirements across developmental time. Embryonic loss of CDX2 triggers homeotic transformation of intestinal fate, whereas adult-onset loss compromises crucial physiological functions but preserves intestinal identity. It is unclear how such diverse requirements are executed across the developmental continuum. Using primary and engineered human tissues, mouse genetics, and a multi-omics approach, we demonstrate that divergent CDX2 loss-of-function phenotypes in embryonic versus adult intestines correspond to divergent CDX2 chromatin-binding profiles in embryonic versus adult stages. CDX2 binds and activates distinct target genes in developing versus adult mouse and human intestinal cells. We find that temporal shifts in chromatin accessibility correspond to these context-specific CDX2 activities. Thus, CDX2 is not sufficient to activate a mature intestinal program; rather, CDX2 responds to its environment, targeting stage-specific genes to contribute to either intestinal patterning or mature intestinal function. This study provides insights into the mechanisms through which lineage-specific regulatory factors achieve divergent functions over developmental time.

A YY1-dependent increase in aerobic metabolism is indispensable for intestinal organogenesis.

During late gestation, villi extend into the intestinal lumen to dramatically increase the surface area of the intestinal epithelium, preparing the gut for the neonatal diet. Incomplete development of the intestine is the most common gastrointestinal complication in neonates, but the causes are unclear. We provide evidence in mice that Yin Yang 1 (Yy1) is crucial for intestinal villus development. YY1 loss in the developing endoderm had no apparent consequences until late gestation, after which the intestine differentiated poorly and exhibited severely stunted villi. Transcriptome analysis revealed that YY1 is required for mitochondrial gene expression, and ultrastructural analysis confirmed compromised mitochondrial integrity in the mutant intestine. We found increased oxidative phosphorylation gene expression at the onset of villus elongation, suggesting that aerobic respiration might function as a regulator of villus growth. Mitochondrial inhibitors blocked villus growth in a fashion similar to Yy1 loss, thus further linking oxidative phosphorylation with late-gestation intestinal development. Interestingly, we find that necrotizing enterocolitis patients also exhibit decreased expression of oxidative phosphorylation genes. Our study highlights the still unappreciated role of metabolic regulation during organogenesis, and suggests that it might contribute to neonatal gastrointestinal disorders.

PipelineDog: a simple and flexible graphic pipeline construction and maintenance tool.

Summary: Analysis pipelines are an essential part of bioinformatics research, and ad hoc pipelines are frequently created by researchers for prototyping and proof-of-concept purposes. However, most existing pipeline management system or workflow engines are too complex for rapid prototyping or learning the pipeline concept. A lightweight, user-friendly and flexible solution is thus desirable. In this study, we developed a new pipeline construction and maintenance tool, PipelineDog. This is a web-based integrated development environment with a modern web graphical user interface. It offers cross-platform compatibility, project management capabilities, code formatting and error checking functions and an online repository. It uses an easy-to-read/write script system that encourages code reuse. With the online repository, it also encourages sharing of pipelines, which enhances analysis reproducibility and accountability. For most users, PipelineDog requires no software installation. Overall, this web application provides a way to rapidly create and easily manage pipelines. Availability and implementation: PipelineDog web app is freely available at http://web.pipeline.dog. The command line version is available at http://www.npmjs.com/package/pipelinedog and online repository at http://repo.pipeline.dog. Contact: ysun@kean.edu or xing@biology.rutgers.edu or ysun@diagnoa.com. Supplementary information: Supplementary data are available at Bioinformatics online.

Ancient hybridization and admixture in macaques (genus Macaca) inferred from whole genome sequences.

The evolutionary history of the stump-tailed macaque (Macaca arctoides) and its genetic relationship to other macaques is a subject of continuing controversy. Here, we have reported the first genome sequences of two stump-tailed macaques and one Assamese macaque (M. assamensis). Additionally, we have investigated the genetic diversity between macaque species and analyzed ancient hybridization events. Genome-wide analyses demonstrated that the stump-tailed macaque is more closely related to sinica species than to fascicularis/mulatta species. This topology contradicts the mitochondrial sequence-based phylogeny that places the stump-tailed macaque into the fascicularis/mulatta group. However, our results further show that stump-tailed macaques have genetic backgrounds distinct from sinica species, and present evidence of gene flows with rhesus macaques. We suggest that an ancient introgression occurred after stump-tailed macaques diverged from sinica species. The distinct gene flow between proto-arctoides and proto-mulatta resulted in the transfer of rhesus macaque-type mitochondria into proto-arctoides. The rhesus macaque-type mitochondria remained in populations because of genetic drift during the bottleneck. The PSMC results and morphological and geographic evidence are consistent with the mitochondria capture pattern in the stump-tailed macaque. The molecular clock estimates suggest that the mitochondrial transference into stump-tailed macaques occurred 0.4-1.4 million years ago. Furthermore, we detected extensive admixtures between different macaque species, indicating that gene flow has played an important role in the evolutionary history of the genus Macaca.

Extremely low-coverage whole genome sequencing in South Asians captures population genomics information.

BACKGROUND: The cost of Whole Genome Sequencing (WGS) has decreased tremendously in recent years due to advances in next-generation sequencing technologies. Nevertheless, the cost of carrying out large-scale cohort studies using WGS is still daunting. Past simulation studies with coverage at ~2x have shown promise for using low coverage WGS in studies focused on variant discovery, association study replications, and population genomics characterization. However, the performance of low coverage WGS in populations with a complex history and no reference panel remains to be determined. RESULTS: South Indian populations are known to have a complex population structure and are an example of a major population group that lacks adequate reference panels. To test the performance of extremely low-coverage WGS (EXL-WGS) in populations with a complex history and to provide a reference resource for South Indian populations, we performed EXL-WGS on 185 South Indian individuals from eight populations to ~1.6x coverage. Using two variant discovery pipelines, SNPTools and GATK, we generated a consensus call set that has ~90% sensitivity for identifying common variants (minor allele frequency ≥ 10%). Imputation further improves the sensitivity of our call set. In addition, we obtained high-coverage for the whole mitochondrial genome to infer the maternal lineage evolutionary history of the Indian samples. CONCLUSIONS: Overall, we demonstrate that EXL-WGS with imputation can be a valuable study design for variant discovery with a dramatically lower cost than standard WGS, even in populations with a complex history and without available reference data. In addition, the South Indian EXL-WGS data generated in this study will provide a valuable resource for future Indian genomic studies.

Identification and characterization of Aurora kinase B and C variants associated with maternal aneuploidy.

STUDY QUESTION: Are single nucleotide variants (SNVs) in Aurora kinases B and C (AURKB, AURKC) associated with risk of aneuploid conception? SUMMARY ANSWER: Two SNVs were found in patients with extreme aneuploid concepti rates with respect to their age; one variant, AURKC p.I79V, is benign, while another, AURKB p.L39P, is a potential gain-of-function mutant with increased efficiency in promoting chromosome alignment. WHAT IS KNOWN ALREADY: Maternal age does not always predict aneuploidy risk, and rare gene variants can be drivers of disease. The AURKB and AURKC regulate chromosome segregation, and are associated with reproductive impairments in mouse and human. STUDY DESIGN, SIZE, DURATION: An extreme phenotype sample selection scheme was performed for variant discovery. Ninety-six DNA samples were from young patients with higher than average embryonic aneuploidy rates and an additional 96 DNA samples were from older patients with lower than average aneuploidy rates. PARTICIPANTS/MATERIALS, SETTING, METHODS: Using the192 DNA samples, the coding regions of AURKB and AURKC were sequenced using next generation sequencing. To assess biological significance, we expressed complementary RNA encoding the human variants in mouse oocytes. Assays such as determining subcellular localization and assessing catalytic activity were performed to determine alterations in protein function during meiosis. MAIN RESULTS AND THE ROLE OF CHANCE: Ten SNVs were identified using three independent variant-calling methods. Two of the SNVs (AURKB p.L39P and AURKC p.I79V) were non-synonymous and identified by at least two variant-identification methods. The variant encoding AURKC p.I79V, identified in a young woman with a higher than average rate of aneuploid embryos, showed wild-type localization pattern and catalytic activity. On the other hand, the variant encoding AURKB p.L39P, identified in an older woman with lower than average rates of aneuploid embryos, increased the protein's ability to regulate alignment of chromosomes at the metaphase plate. These experiments were repeated three independent times using 2-3 mice for each trial. LARGE SCALE DATA: N/A. LIMITATIONS, REASONS FOR CAUTION: Biological significance of the human variants was assessed in an in vitro mouse oocyte model where the variants are over-expressed. Therefore, the human protein may not function identically to the mouse homolog, or the same in mouse oocytes as in human oocytes. Furthermore, supraphysiological expression levels may not accurately reflect endogenous activity. Moreover, the evaluated variants were identified in one patient each, and no trial linking the SNV to pregnancy outcomes was conducted. Finally, the patient aneuploidy rates were established by performing comprehensive chromosome screening in blastocysts, and because of the link between female gamete aneuploidy giving rise to aneuploid embryos, we evaluate the role of the variants in Meiosis I. However, it is possible that the chromosome segregation mistake arose during Meiosis II or in mitosis in the preimplantation embryo. Their implications in human female meiosis and aneuploidy risk remain to be determined. WIDER IMPLICATIONS OF THE FINDINGS: The data provide evidence that gene variants exist in reproductively younger or advanced aged women that are predictive of the risk of producing aneuploid concepti in humans. Furthermore, a single amino acid in the N-terminus of AURKB is a gain-of-function mutant that could be protective of euploidy. STUDY FUNDING/COMPETING INTERESTS: This work was supported by a Research Grant from the American Society of Reproductive Medicine and support from the Charles and Johanna Busch Memorial Fund at Rutgers, the State University of NJ to K.S. and the Foundation for Embryonic Competence, Inc to N.T. The authors declare no conflicts of interest.

MicroRNA and MicroRNA-Target Variants Associated with Autism Spectrum Disorder and Related Disorders.

Autism spectrum disorder (ASD) is a childhood neurodevelopmental disorder with a complex and heterogeneous genetic etiology. MicroRNA (miRNA), a class of small non-coding RNAs, could regulate ASD risk genes post-transcriptionally and affect broad molecular pathways related to ASD and associated disorders. Using whole-genome sequencing, we analyzed 272 samples in 73 families in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Families with at least one ASD patient were recruited and were further assessed for language impairment, reading impairment, and other associated phenotypes. A total of 5104 miRNA variants and 1,181,148 3' untranslated region (3' UTR) variants were identified in the dataset. After applying several filtering criteria, including population allele frequency, brain expression, miRNA functional regions, and inheritance patterns, we identified high-confidence variants in five brain-expressed miRNAs (targeting 326 genes) and 3' UTR miRNA target regions of 152 genes. Some genes, such as SCP2 and UCGC, were identified in multiple families. Using Gene Ontology overrepresentation analysis and protein-protein interaction network analysis, we identified clusters of genes and pathways that are important for neurodevelopment. The miRNAs and miRNA target genes identified in this study are potentially involved in neurodevelopmental disorders and should be considered for further functional studies.

Proteogenomic identification of Hepatitis B virus (HBV) genotype-specific HLA-I restricted peptides from HBV-positive patient liver tissues.

The presentation of virus-derived peptides by HLA class I molecules on the surface of an infected cell and the recognition of these HLA-peptide complexes by, and subsequent activation of, CD8+ cytotoxic T cells provides an important mechanism for immune protection against viruses. Recent advances in proteogenomics have allowed researchers to discover a growing number of unique HLA-restricted viral peptides, resulting in a rapidly expanding repertoire of targets for immunotherapeutics (i.e. bispecific antibodies, engineered T-cell receptors (TCRs), chimeric antigen receptor T-cells (CAR-Ts)) to infected tissues. However, genomic variability between viral strains, such as Hepatitis-B virus (HBV), in combination with differences in patient HLA alleles, make it difficult to develop therapeutics against these targets. To address this challenge, we developed a novel proteogenomics approach for generating patient-specific databases that enable the identification of viral peptides based on the viral transcriptomes sequenced from individual patient liver samples. We also utilized DNA sequencing of patient samples to identify HLA genotypes and assist in target selection. Liver samples from 48 HBV infected patients, primarily from Asia, were examined to reconstruct patient-specific HBV genomes, identify regions within the human chromosomes targeted by HBV integrations and obtain a comprehensive view of HBV peptide epitopes using our HLA class-I (HLA-I) immunopeptidomics discovery platform. Two previously reported HLA associated HBV-derived peptides, HLA-A02 binder FLLTRILTI (S194-202) from the large surface antigen and HLA-A11 binder STLPETTVVRR (C141-151) from the capsid protein were validated by our discovery platform, but both were detected at very low frequencies. In addition, we identified and validated, using heavy peptide analogues, novel strain-specific HBV-HLA associated peptides, such as GSLPQEHIVQK (P606-616) and variants. Overall, our novel approach can guide the development of bispecific antibody, TCR-T, or CAR-T based therapeutics for the treatment of HBV-related HCC and inform vaccine development.

REGN-COV2 antibodies prevent and treat SARS-CoV-2 infection in rhesus macaques and hamsters.

An urgent global quest for effective therapies to prevent and treat coronavirus disease 2019 (COVID-19) is ongoing. We previously described REGN-COV2, a cocktail of two potent neutralizing antibodies (REGN10987 and REGN10933) that targets nonoverlapping epitopes on the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein. In this report, we evaluate the in vivo efficacy of this antibody cocktail in both rhesus macaques, which may model mild disease, and golden hamsters, which may model more severe disease. We demonstrate that REGN-COV-2 can greatly reduce virus load in the lower and upper airways and decrease virus-induced pathological sequelae when administered prophylactically or therapeutically in rhesus macaques. Similarly, administration in hamsters limits weight loss and decreases lung titers and evidence of pneumonia in the lungs. Our results provide evidence of the therapeutic potential of this antibody cocktail.

Evaluating nanopore sequencing data processing pipelines for structural variation identification.

BACKGROUND: Structural variations (SVs) account for about 1% of the differences among human genomes and play a significant role in phenotypic variation and disease susceptibility. The emerging nanopore sequencing technology can generate long sequence reads and can potentially provide accurate SV identification. However, the tools for aligning long-read data and detecting SVs have not been thoroughly evaluated. RESULTS: Using four nanopore datasets, including both empirical and simulated reads, we evaluate four alignment tools and three SV detection tools. We also evaluate the impact of sequencing depth on SV detection. Finally, we develop a machine learning approach to integrate call sets from multiple pipelines. Overall SV callers' performance varies depending on the SV types. For an initial data assessment, we recommend using aligner minimap2 in combination with SV caller Sniffles because of their speed and relatively balanced performance. For detailed analysis, we recommend incorporating information from multiple call sets to improve the SV call performance. CONCLUSIONS: We present a workflow for evaluating aligners and SV callers for nanopore sequencing data and approaches for integrating multiple call sets. Our results indicate that additional optimizations are needed to improve SV detection accuracy and sensitivity, and an integrated call set can provide enhanced performance. The nanopore technology is improving, and the sequencing community is likely to grow accordingly. In turn, better benchmark call sets will be available to more accurately assess the performance of available tools and facilitate further tool development.

SMAD4 Suppresses WNT-Driven Dedifferentiation and Oncogenesis in the Differentiated Gut Epithelium.

The cell of origin of colon cancer is typically thought to be the resident somatic stem cells, which are immortal and escape the continual cellular turnover characteristic of the intestinal epithelium. However, recent studies have identified certain conditions in which differentiated cells can acquire stem-like properties and give rise to tumors. Defining the origins of tumors will inform cancer prevention efforts as well as cancer therapies, as cancers with distinct origins often respond differently to treatments. We report here a new condition in which tumors arise from the differentiated intestinal epithelium. Inactivation of the differentiation-promoting transcription factor SMAD4 in the intestinal epithelium was surprisingly well tolerated in the short term. However, after several months, adenomas developed with characteristics of activated WNT signaling. Simultaneous loss of SMAD4 and activation of the WNT pathway led to dedifferentiation and rapid adenoma formation in differentiated tissue. Transcriptional profiling revealed acquisition of stem cell characteristics, and colabeling indicated that cells expressing differentiated enterocyte markers entered the cell cycle and reexpressed stem cell genes upon simultaneous loss of SMAD4 and activation of the WNT pathway. These results indicate that SMAD4 functions to maintain differentiated enterocytes in the presence of oncogenic WNT signaling, thus preventing dedifferentiation and tumor formation in the differentiated intestinal epithelium.Significance: This work identifies a mechanism through which differentiated cells prevent tumor formation by suppressing oncogenic plasticity. Cancer Res; 78(17); 4878-90. ©2018 AACR.

Degree of Tissue Differentiation Dictates Susceptibility to BRAF-Driven Colorectal Cancer.

Oncogenic mutations in BRAF are believed to initiate serrated colorectal cancers; however, the mechanisms of BRAF-driven colon cancer are unclear. We find that oncogenic BRAF paradoxically suppresses stem cell renewal and instead promotes differentiation. Correspondingly, tumor formation is inefficient in BRAF-driven mouse models of colon cancer. By reducing levels of differentiation via genetic manipulation of either of two distinct differentiation-promoting factors (Smad4 or Cdx2), stem cell activity is restored in BRAFV600E intestines, and the oncogenic capacity of BRAFV600E is amplified. In human patients, we observe that reduced levels of differentiation in normal tissue is associated with increased susceptibility to serrated colon tumors. Together, these findings help resolve the conditions necessary for BRAF-driven colon cancer initiation. Additionally, our results predict that genetic and/or environmental factors that reduce tissue differentiation will increase susceptibility to serrated colon cancer. These findings offer an opportunity to identify susceptible individuals by assessing their tissue-differentiation status.