Novel insight into the genetic basis of high-altitude pulmonary hypertension in Kyrgyz highlanders.

The Central Asian Kyrgyz highland population provides a unique opportunity to address genetic diversity and understand the genetic mechanisms underlying high-altitude pulmonary hypertension (HAPH). Although a significant fraction of the population is unaffected, there are susceptible individuals who display HAPH in the absence of any lung, cardiac or hematologic disease. We report herein the analysis of the whole-genome sequencing of healthy individuals compared with HAPH patients and other controls (total n = 33). Genome scans reveal selection signals in various regions, encompassing multiple genes from the first whole-genome sequences focusing on HAPH. We show here evidence of three candidate genes MTMR4, TMOD3 and VCAM1 that are functionally associated with well-known molecular and pathophysiological processes and which likely lead to HAPH in this population. These processes are (a) dysfunctional BMP signaling, (b) disrupted tissue repair processes and (c) abnormal endothelial cell function. Whole-genome sequence of well-characterized patients and controls and using multiple statistical tools uncovered novel candidate genes that belong to pathways central to the pathogenesis of HAPH. These studies on high-altitude human populations are pertinent to the understanding of sea level diseases involving hypoxia as a main element of their pathophysiology.

The human noncoding genome defined by genetic diversity.

Understanding the significance of genetic variants in the noncoding genome is emerging as the next challenge in human genomics. We used the power of 11,257 whole-genome sequences and 16,384 heptamers (7-nt motifs) to build a map of sequence constraint for the human species. This build differed substantially from traditional maps of interspecies conservation and identified regulatory elements among the most constrained regions of the genome. Using new Hi-C experimental data, we describe a strong pattern of coordination over 2 Mb where the most constrained regulatory elements associate with the most essential genes. Constrained regions of the noncoding genome are up to 52-fold enriched for known pathogenic variants as compared to unconstrained regions (21-fold when compared to the genome average). This map of sequence constraint across thousands of individuals is an asset to help interpret noncoding elements in the human genome, prioritize variants and reconsider gene units at a larger scale.

Genetic risk, dysbiosis, and treatment stratification using host genome and gut microbiome in inflammatory bowel disease.

OBJECTIVES: Inflammatory bowel diseases (IBD), comprised of Crohn's disease (CD) and ulcerative colitis (UC), are characterized by a complex pathophysiology that is thought to result from an aberrant immune response to a dysbiotic luminal microbiota in genetically susceptible individuals. New technologies support the joint assessment of host-microbiome interaction. METHODS: Using whole genome sequencing and shotgun metagenomics, we studied the clinical features, host genome, and stool microbial metagenome of 85 IBD patients, and compared the results to 146 control individuals. Genetic risk scores, computed on 159 single nucleotide variants, and human leukocyte antigen (HLA) types differentiated IBD patients from healthy controls. RESULTS: Genetic risk was associated with the need for use of biologics in IBD and, modestly, with the composition of the gut microbiome. As compared with healthy controls, IBD patients had hallmarks of stool microbiome dysbiosis, with loss of a diversified core microbiome, enrichment and depletion of specific bacteria, and enrichment of bacterial virulence factors. CONCLUSIONS: We show that genetic risk may have a role in early risk stratification in the care of IBD patients and propose that expression of virulence factors in a dysbiotic microbiome may contribute to pathogenesis in IBD.

New Insights into the Genetic Basis of Monge's Disease and Adaptation to High-Altitude.

Human high-altitude (HA) adaptation or mal-adaptation is explored to understand the physiology, pathophysiology, and molecular mechanisms that underlie long-term exposure to hypoxia. Here, we report the results of an analysis of the largest whole-genome-sequencing of Chronic Mountain Sickness (CMS) and nonCMS individuals, identified candidate genes and functionally validated these candidates in a genetic model system (Drosophila). We used PreCIOSS algorithm that uses Haplotype Allele Frequency score to separate haplotypes carrying the favored allele from the noncarriers and accordingly, prioritize genes associated with the CMS or nonCMS phenotype. Haplotypes in eleven candidate regions, with SNPs mostly in nonexonic regions, were significantly different between CMS and nonCMS subjects. Closer examination of individual genes in these regions revealed the involvement of previously identified candidates (e.g., SENP1) and also unreported ones SGK3, COPS5, PRDM1, and IFT122 in CMS. Remarkably, in addition to genes like SENP1, SGK3, and COPS5 which are HIF-dependent, our study reveals for the first time HIF-independent gene PRDM1, indicating an involvement of wider, nonHIF pathways in HA adaptation. Finally, we observed that down-regulating orthologs of these genes in Drosophila significantly enhanced their hypoxia tolerance. Taken together, the PreCIOSS algorithm, applied on a large number of genomes, identifies the involvement of both new and previously reported genes in selection sweeps, highlighting the involvement of multiple hypoxia response systems. Since the overwhelming majority of SNPs are in nonexonic (and possibly regulatory) regions, we speculate that adaptation to HA necessitates greater genetic flexibility allowing for transcript variability in response to graded levels of hypoxia.

Fine population structure analysis method for genomes of many.

Fine population structure can be examined through the clustering of individuals into subpopulations. The clustering of individuals in large sequence datasets into subpopulations makes the calculation of subpopulation specific allele frequency possible, which may shed light on selection of candidate variants for rare diseases. However, as the magnitude of the data increases, computational burden becomes a challenge in fine population structure analysis. To address this issue, we propose fine population structure analysis (FIPSA), which is an individual-based non-parametric method for dissecting fine population structure. FIPSA maximizes the likelihood ratio of the contingency table of the allele counts multiplied by the group. We demonstrated that its speed and accuracy were superior to existing non-parametric methods when the simulated sample size was up to 5,000 individuals. When applied to real data, the method showed high resolution on the Human Genome Diversity Project (HGDP) East Asian dataset. FIPSA was independently validated on 11,257 human genomes. The group assignment given by FIPSA was 99.1% similar to those assigned based on supervised learning. Thus, FIPSA provides high resolution and is compatible with a real dataset of more than ten thousand individuals.

Reconstructing genetic history of Siberian and Northeastern European populations.

Siberia and Northwestern Russia are home to over 40 culturally and linguistically diverse indigenous ethnic groups, yet genetic variation and histories of peoples from this region are largely uncharacterized. We present deep whole-genome sequencing data (∼38×) from 28 individuals belonging to 14 distinct indigenous populations from that region. We combined these data sets with additional 32 modern-day and 46 ancient human genomes to reconstruct genetic histories of several indigenous Northern Eurasian populations. We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ust'-Ishim individual who was previously believed to have no modern-day descendants. Western Siberians trace 57% of their ancestry to ancient North Eurasians, represented by the 24,000-yr-old Siberian Mal'ta boy MA-1. Eastern Siberian populations formed a distinct sublineage that separated from other East Asian populations ∼10,000 yr ago. In addition, we uncovered admixtures between Siberians and Eastern European hunter-gatherers from Samara, Karelia, Hungary, and Sweden (from 8000-6600 yr ago); Yamnaya people (5300-4700 yr ago); and modern-day Northeastern Europeans. Our results provide new insights into genetic histories of Siberian and Northeastern European populations and evidence of ancient gene flow from Siberia into Europe.

Deep sequencing of 10,000 human genomes.

We report on the sequencing of 10,545 human genomes at 30×-40× coverage with an emphasis on quality metrics and novel variant and sequence discovery. We find that 84% of an individual human genome can be sequenced confidently. This high-confidence region includes 91.5% of exon sequence and 95.2% of known pathogenic variant positions. We present the distribution of over 150 million single-nucleotide variants in the coding and noncoding genome. Each newly sequenced genome contributes an average of 8,579 novel variants. In addition, each genome carries on average 0.7 Mb of sequence that is not found in the main build of the hg38 reference genome. The density of this catalog of variation allowed us to construct high-resolution profiles that define genomic sites that are highly intolerant of genetic variation. These results indicate that the data generated by deep genome sequencing is of the quality necessary for clinical use.

Chinese family with diffuse oesophageal leiomyomatosis: a new COL4A5/COL4A6 deletion and a case of gonosomal mosaicism.

BACKGROUND: Diffuse oesophageal leiomyomatosis (DOL) is a rare disorder characterized by tumorous overgrowth of the muscular wall of the oesophagus. DOL is present in 5 % of Alport syndrome (AS) patients. AS is a rare hereditary disease that involves varying degrees of hearing impairment, ocular changes and progressive glomerulonephritis leading to renal failure. In DOL-AS patients, the genetic defect consists of a deletion involving the COL4A5 and COL4A6 genes on the X chromosome. CASE PRESENTATION: We report a two-generation family (4 individuals; parents and two children, one male and one female) with two members (mother and son) affected with oesophageal leiomyomatosis. Signs of potential renal failure, which characterizes AS, were only apparent in the index patient (son) 2 years and three months after the initial diagnosis of DOL. Blood DNA from the four family members were submitted to exome sequencing and array genotyping to perform a genome wide screening for disease causal single nucleotide (SN) and copy number (CN) variations. Analyses revealed a new 40kb deletion encompassing from intron 2 of COL4A5 to intron 1 of COL4A6 at Xq22.3. The breakpoints were also identified. Possible confounding pathogenic exonic variants in genes known to be involved in other extracellular matrices disorders were also shared by the two affected individuals. Meticulous analysis of the maternal DNA revealed a case of gonosomal mosaicism. CONCLUSIONS: This is the first report of gonadosomal mosaicism associated to DOL-AS.

Common variants on Xq28 conferring risk of schizophrenia in Han Chinese.

Schizophrenia is a highly heritable, severe psychiatric disorder affecting approximately 1% of the world population. A substantial portion of heritability is still unexplained and the pathophysiology of schizophrenia remains to be elucidated. To identify more schizophrenia susceptibility loci, we performed a genome-wide association study (GWAS) on 498 patients with schizophrenia and 2025 controls from the Han Chinese population, and a follow-up study on 1027 cases and 1005 controls. In the follow-up study, we included 384 single nucleotide polymorphisms (SNPs) which were selected from the top hits in our GWAS (130 SNPs) and from previously implicated loci for schizophrenia based on the SZGene database, NHGRI GWAS Catalog, copy number variation studies, GWAS meta-analysis results from the international Psychiatric Genomics Consortium (PGC) and candidate genes from plausible biological pathways (254 SNPs). Within the chromosomal region Xq28, SNP rs2269372 in RENBP achieved genome-wide significance with a combined P value of 3.98 × 10(-8) (OR of allele A = 1.31). SNPs with suggestive P values were identified within 2 genes that have been previously implicated in schizophrenia, MECP2 (rs2734647, P combined = 8.78 × 10(-7), OR = 1.28; rs2239464, P combined = 6.71 × 10(-6), OR = 1.26) and ARHGAP4 (rs2269368, P combined = 4.74 × 10(-7), OR = 1.25). In addition, the patient sample in our follow-up study showed a significantly greater burden for pre-defined risk alleles based on the SNPs selected than the controls. This indicates the existence of schizophrenia susceptibility loci among the SNPs we selected. This also further supports multigenic inheritance in schizophrenia. Our findings identified a new schizophrenia susceptibility locus on Xq28, which harbor the genes RENBP, MECP2, and ARHGAP4.

Gene network analysis of candidate loci for human anorectal malformations.

Anorectal malformations (ARMs) are birth defects that require surgery and carry significant chronic morbidity. Our earlier genome-wide copy number variation (CNV) study had provided a wealth of candidate loci. To find out whether these candidate loci are related to important developmental pathways, we have performed an extensive literature search coupled with the currently available bioinformatics tools. This has allowed us to assign both genic and non-genic CNVs to interrelated pathways known to govern the development of the anorectal region. We have linked 11 candidate genes to the WNT signalling pathway and 17 genes to the cytoskeletal network. Interestingly, candidate genes with similar functions are disrupted by the same type of CNV. The gene network we discovered provides evidence that rare mutations in different interrelated genes may lead to similar phenotypes, accounting for genetic heterogeneity in ARMs. Classification of patients according to the affected pathway and lesion type should eventually improve the diagnosis and the identification of common genes/molecules as therapeutic targets.

Genome-wide copy number variation study in anorectal malformations.

Anorectal malformations (ARMs, congenital obstruction of the anal opening) are among the most common birth defects requiring surgical treatment (2-5/10 000 live-births) and carry significant chronic morbidity. ARMs present either as isolated or as part of the phenotypic spectrum of some chromosomal abnormalities or monogenic syndromes. The etiology is unknown. To assess the genetic contribution to ARMs, we investigated single-nucleotide polymorphisms and copy number variations (CNVs) at genome-wide scale. A total of 363 Han Chinese sporadic ARM patients and 4006 Han Chinese controls were included. Overall, we detected a 1.3-fold significant excess of rare CNVs in patients. Stratification of patients by presence/absence of other congenital anomalies showed that while syndromic ARM patients carried significantly longer rare duplications than controls (P = 0.049), non-syndromic patients were enriched with both rare deletions and duplications when compared with controls (P = 0.00031). Twelve chromosomal aberrations and 114 rare CNVs were observed in patients but not in 868 controls nor 11 943 healthy individuals from the Database of Genomic Variants. Importantly, these aberrations were observed in isolated ARM patients. Gene-based analysis revealed 79 genes interfered by CNVs in patients only. In particular, we identified a de novo DKK4 duplication. DKK4 is a member of the WNT signaling pathway which is involved in the development of the anorectal region. In mice, Wnt disruption results in ARMs. Our data suggest a role for rare CNVs not only in syndromic but also in isolated ARM patients and provide a list of plausible candidate genes for the disorder.

Codon usage bias and the evolution of influenza A viruses. Codon Usage Biases of Influenza Virus.

BACKGROUND: The influenza A virus is an important infectious cause of morbidity and mortality in humans and was responsible for 3 pandemics in the 20th century. As the replication of the influenza virus is based on its host's machinery, codon usage of its viral genes might be subject to host selection pressures, especially after interspecies transmission. A better understanding of viral evolution and host adaptive responses might help control this disease. RESULTS: Relative Synonymous Codon Usage (RSCU) values of the genes from segment 1 to segment 6 of avian and human influenza viruses, including pandemic H1N1, were studied via Correspondence Analysis (CA). The codon usage patterns of seasonal human influenza viruses were distinct among their subtypes and different from those of avian viruses. Newly isolated viruses could be added to the CA results, creating a tool to investigate the host origin and evolution of viral genes. It was found that the 1918 pandemic H1N1 virus contained genes with mammalian-like viral codon usage patterns, indicating that the introduction of this virus to humans was not through in toto transfer of an avian influenza virus.Many human viral genes had directional changes in codon usage over time of viral isolation, indicating the effect of host selection pressures. These changes reduced the overall GC content and the usage of G at the third codon position in the viral genome. Limited evidence of translational selection pressure was found in a few viral genes. CONCLUSIONS: Codon usage patterns from CA allowed identification of host origin and evolutionary trends in influenza viruses, providing an alternative method and a tool to understand the evolution of influenza viruses. Human influenza viruses are subject to selection pressure on codon usage which might assist in understanding the characteristics of newly emerging viruses.