Identification of additional body weight QTLs in the Berlin Fat Mouse BFMI861 lines using time series data

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and metabolic syndrome. The sublines BFMI861-S1 and BFMI861-S2 differ in weight despite high genetic similarity and a shared obesity-related locus. This study focused on identifying additional body weight quantitative trait loci (QTLs) by analyzing weekly weight measurements in a male population of the advanced intercross line BFMI861-S1 x BFMI861-S2. QTL analysis, utilizing 200 selectively genotyped mice (GigaMUGA) and 197 males genotyped for top SNPs, revealed a genome-wide significant QTL on Chr 15 (68.46 to 81.40 Mb) for body weight between weeks 9 to 20. Notably, this QTL disappeared (weeks 21 to 23) and reappeared (weeks 24 and 25) coinciding with a diet change. Additionally, a significant body weight QTL on Chr 16 (3.89 to 22.79 Mb) was identified from weeks 6 to 25. Candidate genes, including Gpt, Cbx6, Apol6, Apol8, Sun2 (Chr 15) and Trap1, Rrn3, Mapk1 (Chr 16), were prioritized. This study unveiled two additional body weight QTLs, one of which is novel and responsive to diet changes. These findings illuminate genomic regions influencing weight in BFMI and emphasize the utility of time series data in uncovering novel genetic factors.

Examination of a novel expression-based gene-SNP annotation strategy to identify tissue-specific contributions to heritability in multiple traits

Complex traits show clear patterns of tissue-specific expression influenced by single nucleotide polymorphisms (SNPs), yet current strategies aggregate SNP effects to genes by employing simple physical proximity-based windows. Here, we examined whether incorporating SNPs with effects on tissue-specific cis-expression would improve our ability to detect trait-relevant tissues across 31 complex traits using stratified linkage disequilibrium score regression (S-LDSC). We found that a physical proximity annotation produced more significant tissue enrichments and larger S-LDSC regression coefficients, as compared to an expression-based annotation. Furthermore, we showed that our expression-based annotation did not outperform an annotation strategy in which an equal number of randomly chosen SNPs were annotated to genes within the same genomic window, suggesting extensive redundancy among SNP effect estimates due to linkage disequilibrium. That said, current sample sizes limit estimation of cis-genetic SNP effects; therefore, we recommend reexamination of the expression-based annotation when larger tissue-specific expression datasets become available. To examine the influence of sample size, we used a large whole blood eQTL reference panel (N = 31,684) applying a similar expression-based annotation strategy. We found that significant cis-expression QTLs in whole blood did not outperform the physical proximity annotation when estimating tissue-specific SNP heritability enrichment for either high- or low-density lipoprotein phenotypes but performed similarly for inflammatory bowel disease. Finally, we report new and updated tissue enrichment estimates across 31 complex traits, such as significant heritability enrichment of the frontal cortex for cognitive performance, educational attainment, and intelligence, providing further evidence of this structure's importance in higher cognitive function.

Mendelian randomization reveals apolipoprotein B shortens healthspan and possibly increases risk for Alzheimer's disease

Apolipoprotein B-100 (APOB) is a component of fat- and cholesterol-transporting molecules in the bloodstream. It is the main lipoprotein in low-density lipoprotein cholesterol (LDL) and has been implicated in conditions that end healthspan (the interval between birth and onset of chronic disease). However, APOB's direct relationship with healthspan remains uncertain. With Mendelian randomization, we show that higher levels of APOB and LDL shorten healthspan in humans. Multivariable Mendelian randomization of APOB and LDL on healthspan suggests that the predominant trait accounting for the relationship is APOB. In addition, we provide preliminary evidence that APOB increases risk for Alzheimer's disease, a condition that ends healthspan. If these relationships are causal, they suggest that interventions to improve healthspan in aging populations could include strategies targeting APOB. Ultimately, given that more than 44 million people currently suffer from Alzheimer's disease worldwide, such interventions are needed.

Mean performances, character associations and multi-environmental evaluation of chilli landraces in north western Himalayas

Even though many varieties have been recommended across agro-climate zones of Himachal Pradesh, yet the information on stability is lacking in this State. Hence, the present investigation was carried out to identify high yielding stable genotypes among various pre-adapted landraces. The material consists of 20 chilli landraces including check i.e. DKC-8. The experiment was laid out in a RCBD. The data were recorded and analyzed to work out mean performances and the inferences were drawn for parameters of variability, correlation coefficients, path coefficients and stability analysis. As per mean performances, CS7 and CS9 were earliest in flowering, CS13 is earliest in days to ripe maturity, CS10 had highest plant height and CS9 had highest average fruit weight and ripe fruit yield plant-1. High PCV and GCV were recorded for ripe fruit yield plant-1. Heritability and genetic advance were recorded maximum for plant height in summer seasons and were recorded maximum for number of ripe fruits plant-1 in winter season. Correlation coefficients showed that number of ripe fruits plant-1 and average ripe fruit weight were positively and significantly correlated with ripe fruit yield plant-1. Path coefficient analysis in summer and winter seasons showed that average ripe fruit weight had the highest positive direct effect on ripe fruit yield plant-1. The pooled data over environments were analyzed to estimate the interaction effects between genotypes × environment. The mean sum of squares due to genotypes, environments and genotypes × environment interaction were significant for all the characteristics. CS1, CS3, CS6, CS10, CS13, CS15 were adapted to all environments, CS7 and CS9 were specifically adapted to favourable environment and CS2 was specifically adapted to unfavorable environment for 50% flowering, landraces CS1, CS2 and CS3were well adapted to all environments for ripe maturity whereas landraces CS6, CS10 and CS19 were well adapted to all environment for number of ripe fruit and ripe fruit yield.

Genome-wide association in Drosophila identifies a role for Piezo and Proc-R in sleep latency

Sleep latency, the amount of time that it takes an individual to fall asleep, is a key indicator of sleep need. Sleep latency varies considerably both among and within species and is heritable, but lacks a comprehensive description of its underlying genetic network. Here we conduct a genome-wide association study of sleep latency. Using previously collected sleep and activity data on a wild-derived population of flies, we calculate sleep latency, confirming significant, heritable genetic variation for this complex trait. We identify 520 polymorphisms in 248 genes contributing to variability in sleep latency. Tests of mutations in 23 candidate genes and additional putative pan-neuronal knockdown of 9 of them implicated CG44153, Piezo, Proc-R and Rbp6 in sleep latency. Two large-effect mutations in the genes Proc-R and Piezo were further confirmed via genetic rescue. This work greatly enhances our understanding of the genetic factors that influence variation in sleep latency.

Genetic architecture of individual meiotic crossover rate and distribution in Atlantic Salmon

Meiotic recombination through chromosomal crossovers ensures proper segregation of homologous chromosomes during meiosis, while also breaking down linkage disequilibrium and shuffling alleles at loci located on the same chromosome. Rates of recombination can vary between species, but also between and within individuals, sex and chromosomes within species. Indeed, the Atlantic salmon genome is known to have clear sex differences in recombination with female biased heterochiasmy and markedly different landscapes of crossovers between males and females. In male meiosis, crossovers occur strictly in the telomeric regions, whereas in female meiosis crossovers tend to occur closer to the centromeres. However, little is known about the genetic control of these patterns and how this differs at the individual level. Here, we investigate genetic variation in individual measures of recombination in > 5000 large full-sib families of a Norwegian Atlantic salmon breeding population with high-density SNP genotypes. We show that females had 1.6 × higher crossover counts (CC) than males, with autosomal linkage maps spanning a total of 2174 cM in females and 1483 cM in males. However, because of the extreme telomeric bias of male crossovers, female recombination is much more important for generation of new haplotypes with 8 × higher intra-chromosomal genetic shuffling than males. CC was heritable in females (h2 = 0.11) and males (h2 = 0.10), and shuffling was also heritable in both sex but with a lower heritability in females (h2 = 0.06) than in males (h2 = 0.11). Inter-sex genetic correlations for both traits were close to zero, suggesting that rates and distribution of crossovers are genetically distinct traits in males and females, and that there is a potential for independent genetic change in both sexes in the Atlantic Salmon. Together, these findings give novel insights into the genetic architecture of recombination in salmonids and contribute to a better understanding of how rates and distribution of recombination may evolve in eukaryotes more broadly.

New insights into the genetic predisposition of brucellosis and its effect on the gut and vaginal microbiota in goats

Goats contribute significantly to the global food security and industry. They constitute a main supplier of meat and milk for large proportions of people in Egypt and worldwide. Brucellosis is a zoonotic infectious disease that causes a significant economic loss in animal production. A case-control genome-wide association analysis (GWAS) was conducted using the infectious status of the animal as a phenotype. The does that showed abortion during the last third period of pregnancy and which were positive to both rose bengal plate and serum tube agglutination tests, were considered as cases. Otherwise, they were considered as controls. All animals were genotyped using the Illumina 65KSNP BeadChip. Additionally, the diversity and composition of vaginal and fecal microbiota in cases and controls were investigated using PCR-amplicone sequencing of the V4 region of 16S rDNA. After applying quality control criteria, 35,818 markers and 66 does were available for the GWAS test. The GWAS revealed a significantly associated SNP (P = 5.01 × 10-7) located on Caprine chromosome 15 at 29 megabases. Four other markers surpassed the proposed threshold (P = 2.5 × 10-5). Additionally, fourteen genomic regions accounted for more than 0.1% of the variance explained by all genome windows. Corresponding markers were located within or in close vicinity to several candidate genes, such as ARRB1, RELT, ATG16L2, IGSF21, UBR4, ULK1, DCN, MAPB1, NAIP, CD26, IFIH1, NDFIP2, DOK4, MAF, IL2RB, USP18, ARID5A, ZAP70, CNTN5, PIK3AP1, DNTT, BLNK, and NHLRC3. These genes play important roles in the regulation of immune responses to the infections through several biological pathways. Similar vaginal bacterial community was observed in both cases and controls while the fecal bacterial composition and diversity differed between the groups (P < 0.05). Faeces from the control does showed a higher relative abundance of the phylum Bacteroidota compared to cases (P < 0.05), while the latter showed more Firmicutes, Spirochaetota, Planctomycetota, and Proteobacteria. On the genus level, the control does exhibited higher abundances of Rikenellaceae RC9 gut group and Christensenellaceae R-7 group (P < 0.05), while the infected does revealed higher Bacteroides, Alistipes, and Prevotellaceae UCG-003 (P < 0.05). This information increases our understanding of the genetics of the susceptibility to Brucella in goats and may be useful in breeding programs and selection schemes that aim at controlling the disease in livestock.

Dynamic QTL mapping revealed primarily the genetic structure of photosynthetic traits in castor (Ricinus communis L.)

High photosynthetic efficiency is the basis of high biomass and high harvest index in castor (Ricinus communis L.). Understanding the genetic law of photosynthetic traits will facilitate the breeding for high photosynthetic efficiency. In this study, the dynamic QTL mapping was performed with the populations F2 and BC1 derived from 2 parents with significant difference in net photosynthetic rate (Pn) at 3 stages, in order to reveal the genetic structure of photosynthetic traits. In F2 population, 26 single-locus QTLs were identified, including 3/3/1 (the QTL number at stage I/II/III, the same below), 1/2/0, 1/2/2, 1/3/1, 0/1/1, and 1/1/2 QTLs conferring Pn, water use efficiency (Wue), transpiration rate (Tr), stomatal conductance (Gs), intercellular CO2 concentration (Ci) and chlorophyll content (Cc), with a phenotypic variation explained (PVE) of 8.40%/8.91%/6.17%, 5.36%/31.74%/0, 7.31%/12.80%/15.15%, 1.60%/6.44%/0.02%, 0/1.10%/0.70% and 2.77%/3.96%/6.50% respectively. And 53 epistatic QTLs (31 pairs) were identified, including 2/2/5, 5/6/3, 4/4/2, 6/3/2, 3/2/0 and 4/0/0 ones conferring the above 6 traits, with a PVE of 6.52%/6.47%/19.04%, 16.72%/15.67%/14.12%, 18.57%/15.58%/7.34%, 21.72%/8.52%/7.13%, 13.33%/4.94%/0 and 7.84%/0/0 respectively. The QTL mapping results in BC1 population were consistent with those in F2 population, except fewer QTLs detected. Most QTLs identified were minor-effect ones, only a few were main-effect ones (PVE > 10%), focused on 2 traits, Wue and Tr, such as qWue1.1, qWue1.2, FqTr1.1, FqTr6, BqWue1.1 and BqTr3; The epistatic effects, especially those related to the dominance effects were the main genetic component of photosynthetic traits, and all the epistatic QTLs had no single-locus effects except qPn1.2, FqGs1.2, FqCi1.2 and qCc3.2; The detected QTLs underlying each trait varied at different stages except stable QTLs qGs1.1, detected at 3 stages, qWue2, qTr1.2 and qCc3.2, detected at 2 stages; 6 co-located QTLs were identified, each of which conferring 2-5 different traits, demonstrated the gene pleiotropy between photosynthetic traits; 2 QTL clusters, located within the marker intervals RCM1842-RCM1335 and RCM523-RCM83, contained 15/5 (F2/BC1) and 4/4 (F2/BC1) QTLs conferring multiple traits, including co-located QTLs and main-effect QTLs. The above results provided new insights into the genetic structure of photosynthetic traits and important references for the high photosynthetic efficiency breeding in castor plant.

Biobank-scale methods and projections for sparse polygenic prediction from machine learning

In this paper we characterize the performance of linear models trained via widely-used sparse machine learning algorithms. We build polygenic scores and examine performance as a function of training set size, genetic ancestral background, and training method. We show that predictor performance is most strongly dependent on size of training data, with smaller gains from algorithmic improvements. We find that LASSO generally performs as well as the best methods, judged by a variety of metrics. We also investigate performance characteristics of predictors trained on one genetic ancestry group when applied to another. Using LASSO, we develop a novel method for projecting AUC and correlation as a function of data size (i.e., for new biobanks) and characterize the asymptotic limit of performance. Additionally, for LASSO (compressed sensing) we show that performance metrics and predictor sparsity are in agreement with theoretical predictions from the Donoho-Tanner phase transition. Specifically, a future predictor trained in the Taiwan Precision Medicine Initiative for asthma can achieve an AUC of [Formula: see text] and for height a correlation of [Formula: see text] for a Taiwanese population. This is above the measured values of [Formula: see text] and [Formula: see text], respectively, for UK Biobank trained predictors applied to a European population.

Longitudinal telomere dynamics within natural lifespans of a wild bird

Telomeres, the nucleotide sequences that protect the ends of eukaryotic chromosomes, shorten with each cell division and telomere loss may be influenced by environmental factors. Telomere length (TL) decreases with age in several species, but little is known about the sources of genetic and environmental variation in the change in TL (∆TL) in wild animals. In this study, we tracked changes in TL throughout the natural lifespan (from a few months to almost 9 years) of free-living house sparrows (Passer domesticus) in two different island populations. TL was measured in nestlings and subsequently up to four times during their lifetime. TL generally decreased with age (senescence), but we also observed instances of telomere lengthening within individuals. We found some evidence for selective disappearance of individuals with shorter telomeres through life. Early-life TL positively predicted later-life TL, but the within-individual repeatability in TL was low (9.2%). Using genetic pedigrees, we found a moderate heritability of ∆TL (h² = 0.21), which was higher than the heritabilities of early-life TL (h² = 0.14) and later-life TL measurements (h² = 0.15). Cohort effects explained considerable proportions of variation in early-life TL (60%), later-life TL (53%), and ∆TL (37%), which suggests persistent impacts of the early-life environment on lifelong telomere dynamics. Individual changes in TL were independent of early-life TL. Finally, there was weak evidence for population differences in ∆TL that may be linked to ecological differences in habitat types. Combined, our results show that individual telomere biology is highly dynamic and influenced by both genetic and environmental variation in natural conditions.

Identifying polymorphic cis-regulatory variants as risk markers for lung carcinogenesis and chemotherapy responses in tobacco smokers from eastern India

Aberrant expression of xenobiotic metabolism and DNA repair genes is critical to lung cancer pathogenesis. This study aims to identify the cis-regulatory variants of the genes modulating lung cancer risk among tobacco smokers and altering their chemotherapy responses. From a list of 2984 SNVs, prioritization and functional annotation revealed 22 cis-eQTLs of 14 genes within the gene expression-correlated DNase I hypersensitive sites using lung tissue-specific ENCODE, GTEx, Roadmap Epigenomics, and TCGA datasets. The 22 cis-regulatory variants predictably alter the binding of 44 transcription factors (TFs) expressed in lung tissue. Interestingly, 6 reported lung cancer-associated variants were found in linkage disequilibrium (LD) with 5 prioritized cis-eQTLs from our study. A case-control study with 3 promoter cis-eQTLs (p < 0.01) on 101 lung cancer patients and 401 healthy controls from eastern India with confirmed smoking history revealed an association of rs3764821 (ALDH3B1) (OR = 2.53, 95% CI = 1.57-4.07, p = 0.00014) and rs3748523 (RAD52) (OR = 1.69, 95% CI = 1.17-2.47, p = 0.006) with lung cancer risk. The effect of different chemotherapy regimens on the overall survival of lung cancer patients to the associated variants showed that the risk alleles of both variants significantly decreased (p < 0.05) patient survival.

A method for an unbiased estimate of cross-ancestry genetic correlation using individual-level data

Cross-ancestry genetic correlation is an important parameter to understand the genetic relationship between two ancestry groups. However, existing methods cannot properly account for ancestry-specific genetic architecture, which is diverse across ancestries, producing biased estimates of cross-ancestry genetic correlation. Here, we present a method to construct a genomic relationship matrix (GRM) that can correctly account for the relationship between ancestry-specific allele frequencies and ancestry-specific allelic effects. Through comprehensive simulations, we show that the proposed method outperforms existing methods in the estimations of SNP-based heritability and cross-ancestry genetic correlation. The proposed method is further applied to anthropometric and other complex traits from the UK Biobank data across ancestry groups. For obesity, the estimated genetic correlation between African and European ancestry cohorts is significantly different from unity, suggesting that obesity is genetically heterogenous between these two ancestries.

The effect of the 7R allele at the DRD4 locus on risk tolerance is independent of background risk in Senegalese fishermen

It has been shown that living in risky environments, as well as having a risky occupation, can moderate risk-tolerance. Despite the involvement of dopamine in the expectation of reward described by neurobiologists, a GWAS study was not able to demonstrate a genetic contribution of genes involved in the dopaminergic pathway in risk attitudes and gene candidate studies gave contrasting results. We test the possibility that a genetic effect of the DRD4-7R allele in risk-taking behavior could be modulated by environmental factors. We show that the increase in risk-tolerance due to the 7R allele is independent of the environmental risk in two populations in Northern Senegal, one of which is exposed to a very high risk due to dangerous fishing.

How people wake up is associated with previous night's sleep together with physical activity and food intake

How people wake up and regain alertness in the hours after sleep is related to how they are sleeping, eating, and exercising. Here, in a prospective longitudinal study of 833 twins and genetically unrelated adults, we demonstrate that how effectively an individual awakens in the hours following sleep is not associated with their genetics, but instead, four independent factors: sleep quantity/quality the night before, physical activity the day prior, a breakfast rich in carbohydrate, and a lower blood glucose response following breakfast. Furthermore, an individual's set-point of daily alertness is related to the quality of their sleep, their positive emotional state, and their age. Together, these findings reveal a set of non-genetic (i.e., not fixed) factors associated with daily alertness that are modifiable.

Identifying nootropic drug targets via large-scale cognitive GWAS and transcriptomics

Broad-based cognitive deficits are an enduring and disabling symptom for many patients with severe mental illness, and these impairments are inadequately addressed by current medications. While novel drug targets for schizophrenia and depression have emerged from recent large-scale genome-wide association studies (GWAS) of these psychiatric disorders, GWAS of general cognitive ability can suggest potential targets for nootropic drug repurposing. Here, we (1) meta-analyze results from two recent cognitive GWAS to further enhance power for locus discovery; (2) employ several complementary transcriptomic methods to identify genes in these loci that are credibly associated with cognition; and (3) further annotate the resulting genes using multiple chemoinformatic databases to identify "druggable" targets. Using our meta-analytic data set (N = 373,617), we identified 241 independent cognition-associated loci (29 novel), and 76 genes were identified by 2 or more methods of gene identification. Actin and chromatin binding gene sets were identified as novel pathways that could be targeted via drug repurposing. Leveraging our transcriptomic and chemoinformatic databases, we identified 16 putative genes targeted by existing drugs potentially available for cognitive repurposing.

A combined GWAS approach reveals key loci for socially-affected traits in Yorkshire pigs

Socially affected traits in pigs are controlled by direct genetic effects and social genetic effects, which can make elucidation of their genetic architecture challenging. We evaluated the genetic basis of direct genetic effects and social genetic effects by combining single-locus and haplotype-based GWAS on imputed whole-genome sequences. Nineteen SNPs and 25 haplotype loci are identified for direct genetic effects on four traits: average daily feed intake, average daily gain, days to 100 kg and time in feeder per day. Nineteen SNPs and 11 haplotype loci are identified for social genetic effects on average daily feed intake, average daily gain, days to 100 kg and feeding speed. Two significant SNPs from single-locus GWAS (SSC6:18,635,874 and SSC6:18,635,895) are shared by a significant haplotype locus with haplotype alleles 'GGG' for both direct genetic effects and social genetic effects in average daily feed intake. A candidate gene, MT3, which is involved in growth, nervous, and immune processes, is identified. We demonstrate the genetic differences between direct genetic effects and social genetic effects and provide an anchor for investigating the genetic architecture underlying direct genetic effects and social genetic effects on socially affected traits in pigs.

The limits of normal approximation for adult height

Adult height inspired the first biometrical and quantitative genetic studies and is a test-case trait for understanding heritability. The studies of height led to formulation of the classical polygenic model, that has a profound influence on the way we view and analyse complex traits. An essential part of the classical model is an assumption of additivity of effects and normality of the distribution of the residuals. However, it may be expected that the normal approximation will become insufficient in bigger studies. Here, we demonstrate that when the height of hundreds of thousands of individuals is analysed, the model complexity needs to be increased to include non-additive interactions between sex, environment and genes. Alternatively, the use of log-normal approximation allowed us to still use the additive effects model. These findings are important for future genetic and methodologic studies that make use of adult height as an exemplar trait.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Selection on heritable social network positions is context-dependent in Drosophila melanogaster

Social group structure is highly variable and can be important for nearly every aspect of behavior and its fitness consequences. Group structure can be modeled using social network analysis, but we know little about the evolutionary factors shaping and maintaining variation in how individuals are embedded within their networks (i.e., network position). While network position is a pervasive target of selection, it remains unclear whether network position is heritable and can respond to selection. Furthermore, it is unclear how environmental factors interact with genotypic effects on network positions, or how environmental factors shape selection on heritable network structure. Here we show multiple measures of social network position are heritable, using replicate genotypes and replicate social groups of Drosophila melanogaster flies. Our results indicate genotypic differences in network position are largely robust to changes in the environment flies experience, though some measures of network position do vary across environments. We also show selection on multiple network position metrics depends on the environmental context they are expressed in, laying the groundwork for better understanding how spatio-temporal variation in selection contributes to the evolution of variable social group structure.

The sources of variation for individual prey-to-predator size ratios

The relative body size at which predators are willing to attack prey, a key trait for predator-prey interactions, is usually considered invariant. However, this ratio can vary widely among individuals or populations. Identifying the range and origin of such variation is key to understanding the strength and constraints on selection in both predators and prey. Still, these sources of variation remain largely unknown. We filled this gap by measuring the genetic, maternal and environmental variation of the maximum prey-to-predator size ratio (PPSR_max) in juveniles of the wolf spider Lycosa fasciiventris using a paternal half-sib split-brood design, in which each male was paired with two females and the offspring reared in two food environments: poor and rich. Each juvenile spider was then sequentially offered crickets of decreasing size and the maximum prey size killed was determined. We also measured body size and body condition of spiders upon emergence and just before the trial. We found low, but significant heritability (h² = 0.069) and dominance and common environmental variance (d² + 4c² = 0.056). PPSR_max was also partially explained by body condition (during trial) but there was no effect of the rearing food environment. Finally, a maternal correlation between body size early in life and PPSR_max indicated that offspring born larger were less predisposed to feed on larger prey later in life. Therefore, PPSR_max, a central trait in ecosystems, can vary widely and this variation is due to different sources, with important consequences for changes in this trait in the short and long terms.

Endocranial volume is variable and heritable, but not related to fitness, in a free-ranging primate

Large relative brain size is a defining characteristic of the order Primates. Arguably, this can be attributed to selection for behavioral aptitudes linked to a larger brain size. In order for selection of a trait to occur, the trait must vary, that variation must be heritable, and enhance fitness. In this study, we use a quantitative genetic approach to investigate the production and maintenance of variation in endocranial volume in a population of free-ranging rhesus macaques. We measured the endocranial volume and body mass proxies of 542 rhesus macaques from Cayo Santiago. We investigated variation in endocranial volume within and between sexes. Using a genetic pedigree, we estimated heritability of absolute and relative endocranial volume, and selection gradients of both traits as well as estimated body mass in the sample. Within this population, both absolute and relative endocranial volume display variation and sexual dimorphism. Both absolute and relative endocranial volume are highly heritable, but we found no evidence of selection on absolute or relative endocranial volume. These findings suggest that endocranial volume is not undergoing selection, or that we did not detect it because selection is neither linear nor quadratic, or that we lacked sufficient sample sizes to detect it.

Population-specific causal disease effect sizes in functionally important regions impacted by selection

Many diseases exhibit population-specific causal effect sizes with trans-ethnic genetic correlations significantly less than 1, limiting trans-ethnic polygenic risk prediction. We develop a new method, S-LDXR, for stratifying squared trans-ethnic genetic correlation across genomic annotations, and apply S-LDXR to genome-wide summary statistics for 31 diseases and complex traits in East Asians (average N = 90K) and Europeans (average N = 267K) with an average trans-ethnic genetic correlation of 0.85. We determine that squared trans-ethnic genetic correlation is 0.82× (s.e. 0.01) depleted in the top quintile of background selection statistic, implying more population-specific causal effect sizes. Accordingly, causal effect sizes are more population-specific in functionally important regions, including conserved and regulatory regions. In regions surrounding specifically expressed genes, causal effect sizes are most population-specific for skin and immune genes, and least population-specific for brain genes. Our results could potentially be explained by stronger gene-environment interaction at loci impacted by selection, particularly positive selection.

A quantitative trait rare variant nonparametric linkage method with application to age-at-onset of Alzheimer's disease

To analyze pedigrees with quantitative trait (QT) and sequence data, we developed a rare variant (RV) quantitative nonparametric linkage (QNPL) method, which evaluates sharing of minor alleles. RV-QNPL has greater power than the traditional QNPL that tests for excess sharing of minor and major alleles. RV-QNPL is robust to population substructure and admixture, locus heterogeneity, and inclusion of nonpathogenic variants and can be readily applied outside of coding regions. When QNPL was used to analyze common variants, it often led to loci mapping to large intervals, e.g., >40 Mb. In contrast, when RVs are analyzed, regions are well defined, e.g., a gene. Using simulation studies, we demonstrate that RV-QNPL is substantially more powerful than applying traditional QNPL methods to analyze RVs. RV-QNPL was also applied to analyze age-at-onset (AAO) data for 107 late-onset Alzheimer's disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with whole-genome sequence data. When AAO of AD was analyzed regardless of APOE ε4 status, suggestive linkage (LOD = 2.4) was observed with RVs in KNDC1 and nominally significant linkage (p < 0.05) was observed with RVs in LOAD genes ABCA7 and IQCK. When AAO of AD was analyzed for APOE ε4 positive family members, nominally significant linkage was observed with RVs in APOE, while when AAO of AD was analyzed for APOE ε4 negative family members, nominal significance was observed for IQCK and ADAMTS1. RV-QNPL provides a powerful resource to analyze QTs in families to elucidate their genetic etiology.

Identity-by-descent detection across 487,409 British samples reveals fine scale population structure and ultra-rare variant associations

Detection of Identical-By-Descent (IBD) segments provides a fundamental measure of genetic relatedness and plays a key role in a wide range of analyses. We develop FastSMC, an IBD detection algorithm that combines a fast heuristic search with accurate coalescent-based likelihood calculations. FastSMC enables biobank-scale detection and dating of IBD segments within several thousands of years in the past. We apply FastSMC to 487,409 UK Biobank samples and detect ~214 billion IBD segments transmitted by shared ancestors within the past 1500 years, obtaining a fine-grained picture of genetic relatedness in the UK. Sharing of common ancestors strongly correlates with geographic distance, enabling the use of genomic data to localize a sample's birth coordinates with a median error of 45 km. We seek evidence of recent positive selection by identifying loci with unusually strong shared ancestry and detect 12 genome-wide significant signals. We devise an IBD-based test for association between phenotype and ultra-rare loss-of-function variation, identifying 29 association signals in 7 blood-related traits.

Genome-wide association studies and heritability analysis reveal the involvement of host genetics in the Japanese gut microbiota

Numerous host extrinsic and intrinsic factors affect the gut microbiota composition, but their cumulative effects do not sufficiently explain the variation in the microbiota, suggesting contributions of missing factors. The Japanese population possesses homogeneous genetic features suitable for genome-wide association study (GWAS). Here, we performed GWASs for human gut microbiota using 1068 healthy Japanese adults. To precisely evaluate genetic effects, we corrected for the impacts of numerous host extrinsic and demographic factors by introducing them as covariates, enabling us to discover five loci significantly associated with microbiome diversity measures: HS3ST4, C2CD2, 2p16.1, 10p15.1, and 18q12.2. Nevertheless, these five variants explain only a small fraction of the variation in the gut microbiota. We subsequently investigated the heritability of each of the 21 core genera and found that the abundances of six genera are heritable. We propose that the gut microbiota composition is affected by a highly polygenic architecture rather than several strongly associated variants in the Japanese population.

Genomic selection for heterobothriosis resistance concurrent with body size in the tiger pufferfish, Takifugu rubripes

Parasite resistance traits in aquaculture species often have moderate heritability, indicating the potential for genetic improvements by selective breeding. However, parasite resistance is often synonymous with an undesirable negative correlation with body size. In this study, we first tested the feasibility of genomic selection (GS) on resistance to heterobothriosis, caused by the monogenean parasite Heterobothrium okamotoi, which leads to huge economic losses in aquaculture of the tiger pufferfish Takifugu rubripes. Then, using a simulation study, we tested the possibility of simultaneous improvement of parasite resistance, assessed by parasite counts on host fish (HC), and standard length (SL). Each trait showed moderate heritability (square-root transformed HC: h2 = 0.308 ± 0.123, S.E.; SL: h2 = 0.405 ± 0.131). The predictive abilities of genomic prediction among 12 models, including genomic Best Linear Unbiased Predictor (GBLUP), Bayesian regressions, and machine learning procedures, were also moderate for both transformed HC (0.248‒0.344) and SL (0.340‒0.481). These results confirmed the feasibility of GS for this trait. Although an undesirable genetic correlation was suggested between transformed HC and SL (rg = 0.228), the simulation study suggested the desired gains index can help achieve simultaneous genetic improvements in both traits.

The reliability and heritability of cortical folds and their genetic correlations across hemispheres

Cortical folds help drive the parcellation of the human cortex into functionally specific regions. Variations in the length, depth, width, and surface area of these sulcal landmarks have been associated with disease, and may be genetically mediated. Before estimating the heritability of sulcal variation, the extent to which these metrics can be reliably extracted from in-vivo MRI must be established. Using four independent test-retest datasets, we found high reliability across the brain (intraclass correlation interquartile range: 0.65-0.85). Heritability estimates were derived for three family-based cohorts using variance components analysis and pooled (total N > 3000); the overall sulcal heritability pattern was correlated to that derived for a large population cohort (N > 9000) calculated using genomic complex trait analysis. Overall, sulcal width was the most heritable metric, and earlier forming sulci showed higher heritability. The inter-hemispheric genetic correlations were high, yet select sulci showed incomplete pleiotropy, suggesting hemisphere-specific genetic influences.

Heritability and genome-wide association of swine gut microbiome features with growth and fatness parameters

Despite recent efforts to characterize longitudinal variation in the swine gut microbiome, the extent to which a host's genome impacts the composition of its gut microbiome is not yet well understood in pigs. The objectives of this study were: i) to identify pig gut microbiome features associated with growth and fatness, ii) to estimate the heritability of those features, and, iii) to conduct a genome-wide association study exploring the relationship between those features and single nucleotide polymorphisms (SNP) in the pig genome. A total of 1,028 pigs were characterized. Animals were genotyped with the Illumina PorcineSNP60 Beadchip. Microbiome samples from fecal swabs were obtained at weaning (Wean), at mid-test during the growth trial (MidTest), and at the end of the growth trial (OffTest). Average daily gain was calculated from birth to week 14 of the growth trial, from weaning to week 14, from week 14 to week 22, and from week 14 to harvest. Backfat and loin depth were also measured at weeks 14 and 22. Heritability estimates (±SE) of Operational Taxonomic Units ranged from 0.025 (±0.0002) to 0.139 (±0.003), from 0.029 (±0.003) to 0.289 (±0.004), and from 0.025 (±0.003) to 0.545 (±0.034) at Wean, MidTest, and OffTest, respectively. Several SNP were significantly associated with taxa at the three time points. These SNP were located in genomic regions containing a total of 68 genes. This study provides new evidence linking gut microbiome composition with growth and carcass traits in swine, while also identifying putative host genetic markers associated with significant differences in the abundance of several prevalent microbiome features.

Genetics of the human circadian clock and sleep homeostat

Timing and duration of sleep are controlled by the circadian system, which keeps an ~24-h internal rhythm that entrains to environmental stimuli, and the sleep homeostat, which rises as a function of time awake. There is a normal distribution across the population in how the circadian system aligns with typical day and night resulting in varying circadian preferences called chronotypes. A portion of the variation in the population is controlled by genetics as shown by the single-gene mutations that confer extreme early or late chronotypes. Similarly, there is a normal distribution across the population in sleep duration. Genetic variations have been identified that lead to a short sleep phenotype in which individuals sleep only 4-6.5 h nightly. Negative health consequences have been identified when individuals do not sleep at their ideal circadian timing or are sleep deprived relative to intrinsic sleep need. Whether familial natural short sleepers are at risk of the health consequences associated with a short sleep duration based on population data is not known. More work needs to be done to better assess for an individual's chronotype and degree of sleep deprivation to answer these questions.

The transferability of lipid loci across African, Asian and European cohorts

Most genome-wide association studies are based on samples of European descent. We assess whether the genetic determinants of blood lipids, a major cardiovascular risk factor, are shared across populations. Genetic correlations for lipids between European-ancestry and Asian cohorts are not significantly different from 1. A genetic risk score based on LDL-cholesterol-associated loci has consistent effects on serum levels in samples from the UK, Uganda and Greece (r = 0.23-0.28, p < 1.9 × 10^-14). Overall, there is evidence of reproducibility for ~75% of the major lipid loci from European discovery studies, except triglyceride loci in the Ugandan samples (10% of loci). Individual transferable loci are identified using trans-ethnic colocalization. Ten of fourteen loci not transferable to the Ugandan population have pleiotropic associations with BMI in Europeans; none of the transferable loci do. The non-transferable loci might affect lipids by modifying food intake in environments rich in certain nutrients, which suggests a potential role for gene-environment interactions.

White Matter Hyperintensities Are Under Strong Genetic Influence

BACKGROUND AND PURPOSE

The genetic basis of white matter hyperintensities (WMH) is still unknown. This study examines the heritability of WMH in both sexes and in different brain regions, and the influence of age.

METHODS

Participants from the Older Australian Twins Study were recruited (n=320; 92 monozygotic and 68 dizygotic pairs) who volunteered for magnetic resonance imaging scans and medical assessments. Heritability, that is, the ratio of the additive genetic variance to the total phenotypic variance, was estimated using the twin design.

RESULTS

Heritability was high for total WMH volume (0.76), and for periventricular WMH (0.64) and deep WMH (0.77), and varied from 0.18 for the cerebellum to 0.76 for the occipital lobe. The genetic correlation between deep and periventricular WMH regions was 0.85, with one additive genetics factor accounting for most of the shared variance. Heritability was consistently higher in women in the cerebral regions. Heritability in deep but not periventricular WMH declined with age, in particular after the age of 75.

CONCLUSIONS

WMH have a strong genetic influence but this is not uniform through the brain, being higher for deep than periventricular WMH and in the cerebral regions. The genetic influence is higher in women, and there is an age-related decline, most markedly for deep WMH. The data suggest some heterogeneity in the pathogenesis of WMH for different brain regions and for men and women.