Evolutionary history of Tibetans inferred from whole-genome sequencing.

The indigenous people of the Tibetan Plateau have been the subject of much recent interest because of their unique genetic adaptations to high altitude. Recent studies have demonstrated that the Tibetan EPAS1 haplotype is involved in high altitude-adaptation and originated in an archaic Denisovan-related population. We sequenced the whole-genomes of 27 Tibetans and conducted analyses to infer a detailed history of demography and natural selection of this population. We detected evidence of population structure between the ancestral Han and Tibetan subpopulations as early as 44 to 58 thousand years ago, but with high rates of gene flow until approximately 9 thousand years ago. The CMS test ranked EPAS1 and EGLN1 as the top two positive selection candidates, and in addition identified PTGIS, VDR, and KCTD12 as new candidate genes. The advantageous Tibetan EPAS1 haplotype shared many variants with the Denisovan genome, with an ancient gene tree divergence between the Tibetan and Denisovan haplotypes of about 1 million years ago. With the exception of EPAS1, we observed no evidence of positive selection on Denisovan-like haplotypes.

Genomic characterization of pediatric B-lymphoblastic lymphoma and B-lymphoblastic leukemia using formalin-fixed tissues.

BACKGROUND: Recurrent genomic changes in B-lymphoblastic leukemia (B-ALL) identified by genome-wide single-nucleotide polymorphism (SNP) microarray analysis provide important prognostic information, but gene copy number analysis of its rare lymphoma counterpart, B-lymphoblastic lymphoma (B-LBL), is limited by the low incidence and lack of fresh tissue for genomic testing. PROCEDURE: We used molecular inversion probe (MIP) technology to analyze and compare copy number alterations (CNAs) in archival formalin-fixed paraffin-embedded pediatric B-LBL (n = 23) and B-ALL (n = 55). RESULTS: Similar to B-ALL, CDKN2A/B deletions were the most common alteration identified in 6/23 (26%) B-LBL cases. Eleven of 23 (48%) B-LBL patients were hyperdiploid, but none showed triple trisomies (chromosomes 4, 10, and 17) characteristic of B-ALL. IKZF1 and PAX5 deletions were observed in 13 and 17% of B-LBL, respectively, which was similar to the reported frequency in B-ALL. Immunoglobulin light chain lambda (IGL) locus deletions consistent with normal light chain rearrangement were observed in 5/23 (22%) B-LBL cases, compared with only 1% in B-ALL samples. None of the B-LBL cases showed abnormal, isolated VPREB1 deletion adjacent to IGL locus, which we identified in 25% of B-ALL. CONCLUSIONS: Our study demonstrates that the copy number profile of B-LBL is distinct from B-ALL, suggesting possible differences in pathogenesis between these closely related diseases.

Human spermatogenic failure purges deleterious mutation load from the autosomes and both sex chromosomes, including the gene DMRT1.

Gonadal failure, along with early pregnancy loss and perinatal death, may be an important filter that limits the propagation of harmful mutations in the human population. We hypothesized that men with spermatogenic impairment, a disease with unknown genetic architecture and a common cause of male infertility, are enriched for rare deleterious mutations compared to men with normal spermatogenesis. After assaying genomewide SNPs and CNVs in 323 Caucasian men with idiopathic spermatogenic impairment and more than 1,100 controls, we estimate that each rare autosomal deletion detected in our study multiplicatively changes a man's risk of disease by 10% (OR 1.10 [1.04-1.16], p<2 × 10(-3)), rare X-linked CNVs by 29%, (OR 1.29 [1.11-1.50], p<1 × 10(-3)), and rare Y-linked duplications by 88% (OR 1.88 [1.13-3.13], p<0.03). By contrasting the properties of our case-specific CNVs with those of CNV callsets from cases of autism, schizophrenia, bipolar disorder, and intellectual disability, we propose that the CNV burden in spermatogenic impairment is distinct from the burden of large, dominant mutations described for neurodevelopmental disorders. We identified two patients with deletions of DMRT1, a gene on chromosome 9p24.3 orthologous to the putative sex determination locus of the avian ZW chromosome system. In an independent sample of Han Chinese men, we identified 3 more DMRT1 deletions in 979 cases of idiopathic azoospermia and none in 1,734 controls, and found none in an additional 4,519 controls from public databases. The combined results indicate that DMRT1 loss-of-function mutations are a risk factor and potential genetic cause of human spermatogenic failure (frequency of 0.38% in 1306 cases and 0% in 7,754 controls, p = 6.2 × 10(-5)). Our study identifies other recurrent CNVs as potential causes of idiopathic azoospermia and generates hypotheses for directing future studies on the genetic basis of male infertility and IVF outcomes.

Optimizing single-cell RNA sequencing methods for human colon biopsies: droplet-based vs. picowell-based platforms.

Background & Aims: Single-cell RNA sequencing (scRNA) has empowered many insights into gastrointestinal microenvironments. However, profiling human biopsies using droplet-based scRNA (D-scRNA) is challenging since it requires immediate processing to minimize epithelial cell damage. In contrast, picowell-based (P-scRNA) platforms permit short-term frozen storage before sequencing. We compared P- and D-scRNA platforms on cells derived from human colon biopsies. Methods: Endoscopic rectosigmoid mucosal biopsies were obtained from two adults and conducted D-scRNA (10X Chromium) and P-scRNA (Honeycomb HIVE) in parallel using an individual's pool of single cells (> 10,000 cells/participant). Three experiments were performed to evaluate 1) P-scRNA with cells under specific storage conditions (immediately processed [fresh], vs. frozen at -20C vs. -80C [2 weeks]); 2) fresh P-scRNA versus fresh D-scRNA; and 3) P-scRNA stored at -80C with fresh D-scRNA. Results: Significant recovery of loaded cells was achieved for fresh (80.9%) and -80C (48.5%) P-scRNA and D-scRNA (76.6%), but not -20C P-scRNA (3.7%). However, D-scRNA captures more typeable cells among recovered cells (71.5% vs. 15.8% Fresh and 18.4% -80C P-scRNA), and these cells exhibit higher gene coverage at the expense of higher mitochondrial read fractions across most cell types. Cells profiled using D-scRNA demonstrated more consistent gene expression profiles among the same cell type than those profiled using P-scRNA. Significant intra-cell-type differences were observed in profiled gene classes across platforms. Conclusions: Our results highlight non-overlapping advantages of P-scRNA and D-scRNA and underscore the need for innovation to enable high-fidelity capture of colonic epithelial cells. The platform-specific variation highlights the challenges of maintaining rigor and reproducibility across studies that use different platforms.

Incident Cancer Risk and Signatures Among Older MUTYH Carriers: Analysis of Population-Based and Genomic Cohorts.

MUTYH carriers have an increased colorectal cancer risk in case-control studies, with loss of heterozygosity (LOH) as the presumed mechanism. We evaluated cancer risk among carriers in a prospective, population-based cohort of older adults. In addition, we assessed if cancers from carriers demonstrated mutational signatures (G:C>T:A transversions) associated with early LOH. We calculated incident risk of cancer and colorectal cancer among 13,131 sequenced study participants of the ASPirin in Reducing Events in the Elderly cohort, stratified by sex and adjusting for age, smoking, alcohol use, BMI, polyp history, history of cancer, and aspirin use. MUTYH carriers were identified among 13,033 participants in The Cancer Genome Atlas and International Cancer Genome Consortium, and somatic signatures of cancers were analyzed. Male MUTYH carriers demonstrated an increased risk for overall cancer incidence [multivariable HR, 1.66; 95% confidence interval (CI), 1.03-2.68; P = 0.038] driven by increased colorectal cancer incidence (multivariable HR, 3.55; 95% CI, 1.42-8.78; P = 0.007), as opposed to extracolonic cancer incidence (multivariable HR, 1.40; 95% CI, 0.81-2.44; P = 0.229). Female carriers did not demonstrate increased risk of cancer, colorectal cancer, or extracolonic cancers. Analysis of mutation signatures from cancers of MUTYH carriers revealed no significant contribution toward early mutagenesis from widespread G:C>T:A transversions among gastrointestinal epithelial cancers. Among cancers from carriers, somatic transversions associated with base-excision repair deficiency are uncommon, suggestive of diverse mechanisms of carcinogenesis in carriers compared with those who inherit biallelic MUTYH mutations. PREVENTION RELEVANCE: Despite absence of loss of heterozygosity in colorectal cancers, elderly male MUTYH carriers appeared to be at increased of colorectal cancer.

Genome wide copy number analysis of paediatric Burkitt lymphoma using formalin-fixed tissues reveals a subset with gain of chromosome 13q and corresponding miRNA over expression.

The majority of paediatric Burkitt lymphoma (pBL) patients that relapse will die of disease, but markers for this high-risk subset are unknown. MYC translocations characterize pBL, but additional genetic changes may relate to prognosis and serve as potential biomarkers. We utilized a molecular inversion probe single nucleotide polymorphism assay to perform high resolution, genome-wide copy number analysis on archival formalin-fixed, paraffin-embedded pBL and germline tissues. We identified copy number abnormalities (CNAs) in 18/28 patients (64%) with a total of 62 CNAs that included 32 gains and 30 copy number losses. We identified seven recurrent CNAs including 1q gain (7/28, 25%), 13q gain (3/28, 11%), and 17p loss (4/28, 14%). The minimum common amplified region on 13q was at 13q31 and included the MIR17HG (MIR17-92) locus. Samples with this gain had higher levels of MIR17 RNA and showed a tendency for early relapse. Tumour-specific uniparental disomy was identified in 32% of cases and usually was recurrent. These results demonstrate that high-resolution copy number analysis can be performed on archival lymphoma tissue specimens, which has significance for the study of rare diseases.

Pathogenic effect of TP73 Gene Variants in People With Amyotrophic Lateral Sclerosis.

OBJECTIVE: To identify novel disease associated loci for amyotrophic lateral sclerosis (ALS), we utilized sequencing data and performed in vitro and in vivo experiments to demonstrate pathogenicity of mutations identified in TP73. METHODS: We analyzed exome sequences of 87 sporadic ALS patients and 324 controls, with confirmatory sequencing in independent ALS cohorts of >2,800 patients. For the top hit, TP73, a regulator of apoptosis, differentiation, and a binding partner as well as homolog of the tumor suppressor gene TP53, we assayed mutation effects using in vitro and in vivo experiments. C2C12 myoblast differentiation assays, characterization of myotube appearance, and immunoprecipitation of p53-p73 complexes were perform in vitro. In vivo, we used CRISPR/Cas9 targeting of zebrafish tp73 to assay motor neuron number and axon morphology. RESULTS: Five heterozygous rare, nonsynonymous mutations in TP73 were identified in our sporadic ALS cohort. In independent ALS cohorts, we identified an additional 19 rare, deleterious variants in TP73. Patient TP73 mutations caused abnormal differentiation and increased apoptosis in the myoblast differentiation assay, with abnormal myotube appearance. Immunoprecipitation of mutant ΔN-p73 demonstrated that patient mutations hinder ΔN-p73's ability to bind p53. CRISPR/Cas9 knockout of tp73 in zebrafish led to impaired motor neuron development and abnormal axonal morphology, concordant with ALS pathology. CONCLUSION: Together, these results strongly suggest that variants in TP73 correlate with risk for ALS and indicate a novel role for apoptosis in ALS disease pathology.

Association of Combined Focal 22q11.22 Deletion and IKZF1 Alterations With Outcomes in Childhood Acute Lymphoblastic Leukemia.

IMPORTANCE: Alterations in the IKZF1 gene drive B-cell acute lymphoblastic leukemia (B-ALL) but are not routinely used to stratify patients by risk because of inconsistent associations with outcomes. We describe a novel deletion in 22q11.22 that was consistently associated with very poor outcomes in patients with B-ALL with IKZF1 alterations. OBJECTIVE: To determine whether focal deletions within the λ variable chain region in chromosome 22q11.22 were associated with patients with B-ALL with IKZF1 alterations with the highest risk of relapse and/or death. DESIGN, SETTING, AND PARTICIPANTS: This cohort study included 1310 primarily high-risk pediatric patients with B-ALL who were taken from 6 independent clinical cohorts, consisting of 3 multicenter cohorts (AALL0232 [2004-2011], P9906 [2000-2003], and patients with Down syndrome who were pooled from national and international studies) and 3 single-institution cohorts (University of Utah [Salt Lake City], Children's Hospital of Philadelphia [Philadelphia, Pennsylvania], and St. Jude Children's Hospital [Memphis, Tennessee]). Data analysis began in 2011 using patients from the older studies first, and data analysis concluded in 2021. EXPOSURES: Focal 22q11.22 deletions. MAIN OUTCOMES AND MEASURES: Event-free and overall survival was investigated. The hypothesis that 22q11.22 deletions stratified the prognostic effect of IKZF1 alterations was formulated while investigating nearby deletions in VPREB1 in 2 initial cohorts (n = 270). Four additional cohorts were then obtained to further study this association (n = 1040). RESULTS: This study of 1310 patients with B-ALL (717 male [56.1%] and 562 female patients [43.9%]) found that focal 22q11.22 deletions are frequent (518 of 1310 [39.5%]) in B-ALL and inconsistent with physiologic V(D)J recombination. A total of 299 of 1310 patients with B-ALL had IKZF1 alterations. Among patients with IKZF1 alterations, more than half shared concomitant focal 22q11.22 deletions (159 of 299 [53.0%]). Patients with combined IKZF1 alterations and 22q11.22 deletions had worse outcomes compared with patients with IKZF1 alterations and wild-type 22q11.22 alleles in every cohort examined (combined cohorts: 5-year event-free survival rates, 43.3% vs 68.5%; hazard ratio [HR], 2.18; 95% CI, 1.54-3.07; P < .001; 5-year overall survival rates, 66.9% vs 83.9%; HR, 2.05; 95% CI, 1.32-3.21; P = .001). While 22q11.22 deletions were not prognostic in patients with wild-type IKZF1 , concomitant 22q11.22 deletions in patients with IKZF1 alterations stratified outcomes across additional risk groups, including patients who met the IKZF1plus criteria, and maintained independent significance in multivariate analysis for event-free survival (HR, 2.05; 95% CI, 1.27-3.29; P = .003) and overall survival (HR, 1.83; 95% CI, 1.01-3.34; P = .05). CONCLUSIONS AND RELEVANCE: This cohort study suggests that 22q11.22 deletions identify patients with B-ALL and IKZF1 alterations who have very poor outcomes and may offer a new genetic biomarker to further refine B-ALL risk stratification and treatment strategies.

The evolving genetic risk for sporadic ALS.

OBJECTIVE: To estimate the genetic risk conferred by known amyotrophic lateral sclerosis (ALS)-associated genes to the pathogenesis of sporadic ALS (SALS) using variant allele frequencies combined with predicted variant pathogenicity. METHODS: Whole exome sequencing and repeat expansion PCR of C9orf72 and ATXN2 were performed on 87 patients of European ancestry with SALS seen at the University of Utah. DNA variants that change the protein coding sequence of 31 ALS-associated genes were annotated to determine which were rare and deleterious as predicted by MetaSVM. The percentage of patients with SALS with a rare and deleterious variant or repeat expansion in an ALS-associated gene was calculated. An odds ratio analysis was performed comparing the burden of ALS-associated genes in patients with SALS vs 324 normal controls. RESULTS: Nineteen rare nonsynonymous variants in an ALS-associated gene, 2 of which were found in 2 different individuals, were identified in 21 patients with SALS. Further, 5 deleterious C9orf72 and 2 ATXN2 repeat expansions were identified. A total of 17.2% of patients with SALS had a rare and deleterious variant or repeat expansion in an ALS-associated gene. The genetic burden of ALS-associated genes in patients with SALS as predicted by MetaSVM was significantly higher than in normal controls. CONCLUSIONS: Previous analyses have identified SALS-predisposing variants only in terms of their rarity in normal control populations. By incorporating variant pathogenicity as well as variant frequency, we demonstrated that the genetic risk contributed by these genes for SALS is substantially lower than previous estimates.

A Genome-Wide Search for Greek and Jewish Admixture in the Kashmiri Population.

The Kashmiri population is an ethno-linguistic group that resides in the Kashmir Valley in northern India. A longstanding hypothesis is that this population derives ancestry from Jewish and/or Greek sources. There is historical and archaeological evidence of ancient Greek presence in India and Kashmir. Further, some historical accounts suggest ancient Hebrew ancestry as well. To date, it has not been determined whether signatures of Greek or Jewish admixture can be detected in the Kashmiri population. Using genome-wide genotyping and admixture detection methods, we determined there are no significant or substantial signs of Greek or Jewish admixture in modern-day Kashmiris. The ancestry of Kashmiri Tibetans was also determined, which showed signs of admixture with populations from northern India and west Eurasia. These results contribute to our understanding of the existing population structure in northern India and its surrounding geographical areas.