Toward better understanding of artifacts in variant calling from high-coverage samples

Variant calling and benchmarking in an era of complete human genome sequences

Genetic variant calling from DNA sequencing has enabled understanding of germline variation in hundreds of thousands of humans. Sequencing technologies and variant-calling methods have advanced rapidly, routinely providing reliable variant calls in most of the human genome. We describe how advances in long reads, deep learning, de novo assembly and pangenomes have expanded access to variant calls in increasingly challenging, repetitive genomic regions, including medically relevant regions, and how new benchmark sets and benchmarking methods illuminate their strengths and limitations. Finally, we explore the possible future of more complete characterization of human genome variation in light of the recent completion of a telomere-to-telomere human genome reference assembly and human pangenomes, and we consider the innovations needed to benchmark their newly accessible repetitive regions and complex variants.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Large-scale analysis of de novo mutations identifies risk genes for female infertility characterized by oocyte and early embryo defects

BACKGROUND

Oocyte maturation arrest and early embryonic arrest are important reproductive phenotypes resulting in female infertility and cause the recurrent failure of assisted reproductive technology (ART). However, the genetic etiologies of these female infertility-related phenotypes are poorly understood. Previous studies have mainly focused on inherited mutations based on large pedigrees or consanguineous patients. However, the role of de novo mutations (DNMs) in these phenotypes remains to be elucidated.

RESULTS

To decipher the role of DNMs in ART failure and female infertility with oocyte and embryo defects, we explore the landscape of DNMs in 473 infertile parent-child trios and identify a set of 481 confident DNMs distributed in 474 genes. Gene ontology analysis reveals that the identified genes with DNMs are enriched in signaling pathways associated with female reproductive processes such as meiosis, embryonic development, and reproductive structure development. We perform functional assays on the effects of DNMs in a representative gene Tubulin Alpha 4a (TUBA4A), which shows the most significant enrichment of DNMs in the infertile parent-child trios. DNMs in TUBA4A disrupt the normal assembly of the microtubule network in HeLa cells, and microinjection of DNM TUBA4A cRNAs causes abnormalities in mouse oocyte maturation or embryo development, suggesting the pathogenic role of these DNMs in TUBA4A.

CONCLUSIONS

Our findings suggest novel genetic insights that DNMs contribute to female infertility with oocyte and embryo defects. This study also provides potential genetic markers and facilitates the genetic diagnosis of recurrent ART failure and female infertility.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Genetic dysregulation of an endothelial Ras signaling network in vein of Galen malformations

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and severe congenital brain arteriovenous malformation, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP ( RASA1 ) harbored a genome-wide significant burden of loss-of-function de novo variants (p=4.79×10 ^-7 ). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 ( EPHB4 ) (p=1.22×10 ^-5 ), which cooperates with p120 RasGAP to limit Ras activation. Other probands had pathogenic variants in ACVRL1 , NOTCH1 , ITGB1 , and PTPN11 . ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomics defined developing endothelial cells as a key spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant exhibited constitutive endothelial Ras/ERK/MAPK activation and impaired hierarchical development of angiogenesis-regulated arterial-capillary-venous networks, but only when carrying a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have clinical implications.

Genome-wide identification of copy neutral loss of heterozygosity reveals its possible association with spatial positioning of chromosomes

Loss of heterozygosity (LOH) is a genetic alteration that results from the loss of one allele at a heterozygous locus. In particular, copy neutral LOH (CN-LOH) events are generated, for example, by mitotic homologous recombination after monoallelic defection or gene conversion, resulting in novel homozygous locus having two copies of the normal counterpart allele. This phenomenon can serve as a source of genome diversity and is associated with various diseases. To clarify the nature of the CN-LOH such as the frequency, genomic distribution and inheritance pattern, we made use of whole-genome sequencing data of the three-generation CEPH/Utah family cohort, with the pedigree consisting of grandparents, parents and offspring. We identified an average of 40.7 CN-LOH events per individual taking advantage of 285 healthy individuals from 33 families in the cohort. On average 65% of them were classified as gonosomal-mosaicism-associated CN-LOH, which exists in both germline and somatic cells. We also confirmed that the incidence of the CN-LOH has little to do with the parents' age and sex. Furthermore, through the analysis of the genomic region including the CN-LOH, we found that the chance of the occurrence of the CN-LOH tends to increase at the GC-rich locus and/or on the chromosome having a relatively close inter-homolog distance. We expect that these results provide significant insights into the association between genetic alteration and spatial position of chromosomes as well as the intrinsic genetic property of the CN-LOH.

Performance comparisons between clustering models for reconstructing NGS results from technical replicates

To improve the performance of individual DNA sequencing results, researchers often use replicates from the same individual and various statistical clustering models to reconstruct a high-performance callset. Here, three technical replicates of genome NA12878 were considered and five model types were compared (consensus, latent class, Gaussian mixture, Kamila–adapted k-means, and random forest) regarding four performance indicators: sensitivity, precision, accuracy, and F1-score. In comparison with no use of a combination model, i) the consensus model improved precision by 0.1%; ii) the latent class model brought 1% precision improvement (97%–98%) without compromising sensitivity (= 98.9%); iii) the Gaussian mixture model and random forest provided callsets with higher precisions (both >99%) but lower sensitivities; iv) Kamila increased precision (>99%) and kept a high sensitivity (98.8%); it showed the best overall performance. According to precision and F1-score indicators, the compared non-supervised clustering models that combine multiple callsets are able to improve sequencing performance vs. previously used supervised models. Among the models compared, the Gaussian mixture model and Kamila offered non-negligible precision and F1-score improvements. These models may be thus recommended for callset reconstruction (from either biological or technical replicates) for diagnostic or precision medicine purposes.

Population Structure and Genetic Diversity Analysis of “Yufen 1” H Line Chickens Using Whole-Genome Resequencing

The effective protection and utilization of poultry resources depend on an accurate understanding of the genetic diversity and population structure. The breeding of the specialized poultry lineage “Yufen 1”, with its defined characteristics, was approved by the China Poultry Genetic Resource Committee in 2015. Thus, to investigate the relationship between the progenitor H line and other poultry breeds, the genetic diversity and population structure of “Yufen 1” H line (YF) were investigated and compared with those of 2 commercial chicken breeds, the ancestor breed Red Jungle Fowls, and 11 Chinese indigenous chicken breeds based on a whole-genome resequencing approach using 8,112,424 SNPs. The genetic diversity of YF was low, and the rate of linkage disequilibrium decay was significantly slower than that of the other Chinese indigenous breeds. In addition, it was shown that the YF population was strongly selected during intensive breeding and that genetic resources have been seriously threatened, which highlights the need to establish a systematic conservation strategy as well as utilization techniques to maintain genetic diversity within YF. Moreover, a principal component analysis, a neighbor-joining tree analysis, a structure analysis, and genetic differentiation indices indicated that YF harbors a distinctive genetic resource with a unique genetic structure separate from that of Chinese indigenous breeds at the genome level. The findings provide a valuable resource and the theoretical basis for the further conservation and utilization of YF.

Revisiting mutagenesis at non-B DNA motifs in the human genome

Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.

Horizontal Gene Transfer, Fitness Costs and Mobility Shape the Spread of Antibiotic Resistance Genes into Experimental Populations of Acinetobacter Baylyi

Horizontal gene transfer (HGT) is important for microbial evolution, but how evolutionary forces shape the frequencies of horizontally transferred genetic variants in the absence of strong selection remains an open question. In this study, we evolve laboratory populations of Acinetobacter baylyi (ADP1) with HGT from two clinically relevant strains of multidrug-resistant Acinetobacter baumannii (AB5075 and A9844). We find that DNA can cross the species barrier, even without strong selection, and despite substantial DNA sequence divergence between the two species. Our results confirm previous findings that HGT can drive the spread of antibiotic resistance genes (ARGs) without selection for that antibiotic, but not for all of the resistance genes present in the donor genome. We quantify the costs and benefits of horizontally transferred variants and use whole population sequencing to track the spread of ARGs from HGT donors into antibiotic-sensitive recipients. We find that even though most ARGs are taken up by populations of A. baylyi, the long-term fate of an individual gene depends both on its fitness cost and on the type of genetic element that carries the gene. Interestingly, we also found that an integron, but not its host plasmid, is able to spread in A. baylyi populations despite its strong deleterious effect. Altogether, our results show how HGT provides an evolutionary advantage to evolving populations by facilitating the spread of non-selected genetic variation including costly ARGs.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

Introgressive hybridization in the west Pacific pen shells (genus Atrina): Restricted interspecies gene flow within the genome

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.

Molecular insights into the dynamics of species invasion by hybridisation in Tasmanian eucalypts

In plants where seed dispersal is limited compared with pollen dispersal, hybridisation may enhance gene exchange and species dispersal. We provide genetic evidence of hybridisation contributing to the expansion of the rare Eucalyptus risdonii into the range of the widespread Eucalyptus amygdalina. These closely related tree species are morphologically distinct, and observations suggest that natural hybrids occur along their distribution boundaries and as isolated trees or in small patches within the range of E. amygdalina. Hybrid phenotypes occur outside the range of normal dispersal for E. risdonii seed, yet in some hybrid patches small individuals resembling E. risdonii occur and are hypothesised to be a result of backcrossing. Using 3362 genome-wide SNPs assessed from 97 individuals of E. risdonii and E. amygdalina and 171 hybrid trees, we show that (i) isolated hybrids match the genotypes expected of F₁ /F₂ hybrids, (ii) there is a continuum in the genetic composition among the isolated hybrid patches from patches dominated by F₁ /F₂ -like genotypes to those dominated by E. risdonii-backcross genotypes, and (iii) the E. risdonii-like phenotypes in the isolated hybrid patches are most-closely related to proximal larger hybrids. These results suggest that the E. risdonii phenotype has been resurrected in isolated hybrid patches established from pollen dispersal, providing the first steps in its invasion of suitable habitat by long-distance pollen dispersal and complete introgressive displacement of E. amygdalina. Such expansion accords with the population demographics, common garden performance data, and climate modelling which favours E. risdonii and highlights a role of interspecific hybridisation in climate change adaptation and species expansion.

Advantages of long- and short-reads sequencing for the hybrid investigation of the Mycobacterium tuberculosis genome

Introduction

In the fight to limit the global spread of antibiotic resistance, computational challenges associated with sequencing technology can impact the accuracy of downstream analysis, including drug resistance identification, transmission, and genome resolution. About 10% of Mycobacterium tuberculosis (MTB) genome is constituted by the PE/PPE family, a GC-rich repetitive genome region. Although sequencing using short read technology is widely used, it is well recognized its limit in the PE/PPE regions due to the unambiguously mapping process onto the reference genome. The aim of this study was to compare the performances of short-reads (SRS), long-reads (LRS) and hybrid-reads (HYBR) based analysis over different common investigative tasks: genome coverage estimation, variant calling and cluster analysis, drug resistance detection and de novo assembly.

Methods

For the study 13 model MTB clinical isolates were sequenced with both SRS and LRS. HYBR were produced correcting the long reads with the short reads. The fastq from the three approaches were then processed using a customized version of MTBseq for genome coverage estimation and variant calling and using two different assemblers for de novo assembly evaluation.

Results

Estimation of genome coverage performances showed lower 8X breadth coverage for SRS respect to LRS and HYBR: considering the PE/PPE genes, SRS showed low results for the PE_PGRS family, while obtained acceptable coverage in PE and PPE genes; LRS and HYBR reached optimal coverages in PE/PPE genes. For variant calling HYBR showed the highest resolution, detecting the highest percentage of uniquely identified mutations compared to LRS and SRS. All three approaches agreed on the identification of two major clusters, with HYBR identifying an higher number of SNPs between the two clusters. Comparing the quality of the assemblies, HYBR and LRS obtained better results than SRS.

Discussion

In conclusion, depending on the aim of the investigation, both SRS and LRS present complementary advantages and limitations implying that for a full resolution of MTB genomes, where all the mentioned analyses and both technologies are needed, the use of the HYBR approach represents a valid option and a well-rounded strategy.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

An efficient genotyper and star-allele caller for pharmacogenomics

High-throughput sequencing provides sufficient means for determining genotypes of clinically important pharmacogenes that can be used to tailor medical decisions to individual patients. However, pharmacogene genotyping, also known as star-allele calling, is a challenging problem that requires accurate copy number calling, structural variation identification, variant calling, and phasing within each pharmacogene copy present in the sample. Here we introduce Aldy 4, a fast and efficient tool for genotyping pharmacogenes that uses combinatorial optimization for accurate star-allele calling across different sequencing technologies. Aldy 4 adds support for long reads and uses a novel phasing model and improved copy number and variant calling models. We compare Aldy 4 against the current state-of-the-art star-allele callers on a large and diverse set of samples and genes sequenced by various sequencing technologies, such as whole-genome and targeted Illumina sequencing, barcoded 10x Genomics, and Pacific Biosciences (PacBio) HiFi. We show that Aldy 4 is the most accurate star-allele caller with near-perfect accuracy in all evaluated contexts, and hope that Aldy remains an invaluable tool in the clinical toolbox even with the advent of long-read sequencing technologies.

Pharmacogenomic and Statistical Analysis

Genetic variants can alter response to drugs and other therapeutic interventions. The study of this phenomenon, called pharmacogenomics, is similar in many ways to other types of genetic studies but has distinct methodological and statistical considerations. Genetic variants involved in the processing of exogenous compounds exhibit great diversity and complexity, and the phenotypes studied in pharmacogenomics are also more complex than typical genetic studies. In this chapter, we review basic concepts in pharmacogenomic study designs, data generation techniques, statistical analysis approaches, and commonly used methods and briefly discuss the ultimate translation of findings to clinical care.

Multi-omics profiling visualizes dynamics of cardiac development and functions

Cardiogenesis is a tightly regulated dynamic process through a continuum of differentiation and proliferation events. Key factors and pathways governing this process remain incompletely understood. Here, we investigate mice hearts from embryonic day 10.5 to postnatal week 8 and dissect developmental changes in phosphoproteome-, proteome-, metabolome-, and transcriptome-encompassing cardiogenesis and cardiac maturation. We identify mitogen-activated protein kinases as core kinases involved in transcriptional regulation by mediating the phosphorylation of chromatin remodeling proteins during early cardiogenesis. We construct the reciprocal regulatory network of transcription factors (TFs) and identify a series of TFs controlling early cardiogenesis involved in cycling-dependent proliferation. After birth, we identify cardiac resident macrophages with high arachidonic acid metabolism activities likely involved in the clearance of injured apoptotic cardiomyocytes. Together, our comprehensive multi-omics data offer a panoramic view of cardiac development and maturation that provides a resource for further in-depth functional exploration.

Validation of Polygenic Risk Scores for Coronary Heart Disease in a Middle Eastern Cohort Using Whole Genome Sequencing

BACKGROUND

Enthusiasm for using polygenic risk scores (PRSs) in clinical practice is tempered by concerns about their portability to diverse ancestry groups, thus motivating genome-wide association studies in non-European ancestry cohorts.

METHODS

We conducted a genome-wide association study for coronary heart disease in a Middle Eastern cohort using whole genome sequencing and assessed the performance of 6 PRSs developed with methods including LDpred (PGS000296), metaGRS (PGS000018), Pruning and Thresholding (PGS000337), and an EnsemblePRS we developed. Additionally, we evaluated the burden of rare variants in lipid genes in cases and controls. Whole genome sequencing at 30× coverage was performed in 1067 coronary heart disease cases (mean age=59 years; 70.3% males) and 6170 controls (mean age=40 years; 43.5% males).

RESULTS

The majority of PRSs performed well; odds ratio (OR) per 1 SD increase (OR1sd) was highest for PGS000337 (OR1sd=1.81, 95% CI [1.66-1.98], P=3.07×10-41). EnsemblePRS performed better than individual PRSs (OR1sd=1.8, 95% CI [1.66-1.96], P=5.89×10-44). The OR for the 10th decile versus the remaining deciles was >3.2 for PGS000337, PGS000296, PGS000018, and reached 4.58 for EnsemblePRS. Of 400 known genome-wide significant loci, 33 replicated at P<10-4. However, the 9p21 locus did not replicate. Six suggestive (P<10-5) new loci/genes with plausible biological function were identified (eg, CORO7, RBM47, PDE4D). The burden of rare functional variants in LDLR, APOB, PCSK9, and ANGPTL4 was greater in cases than controls.

CONCLUSIONS

Overall, we demonstrate that PRSs derived from European ancestry genome-wide association studies performed well in a Middle Eastern cohort, suggesting these could be used in the clinical setting while ancestry-specific PRSs are developed.

Gene Sequencing Identifies Perturbation in Nitric Oxide Signaling as a Nonlipid Molecular Subtype of Coronary Artery Disease

BACKGROUND

A key goal of precision medicine is to disaggregate common, complex diseases into discrete molecular subtypes. Rare coding variants in the low-density lipoprotein receptor gene (LDLR) are identified in 1% to 2% of coronary artery disease (CAD) patients, defining a molecular subtype with risk driven by hypercholesterolemia.

METHODS

To search for additional subtypes, we compared the frequency of rare, predicted loss-of-function and damaging missense variants aggregated within a given gene in 41 081 CAD cases versus 217 115 controls.

RESULTS

Rare variants in LDLR were most strongly associated with CAD, present in 1% of cases and associated with 4.4-fold increased CAD risk. A second subtype was characterized by variants in endothelial nitric oxide synthase gene (NOS3), a key enzyme regulating vascular tone, endothelial function, and platelet aggregation. A rare predicted loss-of-function or damaging missense variants in NOS3 was present in 0.6% of cases and associated with 2.42-fold increased risk of CAD (95% CI, 1.80-3.26; P=5.50×10-9). These variants were associated with higher systolic blood pressure (+3.25 mm Hg; [95% CI, 1.86-4.65]; P=5.00×10-6) and increased risk of hypertension (adjusted odds ratio 1.31; [95% CI, 1.14-1.51]; P=2.00×10-4) but not circulating cholesterol concentrations, suggesting that, beyond lipid pathways, nitric oxide synthesis is a key nonlipid driver of CAD risk.

CONCLUSIONS

Beyond LDLR, we identified an additional nonlipid molecular subtype of CAD characterized by rare variants in the NOS3 gene.

DNA Methylation Analysis Reveals Distinct Patterns in Satellite Cell-Derived Myogenic Progenitor Cells of Subjects with Spastic Cerebral Palsy

Spastic type cerebral palsy (CP) is a complex neuromuscular disorder that involves altered skeletal muscle microanatomy and growth, but little is known about the mechanisms contributing to muscle pathophysiology and dysfunction. Traditional genomic approaches have provided limited insight regarding disease onset and severity, but recent epigenomic studies indicate that DNA methylation patterns can be altered in CP. Here, we examined whether a diagnosis of spastic CP is associated with intrinsic DNA methylation differences in myoblasts and myotubes derived from muscle resident stem cell populations (satellite cells; SCs). Twelve subjects were enrolled (6 CP; 6 control) with informed consent/assent. Skeletal muscle biopsies were obtained during orthopedic surgeries, and SCs were isolated and cultured to establish patient-specific myoblast cell lines capable of proliferation and differentiation in culture. DNA methylation analyses indicated significant differences at 525 individual CpG sites in proliferating SC-derived myoblasts (MB) and 1774 CpG sites in differentiating SC-derived myotubes (MT). Of these, 79 CpG sites were common in both culture types. The distribution of differentially methylated 1 Mbp chromosomal segments indicated distinct regional hypo- and hyper-methylation patterns, and significant enrichment of differentially methylated sites on chromosomes 12, 13, 14, 15, 18, and 20. Average methylation load across 2000 bp regions flanking transcriptional start sites was significantly different in 3 genes in MBs, and 10 genes in MTs. SC derived MBs isolated from study participants with spastic CP exhibited fundamental differences in DNA methylation compared to controls at multiple levels of organization that may reveal new targets for studies of mechanisms contributing to muscle dysregulation in spastic CP.

African Suid Genomes Provide Insights into the Local Adaptation to Diverse African Environments

African wild suids consist of several endemic species that represent ancient members of the family Suidae and have colonized diverse habitats on the African continent. However, limited genomic resources for African wild suids hinder our understanding of their evolution and genetic diversity. In this study, we assembled high-quality genomes of a common warthog (Phacochoerus africanus), a red river hog (Potamochoerus porcus), as well as an East Asian Diannan small-ear pig (Sus scrofa). Phylogenetic analysis showed that common warthog and red river hog diverged from their common ancestor around the Miocene/Pliocene boundary, putatively predating their entry into Africa. We detected species-specific selective signals associated with sensory perception and interferon signaling pathways in common warthog and red river hog, respectively, which contributed to their local adaptation to savannah and tropical rainforest environments, respectively. The structural variation and evolving signals in genes involved in T-cell immunity, viral infection, and lymphoid development were identified in their ancestral lineage. Our results provide new insights into the evolutionary histories and divergent genetic adaptations of African suids.

Understanding Mendelian errors in SNP arrays data using a Gochu Asturcelta pig pedigree: genomic alterations, family size and calling errors

Up to 478 Gochu Asturcelta pig parents-offspring trios (61 different families) were genotyped using the Axiom_PigHDv1 Array to identify the causes of Mendelian errors (ME). Up to 545,364 SNPs were retained. Up to 40,540 SNPs gathering 292,297 allelic mismatches were identified and were overlapped with SINEs and LINEs (Sscrofa genome 11.1). Copy number variations (CNV) were called using PennCNV. ME were classified into eight different classes according to the trio member ("Trio" meaning no assignment) and the allele on which ME was identified: TrioA/B, FatherA/B, MotherA/B, OffspringA/B. Most ME occurred due to systematic causes: (a) those assigned to the Father, Mother or Offspring occurred by null or partial null alleles characterized by heterozygote deficiency, varied with family size, involved a low number of loci (6506), and gathered most mismatches (228,145); (b) TrioB errors varied with family size, covaried with SINEs, LINEs and CNV, and involved most ME loci (33,483) and mismatches (65,682); and (c) TrioA errors were non-systematic ME with no sampling bias involving 1.2% of mismatches only and a low number of loci (1939). The influence of TrioB errors on the overall genotyping quality may be low and, since CNV vary among populations, their removal should be considered in each particular dataset. ME assignable to the Father, Mother or Offspring may be consistent within technological platforms and may bias severely linkage or association studies. Most ME caused by null or partial null alleles can be removed using heterozygote deficiency without affecting the size of the datasets.

Genomic trajectories of a near-extinction event in the Chatham Island black robin

BACKGROUND

Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation.

RESULTS

We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations.

CONCLUSION

Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

A database of 5305 healthy Korean individuals reveals genetic and clinical implications for an East Asian population

Despite substantial advances in disease genetics, studies to date have largely focused on individuals of European descent. This limits further discoveries of novel functional genetic variants in other ethnic groups. To alleviate the paucity of East Asian population genome resources, we established the Korean Variant Archive 2 (KOVA 2), which is composed of 1896 whole-genome sequences and 3409 whole-exome sequences from healthy individuals of Korean ethnicity. This is the largest genome database from the ethnic Korean population to date, surpassing the 1909 Korean individuals deposited in gnomAD. The variants in KOVA 2 displayed all the known genetic features of those from previous genome databases, and we compiled data from Korean-specific runs of homozygosity, positively selected intervals, and structural variants. In doing so, we found loci, such as the loci of ADH1A/1B and UHRF1BP1, that are strongly selected in the Korean population relative to other East Asian populations. Our analysis of allele ages revealed a correlation between variant functionality and evolutionary age. The data can be browsed and downloaded from a public website ( https://www.kobic.re.kr/kova/ ). We anticipate that KOVA 2 will serve as a valuable resource for genetic studies involving East Asian populations.

Analytical and Clinical Validation of Cell-Free Circulating Tumor DNA Assay for the Estimation of Tumor Mutational Burden

Background

Ultra-deep sequencing to detect low-frequency mutations in circulating tumor-derived DNA (ctDNA) increases the diagnostic value of liquid biopsy. The demand for large ctDNA panels for comprehensive genomic profiling and tumor mutational burden (TMB) estimation is increasing; however, few ctDNA panels for TMB have been validated. Here, we designed a ctDNA panel with 531 genes, named TMB500, along with a technical and clinical validation.

Methods

Synthetic reference cell-free DNA materials with predefined allele frequencies were sequenced in a total of 92 tests in 6 batches to evaluate the precision, linearity, and limit of detection of the assay. We used clinical samples from 50 patients with various cancers, 11 healthy individuals, and paired tissue samples. Molecular barcoding and data analysis were performed using customized pipelines.

Results

The assay showed high precision and linearity (coefficient of determination, r2 = 0.87) for all single nucleotide variants, with a limit of detection of 0.24%. In clinical samples, the TMB500 ctDNA assay detected most variants present and absent in tissues, showing that ctDNA could assess tumor heterogeneity in different tissues and metastasis sites. The estimated TMBs correlated well between tissue and blood, except in 4 cases with extreme heterogeneity that showed very high blood TMBs compared to tissue TMBs. A pilot evaluation showed that the TMB500 assay could be used for disease monitoring.

Conclusions

The TMB500 assay is an accurate and reliable ctDNA assay for many clinical purposes. It may be useful for guiding the treatment of cancers with diverse genomic profiles, estimating TMB in immune therapy, and disease monitoring.

Whole-exome sequencing of epithelial ovarian carcinomas differing in resistance to platinum therapy

Exploration of the prognostic and predictive significance of exome variation in epithelial ovarian carcinoma patients, with TP53, Hippo, homologous recombination genes, and the SBS6 signature as the most interesting results.

Epithelial ovarian carcinoma (EOC) is highly fatal because of the risk of resistance to therapy and recurrence. We performed whole-exome sequencing of blood and tumor tissue pairs of 50 patients with surgically resected EOC. Compared with sensitive patients, platinum-resistant patients had a significantly higher somatic mutational rate in TP53 and lower in several genes from the Hippo pathway. We confirmed the pivotal role of somatic mutations in homologous recombination repair genes in platinum sensitivity and favorable prognosis of EOC patients. Implementing the germline homologous recombination repair profile significantly improved the prediction. In addition, distinct mutational signatures, for example, SBS6, and overall mutational load, somatic mutations in PABPC1, PABPC3, and TFAM co-segregated with the resistance status, high-grade serous carcinoma subtype, or overall survival of patients. We generated germline and somatic genetic landscapes of prognostically different subgroups of EOC patients for further follow-up studies focused on utilizing the observed associations in precision oncology.

Novel trends of genome evolution in highly complex tropical sponge microbiomes

BACKGROUND

Tropical members of the sponge genus Ircinia possess highly complex microbiomes that perform a broad spectrum of chemical processes that influence host fitness. Despite the pervasive role of microbiomes in Ircinia biology, it is still unknown how they remain in stable association across tropical species. To address this question, we performed a comparative analysis of the microbiomes of 11 Ircinia species using whole-metagenomic shotgun sequencing data to investigate three aspects of bacterial symbiont genomes-the redundancy in metabolic pathways across taxa, the evolution of genes involved in pathogenesis, and the nature of selection acting on genes relevant to secondary metabolism.

RESULTS

A total of 424 new, high-quality bacterial metagenome-assembled genomes (MAGs) were produced for 10 Caribbean Ircinia species, which were evaluated alongside 113 publicly available MAGs sourced from the Pacific species Ircinia ramosa. Evidence of redundancy was discovered in that the core genes of several primary metabolic pathways could be found in the genomes of multiple bacterial taxa. Across hosts, the metagenomes were depleted in genes relevant to pathogenicity and enriched in eukaryotic-like proteins (ELPs) that likely mimic the hosts' molecular patterning. Finally, clusters of steroid biosynthesis genes (CSGs), which appear to be under purifying selection and undergo horizontal gene transfer, were found to be a defining feature of Ircinia metagenomes.

CONCLUSIONS

These results illustrate patterns of genome evolution within highly complex microbiomes that illuminate how associations with hosts are maintained. The metabolic redundancy within the microbiomes could help buffer the hosts from changes in the ambient chemical and physical regimes and from fluctuations in the population sizes of the individual microbial strains that make up the microbiome. Additionally, the enrichment of ELPs and depletion of LPS and cellular motility genes provide a model for how alternative strategies to virulence can evolve in microbiomes undergoing mixed-mode transmission that do not ultimately result in higher levels of damage (i.e., pathogenicity) to the host. Our last set of results provides evidence that sterol biosynthesis in Ircinia-associated bacteria is widespread and that these molecules are important for the survival of bacteria in highly complex Ircinia microbiomes. Video Abstract.

Integrating de novo and inherited variants in 42,607 autism cases identifies mutations in new moderate-risk genes

To capture the full spectrum of genetic risk for autism, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 autism cases, including 35,130 new cases recruited online by SPARK. We identified 60 genes with exome-wide significance (P < 2.5 × 10-6), including five new risk genes (NAV3, ITSN1, MARK2, SCAF1 and HNRNPUL2). The association of NAV3 with autism risk is primarily driven by rare inherited loss-of-function (LoF) variants, with an estimated relative risk of 4, consistent with moderate effect. Autistic individuals with LoF variants in the four moderate-risk genes (NAV3, ITSN1, SCAF1 and HNRNPUL2; n = 95) have less cognitive impairment than 129 autistic individuals with LoF variants in highly penetrant genes (CHD8, SCN2A, ADNP, FOXP1 and SHANK3) (59% vs 88%, P = 1.9 × 10-6). Power calculations suggest that much larger numbers of autism cases are needed to identify additional moderate-risk genes.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

AnFiSA: An open-source computational platform for the analysis of sequencing data for rare genetic disease

Despite genomic sequencing rapidly transforming from being a bench-side tool to a routine procedure in a hospital, there is a noticeable lack of genomic analysis software that supports both clinical and research workflows as well as crowdsourcing. Furthermore, most existing software packages are not forward-compatible in regards to supporting ever-changing diagnostic rules adopted by the genetics community. Regular updates of genomics databases pose challenges for reproducible and traceable automated genetic diagnostics tools. Lastly, most of the software tools score low on explainability amongst clinicians. We have created a fully open-source variant curation tool, AnFiSA, with the intention to invite and accept contributions from clinicians, researchers, and professional software developers. The design of AnFiSA addresses the aforementioned issues via the following architectural principles: using a multidimensional database management system (DBMS) for genomic data to address reproducibility, curated decision trees adaptable to changing clinical rules, and a crowdsourcing-friendly interface to address difficult-to-diagnose cases. We discuss how we have chosen our technology stack and describe the design and implementation of the software. Finally, we show in detail how selected workflows can be implemented using the current version of AnFiSA by a medical geneticist.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Chronology of natural selection in Oceanian genomes

As human populations left Asia to first settle in Oceania around 50,000 years ago, they entered a territory ecologically separated from the Old World for millions of years. We analyzed genomic data of 239 modern Oceanian individuals to detect and date signals of selection specific to this region. Combining both relative and absolute dating approaches, we identified a strong selection pattern between 52,000 and 54,000 years ago in the genomes of descendants of the first settlers of Sahul. This strikingly corresponds to the dates of initial settlement as inferred from archaeological evidence. Loci under selection during this period, some showing enrichment in Denisovan ancestry, overlap genes involved in the immune response and diet, especially based on plants. Pathogens and natural resources, especially from endemic plants, therefore appear to have acted as strong selective pressures on the genomes of the first settlers of Sahul.

Reduced representation sequencing to understand the evolutionary history of Torrey pine (Pinus torreyana Parry) with implications for rare species conservation

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from inter-population genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has likely contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter-population admixture is necessary to empirically evaluate the trade-offs associated with genetic rescue in Torrey pine.

Exome sequencing analysis of Japanese autism spectrum disorder case-control sample supports an increased burden of synaptic function-related genes

Autism spectrum disorder (ASD) is a highly heritable, complex disorder in which rare variants contribute significantly to disease risk. Although many genes have been associated with ASD, there have been few genetic studies of ASD in the Japanese population. In whole exomes from a Japanese ASD sample of 309 cases and 299 controls, rare variants were associated with ASD within specific neurodevelopmental gene sets, including highly constrained genes, fragile X mental retardation protein target genes, and genes involved in synaptic function, with the strongest enrichment in trans-synaptic signaling (p = 4.4 × 10^-4, Q-value = 0.06). In particular, we strengthen the evidence regarding the role of ABCA13, a synaptic function-related gene, in Japanese ASD. The overall results of this case-control exome study showed that rare variants related to synaptic function are associated with ASD susceptibility in the Japanese population.

Software Choice and Sequencing Coverage Can Impact Plastid Genome Assembly-A Case Study in the Narrow Endemic Calligonum bakuense

Most plastid genome sequences are assembled from short-read whole-genome sequencing data, yet the impact that sequencing coverage and the choice of assembly software can have on the accuracy of the resulting assemblies is poorly understood. In this study, we test the impact of both factors on plastid genome assembly in the threatened and rare endemic shrub Calligonum bakuense. We aim to characterize the differences across plastid genome assemblies generated by different assembly software tools and levels of sequencing coverage and to determine if these differences are large enough to affect the phylogenetic position inferred for C. bakuense compared to congeners. Four assembly software tools (FastPlast, GetOrganelle, IOGA, and NOVOPlasty) and seven levels of sequencing coverage across the plastid genome (original sequencing depth, 2,000x, 1,000x, 500x, 250x, 100x, and 50x) are compared in our analyses. The resulting assemblies are evaluated with regard to reproducibility, contig number, gene complement, inverted repeat length, and computation time; the impact of sequence differences on phylogenetic reconstruction is assessed. Our results show that software choice can have a considerable impact on the accuracy and reproducibility of plastid genome assembly and that GetOrganelle produces the most consistent assemblies for C. bakuense. Moreover, we demonstrate that a sequencing coverage between 500x and 100x can reduce both the sequence variability across assembly contigs and computation time. When comparing the most reliable plastid genome assemblies of C. bakuense, a sequence difference in only three nucleotide positions is detected, which is less than the difference potentially introduced through software choice.

The evaluation of Bcftools mpileup and GATK HaplotypeCaller for variant calling in non-human species

Identification of genetic variations is a central part of population and quantitative genomics studies based on high-throughput sequencing data. Even though popular variant callers such as Bcftools mpileup and GATK HaplotypeCaller were developed nearly 10 years ago, their performance is still largely unknown for non-human species. Here, we showed by benchmark analyses with a simulated insect population that Bcftools mpileup performs better than GATK HaplotypeCaller in terms of recovery rate and accuracy regardless of mapping software. The vast majority of false positives were observed from repeats, especially for GATK HaplotypeCaller. Variant scores calculated by GATK did not clearly distinguish true positives from false positives in the vast majority of cases, implying that hard-filtering with GATK could be challenging. These results suggest that Bcftools mpileup may be the first choice for non-human studies and that variants within repeats might have to be excluded for downstream analyses.

Genetic Analysis of Six Transmembrane Protein Family Genes in Parkinson's Disease in a Large Chinese Cohort

Objectives

Parkinson's disease (PD) is a neurodegenerative disorder with the manifestation of motor symptoms and non-motor symptoms. Previous studies have indicated the role of several transmembrane (TMEM) protein family genes in PD pathogenesis.

Materials and Methods

In order to better investigate the genetic role of PD-related TMEM protein family genes in PD, including TMEM230, TMEM59, TMEM108, TMEM163, TMEM175, and TMEM229B, 1,917 sporadic early onset PD (sEOPD) or familial PD (FPD) patients and 1,652 healthy controls were analyzed by whole-exome sequencing (WES) while 1,962 sporadic late-onset PD (sLOPD) and 1,279 healthy controls were analyzed by whole-genome sequencing (WGS). Rare and common variants for each gene were included in the analysis.

Results

One hundred rare damaging or loss of function variants of six genes were found at the threshold of MAF < 0.1%. Three rare Dmis variants of TMEM230 were specifically identified in PD. Rare missense variants of TMEM59 were statistically significantly associated with PD in the WES cohort, indicating the role of TMEM59 in FPD and sEOPD. Rare missense variants of TMEM108 were suggestively associated with PD in the WGS cohort, indicating the potential role of TMEM108 in sLOPD. The rare variant of the other three genes and common variants of six genes were not significantly associated with PD.

Conclusion

We performed a large case-control study to systematically investigate the role of several PD-related TMEM protein family genes in PD. We identified three PD-specific variants in TMEM230, the significant association of TMEM59 with FPD, and sEOPD and the suggestive association of TMEM108 with sLOPD.

De novo Mutations in Domestic Cat are Consistent with an Effect of Reproductive Longevity on Both the Rate and Spectrum of Mutations

The mutation rate is a fundamental evolutionary parameter with direct and appreciable effects on the health and function of individuals. Here, we examine this important parameter in the domestic cat, a beloved companion animal as well as a valuable biomedical model. We estimate a mutation rate of 0.86 × 10-8 per bp per generation for the domestic cat (at an average parental age of 3.8 years). We find evidence for a significant paternal age effect, with more mutations transmitted by older sires. Our analyses suggest that the cat and the human have accrued similar numbers of mutations in the germline before reaching sexual maturity. The per-generation mutation rate in the cat is 28% lower than what has been observed in humans, but is consistent with the shorter generation time in the cat. Using a model of reproductive longevity, which takes into account differences in the reproductive age and time to sexual maturity, we are able to explain much of the difference in per-generation rates between species. We further apply our reproductive longevity model in a novel analysis of mutation spectra and find that the spectrum for the cat resembles the human mutation spectrum at a younger age of reproduction. Together, these results implicate changes in life-history as a driver of mutation rate evolution between species. As the first direct observation of the paternal age effect outside of rodents and primates, our results also suggest a phenomenon that may be universal among mammals.

Association of Pathogenic DNA Variants Predisposing to Cardiomyopathy With Cardiovascular Disease Outcomes and All-Cause Mortality

Importance

Pathogenic variants associated with inherited cardiomyopathy are recognized as important and clinically actionable when identified, leading some clinicians to recommend population-wide genomic screening.

Objective

To determine the prevalence and clinical importance of pathogenic variants associated with inherited cardiomyopathy within the context of contemporary clinical care.

Design, Setting, and Participants

This was a genetic association study of participants in Atherosclerosis in Risk Communities (ARIC), recruited from 1987 to 1989, with median follow-up of 27 years, and the UK Biobank, recruited from 2006 to 2010, with median follow-up of 10 years. ARIC participants were recruited from 4 sites across the US. UK Biobank participants were recruited from 22 sites across the UK. Participants in the US were of African and European ancestry; those in the UK were of African, East Asian, South Asian, and European ancestry. Statistical analyses were performed between August 1, 2021, and February 9, 2022.

Exposures

Rare genetic variants predisposing to inherited cardiomyopathy.

Main Outcomes and Measures

Pathogenicity of observed DNA sequence variants in sequenced exomes of 13 genes (ACTC1, FLNC, GLA, LMNA, MYBPC3, MYH7, MYL2, MYL3, PRKAG2, TNNI3, TNNT2, TPM1, and TTN) associated with inherited cardiomyopathies were classified by a blinded clinical geneticist per American College of Medical Genetics recommendations. Incidence of all-cause mortality, heart failure, and atrial fibrillation were determined. Cardiac magnetic resonance imaging, echocardiography, and electrocardiogram measures were assessed in a subset of participants.

Results

A total of 9667 ARIC participants (mean [SD] age, 54.0 [5.7] years; 4232 women [43.8%]; 2658 African [27.5%] and 7009 European [72.5%] ancestry) and 49 744 UK Biobank participants (mean [SD] age, 57.1 [8.0] years; 27 142 women [54.5%]; 1006 African [2.0%], 173 East Asian [0.3%], 939 South Asian [1.9%], and 46 449 European [93.4%] European ancestry) were included in the study. Of those, 59 participants (0.61%) in ARIC and 364 participants (0.73%) in UK Biobank harbored an actionable pathogenic or likely pathogenic variant associated with dilated or hypertrophic cardiomyopathy. Carriers of these variants were not reliably identifiable by imaging. However, the presence of these variants was associated with increased risk of heart failure (hazard ratio [HR], 1.7; 95% CI, 1.1-2.8), atrial fibrillation (HR, 2.9; 95% CI, 1.9-4.5), and all-cause mortality (HR, 1.5; 95% CI, 1.1-2.2) in ARIC. Similar risk patterns were observed in the UK Biobank.

Conclusions and Relevance

Results of this genetic association study suggest that approximately 0.7% of study participants harbored a pathogenic variant associated with inherited cardiomyopathy. These variant carriers would be challenging to identify within clinical practice without genetic testing but are at increased risk for cardiovascular disease and all-cause mortality.

Inherited basis of visceral, abdominal subcutaneous and gluteofemoral fat depots

For any given level of overall adiposity, individuals vary considerably in fat distribution. The inherited basis of fat distribution in the general population is not fully understood. Here, we study up to 38,965 UK Biobank participants with MRI-derived visceral (VAT), abdominal subcutaneous (ASAT), and gluteofemoral (GFAT) adipose tissue volumes. Because these fat depot volumes are highly correlated with BMI, we additionally study six local adiposity traits: VAT adjusted for BMI and height (VATadj), ASATadj, GFATadj, VAT/ASAT, VAT/GFAT, and ASAT/GFAT. We identify 250 independent common variants (39 newly-identified) associated with at least one trait, with many associations more pronounced in female participants. Rare variant association studies extend prior evidence for PDE3B as an important modulator of fat distribution. Local adiposity traits (1) highlight depot-specific genetic architecture and (2) enable construction of depot-specific polygenic scores that have divergent associations with type 2 diabetes and coronary artery disease. These results - using MRI-derived, BMI-independent measures of local adiposity - confirm fat distribution as a highly heritable trait with important implications for cardiometabolic health outcomes.

Patient-reported outcomes in a Chinese cohort of osteogenesis imperfecta unveil psycho-physical stratifications associated with clinical manifestations

Background

Osteogenesis imperfecta (OI) is a rare congenital disorder of the skeletal system, inflicting debilitating physical and psychological distress on patients and caregivers. Over the decades, much effort has been channeled towards understanding molecular mechanisms and developing new treatments. It has recently become more apparent that patient-reported outcome measurements (PROM) during treatment, healing and rehabilitation are helpful in facilitating smoother communication, refining intervention strategies and achieving higher quality of life. To date, systematic analyses of PROM in OI patients remain scarce.

Results

Here, utilizing a PROM Information System, we report a cross-sectional and longitudinal study in a southern Chinese cohort of 90 OI patients, covering both the child and adult age-groups. In the child group where both self and parental surveys were obtained, we identified two clusters of comparable sizes showing different outlooks in physical mobility and emotional experiences. One cluster (Cluster 1) is more negative about themselves than the other (Cluster 2). A concordance of 84.7% between self and parental assessments was recorded, suggesting the stability and validity of PROM-based stratification. Clinical subtyping, deformity, leg length discrepancy, and limited joint mobility were significantly associated with this stratification, with Cluster 1 showing higher percentages of severe phenotypes than Cluster 2. Since OI is a genetic disorder, we performed genetic testing on 72 of the 90 patients, but found no obvious association between genotypes and the PROM stratification. Analyses of longitudinal data suggested that patients tended to stay in the same psychological state, in both clusters. Adult patients also showed a continuous spectrum of self-evaluation that matches their clinical manifestations.

Conclusion

By systematically analyzing patient-reported outcomes, our study demonstrated the link between the sociopsychological wellbeing of OI patients, and their clinical manifestations, which may serve as the basis for evaluating clinical interventions and help achieve better patient-centric medical practices. The lack of genotype-PROM association may be due to the diverse mutational spectrum in OI, which warrants further investigation when a larger sample size is available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02394-7.

Moment estimators of relatedness from low-depth whole-genome sequencing data

Background

Estimating relatedness is an important step for many genetic study designs. A variety of methods for estimating coefficients of pairwise relatedness from genotype data have been proposed. Both the kinship coefficient \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ and the fraternity coefficient \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi$$\end{document}ψ for all pairs of individuals are of interest. However, when dealing with low-depth sequencing or imputation data, individual level genotypes cannot be confidently called. To ignore such uncertainty is known to result in biased estimates. Accordingly, methods have recently been developed to estimate kinship from uncertain genotypes.

Results

We present new method-of-moment estimators of both the coefficients \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi$$\end{document}ψ calculated directly from genotype likelihoods. We have simulated low-depth genetic data for a sample of individuals with extensive relatedness by using the complex pedigree of the known genetic isolates of Cilento in South Italy. Through this simulation, we explore the behaviour of our estimators, demonstrate their properties, and show advantages over alternative methods. A demonstration of our method is given for a sample of 150 French individuals with down-sampled sequencing data.

Conclusions

We find that our method can provide accurate relatedness estimates whilst holding advantages over existing methods in terms of robustness, independence from external software, and required computation time. The method presented in this paper is referred to as LowKi (Low-depth Kinship) and has been made available in an R package (https://github.com/genostats/LowKi).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04795-8.

RDscan: A New Method for Improving Germline and Somatic Variant Calling Based on Read Depth Distribution

Several tools have been developed for calling variants from next-generation sequencing (NGS) data. Although they are generally accurate and reliable, most of them have room for improvement, especially regarding calling variants in datasets with low read depth. In addition, the somatic variants predicted by several somatic variant callers tend to have very low concordance rates. In this study, we developed a new method (RDscan) for improving germline and somatic variant calling in NGS data. RDscan removes misaligned reads, repositions reads, and calculates RDscore based on the read depth distribution. With RDscore, RDscan improves the precision of variant callers by removing false-positive variant calls. When we tested our new tool using the latest variant calling algorithms and data from the 1000 Genomes Project and Illumina's public datasets, accuracy was improved for most of the algorithms. After screening variants with RDscan, calling accuracies increased for germline variants in 11 of 12 cases and for somatic variants in 21 of 24 cases. RDscan is simple to use and can effectively remove false-positive variants while maintaining a low computation load. Therefore, RDscan, along with existing variant callers, should contribute to improvements in genome analysis.

Identifying Rare Genetic Variants of Immune Mediators as Risk Factors for Autism Spectrum Disorder

Autism spectrum disorder (ASD) affects more than 1% of children, and there is no viable pharmacotherapeutic agent to treat the core symptoms of ASD. Studies have shown that children with ASD show changes in their levels of immune response molecules. Our previous studies have shown that ASD is more common in children with folate receptor autoantibodies. We also found that children with ASD have abnormal gut immune function, which was characterized by a significant increase in the content of immunoglobulin A and an increase in gut-microbiota-associated epitope diversity. These studies suggest that the immune mechanism plays an important role in the occurrence of ASD. The present study aims to systematically assess gene mutations in immune mediators in patients with ASD. We collected genetic samples from 72 children with ASD (2−12 years old) and 107 healthy controls without ASD (20−78 years old). We used our previously-designed immune gene panel, which can capture cytokine and receptor genes, the coding regions of MHC genes, and genes of innate immunity. Target region sequencing (500×) and bioinformatics analytical methods were used to identify variants in immune response genes associated with patients with ASD. A total of 4 rare variants were found to be associated with ASD, including HLA-B: p.A93G, HLA-DQB1: p.S229N, LILRB2: p.R322H, and LILRB2: c.956-4C>T. These variants were present in 44.44% (32/72) of the ASD patients and were detected in 3.74% (4/107) of the healthy controls. We expect these genetic variants will serve as new targets for the clinical genetic assessment of ASD, and our findings suggest that immune abnormalities in children with ASD may have a genetic basis.

Protocol for unbiased, consolidated variant calling from whole exome sequencing data

Whole Exome Sequencing (WES) is used for querying DNA variants using the protein coding parts of genomes (exomes). However, WES analysis can be challenging because of the complexity of the data. Here, we describe a consolidated protocol for unbiased WES analysis. The protocol uses three variant callers (HaplotypeCaller, FreeBayes, and DeepVariant), which have different underlying models. We provide detailed execution steps, as well as basic variant filtering, annotation, visualization, and consolidation aspects.

•

Protocol to enable whole exome data analysis in an unbiased approach

•

A protocol for unbiased analysis using 3 variant callers with different underlying models

•

From raw data to filtered, consolidated, and annotated DNA variant calls

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.