Cultural, demographic, and other non-demographic factors associated with cancer genetic counseling patients' appointment accompaniment preferences in the United States.

Although the presence of companion(s) in a genetic counseling session can positively influence session dynamics, research has found that some patients prefer to attend their appointments alone. To date, no studies have examined patient accompaniment preferences across different cultural groups in the context of genetic counseling. This quantitative study aimed to identify factors associated with individual preferences in accompaniment at cancer genetic counseling appointments in a sample (N = 130) of Hispanic/Latine (n = 29) and non-Hispanic/Latine White (n = 101) participants at a large academic medical institution. Variables examined included demographics, horizontal and vertical collectivism, and Hispanic and American acculturation. A link to an online questionnaire was emailed to patients who met four criteria: (1) identified as either Hispanic/Latine or non-Hispanic/Latine White; (2) had attended a cancer genetic counseling appointment at UCLA Health to discuss genetic testing options between October 2020 and December 2022; (3) were at least 18 years of age at the time of their appointment; and (4) indicated they were comfortable reading in Spanish or English; responses were anonymous. Logistic regression analyses identified four significant variables in the model associated with accompaniment preferences: individuals with at least one parent born outside of the US, those who attended their appointment in-person, and those with a higher horizontal collectivism score were less likely to want to attend their cancer genetic counseling appointment alone, while the converse was true among those with a higher American acculturation score. These findings highlight cultural and demographic factors that are associated with patient accompaniment preferences unrelated to ethnicity, indicating genetic counselors should not make assumptions regarding accompaniment preferences based solely on cultural or racial/ethnic background. Genetic counselors should incorporate this understanding when assessing patients' accompaniment preferences.

Genome-wide mapping of brain phenotypes in extended pedigrees with strong genetic loading for bipolar disorder.

Bipolar disorder is a highly heritable illness, associated with alterations of brain structure. As such, identification of genes influencing inter-individual differences in brain morphology may help elucidate the underlying pathophysiology of bipolar disorder (BP). To identify quantitative trait loci (QTL) that contribute to phenotypic variance of brain structure, structural neuroimages were acquired from family members (n = 527) of extended pedigrees heavily loaded for bipolar disorder ascertained from genetically isolated populations in Latin America. Genome-wide linkage and association analysis were conducted on the subset of heritable brain traits that showed significant evidence of association with bipolar disorder (n = 24) to map QTL influencing regional measures of brain volume and cortical thickness. Two chromosomal regions showed significant evidence of linkage; a QTL on chromosome 1p influencing corpus callosum volume and a region on chromosome 7p linked to cortical volume. Association analysis within the two QTLs identified three SNPs correlated with the brain measures.

Genetic analysis of activity, brain and behavioral associations in extended families with heavy genetic loading for bipolar disorder.

BACKGROUND: Disturbed sleep and activity are prominent features of bipolar disorder type I (BP-I). However, the relationship of sleep and activity characteristics to brain structure and behavior in euthymic BP-I patients and their non-BP-I relatives is unknown. Additionally, underlying genetic relationships between these traits have not been investigated. METHODS: Relationships between sleep and activity phenotypes, assessed using actigraphy, with structural neuroimaging (brain) and cognitive and temperament (behavior) phenotypes were investigated in 558 euthymic individuals from multi-generational pedigrees including at least one member with BP-I. Genetic correlations between actigraphy-brain and actigraphy-behavior associations were assessed, and bivariate linkage analysis was conducted for trait pairs with evidence of shared genetic influences. RESULTS: More physical activity and longer awake time were significantly associated with increased brain volumes and cortical thickness, better performance on neurocognitive measures of long-term memory and executive function, and less extreme scores on measures of temperament (impulsivity, cyclothymia). These associations did not differ between BP-I patients and their non-BP-I relatives. For nine activity-brain or activity-behavior pairs there was evidence for shared genetic influence (genetic correlations); of these pairs, a suggestive bivariate quantitative trait locus on chromosome 7 for wake duration and verbal working memory was identified. CONCLUSIONS: Our findings indicate that increased physical activity and more adequate sleep are associated with increased brain size, better cognitive function and more stable temperament in BP-I patients and their non-BP-I relatives. Additionally, we found evidence for pleiotropy of several actigraphy-behavior and actigraphy-brain phenotypes, suggesting a shared genetic basis for these traits.

Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect.

BACKGROUND: Bipolar disorder (BD) is a highly heritable mood disorder with complex genetic architecture and poorly understood etiology. Previous transcriptomic BD studies have had inconsistent findings due to issues such as small sample sizes and difficulty in adequately accounting for confounders like medication use. METHODS: We performed a differential expression analysis in a well-characterized BD case-control sample (Nsubjects = 480) by RNA sequencing of whole blood. We further performed co-expression network analysis, functional enrichment, and cell type decomposition, and integrated differentially expressed genes with genetic risk. RESULTS: While we observed widespread differential gene expression patterns between affected and unaffected individuals, these effects were largely linked to lithium treatment at the time of blood draw (FDR < 0.05, Ngenes = 976) rather than BD diagnosis itself (FDR < 0.05, Ngenes = 6). These lithium-associated genes were enriched for cell signaling and immune response functional annotations, among others, and were associated with neutrophil cell-type proportions, which were elevated in lithium users. Neither genes with altered expression in cases nor in lithium users were enriched for BD, schizophrenia, and depression genetic risk based on information from genome-wide association studies, nor was gene expression associated with polygenic risk scores for BD. CONCLUSIONS: These findings suggest that BD is associated with minimal changes in whole blood gene expression independent of medication use but emphasize the importance of accounting for medication use and cell type heterogeneity in psychiatric transcriptomic studies. The results of this study add to mounting evidence of lithium's cell signaling and immune-related mechanisms.

Genetic Deficiency of Flavin-Containing Monooxygenase 3 ( Fmo3) Protects Against Thrombosis but Has Only a Minor Effect on Plasma Lipid Levels-Brief Report.

Objective- FMO (flavin-containing monooxygenase) 3 converts bacterial-derived trimethylamine to trimethylamine N-oxide (TMAO), an independent risk factor for cardiovascular disease. We generated FMO3 knockout (FMO3KO) mouse to study its effects on plasma TMAO, lipids, glucose/insulin metabolism, thrombosis, and atherosclerosis. Approach and Results- Previous studies with an antisense oligonucleotide (ASO) knockdown strategy targeting FMO3 in LDLRKO (low-density lipoprotein receptor knockout) mice resulted in major reductions in TMAO levels and atherosclerosis, but also showed effects on plasma lipids, insulin, and glucose. Although FMO3KO mice generated via CRISPR/Cas9 technology bred onto the LDLRKO background did exhibit similar effects on TMAO levels, the effects on lipid metabolism were not as pronounced as with the ASO knockdown model. These differences could result from either off-target effects of the ASO or from a developmental adaptation to the FMO3 deficiency. To distinguish these possibilities, we treated wild-type and FMO3KO mice with control or FMO3 ASOs. FMO3-ASO treatment led to the same extent of lipid-lowering effects in the FMO3KO mice as the wild-type mice, indicating off-target effects. The levels of TMAO in LDLRKO mice fed an atherogenic diet are very low in both wild-type and FMO3KO mice, and no significant effect was observed on atherosclerosis. When FMO3KO and wild-type mice were maintained on a 0.5% choline diet, FMO3KO showed a marked reduction in both TMAO and in vivo thrombosis potential. Conclusions- FMO3KO markedly reduces systemic TMAO levels and thrombosis potential. However, the previously observed large effects of an FMO3 ASO on plasma lipid levels appear to be due partly to off-target effects.

Rare Inherited and De Novo CNVs Reveal Complex Contributions to ASD Risk in Multiplex Families.

Rare mutations, including copy-number variants (CNVs), contribute significantly to autism spectrum disorder (ASD) risk. Although their importance has been established in families with only one affected child (simplex families), the contribution of both de novo and inherited CNVs to ASD in families with multiple affected individuals (multiplex families) is less well understood. We analyzed 1,532 families from the Autism Genetic Resource Exchange (AGRE) to assess the impact of de novo and rare CNVs on ASD risk in multiplex families. We observed a higher burden of large, rare CNVs, including inherited events, in individuals with ASD than in their unaffected siblings (odds ratio [OR] = 1.7), but the rate of de novo events was significantly lower than in simplex families. In previously characterized ASD risk loci, we identified 49 CNVs, comprising 24 inherited events, 19 de novo events, and 6 events of unknown inheritance, a significant enrichment in affected versus control individuals (OR = 3.3). In 21 of the 30 families (71%) in whom at least one affected sibling harbored an established ASD major risk CNV, including five families harboring inherited CNVs, the CNV was not shared by all affected siblings, indicating that other risk factors are contributing. We also identified a rare risk locus for ASD and language delay at chromosomal region 2q24 (implicating NR4A2) and another lower-penetrance locus involving inherited deletions and duplications of WWOX. The genetic architecture in multiplex families differs from that in simplex families and is complex, warranting more complete genetic characterization of larger multiplex ASD cohorts.

Characterization of Expression Quantitative Trait Loci in Pedigrees from Colombia and Costa Rica Ascertained for Bipolar Disorder.

The observation that variants regulating gene expression (expression quantitative trait loci, eQTL) are at a high frequency among SNPs associated with complex traits has made the genome-wide characterization of gene expression an important tool in genetic mapping studies of such traits. As part of a study to identify genetic loci contributing to bipolar disorder and other quantitative traits in members of 26 pedigrees from Costa Rica and Colombia, we measured gene expression in lymphoblastoid cell lines derived from 786 pedigree members. The study design enabled us to comprehensively reconstruct the genetic regulatory network in these families, provide estimates of heritability, identify eQTL, evaluate missing heritability for the eQTL, and quantify the number of different alleles contributing to any given locus. In the eQTL analysis, we utilize a recently proposed hierarchical multiple testing strategy which controls error rates regarding the discovery of functional variants. Our results elucidate the heritability and regulation of gene expression in this unique Latin American study population and identify a set of regulatory SNPs which may be relevant in future investigations of complex disease in this population. Since our subjects belong to extended families, we are able to compare traditional kinship-based estimates with those from more recent methods that depend only on genotype information.

Genetic contributions to circadian activity rhythm and sleep pattern phenotypes in pedigrees segregating for severe bipolar disorder.

Abnormalities in sleep and circadian rhythms are central features of bipolar disorder (BP), often persisting between episodes. We report here, to our knowledge, the first systematic analysis of circadian rhythm activity in pedigrees segregating severe BP (BP-I). By analyzing actigraphy data obtained from members of 26 Costa Rican and Colombian pedigrees [136 euthymic (i.e., interepisode) BP-I individuals and 422 non-BP-I relatives], we delineated 73 phenotypes, of which 49 demonstrated significant heritability and 13 showed significant trait-like association with BP-I. All BP-I-associated traits related to activity level, with BP-I individuals consistently demonstrating lower activity levels than their non-BP-I relatives. We analyzed all 49 heritable phenotypes using genetic linkage analysis, with special emphasis on phenotypes judged to have the strongest impact on the biology underlying BP. We identified a locus for interdaily stability of activity, at a threshold exceeding genome-wide significance, on chromosome 12pter, a region that also showed pleiotropic linkage to two additional activity phenotypes.

DMD genotype correlations from the Duchenne Registry: Endogenous exon skipping is a factor in prolonged ambulation for individuals with a defined mutation subtype.

Antisense oligonucleotide (AON)-mediated exon skipping is an emerging therapeutic for individuals with Duchenne muscular dystrophy (DMD). Skipping of exons adjacent to common exon deletions in DMD using AONs can produce in-frame transcripts and functional protein. Targeted skipping of DMD exons 8, 44, 45, 50, 51, 52, 53, and 55 is predicted to benefit 47% of affected individuals. We observed a correlation between mutation subgroups and age at loss of ambulation in the Duchenne Registry, a large database of phenotypic and genetic data for DMD (N = 765). Males amenable to exon 44 (N = 74) and exon 8 skipping (N = 18) showed prolonged ambulation compared to other exon skip groups and nonsense mutations (P = 0.035 and P < 0.01, respectively). In particular, exon 45 deletions were associated with prolonged age at loss of ambulation relative to the rest of the exon 44 skip amenable cohort and other DMD mutations. Exon 3-7 deletions also showed prolonged ambulation relative to all other exon 8 skippable mutations. Cultured myotubes from DMD patients with deletions of exons 3-7 or exon 45 showed higher endogenous skipping than other mutations, providing a potential biological rationale for our observations. These results highlight the utility of aggregating phenotypic and genotypic data for rare pediatric diseases to reveal progression differences, identify potentially confounding factors, and probe molecular mechanisms that may affect disease severity.

Diet1, bile acid diarrhea, and FGF15/19: mouse model and human genetic variants.

Diet1 modulates intestinal production of the hormone, fibroblast growth factor (FGF)15, which signals in liver to regulate bile acid synthesis. C57BL/6ByJ mice with a spontaneous Diet1-null mutation are resistant to hypercholesterolemia compared with wild-type C57BL/6J mice through enhanced cholesterol conversion to bile acids. To further characterize the role of Diet1 in metabolism, we generated Diet1-/- mice on the C57BL/6J genetic background. C57BL/6J Diet1-/- mice had elevated bile acid levels, reduced Fgf15 expression, and increased gastrointestinal motility and intestinal luminal water content, which are symptoms of bile acid diarrhea (BAD) in humans. Natural genetic variation in Diet1 mRNA expression levels across 76 inbred mouse strains correlated positively with Ffg15 mRNA and negatively with serum bile acid levels. This led us to investigate the role of DIET1 genetic variation in primary BAD patients. We identified a DIET1 coding variant (rs12256835) that had skewed prevalence between BAD cases and controls. This variant causes an H1721Q amino acid substitution that increases the levels of FGF19 protein secreted from cultured cells. We propose that genetic variation in DIET1 may be a determinant of FGF19 secretion levels, and may affect bile acid metabolism in both physiological and pathological conditions.

Detecting responses to treatment with fenofibrate in pedigrees.

BACKGROUND: Fenofibrate (Fb) is a known treatment for elevated triglyceride (TG) levels. The Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study was designed to investigate potential contributors to the effects of Fb on TG levels. Here, we summarize the analyses of 8 papers whose authors had access to the GOLDN data and were grouped together because they pursued investigations into Fb treatment responses as part of GAW20. These papers report explorations of a variety of genetics, epigenetics, and study design questions. Data regarding treatment with 160 mg of micronized Fb per day for 3 weeks included pretreatment and posttreatment TG and methylation levels (ML) at approximately 450,000 epigenetic markers (cytosine-phosphate-guanine [CpG] sites). In addition, approximately 1 million single-nucleotide polymorphisms (SNPs) were genotyped or imputed in each of the study participants, drawn from 188 pedigrees. RESULTS: The analyses of a variety of subsets of the GOLDN data used a number of analytic approaches such as linear mixed models, a kernel score test, penalized regression, and artificial neural networks. CONCLUSIONS: Results indicate that (a) CpG ML are responsive to Fb; (b) CpG ML should be included in models predicting the TG level responses to Fb;

Genetic and hormonal control of hepatic steatosis in female and male mice.

The etiology of nonalcoholic fatty liver disease is complex and influenced by factors such as obesity, insulin resistance, hyperlipidemia, and sex. We now report a study on sex difference in hepatic steatosis in the context of genetic variation using a population of inbred strains of mice. While male mice generally exhibited higher concentration of hepatic TG levels on a high-fat high-sucrose diet, sex differences showed extensive interaction with genetic variation. Differences in percentage body fat were the best predictor of hepatic steatosis among the strains and explained about 30% of the variation in both sexes. The difference in percent gonadal fat and HDL explained 9.6% and 6.7% of the difference in hepatic TGs between the sexes, respectively. Genome-wide association mapping of hepatic TG revealed some striking differences in genetic control of hepatic steatosis between females and males. Gonadectomy increased the hepatic TG to body fat percentage ratio among male, but not female, mice. Our data suggest that the difference between the sexes in hepatic TG can be partly explained by differences in body fat distribution, plasma HDL, and genetic regulation. Future studies are required to understand the molecular interactions between sex, genetics, and the environment.

Diet, gonadal sex, and sex chromosome complement influence white adipose tissue miRNA expression.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNA molecules that regulate gene expression by targeting specific mRNA species for degradation or interfering with translation. Specific miRNAs are key regulators of adipogenesis, and are expressed at different levels in adipose tissue from lean and obese mice. The degree of lipid accumulation and distribution of white adipose tissue differs between males and females, and it is unknown whether sex differences in adipose tissue-specific miRNA expression may contribute to this dimorphism. Typically, sex differences are attributed to hormones secreted from ovaries or testes. However, the sex chromosome complement (XX versus XY) is also a determinant of sex differences and may regulate miRNA expression in adipocytes. RESULTS: To identify sex differences in adipose tissue miRNA expression and to understand the underlying mechanisms, we performed high-throughput miRNA sequencing in gonadal fat depots of the Four Core Genotypes mouse model. This model, which consists of XX female, XX male, XY female, and XY male mice, allowed us to assess independent effects of gonadal type (male vs. female) and sex chromosome complement (XX vs. XY) on miRNA expression profiles. We have also assessed the effects of a high fat diet on sex differences in adipose tissue miRNA profiles. We identified a male-female effect on the overall miRNA expression profile in mice fed a chow diet, with a bias toward higher expression in male compared to female gonadal adipose tissue. This sex bias disappeared after gonadectomy, suggesting that circulating levels of gonadal secretions modulate the miRNA expression profile. After 16 weeks of high fat diet, the miRNA expression distribution was shifted toward higher expression in XY vs. XX adipose tissue. Principal component analysis revealed that high fat diet has a substantial effect on miRNA profile variance, while gonadal secretions and sex chromosome complement each have milder effects. CONCLUSIONS: Our results demonstrate that the overall miRNA expression profile in adipose tissue is influenced by gonadal hormones and the sex chromosome complement, and that expression profiles change in response to gonadectomy and high fat diet. Differential miRNA expression profiles may contribute to sex differences in adipose tissue gene expression, adipose tissue development, and diet-induced obesity.

Transcriptomic analysis of autistic brain reveals convergent molecular pathology.

Autism spectrum disorder (ASD) is a common, highly heritable neurodevelopmental condition characterized by marked genetic heterogeneity. Thus, a fundamental question is whether autism represents an aetiologically heterogeneous disorder in which the myriad genetic or environmental risk factors perturb common underlying molecular pathways in the brain. Here, we demonstrate consistent differences in transcriptome organization between autistic and normal brain by gene co-expression network analysis. Remarkably, regional patterns of gene expression that typically distinguish frontal and temporal cortex are significantly attenuated in the ASD brain, suggesting abnormalities in cortical patterning. We further identify discrete modules of co-expressed genes associated with autism: a neuronal module enriched for known autism susceptibility genes, including the neuronal specific splicing factor A2BP1 (also known as FOX1), and a module enriched for immune genes and glial markers. Using high-throughput RNA sequencing we demonstrate dysregulated splicing of A2BP1-dependent alternative exons in the ASD brain. Moreover, using a published autism genome-wide association study (GWAS) data set, we show that the neuronal module is enriched for genetically associated variants, providing independent support for the causal involvement of these genes in autism. In contrast, the immune-glial module showed no enrichment for autism GWAS signals, indicating a non-genetic aetiology for this process. Collectively, our results provide strong evidence for convergent molecular abnormalities in ASD, and implicate transcriptional and splicing dysregulation as underlying mechanisms of neuronal dysfunction in this disorder.

Linkage analysis in a Dutch population isolate shows no major gene for left-handedness or atypical language lateralization.

Cerebral dominance of language function and hand preference are suggested to be heritable traits with possible shared genetic background. However, joined genetic studies of these traits have never been conducted. We performed a genetic linkage study in 37 multigenerational human pedigrees of both sexes (consisting of 355 subjects) enriched with left-handedness in which we also measured language lateralization. Hand preference was measured with the Edinburgh Handedness Inventory, and language lateralization was measured with functional transcranial Doppler during language production. The estimated heritability of left-handedness and language lateralization in these pedigrees is 0.24 and 0.31, respectively. A parametric major gene model was tested for left-handedness. Nonparametric analyses were performed for left-handedness, atypical lateralization, and degree of language lateralization. We did not observe genome-wide evidence for linkage in the parametric or nonparametric analyses for any of the phenotypes tested. However, multiple regions showed suggestive evidence of linkage. The parametric model showed suggestive linkage for left-handedness in the 22q13 region [heterogeneity logarithm of odds (HLOD) = 2.18]. Nonparametric multipoint analysis of left-handedness showed suggestive linkage in the same region [logarithm of odds (LOD) = 2.80]. Atypical language lateralization showed suggestive linkage in the 7q34 region (LODMax = 2.35). For strength of language lateralization, we observed suggestive linkage in the 6p22 (LODMax = 2.54), 7q32 (LODMax = 1.93), and 9q33 (LODMax = 2.10) regions. We did not observe any overlap of suggestive genetic signal between handedness and the extent of language lateralization. The absence of significant linkage argues against the presence of a major gene coding for both traits; rather, our results are suggestive of these traits being two independent polygenic complex traits.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production.

OBJECTIVES: Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants. METHODS: We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA. RESULTS: We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10-3 and 6.8×10-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10-5 and 2.0×10-4, respectively). CONCLUSION: We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Gene expression in large pedigrees: analytic approaches.

BACKGROUND: We currently have the ability to quantify transcript abundance of messenger RNA (mRNA), genome-wide, using microarray technologies. Analyzing genotype, phenotype and expression data from 20 pedigrees, the members of our Genetic Analysis Workshop (GAW) 19 gene expression group published 9 papers, tackling some timely and important problems and questions. To study the complexity and interrelationships of genetics and gene expression, we used established statistical tools, developed newer statistical tools, and developed and applied extensions to these tools. METHODS: To study gene expression correlations in the pedigree members (without incorporating genotype or trait data into the analysis), 2 papers used principal components analysis, weighted gene coexpression network analysis, meta-analyses, gene enrichment analyses, and linear mixed models. To explore the relationship between genetics and gene expression, 2 papers studied expression quantitative trait locus allelic heterogeneity through conditional association analyses, and epistasis through interaction analyses. A third paper assessed the feasibility of applying allele-specific binding to filter potential regulatory single-nucleotide polymorphisms (SNPs). Analytic approaches included linear mixed models based on measured genotypes in pedigrees, permutation tests, and covariance kernels. To incorporate both genotype and phenotype data with gene expression, 4 groups employed linear mixed models, nonparametric weighted U statistics, structural equation modeling, Bayesian unified frameworks, and multiple regression. RESULTS AND DISCUSSION: Regarding the analysis of pedigree data, we found that gene expression is familial, indicating that at least 1 factor for pedigree membership or multiple factors for the degree of relationship should be included in analyses, and we developed a method to adjust for familiality prior to conducting weighted co-expression gene network analysis. For SNP association and conditional analyses, we found FaST-LMM (Factored Spectrally Transformed Linear Mixed Model) and SOLAR-MGA (Sequential Oligogenic Linkage Analysis Routines -Major Gene Analysis) have similar type 1 and type 2 errors and can be used almost interchangeably. To improve the power and precision of association tests, prior knowledge of DNase-I hypersensitivity sites or other relevant biological annotations can be incorporated into the analyses. On a biological level, eQTL (expression quantitative trait loci) are genetically complex, exhibiting both allelic heterogeneity and epistasis. Including both genotype and phenotype data together with measurements of gene expression was found to be generally advantageous in terms of generating improved levels of significance and in providing more interpretable biological models. CONCLUSIONS: Pedigrees can be used to conduct analyses of and enhance gene expression studies.

Brain structure-function associations in multi-generational families genetically enriched for bipolar disorder.

Recent theories regarding the pathophysiology of bipolar disorder suggest contributions of both neurodevelopmental and neurodegenerative processes. While structural neuroimaging studies indicate disease-associated neuroanatomical alterations, the behavioural correlates of these alterations have not been well characterized. Here, we investigated multi-generational families genetically enriched for bipolar disorder to: (i) characterize neurobehavioural correlates of neuroanatomical measures implicated in the pathophysiology of bipolar disorder; (ii) identify brain-behaviour associations that differ between diagnostic groups; (iii) identify neurocognitive traits that show evidence of accelerated ageing specifically in subjects with bipolar disorder; and (iv) identify brain-behaviour correlations that differ across the age span. Structural neuroimages and multi-dimensional assessments of temperament and neurocognition were acquired from 527 (153 bipolar disorder and 374 non-bipolar disorder) adults aged 18-87 years in 26 families with heavy genetic loading for bipolar disorder. We used linear regression models to identify significant brain-behaviour associations and test whether brain-behaviour relationships differed: (i) between diagnostic groups; and (ii) as a function of age. We found that total cortical and ventricular volume had the greatest number of significant behavioural associations, and included correlations with measures from multiple cognitive domains, particularly declarative and working memory and executive function. Cortical thickness measures, in contrast, showed more specific associations with declarative memory, letter fluency and processing speed tasks. While the majority of brain-behaviour relationships were similar across diagnostic groups, increased cortical thickness in ventrolateral prefrontal and parietal cortical regions was associated with better declarative memory only in bipolar disorder subjects, and not in non-bipolar disorder family members. Additionally, while age had a relatively strong impact on all neurocognitive traits, the effects of age on cognition did not differ between diagnostic groups. Most brain-behaviour associations were also similar across the age range, with the exception of cortical and ventricular volume and lingual gyrus thickness, which showed weak correlations with verbal fluency and inhibitory control at younger ages that increased in magnitude in older subjects, regardless of diagnosis. Findings indicate that neuroanatomical traits potentially impacted by bipolar disorder are significantly associated with multiple neurobehavioural domains. Structure-function relationships are generally preserved across diagnostic groups, with the notable exception of ventrolateral prefrontal and parietal association cortex, volumetric increases in which may be associated with cognitive resilience specifically in individuals with bipolar disorder. Although age impacted all neurobehavioural traits, we did not find any evidence of accelerated cognitive decline specific to bipolar disorder subjects. Regardless of diagnosis, greater global brain volume may represent a protective factor for the effects of ageing on executive functioning.

Preferential binding to Elk-1 by SLE-associated IL10 risk allele upregulates IL10 expression.

Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10⁻8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.

MicroRNA-3148 modulates allelic expression of toll-like receptor 7 variant associated with systemic lupus erythematosus.

We previously reported that the G allele of rs3853839 at 3'untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10(-10), odds ratio (OR) (95%CI) = 1.27 (1.17-1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10(-11), OR = 1.24 [1.18-1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3'UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R(2) = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3'UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta  = 2.0×10(-19), OR = 1.25 [1.20-1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.