Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans

Joint models reveal genetic architecture of pubertal stage transitions and their association with BMI in admixed Chilean population

Early or late pubertal onset can lead to disease in adulthood, including cancer, obesity, type 2 diabetes, metabolic disorders, bone fractures, and psychopathologies. Thus, knowing the age at which puberty is attained is crucial as it can serve as a risk factor for future diseases. Pubertal development is divided into five stages of sexual maturation in boys and girls according to the standardized Tanner scale. We performed genome-wide association studies (GWAS) on the "Growth and Obesity Chilean Cohort Study" cohort composed of admixed children with mainly European and Native American ancestry. Using joint models that integrate time-to-event data with longitudinal trajectories of body mass index (BMI), we identified genetic variants associated with phenotypic transitions between pairs of Tanner stages. We identified $42$ novel significant associations, most of them in boys. The GWAS on Tanner $3\rightarrow 4$ transition in boys captured an association peak around the growth-related genes LARS2 and LIMD1 genes, the former of which causes ovarian dysfunction when mutated. The associated variants are expression and splicing Quantitative Trait Loci regulating gene expression and alternative splicing in multiple tissues. Further, higher individual Native American genetic ancestry proportions predicted a significantly earlier puberty onset in boys but not in girls. Finally, the joint models identified a longitudinal BMI parameter significantly associated with several Tanner stages' transitions, confirming the association of BMI with pubertal timing.

Ancestral Origins and Admixture History of Kazakhs

Kazakh people, like many other populations that settled in Central Asia, demonstrate an array of mixed anthropological features of East Eurasian (EEA) and West Eurasian (WEA) populations, indicating a possible scenario of biological admixture between already differentiated EEA and WEA populations. However, their complex biological origin, genomic makeup, and genetic interaction with surrounding populations are not well understood. To decipher their genetic structure and population history, we conducted, to our knowledge, the first whole-genome sequencing study of Kazakhs residing in Xinjiang (KZK). We demonstrated that KZK derived their ancestries from 4 ancestral source populations: East Asian (∼39.7%), West Asian (∼28.6%), Siberian (∼23.6%), and South Asian (∼8.1%). The recognizable interactions of EEA and WEA ancestries in Kazakhs were dated back to the 15th century BCE. Kazakhs were genetically distinctive from the Uyghurs in terms of their overall genomic makeup, although the 2 populations were closely related in genetics, and both showed a substantial admixture of western and eastern peoples. Notably, we identified a considerable sex-biased admixture, with an excess of western males and eastern females contributing to the KZK gene pool. We further identified a set of genes that showed remarkable differentiation in KZK from the surrounding populations, including those associated with skin color (SLC24A5, OCA2), essential hypertension (HLA-DQB1), hypertension (MTHFR, SLC35F3), and neuron development (CNTNAP2). These results advance our understanding of the complex history of contacts between Western and Eastern Eurasians, especially those living or along the old Silk Road.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Recent advances in Forensic DNA Phenotyping of appearance, ancestry and age

Forensic DNA Phenotyping (FDP) comprises the prediction of a person's externally visible characteristics regarding appearance, biogeographic ancestry and age from DNA of crime scene samples, to provide investigative leads to help find unknown perpetrators that cannot be identified with forensic STR-profiling. In recent years, FDP has advanced considerably in all of its three components, which we summarize in this review article. Appearance prediction from DNA has broadened beyond eye, hair and skin color to additionally comprise other traits such as eyebrow color, freckles, hair structure, hair loss in men, and tall stature. Biogeographic ancestry inference from DNA has progressed from continental ancestry to sub-continental ancestry detection and the resolving of co-ancestry patterns in genetically admixed individuals. Age estimation from DNA has widened beyond blood to more somatic tissues such as saliva and bones as well as new markers and tools for semen. Technological progress has allowed forensically suitable DNA technology with largely increased multiplex capacity for the simultaneous analysis of hundreds of DNA predictors with targeted massively parallel sequencing (MPS). Forensically validated MPS-based FDP tools for predicting from crime scene DNA i) several appearance traits, ii) multi-regional ancestry, iii) several appearance traits together with multi-regional ancestry, and iv) age from different tissue types, are already available. Despite recent advances that will likely increase the impact of FDP in criminal casework in the near future, moving reliable appearance, ancestry and age prediction from crime scene DNA to the level of detail and accuracy police investigators may desire, requires further intensified scientific research together with technical developments and forensic validations as well as the necessary funding.

Investigating the morphology and genetics of scalp and facial hair characteristics for phenotype prediction

Microscopic traits and ultrastructure of hair such as cross-sectional shape, pigmentation, curvature, and internal structure help determine the level of variations between and across human populations. Apart from cosmetics and anthropological applications, such as determining species, somatic origin (body area), and biogeographic ancestry, the evidential value of hair has increased with rapid progression in the area of forensic DNA phenotyping (FDP). Individuals differ in the features of their scalp hair (greying, shape, colour, balding, thickness, and density) and facial hair (eyebrow thickness, monobrow, and beard thickness) features. Scalp and facial hair characteristics are genetically controlled and lead to visible inter-individual variations within and among populations of various ethnic origins. Hence, these characteristics can be exploited and made more inclusive in FDP, thereby leading to more comprehensive, accurate, and robust prediction models for forensic purposes. The present article focuses on understanding the genetics of scalp and facial hair characteristics with the goal to develop a more inclusive approach to better understand hair biology by integrating hair microscopy with genetics for genotype-phenotype correlation research.

Genomic diversity and post-admixture adaptation in the Uyghurs

Population admixture results in genome-wide combinations of genetic variants derived from different ancestral populations of distinct ancestry, thus providing a unique opportunity for understanding the genetic determinants of phenotypic variation in humans. Here, we used whole-genome sequencing of 92 individuals with high coverage (30–60×) to systematically investigate genomic diversity in the Uyghurs living in Xinjiang, China (XJU), an admixed population of both European-like and East-Asian-like ancestry. The XJU population shows greater genetic diversity, especially a higher proportion of rare variants, compared with their ancestral source populations, corresponding to greater phenotypic diversity of XJU. Admixture-induced functional variants in EDAR were associated with the diversity of facial morphology in XJU. Interestingly, the interaction of functional variants between SLC24A5 and OCA2 likely influences the diversity of skin pigmentation. Notably, selection has seemingly been relaxed or canceled in several genes with significantly biased ancestry, such as HERC2–OCA2. Moreover, signatures of post-admixture adaptation in XJU were identified, including genes related to metabolism (e.g. CYP2D6), digestion (e.g. COL11A1), olfactory perception (e.g. ANO2) and immunity (e.g. HLA). Our results demonstrated population admixture as a driving force, locally or globally, in shaping human genetic and phenotypic diversity as well as in adaptive evolution.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

Reproducibility in the UK biobank of genome-wide significant signals discovered in earlier genome-wide association studies

With the establishment of large biobanks, discovery of single nucleotide variants (SNVs, also known as single nucleotide polymorphisms (SNVs)) associated with various phenotypes has accelerated. An open question is whether genome-wide significant SNVs identified in earlier genome-wide association studies (GWAS) are replicated in later GWAS conducted in biobanks. To address this, we examined a publicly available GWAS database and identified two, independent GWAS on the same phenotype (an earlier, “discovery” GWAS and a later, “replication” GWAS done in the UK biobank). The analysis evaluated 136,318,924 SNVs (of which 6289 reached P < 5e−8 in the discovery GWAS) from 4,397,962 participants across nine phenotypes. The overall replication rate was 85.0%; although lower for binary than quantitative phenotypes (58.1% versus 94.8% respectively). There was a 18.0% decrease in SNV effect size for binary phenotypes, but a 12.0% increase for quantitative phenotypes. Using the discovery SNV effect size, phenotype trait (binary or quantitative), and discovery P value, we built and validated a model that predicted SNV replication with area under the Receiver Operator Curve = 0.90. While non-replication may reflect lack of power rather than genuine false-positives, these results provide insights about which discovered associations are likely to be replicated across subsequent GWAS.

From mutation to mechanism: deciphering the molecular function of genetic variants linked to human ageing

Many of the leading causes of death in humans, such as cardiovascular disease, type 2 diabetes and Alzheimer’s disease are influenced by biological mechanisms that become dysregulated with increasing age. Hence, by targeting these ageing-related mechanisms, we may be able to improve health in old age. Ageing is partly heritable and genetic studies have been moderately successful in identifying genetic variants associated with ageing-related phenotypes (lifespan, healthspan and longevity). To decipher the mechanisms by which the identified variants influence ageing, studies that focus on their functional validation are vital. In this perspective, we describe the steps that could be taken in the process of functional validation: (1) in silico characterisation using bioinformatic tools; (2) in vitro characterisation using cell lines or organoids; and (3) in vivo characterisation studies using model organisms. For the in vivo characterisation, it is important to focus on translational phenotypes that are indicative of both healthspan and lifespan, such as the frailty index, to inform subsequent intervention studies. The depth of functional validation of a genetic variant depends on its location in the genome and conservation in model organisms. Moreover, some variants may prove to be hard to characterise due to context-dependent effects related to the experimental environment or genetic background. Future efforts to functionally characterise the (newly) identified genetic variants should shed light on the mechanisms underlying ageing and will help in the design of targeted interventions to improve health in old age.

Long Noncoding RNA SOX2-OT: Regulations, Functions, and Roles on Mental Illnesses, Cancers, and Diabetic Complications

SRY-box transcription factor 2 (SOX2) overlapping transcript (SOX2-OT) is an evolutionarily conserved long noncoding RNA. Its intronic region contains the SOX2 gene, the major regulator of the pluripotency of embryonic stem cells. The human SOX2-OT gene comprises multiple exons and has multiple transcription start sites and generates hundreds of transcripts. Transcription factors (IRF4, AR, and SOX3), transcriptional inhibitors (NSPc1, MTA3, and YY1), and miRNAs (miR-211 and miR-375) have been demonstrated to control certain SOX2-OT transcript level at the transcriptional or posttranscriptional levels. Accumulated evidence indicates its crucial roles in the regulation of the SOX2 gene, miRNAs, and transcriptional process. Restricted expression of SOX2-OT transcripts in the brain results in the association between SOX2-OT single nucleotide polymorphisms and mental illnesses such as schizophrenia and anorexia nervosa. SOX2-OT is notably elevated in tumor tissues, and a high level of SOX2-OT is well correlated with poor clinical outcomes in cancer patients, leading to the establishment of its role as an oncogene and a prognostic or diagnostic biomarker for cancers. The emerging evidence supports that SOX2-OT mediates diabetic complications. In summary, SOX2-OT has diversified functions and could be a therapeutic target for various diseases.

Unravelling the complex genetics of common kidney diseases: from variants to mechanisms

Genome-wide association studies (GWAS) have identified hundreds of loci associated with kidney-related traits such as glomerular filtration rate, albuminuria, hypertension, electrolyte and metabolite levels. However, these impressive, large-scale mapping approaches have not always translated into an improved understanding of disease or development of novel therapeutics. GWAS have several important limitations. Nearly all disease-associated risk loci are located in the non-coding region of the genome and therefore, their target genes, affected cell types and regulatory mechanisms remain unknown. Genome-scale approaches can be used to identify associations between DNA sequence variants and changes in gene expression (quantified through bulk and single-cell methods), gene regulation and other molecular quantitative trait studies, such as chromatin accessibility, DNA methylation, protein expression and metabolite levels. Data obtained through these approaches, used in combination with robust computational methods, can deliver robust mechanistic inferences for translational exploitation. Understanding the genetic basis of common kidney diseases means having a comprehensive picture of the genes that have a causal role in disease development and progression, of the cells, tissues and organs in which these genes act to affect the disease, of the cellular pathways and mechanisms that drive disease, and of potential targets for disease prevention, detection and therapy.

Prefrontal cortex eQTLs/mQTLs enriched in genetic variants associated with alcohol use disorder and other diseases

Aim: This study aimed to investigate the function of genome-wide association study (GWAS)-identified variants associated with alcohol use disorder (AUD)/comorbid psychiatric disorders. Materials & methods: Genome-wide genotype, transcriptome and DNA methylome data were obtained from postmortem prefrontal cortex (PFC) of 48 Caucasians (24 AUD cases/24 controls). Expression/methylation quantitative trait loci (eQTL/mQTL) were identified and their enrichment in GWAS signals for the above disorders were analyzed. Results: PFC cis-eQTLs (923 from cases+controls, 27 from cases and 98 from controls) and cis-mQTLs (9,932 from cases+controls, 264 from cases and 695 from controls) were enriched in GWAS-identified genetic variants for the above disorders. Cis-eQTLs from AUD cases were mapped to morphine addiction-related genes. Conclusion: PFC cis-eQTLs/cis-mQTLs influence gene expression/DNA methylation patterns, thus increasing the disease risk.

Objectives, design and main findings until 2020 from the Rotterdam Study

The Rotterdam Study is an ongoing prospective cohort study that started in 1990 in the city of Rotterdam, The Netherlands. The study aims to unravel etiology, preclinical course, natural history and potential targets for intervention for chronic diseases in mid-life and late-life. The study focuses on cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1700 research articles and reports. This article provides an update on the rationale and design of the study. It also presents a summary of the major findings from the preceding 3 years and outlines developments for the coming period.

CRISPR/Cas9 guided genome and epigenome engineering and its therapeutic applications in immune mediated diseases

Recent developments in the nucleic acid editing technologies have provided a powerful tool to precisely engineer the genome and epigenome for studying many aspects of immune cell differentiation and development as well as several immune mediated diseases (IMDs) including autoimmunity and cancer. Here, we discuss the recent technological achievements of the CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-based RNA-guided genome and epigenome editing toolkit and provide an insight into how CRISPR/Cas9 (CRISPR Associated Protein 9) toolbox could be used to examine genetic and epigenetic mechanisms underlying IMDs. In addition, we will review the progress in CRISPR/Cas9-based genome-wide genome and epigenome screens in various cell types including immune cells. Finally, we will discuss the potential of CRISPR/Cas9 in defining the molecular function of disease associated SNPs overlapping gene regulatory elements.