Hundreds of variants clustered in genomic loci and biological pathways affect human height

The distribution of common-variant effect sizes

The genetic effect-size distribution of a disease describes the number of risk variants, the range of their effect sizes and sample sizes that will be required to discover them. Accurate estimation has been a challenge. Here I propose Fourier Mixture Regression (FMR), validating that it accurately estimates real and simulated effect-size distributions. Applied to summary statistics for ten diseases (average [Formula: see text]), FMR estimates that 100,000-1,000,000 cases will be required for genome-wide significant SNPs to explain 50% of SNP heritability. In such large studies, genome-wide significance becomes increasingly conservative, and less stringent thresholds achieve high true positive rates if confounding is controlled. Across traits, polygenicity varies, but the range of their effect sizes is similar. Compared with effect sizes in the top 10% of heritability, including most discovered thus far, those in the bottom 10-50% are orders of magnitude smaller and more numerous, spanning a large fraction of the genome.

Accelerated deciphering of the genetic architecture of agricultural economic traits in pigs using a low-coverage whole-genome sequencing strategy

BACKGROUND

Uncovering the genetic architecture of economic traits in pigs is important for agricultural breeding. However, high-density haplotype reference panels are unavailable in most agricultural species, limiting accurate genotype imputation in large populations. Moreover, the infinitesimal model of quantitative traits implies that weak association signals tend to be spread across most of the genome, further complicating the genetic analysis. Hence, there is a need to develop new methods for sequencing large cohorts without large reference panels.

RESULTS

We describe a Tn5-based highly accurate, cost- and time-efficient, low-coverage sequencing method to obtain 11.3 million whole-genome single-nucleotide polymorphisms in 2,869 Duroc boars at a mean depth of 0.73×. On the basis of these single-nucleotide polymorphisms, a genome-wide association study was performed, resulting in 14 quantitative trait loci (QTLs) for 7 of 21 important agricultural traits in pigs. These QTLs harbour genes, such as ABCD4 for total teat number and HMGA1 for back fat thickness, and provided a starting point for further investigation. The inheritance models of the different traits varied greatly. Most follow the minor-polygene model, but this can be attributed to different reasons, such as the shaping of genetic architecture by artificial selection for this population and sufficiently interconnected minor gene regulatory networks.

CONCLUSIONS

Genome-wide association study results for 21 important agricultural traits identified 14 QTLs/genes and showed their genetic architectures, providing guidance for genetic improvement harnessing genomic features. The Tn5-based low-coverage sequencing method can be applied to large-scale genome studies for any species without a good reference panel and can be used for agricultural breeding.

Short stature with low insulin-like growth factor 1 availability due to pregnancy-associated plasma protein A2 deficiency in a Saudi family

In 2016 a new syndrome with postnatal short stature and low IGF1 bioavailability caused by biallelic loss-of-function mutations in the gene encoding the metalloproteinase pregnancy-associated plasma protein A2 (PAPP-A2) was described in two families. Here we report two siblings of a third family from Saudi Arabia with postnatal growth retardation and decreased IGF1 availability due to a new homozygous nonsense mutation (p.Glu886* in exon 7) in PAPPA2. The two affected males showed progressively severe short stature starting around 8 years of age, moderate microcephaly, decreased bone mineral density, and high circulating levels of total IGF1, IGFBP3, and the IGF acid-labile subunit (IGFALS), with decreased free IGF1 concentrations. Interestingly, circulating IGF2 and IGFBP5 were not increased. An increase in growth velocity and height was seen in the prepuberal patient in response to rhIGF1. These patients contribute to the confirmation of the clinical picture associated with PAPP-A2 deficiency and that the PAPPA2 gene should be studied in all patients with short stature with this characteristic phenotype. Hence, pediatric endocrinologists should measure circulating PAPP-A2 levels in the study of short stature as very low or undetectable levels of this protein can help to focus the diagnosis and treatment.

Genome-wide association study in the Taiwan biobank identifies four novel genes for human height: NABP2, RASA2, RNF41 and SLC39A5

Numerous genome-wide association studies (GWASs) have been conducted for the identification of genetic variants involved with human height. The vast majority of these studies, however, have been conducted in populations of European ancestry. Here, we report the first GWAS of adult height in the Taiwan Biobank using a discovery sample of 14 571 individuals and an independent replication sample of 20 506 individuals. From our analysis we generalize to the Taiwanese population genome-wide significant associations with height and 18 previously identified genes in European and non-Taiwanese East Asian populations. We also identify and replicate, at the genome-wide significance level, associated variants for height in four novel genes at two loci that have not previously been reported: RASA2 on chromosome 3 and NABP2, RNF41, and SLC39A5 at 12q13.3 on chromosome 12. RASA2 and RNF41 are strong candidates for having a role in height with copy number and loss of function variants in RASA2 previously found to be associated with short stature disorders, and decreased expression of the RNF41 gene resulting in insulin resistance in skeletal muscle. The results from our analysis of the Taiwan Biobank underscore the potential for the identification of novel genetic discoveries in underrepresented worldwide populations, even for traits, such as height, that have been extensively investigated in large-scale studies of European ancestry populations.

Detection of 15-bp Deletion Mutation within PLAG1 Gene and Its Effects on Growth Traits in Goats

Simple Summary

Goats have always served as an important domesticated livestock. PLAG1 is a major gene that affects the stature and growth of animals. Body size traits are very important for goats as it directly affects the economic characteristics of meat and cashmere production. This study showed that the 15-base pair (bp) InDel (rs637141549) can significantly affect growth traits such as body weight, height, height at hip cross, chest circumference, hip width and body index of goats through the detection of large samples (n = 1581) in four indigenous breeds. Accordingly, it is suggested that the deletion mutation can be used as a potential molecular marker that significantly affects goat growth traits. Moreover, the 15bp deletion mutation can be used as a potential molecular marker, which significantly affects the growth traits of goats and plays an important role in animal husbandry production.

Abstract

Stature and weight are important growth and development traits for animals, which also significantly affect the productivity of livestock. Polymorphic adenoma gene 1 (PLAG1) is located in the growth-related quantitative trait nucleotides (QTN), and its variation has been determined to significantly affect the body stature of bovines. This study found that novel 15-bp InDel could significantly influence important growth traits in goats. The frequencies of genotypes of the 15-bp mutation and relationship with core growth traits such as body weight, body height, height at hip cross, chest circumference, hip width and body index were explored in 1581 individuals among 4 Chinese native goat breeds. The most frequent genotypes of Shaanbei white Cashmere goat (SWCG), Inner Mongolia White Cashmere goat (IMCG) and Guanzhong Dairy goat (GZDG) were II genotypes (insertion/insertion), and the frequency of ID genotype (insertion/deletion) was found to be slightly higher than that of II genotype in Hainan Black goat (HNBG), showing that the frequency of the I allele was higher than that of the D allele. In adult goats, there were significant differences between 15-bp variation and body weight, chest circumference and body height traits in SWCG (p < 0.05). Furthermore, the locus was also found to be significantly correlated with the body index of HNBG (p = 0.044) and hip width in GZDG (p = 0.002). In regard to lambs, there were significant differences in height at the hip cross of SWCG (p = 0.036) and hip width in IMWC (p = 0.005). The corresponding results suggest that the 15-bp InDel mutation of PLAG1 is associated with the regulation of important growth characteristics of both adult and lamb of goats, which may serve as efficient molecular markers for goat breeding.

Phenotypic and genomic diversification with isolation by environment along elevational gradients in a neotropical treefrog

Understanding how geographic and environmental heterogeneity drive local patterns of genetic variation is a major goal of ecological genomics and a key question in evolutionary biology. The tropical Andes and inter-Andean valleys are shaped by markedly heterogeneous landscapes, where species experience strong selective processes. We examined genome-wide SNP data together with behavioural and ecological traits (mating calls and body size) known to contribute to genetic isolation in anurans in the banana tree-dwelling frog, Boana platanera, distributed across an environmental gradient in Central Colombia (northern South America). Here, we analysed the relationships between environmentally (temperature and precipitation) associated genetic and phenotypic differentiation and the potential drivers of isolation by environment along an elevation gradient. We identified candidate SNPs associated with temperature and body size, which follow a clinal pattern of genome-wide differentiation tightly coupled with phenotypic variation: as elevation increases, B. platanera exhibits larger body size and longer call duration with more pulses but lower pulse rate and frequency. Thus, the environmental landscape has rendered a scenario where isolation by environment and candidate loci show concordance with phenotypic divergence in this tropical frog along an elevation gradient in the Colombian Andes. Our study sets the basis for evaluating the role of temperature in the genetic structure and local adaptation in tropical treefrogs and its putative effect on life cycle (embryos, tadpoles, adults) along elevation gradients.

Does height and IGF-I determine pubertal timing in girls?

BACKGROUND

Pubertal timing is closely linked to growth regulated by the growth hormone/insulin-like factor (GH/IGF) axis that includes IGF-regulating factors such as pregnancy-associated plasma protein-A/A2 (PAPP-A/PAPP-A2) and stanniocalcin 2 (STC2). We investigated the association between height, IGF-I concentration, and PAPPA, PAPPA2, and STC2 genotypes on the timing of female pubertal milestones.

METHODS

Height, IGF-I, and genotypes were analyzed in 1382 Danish girls from the general population, 67 patients with tall stature (height ≥2 SD), and 124 patients with short stature (height ≤-2 SD). The main outcomes were breast stage and menarche.

RESULTS

Thelarche occurred significantly earlier in patients with tall stature (mean age 9.37 years [95% confidence interval (CI) 8.87-9.87]) and later in patients with short stature (11.07 years [95% CI 10.7-11.43]) compared with girls within the normal range (9.96 years [95% CI 9.85-10.07]) (p = 0.02 and p < 0.01, respectively). Girls with higher IGF-I levels experienced thelarche and menarche earlier compared with the rest of the cohort (p < 0.01). Genotypes were not associated with age at thelarche nor menarche, but the PAPPA2 minor allele carriers were shorter compared with major allele carriers, p = 0.03.

CONCLUSIONS

Height and IGF-I, but not PAPP-A, PAPP-A2, and STC2 genotypes, were negatively associated with age at thelarche and menarche.

IMPACT

Girls with tall and short stature enter puberty earlier and later compared with girls with normal height. Girls with higher insulin-growth factor-I in childhood enter puberty earlier. Pubertal timing is influenced by longitudinal growth and IGF-I levels earlier in childhood. Childhood growth and the levels of IGF-I in childhood may be biomarkers of pubertal timing.

Advancing drug discovery using the power of the human genome

Human genetics plays an increasingly important role in drug development and population health. Here we review the history of human genetics in the context of accelerating the discovery of therapies, present examples of how human genetics evidence supports successful drug targets, and discuss how polygenic risk scores could be beneficial in various clinical settings. We highlight the value of direct-to-consumer platforms in the era of fast-paced big data biotechnology, and how diverse genetic and health data can benefit society. © 2021 23andMe, Inc. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Constrained maximum likelihood-based Mendelian randomization robust to both correlated and uncorrelated pleiotropic effects

With the increasing availability of large-scale GWAS summary data on various complex traits and diseases, there have been tremendous interests in applications of Mendelian randomization (MR) to investigate causal relationships between pairs of traits using SNPs as instrumental variables (IVs) based on observational data. In spite of the potential significance of such applications, the validity of their causal conclusions critically depends on some strong modeling assumptions required by MR, which may be violated due to the widespread (horizontal) pleiotropy. Although many MR methods have been proposed recently to relax the assumptions by mainly dealing with uncorrelated pleiotropy, only a few can handle correlated pleiotropy, in which some SNPs/IVs may be associated with hidden confounders, such as some heritable factors shared by both traits. Here we propose a simple and effective approach based on constrained maximum likelihood and model averaging, called cML-MA, applicable to GWAS summary data. To deal with more challenging situations with many invalid IVs with only weak pleiotropic effects, we modify and improve it with data perturbation. Extensive simulations demonstrated that the proposed methods could control the type I error rate better while achieving higher power than other competitors. Applications to 48 risk factor-disease pairs based on large-scale GWAS summary data of 3 cardio-metabolic diseases (coronary artery disease, stroke, and type 2 diabetes), asthma, and 12 risk factors confirmed its superior performance.

Genome-wide association study of neck circumference identifies sex-specific loci independent of generalized adiposity

BACKGROUND/OBJECTIVES

Neck circumference, an index of upper airway fat, has been suggested to be an important measure of body-fat distribution with unique associations with health outcomes such as obstructive sleep apnea and metabolic disease. This study aims to study the genetic bases of neck circumference.

METHODS

We conducted a multi-ethnic genome-wide association study of neck circumference, adjusted and unadjusted for BMI, in up to 15,090 European Ancestry (EA) and African American (AA) individuals. Because sexually dimorphic associations have been observed for anthropometric traits, we conducted both sex-combined and sex-specific analysis.

RESULTS

We identified rs227724 near the Noggin (NOG) gene as a possible quantitative locus for neck circumference in men (N = 8831, P = 1.74 × 10-9) but not in women (P = 0.08). The association was replicated in men (N = 1554, P = 0.045) in an independent dataset. This locus was previously reported to be associated with human height and with self-reported snoring. We also identified rs13087058 on chromosome 3 as a suggestive locus in sex-combined analysis (N = 15090, P = 2.94 × 10-7; replication P =0.049). This locus was also associated with electrocardiogram-assessed PR interval and is a cis-expression quantitative locus for the PDZ Domain-containing ring finger 2 (PDZRN3) gene. Both NOG and PDZRN3 interact with members of transforming growth factor-beta superfamily signaling proteins.

CONCLUSIONS

Our study suggests that neck circumference may have unique genetic basis independent of BMI.

Genome-wide association study of body fat distribution traits in Hispanics/Latinos from the HCHS/SOL

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs

Background

The Chinese Chenghua pig (CHP) is a typical Chinese domestic fatty pig breed with superior meat quality characteristics, while the Yorkshire pig (YP) has the characteristics of fast growth and a high rate of lean meat. Long term natural selection and artificial selection resulted in great phenotypic differences between the two breeds, including growth, development, production performance, meat quality, and coat color. However, genome-wide DNA methylation differences between CHP and YP remain unclear.

Results

DNA methylation data were generated for muscle tissues of CHP and YP using reduced representation bisulfite sequencing (RRBS). In this study, a total of 2,416,211 CpG sites were identified. Besides, the genome-wide DNA methylation analysis revealed 722 differentially methylated regions (DMRs) and 466 differentially methylated genes (DMGs) in pairwise CHP vs. YP comparison. Six key genomic regions (Sus scrofa chromosome (SSC)1:253.47–274.23 Mb, SSC6:148.71–169.49 Mb, SSC7:0.25–9.86 Mb, SSC12:43.06–61.49 Mb, SSC14:126.43–140.95 Mb, and SSC18:49.17–54.54 Mb) containing multiple DMRs were identified, and differences of methylation patterns in these regions may be related to phenotypic differences between CHP and YP. Based on the functional analysis of DMGs, 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT and MYH8) were considered as important candidate genes associated with muscle development and meat quality traits in pigs.

Conclusions

This study explored the difference in meat quality between CHP and YP from the epigenetic point of view, which has important reference significance for the local pork industry and pork food processing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00977-0.

Dual Characters of GH-IGF1 Signaling Pathways in Radiotherapy and Post-radiotherapy Repair of Cancers

Radiotherapy remains one of the most important cancer treatment modalities. In the course of radiotherapy for tumor treatment, the incidental irradiation of adjacent tissues could not be completely avoided. DNA damage is one of the main factors of cell death caused by ionizing radiation, including single-strand (SSBs) and double-strand breaks (DSBs). The growth hormone-Insulin-like growth factor 1 (GH-IGF1) axis plays numerous roles in various systems by promoting cell proliferation and inhibiting apoptosis, supporting its effects in inducing the development of multiple cancers. Meanwhile, the GH-IGF1 signaling involved in DNA damage response (DDR) and DNA damage repair determines the radio-resistance of cancer cells subjected to radiotherapy and repair of adjacent tissues damaged by radiotherapy. In the present review, we firstly summarized the studies on GH-IGF1 signaling in the development of cancers. Then we discussed the adverse effect of GH-IGF1 signaling in radiotherapy to cancer cells and the favorable impact of GH-IGF1 signaling on radiation damage repair to adjacent tissues after irradiation. This review further summarized recent advances on research into the molecular mechanism of GH-IGF1 signaling pathway in these effects, expecting to specify the dual characters of GH-IGF1 signaling pathways in radiotherapy and post-radiotherapy repair of cancers, subsequently providing theoretical basis of their roles in increasing radiation sensitivity during cancer radiotherapy and repairing damage after radiotherapy.

Translating genetic association of lipid levels for biological and clinical application

PURPOSE OF REVIEW

This review focuses on the foundational evidence from the last two decades of lipid genetics research and describes the current status of data-driven approaches for transethnic GWAS, fine-mapping, transcriptome informed fine-mapping, and disease prediction.

RECENT FINDINGS

Current lipid genetics research aims to understand the association mechanisms and clinical relevance of lipid loci as well as to capture population specific associations found in global ancestries. Recent genome-wide trans-ethnic association meta-analyses have identified 118 novel lipid loci reaching genome-wide significance. Gene-based burden tests of whole exome sequencing data have identified three genes-PCSK9, LDLR, and APOB-with significant rare variant burden associated with familial dyslipidemia. Transcriptome-wide association studies discovered five previously unreported lipid-associated loci. Additionally, the predictive power of genome-wide genetic risk scores amalgamating the polygenic determinants of lipid levels can potentially be used to increase the accuracy of coronary artery disease prediction.

CONCLUSIONS

Lipids are one of the most successful group of traits in the era of genome-wide genetic discovery for identification of novel loci and plausible drug targets. However, a substantial fraction of lipid trait heritability remains unexplained. Further analysis of diverse ancestries and state of the art methods for association locus refinement could potentially reveal some of this missing heritability and increase the clinical application of the genomic association results.

Brain-derived neurotrophic factor Val66Met polymorphism is associated with mild cognitive impairment in elderly patients with type 2 diabetes: a case-controlled study

BACKGROUND

Brain-derived neurotrophic factor (BDNF) Val66Met polymorphism is reported to be associated with cognitive dysfunction, an important comorbidity factor in patients with type 2 diabetes mellitus (T2DM), especially in elderly populations, however, the underlying pathophysiological mechanisms are unclear.

AIM

This study was performed to investigate the association between BDNF Val66Met polymorphism and mild cognitive impairment (MCI) in elderly patients with T2DM.

METHODS

In total, 105 MCI and 105 normal cognition controls of T2DM patients were enrolled; all of the patients underwent neuropsychological assessments. BDNF Val66Met polymorphism was genotyped via TaqMan SNP genotyping assay. Data from clinical and laboratory-based examinations were collected.

RESULTS

The frequency of the BDNF Met allele was significantly higher in the MCI group than in the controls. Multiple regression analysis indicated an association of the Met allele with MCI in patients with T2DM (OR = 2.54; 95% CI 1.33-4.84; p = 0.005). Stratified by educational level, the BDNF Met allele was significantly associated with MCI in elderly T2DM patients (OR = 3.29; 95% CI 1.26-8.57; p = 0.015) among the group of low educational levels (< 12 years); however, the association was insignificant among those with higher educational levels.

DISCUSSION

BDNF Met allele carriers showed a higher frequency of MCI than Val/Val homozygotes in elderly T2DM patients. However, this association was only significant in patients with low education levels.

CONCLUSION

BDNF Val66Met polymorphism may have a potential role in MCI in elderly T2DM patients, especially those with low educational levels.

Dissecting polygenic signals from genome-wide association studies on human behaviour

Genome-wide association studies on human behavioural traits are producing large amounts of polygenic signals with significant predictive power and potentially useful biological clues. Behavioural traits are more distal and are less directly under biological control compared with physical characteristics, which makes the associated genetic effects harder to interpret. The results of genome-wide association studies for human behaviour are likely made up of a composite of signals from different sources. While sample sizes continue to increase, we outline additional steps that need to be taken to better delineate the origin of the increasingly stronger polygenic signals. In addition to genetic effects on the traits themselves, the major sources of polygenic signals are those that are associated with correlated traits, environmental effects and ascertainment bias. Advances in statistical approaches that disentangle polygenic effects from different traits as well as extending data collection to families and social circles with better geographical coverage will probably contribute to filling the gap of knowledge between genetic effects and behavioural outcomes.

Genetics of Body Fat Distribution: Comparative Analyses in Populations with European, Asian and African Ancestries

Preferential fat accumulation in visceral vs. subcutaneous depots makes obese individuals more prone to metabolic complications. Body fat distribution (FD) is regulated by genetics. FD patterns vary across ethnic groups independent of obesity. Asians have more and Africans have less visceral fat compared with Europeans. Consequently, Asians tend to be more susceptible to type 2 diabetes even with lower BMIs when compared with Europeans. To date, genome-wide association studies (GWAS) have identified more than 460 loci related to FD traits. However, the majority of these data were generated in European populations. In this review, we aimed to summarize recent advances in FD genetics with a focus on comparisons between European and non-European populations (Asians and Africans). We therefore not only compared FD-related susceptibility loci identified in three ethnicities but also discussed whether known genetic variants might explain the FD pattern heterogeneity across different ancestries. Moreover, we describe several novel candidate genes potentially regulating FD, including NID2, HECTD4 and GNAS, identified in studies with Asian populations. It is of note that in agreement with current knowledge, most of the proposed FD candidate genes found in Asians belong to the group of developmental genes.

Identification of candidate genes associated with slaughter traits in F2 chicken population using genome-wide association study

Slaughter traits are crucial economic traits of chickens. We performed a GWAS to discover critical loci and candidate genes for 21 slaughter traits in an F2 chicken population resulting from crossing Luxi gamecocks and recessive white feather broilers. We found some SNPs and genes which were significantly associated with keel length, head length, body slope length, bilateral leg weight without shin, bilateral foot weight, subcutaneous fat thickness, heart weight, muscular stomach weight and glandular stomach weight. This study provides references for further investigation of slaughter traits and molecular breeding in chicken.

A Genome-Wide Association Study of Novel Genetic Variants Associated With Anthropometric Traits in Koreans

Most previous genome-wide association studies (GWAS) have identified genetic variants associated with anthropometric traits. However, most of the evidence were reported in European populations. Anthropometric traits such as height and body fat distribution are significantly affected by gender and genetic factors. Here we performed GWAS involving 64,193 Koreans to identify the genetic factors associated with anthropometric phenotypes including height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip ratio. We found nine novel single-nucleotide polymorphisms (SNPs) and 59 independent genetic signals in genomic regions that were reported previously. Of the 19 SNPs reported previously, eight genetic variants at RP11-513I15.6 and one genetic variant at the RP11-977G19.10 region and six Asian-specific genetic variants were newly found. We compared our findings with those of previous studies in other populations. Five overlapping genetic regions (PAN2, ANKRD52, RNF41, HGMA1, and C6orf106) had been reported previously but none of the SNPs were independently identified in the current study. Seven of the nine newly found novel loci associated with height in women revealed a statistically significant skeletal expression of quantitative trait loci. Our study provides additional insight into the genetic effects of anthropometric phenotypes in East Asians.

Long-read sequencing of 3,622 Icelanders provides insight into the role of structural variants in human diseases and other traits

Long-read sequencing (LRS) promises to improve the characterization of structural variants (SVs). We generated LRS data from 3,622 Icelanders and identified a median of 22,636 SVs per individual (a median of 13,353 insertions and 9,474 deletions). We discovered a set of 133,886 reliably genotyped SV alleles and imputed them into 166,281 individuals to explore their effects on diseases and other traits. We discovered an association of a rare deletion in PCSK9 with lower low-density lipoprotein (LDL) cholesterol levels, compared to the population average. We also discovered an association of a multiallelic SV in ACAN with height; we found 11 alleles that differed in the number of a 57-bp-motif repeat and observed a linear relationship between the number of repeats carried and height. These results show that SVs can be accurately characterized at the population scale using LRS data in a genome-wide non-targeted approach and demonstrate how SVs impact phenotypes.

Genetics of Obesity: What We Have Learned Over Decades of Research

There is a genetic component to human obesity that accounts for 40% to 50% of the variability in body weight status but that is lower among normal weight individuals (about 30%) and substantially higher in the subpopulation of individuals with obesity and severe obesity (about 60%-80%). The appreciation that heritability varies across classes of BMI represents an important advance. After controlling for BMI, ectopic fat and fat distribution traits are characterized by heritability levels ranging from 30% to 55%. Defects in at least 15 genes are the cause of monogenic obesity cases, resulting mostly from deficiencies in the leptin-melanocortin signaling pathway. Approximately two-thirds of the BMI heritability can be imputed to common DNA variants, whereas low-frequency and rare variants explain the remaining fraction. Diminishing allele effect size is observed as the number of obesity-associated variants expands, with most BMI-increasing or -decreasing alleles contributing only a few grams or less to body weight. Obesity-promoting alleles exert minimal effects in normal weight individuals but have larger effects in individuals with a proneness to obesity, suggesting a higher penetrance; however, it is not known whether these larger effect sizes precede obesity or are caused by an obese state. The obesity genetic risk is conditioned by thousands of DNA variants that make genetically based obesity prevention and treatment a major challenge.

Joo JW. MARS: leveraging allelic heterogeneity to increase power of association testing

In standard genome-wide association studies (GWAS), the standard association test is underpowered to detect associations between loci with multiple causal variants with small effect sizes. We propose a statistical method, Model-based Association test Reflecting causal Status (MARS), that finds associations between variants in risk loci and a phenotype, considering the causal status of variants, only requiring the existing summary statistics to detect associated risk loci. Utilizing extensive simulated data and real data, we show that MARS increases the power of detecting true associated risk loci compared to previous approaches that consider multiple variants, while controlling the type I error.

Climate-Resilient Dairy Cattle Production: Applications of Genomic Tools and Statistical Models

The current changing climate trend poses a threat to the productive efficacy and welfare of livestock across the globe. This review is an attempt to synthesize information pertaining to the applications of various genomic tools and statistical models that are available to identify climate-resilient dairy cows. The different functional and economical traits which govern milk production play a significant role in determining the cost of milk production. Thus, identification of these traits may revolutionize the breeding programs to develop climate-resilient dairy cattle. Moreover, the genotype–environment interaction also influences the performance of dairy cattle especially during a challenging situation. The recent advancement in molecular biology has led to the development of a few biotechnological tools and statistical models like next-generation sequencing (NGS), microarray technology, whole transcriptome analysis, and genome-wide association studies (GWAS) which can be used to quantify the molecular mechanisms which govern the climate resilience capacity of dairy cows. Among these, the most preferred option for researchers around the globe was GWAS as this approach jointly takes into account all the genotype, phenotype, and pedigree information of farm animals. Furthermore, selection signatures can also help to demarcate functionally important regions in the genome which can be used to detect potential loci and candidate genes that have undergone positive selection in complex milk production traits of dairy cattle. These identified biomarkers can be incorporated in the existing breeding policies using genomic selection to develop climate-resilient dairy cattle.

Genome Reconstruction Attacks Against Genomic Data-Sharing Beacons

Sharing genome data in a privacy-preserving way stands as a major bottleneck in front of the scientific progress promised by the big data era in genomics. A community-driven protocol named genomic data-sharing beacon protocol has been widely adopted for sharing genomic data. The system aims to provide a secure, easy to implement, and standardized interface for data sharing by only allowing yes/no queries on the presence of specific alleles in the dataset. However, beacon protocol was recently shown to be vulnerable against membership inference attacks. In this paper, we show that privacy threats against genomic data sharing beacons are not limited to membership inference. We identify and analyze a novel vulnerability of genomic data-sharing beacons: genome reconstruction. We show that it is possible to successfully reconstruct a substantial part of the genome of a victim when the attacker knows the victim has been added to the beacon in a recent update. In particular, we show how an attacker can use the inherent correlations in the genome and clustering techniques to run such an attack in an efficient and accurate way. We also show that even if multiple individuals are added to the beacon during the same update, it is possible to identify the victim's genome with high confidence using traits that are easily accessible by the attacker (e.g., eye color or hair type). Moreover, we show how a reconstructed genome using a beacon that is not associated with a sensitive phenotype can be used for membership inference attacks to beacons with sensitive phenotypes (e.g., HIV+). The outcome of this work will guide beacon operators on when and how to update the content of the beacon and help them (along with the beacon participants) make informed decisions.

Growth genes are implicated in the evolutionary divergence of sympatric piscivorous and insectivorous rainbow trout (Oncorhynchus mykiss)

Background

Identifying ecologically significant phenotypic traits and the genomic mechanisms that underly them are crucial steps in understanding traits associated with population divergence. We used genome-wide data to identify genomic regions associated with key traits that distinguish two ecomorphs of rainbow trout (Oncorhynchus mykiss)—insectivores and piscivores—that coexist for the non-breeding portion of the year in Kootenay Lake, southeastern British Columbia. “Gerrards” are large-bodied, rapidly growing piscivores with high metabolic rates that spawn north of Kootenay Lake in the Lardeau River, in contrast to the insectivorous populations that are on average smaller in body size, with lower growth and metabolic rates, mainly forage on aquatic insects, and spawn in tributaries immediately surrounding Kootenay Lake. We used pool-seq data representing ~ 60% of the genome and 80 fish per population to assess the level of genomic divergence between ecomorphs and to identify and interrogate loci that may play functional or selective roles in their divergence.

Results

Genomic divergence was high between sympatric insectivores and piscivores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.188), and in fact higher than between insectivorous populations from Kootenay Lake and the Blackwater River (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.159) that are > 500 km apart. A window-based \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST analysis did not reveal “islands” of genomic differentiation; however, the window with highest \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST estimate did include a gene associated with insulin secretion. Although we explored the use of the “Local score” approach to identify genomic outlier regions, this method was ultimately not used because simulations revealed a high false discovery rate (~ 20%). Gene ontology (GO) analysis identified several growth processes as enriched in genes occurring in the ~ 200 most divergent genomic windows, indicating many loci of small effect involved in growth and growth-related metabolic processes are associated with the divergence of these ecomorphs.

Conclusion

Our results reveal a high degree of genomic differentiation between piscivorous and insectivorous populations and indicate that the large body piscivorous phenotype is likely not due to one or a few loci of large effect. Rather, the piscivore phenotype may be controlled by several loci of small effect, thus highlighting the power of whole-genome resequencing in identifying genomic regions underlying population-level phenotypic divergences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01795-9.

The Role of MicroRNAs in Influencing Body Growth and Development

Body growth and development are regulated among others by genetic and epigenetic factors. MicroRNAs (miRNAs) are epigenetic regulators of gene expression that act at the post-transcriptional level, thereby exerting a strong influence on regulatory gene networks. Increasing studies suggest the importance of miRNAs in the regulation of the growth plate and growth hormone (GH)-insulin-like growth factor (IGF) axis during the life course in a broad spectrum of animal species, contributing to longitudinal growth. This review summarizes the role of miRNAs in regulating growth in different in vitro and in vivo models acting on GH, GH receptor (GHR), IGFs, and IGF1R genes besides current knowledge in humans, and highlights that this regulatory system is of importance for growth.

Identifying therapeutic drug targets using bidirectional effect genes

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Pathway analysis for genome-wide genetic variation data: Analytic principles, latest developments, and new opportunities

Pathway analysis, also known as gene-set enrichment analysis, is a multilocus analytic strategy that integrates a priori, biological knowledge into the statistical analysis of high-throughput genetics data. Originally developed for the studies of gene expression data, it has become a powerful analytic procedure for in-depth mining of genome-wide genetic variation data. Astonishing discoveries were made in the past years, uncovering genes and biological mechanisms underlying common and complex disorders. However, as massive amounts of diverse functional genomics data accrue, there is a pressing need for newer generations of pathway analysis methods that can utilize multiple layers of high-throughput genomics data. In this review, we provide an intellectual foundation of this powerful analytic strategy, as well as an update of the state-of-the-art in recent method developments. The goal of this review is threefold: (1) introduce the motivation and basic steps of pathway analysis for genome-wide genetic variation data; (2) review the merits and the shortcomings of classic and newly emerging integrative pathway analysis tools; and (3) discuss remaining challenges and future directions for further method developments.

Functional Antibodies in COVID-19 Convalescent Plasma

In the absence of an effective vaccine or monoclonal therapeutic, transfer of convalescent plasma (CCP) was proposed early in the SARS-CoV-2 pandemic as an easily accessible therapy. However, despite the global excitement around this historically valuable therapeutic approach, results from CCP trials have been mixed and highly debated. Unlike other therapeutic interventions, CCP represents a heterogeneous drug. Each CCP unit is unique and collected from an individual recovered COVID-19 patient, making the interpretation of therapeutic benefit more complicated. While the prevailing view in the field would suggest that it is administration of neutralizing antibodies via CCP that centrally provides therapeutic benefit to newly infected COVID-19 patients, many hospitalized COVID-19 patients already possess neutralizing antibodies. Importantly, the therapeutic benefit of antibodies can extend far beyond their simple ability to bind and block infection, especially related to their ability to interact with the innate immune system. In our work we deeply profiled the SARS-CoV-2-specific Fc-response in CCP donors, along with the recipients prior to and after CCP transfer, revealing striking SARS-CoV-2 specific Fc-heterogeneity across CCP units and their recipients. However, CCP units possessed more functional antibodies than acute COVID-19 patients, that shaped the evolution of COVID-19 patient humoral profiles via distinct immunomodulatory effects that varied by pre-existing SARS-CoV-2 Spike (S)-specific IgG titers in the patients. Our analysis identified surprising influence of both S and Nucleocapsid (N) specific antibody functions not only in direct antiviral activity but also in anti-inflammatory effects. These findings offer insights for more comprehensive interpretation of correlates of immunity in ongoing large scale CCP trials and for the design of next generation therapeutic design.

Genome-wide association study across pediatric central nervous system tumors implicates shared predisposition and points to 1q25.2 (PAPPA2) and 11p12 (LRRC4C) as novel candidate susceptibility loci

INTRODUCTION

Central nervous system (CNS) tumors constitute the most common form of solid neoplasms in children, but knowledge on genetic predisposition is sparse. In particular, whether susceptibility attributable to common variants is shared across CNS tumor types in children has not been investigated. The purpose of this study was to explore potential common genetic risk variants exhibiting pleiotropic effects across pediatric CNS tumors. We also investigated whether such susceptibility differs between early and late onset of disease.

METHOD

A Danish nationwide genome-wide association study (GWAS) of 1,097 consecutive patients (< 15 years of age) with CNS tumors and a cohort of 4,745 population-based controls.

RESULTS

For both the overall cohort and patients diagnosed after the age of four, the strongest association was rs12064625 which maps to PAPPA2 at 1q25.2 (p = 3.400 × 10^-7 and 9.668 × 10^-8, respectively). PAPPA2 regulates local bioavailability of insulin-like growth factor I (IGF-I). IGF-I is fundamental to CNS development and is involved in tumorigenesis across a wide range of different cancers. For the younger children, the strongest association was provided by rs11036373 mapping to LRRC4C at 11p12 (p = 7.620 × 10^-7), which encoded protein acts as an axon guidance molecule during CNS development and has not formerly been associated with brain tumors.

DISCUSSION

This GWAS indicates shared susceptibility attributable to common variants across pediatric CNS tumor types. Variations in genetic loci with roles in CNS development appear to be involved, possibly via altered IGF-I related pathways.

Revisiting the genome-wide significance threshold for common variant GWAS

Over the last decade, GWAS meta-analyses have used a strict P-value threshold of 5 × 10-8 to classify associations as significant. Here, we use our current understanding of frequently studied traits including lipid levels, height, and BMI to revisit this genome-wide significance threshold. We compare the performance of studies using the P = 5 × 10-8 threshold in terms of true and false positive rate to other multiple testing strategies: (1) less stringent P-value thresholds, (2) controlling the FDR with the Benjamini-Hochberg and Benjamini-Yekutieli procedure, and (3) controlling the Bayesian FDR with posterior probabilities. We applied these procedures to re-analyze results from the Global Lipids and GIANT GWAS meta-analysis consortia and supported them with extensive simulation that mimics the empirical data. We observe in simulated studies with sample sizes ∼20,000 and >120,000 that relaxing the P-value threshold to 5 × 10-7 increased discovery at the cost of 18% and 8% of additional loci being false positive results, respectively. FDR and Bayesian FDR are well controlled for both sample sizes with a few exceptions that disappear under a less stringent definition of true positives and the two approaches yield similar results. Our work quantifies the value of using a relaxed P-value threshold in large studies to increase their true positive discovery but also show the excess false positive rates due to such actions in modest-sized studies. These results may guide investigators considering different thresholds in replication studies and downstream work such as gene-set enrichment or pathway analysis. Finally, we demonstrate the viability of FDR-controlling procedures in GWAS.

Whole-genome SNP markers reveal conservation status, signatures of selection, and introgression in Chinese Laiwu pigs

Laiwu pigs are a Chinese indigenous breed that is renowned for its exceptionally high intramuscular fat content (average greater than 6%), providing an excellent genetic resource for the genetic improvement of meat quality of modern commercial pigs. To uncover genetic diversity, population structure, signature of selection, and potential exotic introgression in this breed, we sampled 238 Laiwu pigs from a state-supported conservation population and genotyped these individuals using GeneSeek 80K SNP BeadChip. We then conducted in-depth population genetics analyses for the Laiwu pig in a context of 1,116 pigs from 42 Eurasian diverse breeds. First, we show that the current Laiwu population has more abundant genetic diversity than the population of 18 years ago likely due to gene flow from European commercial breeds. Both neighbor-joining (NJ) and principal component analyses indicate the introgression of European haplotypes into Laiwu pigs. The admixture analysis reveals that an average 26.66% of Laiwu genetic components are of European origin. Then, we assigned the tested individuals to different families according to their clustering patterns in the NJ tree and proposed a family-based conservation strategy to reduce the risk of inbreeding depression in Laiwu pigs. Next, we explored three statistics (ROH and iHS and EigenGWAS) to identify a list of candidate genes for fat deposition, reproduction, and growth in Laiwu pigs. Last, we detected a strong signature of introgression from European pigs into Laiwu pigs at the GPC6 locus that regulates the growth of developing long bones. Further association analyses indicate that the introgressed GPC6 haplotype likely contributed to the improvement of growth performance in Laiwu pigs. Altogether, this study not only benefits the better conservation of the Laiwu pig, but also advances our knowledge of the poorly understood effect of human-mediated introgression on phenotypic traits in Chinese indigenous pigs.

A hypomorphic variant in EYS detected by genome-wide association study contributes toward retinitis pigmentosa

The genetic basis of Japanese autosomal recessive retinitis pigmentosa (ARRP) remains largely unknown. Herein, we applied a 2-step genome-wide association study (GWAS) in 640 Japanese patients. Meta-GWAS identified three independent peaks at P < 5.0 × 10⁻⁸, all within the major ARRP gene EYS. Two of the three were each in linkage disequilibrium with a different low frequency variant (allele frequency < 0.05); a known founder Mendelian mutation (c.4957dupA, p.S1653Kfs*2) and a non-synonymous variant (c.2528 G > A, p.G843E) of unknown significance. mRNA harboring c.2528 G > A failed to restore rhodopsin mislocalization induced by morpholino-mediated knockdown of eys in zebrafish, consistent with the variant being pathogenic. c.2528 G > A solved an additional 7.0% of Japanese ARRP cases. The third peak was in linkage disequilibrium with a common non-synonymous variant (c.7666 A > T, p.S2556C), possibly representing an unreported disease-susceptibility signal. GWAS successfully unraveled genetic causes of a rare monogenic disorder and identified a high frequency variant potentially linked to development of local genome therapeutics.

Koji Nishiguchi et al. identify three genetic variants within the EYS gene that are associated with retinitis pigmentosa using a genome-wide association study. They demonstrate that one of these variants (G843E) causes retinal dysfunction in zebrafish, suggesting a causal role for EYS in retinitis pigmentosa.

Candidate Genes for Age at Menarche Are Associated With Uterine Leiomyoma

Age at menarche (AAM) is an important marker of the pubertal development and function of the hypothalamic-pituitary-ovarian system. It was reported as a possible factor for a risk of uterine leiomyoma (UL). However, while more than 350 loci for AAM have been determined by genome-wide association studies (GWASs) to date, no studies of these loci for their association with UL have been conducted so far. In this study, we analyzed 52 candidate loci for AAM for possible association with UL in a sample of 569 patients and 981 controls. The results of the study suggested that 23 out of the 52 studied polymorphisms had association with UL. Locus rs7759938 LIN28B was individually associated with the disease according to the dominant model. Twenty loci were associated with UL within 11 most significant models of intergenic interactions. Nine loci involved in 16 most significant models of interactions between single-nucleotide polymorphism (SNP), induced abortions, and chronic endometritis were associated with UL. Among the 23 loci associated with UL, 16 manifested association also with either AAM (7 SNPs) or height and/or body mass index (BMI) (13 SNPs). The above 23 SNPs and 514 SNPs linked to them have non-synonymous, regulatory, and expression quantitative trait locus (eQTL) significance for 35 genes, which play roles in the pathways related to development of the female reproductive organs and hormone-mediated signaling [false discovery rate (FDR) ≤ 0.05]. This is the first study reporting associations of candidate genes for AAM with UL.

Pathways involved in pony body size development

BACKGROUND

The mechanism of body growth in mammals is poorly understood. Here, we investigated the regulatory networks involved in body growth through transcriptomic analysis of pituitary and epiphyseal tissues of smaller sized Debao ponies and Mongolian horses at the juvenile and adult stages.

RESULTS

We found that growth hormone receptor (GHR) was expressed at low levels in long bones, although growth hormone (GH) was highly expressed in Debao ponies compared with Mongolian horses. Moreover, significant downregulated of the GHR pathway components m-RAS and ATF3 was found in juvenile ponies, which slowed the proliferation of bone osteocytes. However, WNT2 and PLCβ2 were obviously upregulated in juvenile Debao ponies, which led to premature mineralization of the bone extracellular matrix. Furthermore, we found that the WNT/Ca²⁺ pathway may be responsible for regulating body growth. GHR was demonstrated by q-PCR and Western blot analyses to be expressed at low levels in long bones of Debao ponies. Treatment with WNT antagonistI decreased the expression of WNT pathway components (P < 0.05) in vitro. Transduction of ATDC5 cells with a GHR-RNAi lentiviral vector decreased the expression of the GHR pathway components (P < 0.05). Additionally, the expression of the IGF-1 gene in the liver was lower in Debao ponies than in Mongolian horses at the juvenile and adult stages. Detection of plasma hormone concentrations showed that Debao ponies expressed higher levels of IGF-1 as juveniles and higher levels of GH as adults than Mongolian horses, indicating that the hormone regulation in Debao ponies differs from that in Mongolian horses.

CONCLUSION

Our work provides insights into the genetic regulation of short stature growth in mammals and can provide useful information for the development of therapeutic strategies for small size.

Sporadic Parkinson's Disease Potential Risk Loci Identified in Han Ancestry of Chinese Mainland

Recent investigations demonstrated that genetic factors might play an important role in sporadic Parkinson's disease (sPD). To clarify the specific loci susceptibility to sPD, we analyze the relationship between 30 candidate single nucleotide polymorphisms (SNPs) and sPD in the population of Han ancestry from Chinese mainland (HACM) by using genome-wide association study, sequenom massARRAY, DNA sequence, and biological information analysis. Results showed that the subjects carrying the T allele of rs863108 and rs28499371 exhibited a decreased risk for sPD. The subjects carrying the T allele of rs80315856 exhibited an increased risk for sPD. The A/T genotype of rs863108 and the C/T genotype of rs28499371 were a potential increased risk for sPD, and the G/T genotype of rs80315856 and T/T genotype of rs2270568 were a potential decreased risk for sPD. The minor allele frequency (MAF) of rs80315856 and rs2270568 was higher in sPD. The T allele of rs80315856 and rs2270568 might be a risk locus for sPD. Our data suggested that the alteration of these SNPs might play some roles through changing/affecting LINC01524/LOC105372666, DMRT2/SMARCA2, PLEKHN1, and FLJ23172/FNDC3B genes in the pathogenesis of sPD.

STAB: a spatio-temporal cell atlas of the human brain

The human brain is the most complex organ consisting of billions of neuronal and non-neuronal cells that are organized into distinct anatomical and functional regions. Elucidating the cellular and transcriptome architecture underlying the brain is crucial for understanding brain functions and brain disorders. Thanks to the single-cell RNA sequencing technologies, it is becoming possible to dissect the cellular compositions of the brain. Although great effort has been made to explore the transcriptome architecture of the human brain, a comprehensive database with dynamic cellular compositions and molecular characteristics of the human brain during the lifespan is still not available. Here, we present STAB (a Spatio-Temporal cell Atlas of the human Brain), a database consists of single-cell transcriptomes across multiple brain regions and developmental periods. Right now, STAB contains single-cell gene expression profiling of 42 cell subtypes across 20 brain regions and 11 developmental periods. With STAB, the landscape of cell types and their regional heterogeneity and temporal dynamics across the human brain can be clearly seen, which can help to understand both the development of the normal human brain and the etiology of neuropsychiatric disorders. STAB is available at http://stab.comp-sysbio.org.

Genome-wide identification of major genes and genomic prediction using high-density and text-mined gene-based SNP panels in Hanwoo (Korean cattle)

It was hypothesized that single-nucleotide polymorphisms (SNPs) extracted from text-mined genes could be more tightly related to causal variant for each trait and that differentially weighting of this SNP panel in the GBLUP model could improve the performance of genomic prediction in cattle. Fitting two GRMs constructed by text-mined SNPs and SNPs except text-mined SNPs from 777k SNPs set (exp_777K) as different random effects showed better accuracy than fitting one GRM (Im_777K) for six traits (e.g. backfat thickness: + 0.002, eye muscle area: + 0.014, Warner–Bratzler Shear Force of semimembranosus and longissimus dorsi: + 0.024 and + 0.068, intramuscular fat content of semimembranosus and longissimus dorsi: + 0.008 and + 0.018). These results can suggest that attempts to incorporate text mining into genomic predictions seem valuable, and further study using text mining can be expected to present the significant results.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

Disorders caused by genetic defects associated with GH-dependent genes: PAPPA2 defects

Growth hormone (GH) and its mediator, insulin-like growth factor-1 (IGF-1), have long been recognized as central to human growth physiology. IGF-1 is known to complex with IGF binding proteins as well as with the acid labile subunit (ALS) in order to prolong its half-life in circulation. Factors regulating the bioavailability of IGF-1 (i.e. the balance between free and bound IGF-1) were less well understood. Recently, pregnancy-associated plasma protein-A2 (PAPP-A2) was discovered as a protease which specifically cleaves IGF-binding protein (IGFBP)-3 and -5. PAPP-A2 deficient patients present with characteristic findings including growth failure, elevated total IGF-1 and -2, IGFBPs, and ALS, but decreased percentage of free to total IGF-1. Additionally, patients with PAPP-A2 deficiency have impairments in glucose metabolism and bone mineral density (BMD). Treatment with recombinant human IGF-1 (rhIGF-1) improved height SD scores, growth velocity, body composition, and dysglycemia. Mouse models recapitulate many of the human findings of PAPP-A2 deficiency. This review summarizes the function of PAPP-A2 and its contribution to the GH-IGF axis through an examination of PAPP-A2 deficient patients and mouse models, thereby emphasizing the importance of the regulation of IGF-1 bioavailability in human growth.

Expression of acid-labile subunit (ALS) in developing and adult zebrafish and its role in dorso-ventral patterning during development

Mammalian acid-labile subunit (ALS) is a serum protein that binds binary complexes between Insulin-like growth factors (IGFs) and Insulin-like growth factor-binding proteins (IGFBPs) extending their half-life and keeping them in the vasculature. Human ALS deficiency (ACLSD), due to homozygous or compound heterozygous mutations in IGFALS, leads to moderate short stature with reduced levels of IGF-I and IGFBP-3. There is only one corresponding zebrafish ortholog gene and it has not yet been studied. In this study we elucidate the role of igfals during zebrafish development. In zebrafish embryos igfals mRNA is expressed throughout development, mainly in the brain and subsequently also in the gut and swimbladder. To determine its role during development, we knocked down igfals gene product using morpholinos (MOs). Igfals morphant embryos displayed dorsalization in different degrees of severity, including a shortened trunk and loss of tail. Furthermore, co-injection of human IGFALS (hIGFALS) mRNA was able to rescue the MO-induced phenotype. Finally, overexpression of either hIGFALS or zebrafish igfals (zigfals) mRNA leads to ventralization of embryos including a reduced head and enlarged tail. These findings suggest that als plays an important role in dorso-ventral patterning during zebrafish development.

Growth change in Polish women: Reduction of the secular trends?

The aim of the study was to analyse changes in the average height of adult Polish women born in 1931-2001 in the aspect of dynamically changing economic and socio-economic conditions of the living environment. An ethnically homogeneous group of 6,028 adult women from large Polish cities, born in 1931-2001, living between 1931 and 2020, were examined using the same research methods and research equipment. All women were divided into eight birth cohorts. The Kruskal-Wallis test and multiple regression analyses were used. Root Mean Square Standardized Effect (RMSSE), critical value of the test, and test power were calculated. The average height of women born during 70 years of the study increased by 9.63 cm, from 158.22 cm (SD = 5.57 cm) to 167.85 cm (SD = 6.91 cm) (H = 1084.84, p<0.001). The intensity of the intergenerational trend in subsequent cohorts of years of birth varied strongly between decades, averaging 1.34 cm/decade. The body height in women increased significantly up to the height of those born between 1970 and 1979 and then the trend weakened noticeably, although it remained positive. The observed secular trend confirms positive changes in the standard of living of Polish women between 1931 and 2020. Improving living conditions allow people to fully achieve their genetically determined growth potential.

Whole Transcriptome Analysis Identifies the Taxonomic Status of a New Chinese Native Cattle Breed and Reveals Genes Related to Body Size

Wandong (WD) cattle has recently been identified as a new Chinese native cattle breed by the National Commission for Livestock and Poultry Genetic Resources. The population size of this breed is less than 10,000. WD cattle and Dabieshan (DB) cattle are sympatric but are raised in different ecological environments, on mountains and plains, respectively, and the body sizes of these two breeds are markedly different. Blood samples were obtained from 8 adult female WD cattle and 7 adult female DB cattle (24 months old). The total RNA was extracted from leukocyte cells, and sequencing experiments were conducted on the Illumina HiSeq^TM 4000 platform. After the removal of one outlier sample from the WD cattle breed as determined by principal component analysis (PCA), phylogenetic and population structure analyses indicated that WD and DB cattle formed a distinct Central China cattle group and showed evidence of hybridization between Bos. taurus and Bos. indicus. The immune-regulator CD48 (P = 1.3E-6) was associated with breed-specific traits according to loss-of-function variant enrichment analysis. In addition, 113 differentially expressed genes were identified between the two breeds, many of which are associated with the regulation of body growth, which is the major difference between the two breeds. This study showed that WD cattle belong to the group of hybrids between Bos. Taurus and Bos. indicus, and one novel gene associated with breed traits and multiple differentially expressed genes between these two closely related breeds was identified. The results provide insights into the genetic mechanisms that underlie economically important traits, such as body size, in cattle.

Population-Matched Transcriptome Prediction Increases TWAS Discovery and Replication Rate

Most genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) focus on European populations; however, these results cannot always be accurately applied to non-European populations due to genetic architecture differences. Using GWAS summary statistics in the Population Architecture using Genomics and Epidemiology study, which comprises ∼50,000 Hispanic/Latinos, African Americans, Asians, Native Hawaiians, and Native Americans, we perform TWAS to determine gene-trait associations. We compared results using three transcriptome prediction models derived from Multi-Ethnic Study of Atherosclerosis populations: the African American and Hispanic/Latino (AFHI) model, the European (EUR) model, and the African American, Hispanic/Latino, and European (ALL) model. We identified 240 unique significant trait-associated genes. We found more significant, colocalized genes that replicate in larger cohorts when applying the AFHI model than the EUR or ALL model. Thus, TWAS with population-matched transcriptome models have more power for discovery and replication, demonstrating the need for more transcriptome studies in diverse populations.

•

TWAS mechanistically extends GWAS findings in diverse populations

•

Population-matched transcriptome models detect more replicable associations

•

Colocalization shows GWAS variants likely act through gene expression regulation

•

More GWAS and transcriptome modeling in diverse populations are needed

Multiple-Tissue Integrative Transcriptome-Wide Association Studies Discovered New Genes Associated With Amyotrophic Lateral Sclerosis

Genome-wide association studies (GWAS) have identified multiple causal genes associated with amyotrophic lateral sclerosis (ALS); however, the genetic architecture of ALS remains completely unknown and a large number of causal genes have yet been discovered. To full such gap in part, we implemented an integrative analysis of transcriptome-wide association study (TWAS) for ALS to prioritize causal genes with summary statistics from 80,610 European individuals and employed 13 GTEx brain tissues as reference transcriptome panels. The summary-level TWAS analysis with single brain tissue was first undertaken and then a flexible p-value combination strategy, called summary data-based Cauchy Aggregation TWAS (SCAT), was proposed to pool association signals from single-tissue TWAS analysis while protecting against highly positive correlation among tests. Extensive simulations demonstrated SCAT can produce well-calibrated p-value for the control of type I error and was often much more powerful to identify association signals across various scenarios compared with single-tissue TWAS analysis. Using SCAT, we replicated three ALS-associated genes (i.e., ATXN3, SCFD1, and C9orf72) identified in previous GWASs and discovered additional five genes (i.e., SLC9A8, FAM66D, TRIP11, JUP, and RP11-529H20.6) which were not reported before. Furthermore, we discovered the five associations were largely driven by genes themselves and thus might be new genes which were likely related to the risk of ALS. However, further investigations are warranted to verify these results and untangle the pathophysiological function of the genes in developing ALS.