Increasing diversity of functional genetics studies to advance biological discovery and human health.

In this perspective we discuss the current lack of genetic and environmental diversity in functional genomics datasets. There is a well-described Eurocentric bias in genetic and functional genomic research that has a clear impact on the benefit this research can bring to underrepresented populations. Current research focused on genetic variant-to-function experiments aims to identify molecular QTLs, but the lack of data from genetically diverse individuals has limited analyses to mostly populations of European ancestry. Although some efforts have been established to increase diversity in functional genomic studies, much remains to be done to consistently generate data for underrepresented populations from now on. We discuss the major barriers for this continuity and suggest actionable insights, aiming to empower research and researchers from underserved populations.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries.

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations.

Denisovan introgression has shaped the immune system of present-day Papuans.

Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.

The p.C759F Variant in USH2A Is a Pathogenic Mutation: Systematic Literature Review and Meta-Analysis of 667 Genotypes.

BACKGROUND: Although the p.C759F (c.2276G>T, p.Cys759Phe) variant in the USH2A gene has been identified in association with retinal degeneration by several authors, its pathogenicity has been questioned once by the publication of two unaffected homozygotes from a single family. OBJECTIVES: The objective of the study was to ascertain the role of p.C759F in hereditary retinal disease. METHODS: We examined 87 research articles reporting on patients carrying this variant and then used this information as primary data for a series of meta-analytical tests. RESULTS: Independent statistical analyses showed that p.C759F (i) is highly enriched in patients with respect to healthy individuals, (ii) represents a clear-cut recessive allele causing disease when it is in trans with other mutations, (iii) is pathogenic in homozygotes. CONCLUSIONS: Our results confirm that p.C759F is a bona fide mutation, leading to retinal blindness according to a recessive pattern of inheritance.

PRPH2-Related Retinal Dystrophies: Mutational Spectrum in 103 Families from a Spanish Cohort.

PRPH2, one of the most frequently inherited retinal dystrophy (IRD)-causing genes, implies a high phenotypic variability. This study aims to analyze the PRPH2 mutational spectrum in one of the largest cohorts worldwide, and to describe novel pathogenic variants and genotype-phenotype correlations. A study of 220 patients from 103 families recruited from a database of 5000 families. A molecular diagnosis was performed using classical molecular approaches and next-generation sequencing. Common haplotypes were ascertained by analyzing single-nucleotide polymorphisms. We identified 56 variants, including 11 novel variants. Most of them were missense variants (64%) and were located in the D2-loop protein domain (77%). The most frequently occurring variants were p.Gly167Ser, p.Gly208Asp and p.Pro221_Cys222del. Haplotype analysis revealed a shared region in families carrying p.Leu41Pro or p.Pro221_Cys222del. Patients with retinitis pigmentosa presented an earlier disease onset. We describe the largest cohort of IRD families associated with PRPH2 from a single center. Most variants were located in the D2-loop domain, highlighting its importance in interacting with other proteins. Our work suggests a likely founder effect for the variants p.Leu41Pro and p.Pro221_Cys222del in our Spanish cohort. Phenotypes with a primary rod alteration presented more severe affectation. Finally, the high phenotypic variability in PRPH2 hinders the possibility of drawing genotype-phenotype correlations.

Hypogammaglobulinemia after front-line immunochemotherapy in patients with follicular lymphoma.

INTRODUCTION: Hypogammaglobulinemia after front-line immunochemotherapy for follicular lymphoma is a poorly studied adverse event that could be related to the appearance of severe and/or recurrent non-neutropenic infections which could affect the quality of life of the patients, even motivating a need of long-term replacement therapy with human immunoglobulins. METHODS: Observational, retrospective study aiming to estimate the incidence of hypogammaglobulinemia, as well as its severity and clinical consequences, and to explore possible predictive factors for its development. Specific immunoglobulin deficiencies were also studied. RESULTS: 76.5% of patients had hypogammaglobulinemia during or after front-line treatment, mostly grade 1-2; with 38.8% patients who developed clinically relevant infections and 20% patients requiring human immunoglobulins replacement therapy. A high-risk FLIPI score was identified as a risk factor for hypogammaglobulinemia (ods ratio: 4.51; 95% confidence interval: 1.29-15.68; p < 0.001) and basal gamma globulin level as a protective factor (odds ratio: 0.92; 95% confidence interval: 0.988-0.996; p = 0.018). Any type of immunochemotherapy regimen was associated with different risks of hypogammaglobulinemia in our study. CONCLUSIONS: Hypogammaglobulinemia appears in a high percentage of patients with follicular lymphoma in a real-world population, identifying a high-risk FLIPI score as a risk factor for its development and basal gamma globulins as a protective factor.

Genome-wide DNA methylation and gene expression patterns reflect genetic ancestry and environmental differences across the Indonesian archipelago.

Indonesia is the world's fourth most populous country, host to striking levels of human diversity, regional patterns of admixture, and varying degrees of introgression from both Neanderthals and Denisovans. However, it has been largely excluded from the human genomics sequencing boom of the last decade. To serve as a benchmark dataset of molecular phenotypes across the region, we generated genome-wide CpG methylation and gene expression measurements in over 100 individuals from three locations that capture the major genomic and geographical axes of diversity across the Indonesian archipelago. Investigating between- and within-island differences, we find up to 10.55% of tested genes are differentially expressed between the islands of Sumba and New Guinea. Variation in gene expression is closely associated with DNA methylation, with expression levels of 9.80% of genes correlating with nearby promoter CpG methylation, and many of these genes being differentially expressed between islands. Genes identified in our differential expression and methylation analyses are enriched in pathways involved in immunity, highlighting Indonesia's tropical role as a source of infectious disease diversity and the strong selective pressures these diseases have exerted on humans. Finally, we identify robust within-island variation in DNA methylation and gene expression, likely driven by fine-scale environmental differences across sampling sites. Together, these results strongly suggest complex relationships between DNA methylation, transcription, archaic hominin introgression and immunity, all jointly shaped by the environment. This has implications for the application of genomic medicine, both in critically understudied Indonesia and globally, and will allow a better understanding of the interacting roles of genomic and environmental factors shaping molecular and complex phenotypes.

Prioritization of New Candidate Genes for Rare Genetic Diseases by a Disease-Aware Evaluation of Heterogeneous Molecular Networks.

Screening for pathogenic variants in the diagnosis of rare genetic diseases can now be performed on all genes thanks to the application of whole exome and genome sequencing (WES, WGS). Yet the repertoire of gene-disease associations is not complete. Several computer-based algorithms and databases integrate distinct gene-gene functional networks to accelerate the discovery of gene-disease associations. We hypothesize that the ability of every type of information to extract relevant insights is disease-dependent. We compiled 33 functional networks classified into 13 knowledge categories (KCs) and observed large variability in their ability to recover genes associated with 91 genetic diseases, as measured using efficiency and exclusivity. We developed GLOWgenes, a network-based algorithm that applies random walk with restart to evaluate KCs' ability to recover genes from a given list associated with a phenotype and modulates the prediction of new candidates accordingly. Comparison with other integration strategies and tools shows that our disease-aware approach can boost the discovery of new gene-disease associations, especially for the less obvious ones. KC contribution also varies if obtained using recently discovered genes. Applied to 15 unsolved WES, GLOWgenes proposed three new genes to be involved in the phenotypes of patients with syndromic inherited retinal dystrophies.

Papua New Guinean Genomes Reveal the Complex Settlement of North Sahul.

The settlement of Sahul, the lost continent of Oceania, remains one of the most ancient and debated human migrations. Modern New Guineans inherited a unique genetic diversity tracing back 50,000 years, and yet there is currently no model reconstructing their past population dynamics. We generated 58 new whole-genome sequences from Papua New Guinea, filling geographical gaps in previous sampling, specifically to address alternative scenarios of the initial migration to Sahul and the settlement of New Guinea. Here, we present the first genomic models for the settlement of northeast Sahul considering one or two migrations from Wallacea. Both models fit our data set, reinforcing the idea that ancestral groups to New Guinean and Indigenous Australians split early, potentially during their migration in Wallacea where the northern route could have been favored. The earliest period of human presence in Sahul was an era of interactions and gene flow between related but already differentiated groups, from whom all modern New Guineans, Bismarck islanders, and Indigenous Australians descend. The settlement of New Guinea was probably initiated from its southeast region, where the oldest archaeological sites have been found. This was followed by two migrations into the south and north lowlands that ultimately reached the west and east highlands. We also identify ancient gene flows between populations in New Guinea, Australia, East Indonesia, and the Bismarck Archipelago, emphasizing the fact that the anthropological landscape during the early period of Sahul settlement was highly dynamic rather than the traditional view of extensive isolation.

The Simons Genome Diversity Project: 300 genomes from 142 diverse populations.

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans.

Short-Term Gaseous Treatments Improve Rachis Browning in Red and White Table Grapes Stored at Low Temperature: Molecular Mechanisms Underlying Its Beneficial Effect.

Short-term gaseous treatments improve rachis quality during table grape postharvest, but little is known about the mechanisms involved. In this work, we observed that the application of a 3-day CO2 treatment at 0 °C improved rachis browning of Superior Seedless and Red Globe bunches, affecting the non-enzymatic antioxidant system by reducing the total phenolic content, the antioxidant activity and the expression of different stilbene synthase genes. Lipid peroxidation levels revealed lower oxidative stress in CO2-treated rachis of both cultivars linked to the activation of the enzymatic antioxidant system. Furthermore, whereas a positive correlation was denoted between rachis browning and the accumulation of key ABA regulatory genes in Red Globe bunches, this effect was restricted to ACS1, a key synthetic ethylene gene, in Superior Seedless clusters. This work also corroborated the important role of ethylene-responsive factors in the beneficial effect of the gaseous treatment, not only in the berries but also in the rachis. Finally, the application of the gaseous treatment avoided the induction of cell wall-degrading enzyme-related genes in both cultivars, which could favor the maintenance of rachis quality. This work provides new insight into specific responses modulated by the gaseous treatment focused on mitigating rachis browning independently of the cultivar.

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies.

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort's IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

Impact of regulatory variation across human iPSCs and differentiated cells.

Induced pluripotent stem cells (iPSCs) are an essential tool for studying cellular differentiation and cell types that are otherwise difficult to access. We investigated the use of iPSCs and iPSC-derived cells to study the impact of genetic variation on gene regulation across different cell types and as models for studies of complex disease. To do so, we established a panel of iPSCs from 58 well-studied Yoruba lymphoblastoid cell lines (LCLs); 14 of these lines were further differentiated into cardiomyocytes. We characterized regulatory variation across individuals and cell types by measuring gene expression levels, chromatin accessibility, and DNA methylation. Our analysis focused on a comparison of inter-individual regulatory variation across cell types. While most cell-type-specific regulatory quantitative trait loci (QTLs) lie in chromatin that is open only in the affected cell types, we found that 20% of cell-type-specific regulatory QTLs are in shared open chromatin. This observation motivated us to develop a deep neural network to predict open chromatin regions from DNA sequence alone. Using this approach, we were able to use the sequences of segregating haplotypes to predict the effects of common SNPs on cell-type-specific chromatin accessibility.

Apparent but unconfirmed digenism in an Iranian consanguineous family with syndromic Retinal Disease.

BACKGROUND: Bardet-Biedl syndrome is an autosomal recessive disease characterized by rod-cone dystrophy, postaxial polydactyly, kidney defects, obesity, mental retardation and hypogonadism. Here, we report different genotypes in two Bardet-Biedl syndrome affected sisters with a different clinical phenotype regarding severity. MATERIALS AND METHODS: The proband of the family was examined by Next Generation Sequencing (NGS) using clinical exome and filtering by syndromic and non-syndromic genes associated with retinal dystrophies. RESULTS: Targeted NGS revealed two novel variants in the MKKS and CEP290 genes in homozygosis state in the proband. Segregation analysis revealed the presence of the same MKKS homozygous variant in her younger affected sister but not the CEP290 variant. Both sisters presented different clinical manifestation, at different ages, with a more severe renal and retinal defect in the case of the sister carrying mutations in both genes. Another unaffected sister showed only homozygosity for the CEP290 variant, thus supporting the non-pathogenic role of this mutation in BBS phenotype. CONCLUSIONS: In this study, NGS proved to be a powerful and efficient sequencing method to identify causal variants in different genes. However, it remarks the importance of the segregation analysis and clinical information to establish the pathogenicity of new variants. The two affected sisters present different genotypes and clinical manifestation, suggesting that the novel CEP290 variant could be acting as a modifier, making the phenotype more severe in the sister homozygote for MKKS and CEP290 genes. On the other hand, the difference in the age of both sisters highlight the important role of monitoring disease progression also to confirm the modifier role of genetic variants.

RPE65-related retinal dystrophy: Mutational and phenotypic spectrum in 45 affected patients.

INTRODUCTION: Biallelic pathogenic RPE65 variants are related to a spectrum of clinically overlapping inherited retinal dystrophies (IRD). Most affected individuals progress to severe disease, with 50% of patients becoming legally blind by 20 years of age. Deeper knowledge of the mutational spectrum and the phenotype-genotype correlation in RPE65-related IRD is needed. PATIENTS AND METHODS: Forty-five affected subjects from 27 unrelated families with a clinical diagnosis of RPE65-related IRD were included. Clinical evaluation consisted of self-reported ophthalmological history and objective ophthalmological examination. Patients' genotype was classified according to variant class (truncating or missense) or to variant location at different protein domains. The main phenotypic outcome measure was age at onset (AAO) of symptomatic disease and a Kaplan-Meier analysis of disease symptom event-free survival was performed. RESULTS: Twenty-nine different RPE65 variants were identified in our cohort, 7 of them novel. Patients carrying two missense alleles showed a later disease onset than those with 1 or 2 truncating variants (log-rank test p <0.05). While 60% of patients carrying a missense/missense genotype presented symptoms before or during the first year of life, almost all patients with at least 1 truncating allele (91%) had an AAO ≤1 year (p <0.05). CONCLUSION: Our findings suggest an association between the type of RPE65 variant carried and AAO. These findings provide useful data on RPE65-associated IRD phenotypes and may help improve clinical and therapeutic management of these patients.

The Effect of Ethanol Treatment on the Quality of a New Table Grape Cultivar It 681-30 Stored at Low Temperature and after a 7-Day Shelf-Life Period at 20 °C: A Molecular Approach.

Despite the fact that many studies have examined the effectiveness of different gaseous postharvest treatments applied at low temperature to maintain table grape quality, the use of ethanol vapor has hardly been investigated. Thus, this work has studied the effectiveness of ethanol vapor-generating sachets in the maintenance of It 681-30 table grape quality, a new cultivar, during storage at low temperature and after the shelf-life period at 20 °C. To this end, various quality assessments have been carried out and the effect of the ethanol treatment on the expression of different genes (phenylpropanoids, transcription factors, PRs, and aquaporins) was determined. The results indicated that the application of ethanol vapor reduced the total decay incidence, weight loss, and the rachis browning index in It 681-30 grapes stored at 0 °C and after the shelf-life period at 20 °C, as compared to non-treated samples. Moreover, the modulation of STS7 and the different PR genes analyzed seems to play a part in the molecular mechanisms activated to cope with fungal attacks during the postharvest of It 681-30 grapes, and particularly during the shelf-life period at 20 °C. Furthermore, the expression of aquaporin transcripts was activated in samples showing higher weight loss. Although further work is needed to elucidate the role of ethanol in table grape quality, the results obtained in this work provide new insight into the transcriptional regulation triggered by ethanol treatment.

Ancestry-Specific Analyses Reveal Differential Demographic Histories and Opposite Selective Pressures in Modern South Asian Populations.

Genetic variation in contemporary South Asian populations follows a northwest to southeast decreasing cline of shared West Eurasian ancestry. A growing body of ancient DNA evidence is being used to build increasingly more realistic models of demographic changes in the last few thousand years. Through high-quality modern genomes, these models can be tested for gene and genome level deviations. Using local ancestry deconvolution and masking, we reconstructed population-specific surrogates of the two main ancestral components for more than 500 samples from 25 South Asian populations and showed our approach to be robust via coalescent simulations. Our f3 and f4 statistics-based estimates reveal that the reconstructed haplotypes are good proxies for the source populations that admixed in the area and point to complex interpopulation relationships within the West Eurasian component, compatible with multiple waves of arrival, as opposed to a simpler one wave scenario. Our approach also provides reliable local haplotypes for future downstream analyses. As one such example, the local ancestry deconvolution in South Asians reveals opposite selective pressures on two pigmentation genes (SLC45A2 and SLC24A5) that are common or fixed in West Eurasians, suggesting post-admixture purifying and positive selection signals, respectively.

Table Grapes during Postharvest Storage: A Review of the Mechanisms Implicated in the Beneficial Effects of Treatments Applied for Quality Retention.

Table grape is a fruit with increasing interest due to its attributes and nutritional compounds. During recent years, new cultivars such as those without seeds and with new flavors have reached countries around the world. For this reason, postharvest treatments that retain fruit quality need to be improved. However, little is known to date about the biochemical and molecular mechanisms related with observed quality improvements. This review aims to examine existing literature on the different mechanisms. Special attention will be placed on molecular mechanisms which activate and regulate the different postharvest treatments applied in order to improve table grape quality.

Genetic Variation, Not Cell Type of Origin, Underlies the Majority of Identifiable Regulatory Differences in iPSCs.

The advent of induced pluripotent stem cells (iPSCs) revolutionized human genetics by allowing us to generate pluripotent cells from easily accessible somatic tissues. This technology can have immense implications for regenerative medicine, but iPSCs also represent a paradigm shift in the study of complex human phenotypes, including gene regulation and disease. Yet, an unresolved caveat of the iPSC model system is the extent to which reprogrammed iPSCs retain residual phenotypes from their precursor somatic cells. To directly address this issue, we used an effective study design to compare regulatory phenotypes between iPSCs derived from two types of commonly used somatic precursor cells. We find a remarkably small number of differences in DNA methylation and gene expression levels between iPSCs derived from different somatic precursors. Instead, we demonstrate genetic variation is associated with the majority of identifiable variation in DNA methylation and gene expression levels. We show that the cell type of origin only minimally affects gene expression levels and DNA methylation in iPSCs, and that genetic variation is the main driver of regulatory differences between iPSCs of different donors. Our findings suggest that studies using iPSCs should focus on additional individuals rather than clones from the same individual.