Pooled Sequencing of Candidate Genes Implicates Rare Variants in the Development of Asthma Following Severe RSV Bronchiolitis in Infancy

Severe infection with respiratory syncytial virus (RSV) during infancy is strongly associated with the development of asthma. To identify genetic variation that contributes to asthma following severe RSV bronchiolitis during infancy, we sequenced the coding exons of 131 asthma candidate genes in 182 European and African American children with severe RSV bronchiolitis in infancy using anonymous pools for variant discovery, and then directly genotyped a set of 190 nonsynonymous variants. Association testing was performed for physician-diagnosed asthma before the 7^th birthday (asthma) using genotypes from 6,500 individuals from the Exome Sequencing Project (ESP) as controls to gain statistical power. In addition, among patients with severe RSV bronchiolitis during infancy, we examined genetic associations with asthma, active asthma, persistent wheeze, and bronchial hyperreactivity (methacholine PC₂₀) at age 6 years. We identified four rare nonsynonymous variants that were significantly associated with asthma following severe RSV bronchiolitis, including single variants in ADRB2, FLG and NCAM1 in European Americans (p = 4.6x10^-4, 1.9x10^-13 and 5.0x10^-5, respectively), and NOS1 in African Americans (p = 2.3x10^-11). One of the variants was a highly functional nonsynonymous variant in ADRB2 (rs1800888), which was also nominally associated with asthma (p = 0.027) and active asthma (p = 0.013) among European Americans with severe RSV bronchiolitis without including the ESP. Our results suggest that rare nonsynonymous variants contribute to the development of asthma following severe RSV bronchiolitis in infancy, notably in ADRB2. Additional studies are required to explore the role of rare variants in the etiology of asthma and asthma-related traits following severe RSV bronchiolitis.     

Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.     

Rare Risk Variants Identification by Identity-by-Descent Mapping and Whole-Exome Sequencing Implicates Neuronal Development Pathways in Schizophrenia and Bipolar Disorder

Schizophrenia (SCZ) and bipolar disorder (BPD) are highly heritable disorders with an estimated co-heritability of 68%. Hundreds of common alleles have been implicated, but recently a role for rare, high-penetrant variants has been also suggested in both disorders. This study investigated a familial cohort of SCZ and BPD patients from a closed population sample, where the high recurrence of the disorders and the homogenous genetic background indicate a possible enrichment in rare risk alleles. A total of 230 subjects (161 cases, 22 unaffected relatives, and 47 controls) were genetically investigated through an innovative strategy that integrates identity-by-descent (IBD) mapping and whole-exome sequencing (WES). IBD analysis allowed to track high-risk haplotypes (IBD_risk) shared exclusively by multiple patients from different families and possibly carrying the most penetrant alleles. A total of 444 non-synonymous sequence variants, of which 137 disruptive, were identified in IBD_risk haplotypes by WES. Interestingly, gene sets previously implicated in SCZ (i.e., post-synaptic density (PSD) proteins, voltage-gated calcium channels (VGCCs), and fragile X mental retardation protein (FMRP) targets) were found significantly enriched in genes carrying IBD_risk variants. Further, IBD_risk variants were preferentially affecting genes involved in the extracellular matrix (ECM) biology and axon guidance processes which appeared to be functionally connected in the pathway-derived meta-network analysis. Results thus confirm rare risk variants as key factors in SCZ and BPD pathogenesis and highlight a role for the development of neuronal connectivity in the etiology of both disorders.     

A Fast Method that Uses Polygenic Scores to Estimate the Variance Explained by Genome-wide Marker Panels and the Proportion of Variants Affecting a Trait

Several methods have been proposed to estimate the variance in disease liability explained by large sets of genetic markers. However, current methods do not scale up well to large sample sizes. Linear mixed models require solving high-dimensional matrix equations, and methods that use polygenic scores are very computationally intensive. Here we propose a fast analytic method that uses polygenic scores, based on the formula for the non-centrality parameter of the association test of the score. We estimate model parameters from the results of multiple polygenic score tests based on markers with p values in different intervals. We estimate parameters by maximum likelihood and use profile likelihood to compute confidence intervals. We compare various options for constructing polygenic scores, based on nested or disjoint intervals of p values, weighted or unweighted effect sizes, and different numbers of intervals, in estimating the variance explained by a set of markers, the proportion of markers with effects, and the genetic covariance between a pair of traits. Our method provides nearly unbiased estimates and confidence intervals with good coverage, although estimation of the variance is less reliable when jointly estimated with the covariance. We find that disjoint p value intervals perform better than nested intervals, but the weighting did not affect our results. A particular advantage of our method is that it can be applied to summary statistics from single markers, and so can be quickly applied to large consortium datasets. Our method, named AVENGEME (Additive Variance Explained and Number of Genetic Effects Method of Estimation), is implemented in R software.     

Analysis of Stop-Gain and Frameshift Variants in Human Innate Immunity Genes

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.     

Analysis of Genes Involved in Body Weight Regulation by Targeted Re-Sequencing

Introduction

Genes involved in body weight regulation that were previously investigated in genome-wide association studies (GWAS) and in animal models were target-enriched followed by massive parallel next generation sequencing.

Methods

We enriched and re-sequenced continuous genomic regions comprising FTO, MC4R, TMEM18, SDCCAG8, TKNS, MSRA and TBC1D1 in a screening sample of 196 extremely obese children and adolescents with age and sex specific body mass index (BMI) ≥ 99^th percentile and 176 lean adults (BMI ≤ 15^th percentile). 22 variants were confirmed by Sanger sequencing. Genotyping was performed in up to 705 independent obesity trios (extremely obese child and both parents), 243 extremely obese cases and 261 lean adults.

Results and Conclusion

We detected 20 different non-synonymous variants, one frame shift and one nonsense mutation in the 7 continuous genomic regions in study groups of different weight extremes. For SNP Arg695Cys (rs58983546) in TBC1D1 we detected nominal association with obesity (p_TDT = 0.03 in 705 trios). Eleven of the variants were rare, thus were only detected heterozygously in up to ten individual(s) of the complete screening sample of 372 individuals. Two of them (in FTO and MSRA) were found in lean individuals, nine in extremely obese. In silico analyses of the 11 variants did not reveal functional implications for the mutations. Concordant with our hypothesis we detected a rare variant that potentially leads to loss of FTO function in a lean individual. For TBC1D1, in contrary to our hypothesis, the loss of function variant (Arg443Stop) was found in an obese individual. Functional in vitro studies are warranted.     

Identification of Potential Pleiotropic Genes for Immune and Skeletal Diseases Using Multivariate MetaCCA Analysis

BackgroundImmune and skeletal systems physiologically and pathologically interact with each other. Immune and skeletal diseases may share potential pleiotropic genetics factors, but the shared specific genes are largely unknown.ObjectiveThis study aimed to investigate the overlapping genetic factors between multiple diseases (including rheumatoid arthritis (RA), psoriasis, osteoporosis, osteoarthritis, sarcopenia, and fracture).MethodsThe canonical correlation analysis (metaCCA) approach was used to identify the shared genes for six diseases by integrating genome-wide association study (GWAS)-derived summary statistics. The versatile Gene-based Association Study (VEGAS2) method was further applied to refine and validate the putative pleiotropic genes identified by metaCCA.ResultsAbout 157 (p<8.19E-6), 319 (p<3.90E-6), and 77 (p<9.72E-6) potential pleiotropic genes were identified shared by two immune diseases, four skeletal diseases, and all of the six diseases, respectively. The top three significant putative pleiotropic genes shared by both immune and skeletal diseases, including HLA-B, TSBP1, and TSBP1-AS1 (p<E-300), were located in the major histocompatibility complex (MHC) region. Nineteen of 77 putative pleiotropic genes identified by metaCCA analysis were associated with at least one disease in the VEGAS2 analysis. Specifically, the majority (18) of these 19 putative validated pleiotropic genes were associated with RA.ConclusionThe metaCCA method identified some pleiotropic genes shared by the immune and skeletal diseases. These findings help to improve our understanding of the shared genetic mechanisms and signaling pathways underlying immune and skeletal diseases.     

A Multiplex Human Syndrome Implicates a Key Role for Intestinal Cell Kinase in Development of Central Nervous, Skeletal, and Endocrine Systems

Six infants in an Old Order Amish pedigree were observed to be affected with endocrine-cerebro-osteodysplasia (ECO). ECO is a previously unidentified neonatal lethal recessive disorder with multiple anomalies involving the endocrine, cerebral, and skeletal systems. Autozygosity mapping and sequencing identified a previously unknown missense mutation, R272Q, in ICK, encoding intestinal cell kinase (ICK). Our results established that R272 is conserved across species and among ethnicities, and three-dimensional analysis of the protein structure suggests protein instability due to the R272Q mutation. We also demonstrate that the R272Q mutant fails to localize at the nucleus and has diminished kinase activity. These findings suggest that ICK plays a key role in the development of multiple organ systems.     

甘蓝型油菜ZF-HD基因家族的鉴定与系统进化分析

ZF-HD是一类植物特有的转录因子, 在植物生长发育及胁迫响应过程中发挥重要作用。利用生物信息学方法, 在甘蓝型油菜(Brassica napus)基因组中鉴定到62个ZF-HD基因, 其中83.9%的基因缺乏内含子, 93.5%的BnZF-HD等电点大于7, 预测定位于细胞核的蛋白大多由100个以上氨基酸组成。根据进化关系可将其分为6个亚群, 在每个亚群中, 甘蓝(B. oleracea)和白菜(B. rapa)的ZF-HD基因数量相等或近似相等, 而甘蓝型油菜的ZF-HD基因数量接近或等同于甘蓝和白菜的ZF-HD基因数量之和。同一亚群的motif数量和类型高度相似。共线性分析结果显示, 全基因组三倍体化使ZF-HD基因在二倍体祖先得到扩张, 而异源多倍体化又进一步使甘蓝型油菜ZF-HD基因家族扩张。K_a/K_s值说明大多数ZF-HD基因在进化过程中受到了纯化选择。所有BnZF-HD基因都具有光响应元件, 2/3的基因具有MeJA、ABA和厌氧诱导顺式作用元件, 推测这些基因可能参与相关生物学过程。研究结果为进一步挖掘该家族基因的生物学功能奠定基础, 同时为揭示多基因家族在异源多倍体中的进化式样提供借鉴。     

Identifying the Source of Unknown Microcystin Genes and Predicting Microcystin Variants by Comparing Genes within Uncultured Cyanobacterial Cells

While multiple phylogenetic markers have been used in the culture-independent study of microcystin-producing cyanobacteria, in only a few instances have multiple markers been studied within individual cells, and in all cases these studies have been conducted with cultured isolates. Here, we isolate and evaluate large DNA fragments (>6 kb) encompassing two genes involved in microcystin biosynthesis (mcyA2 and mcyB1) and use them to identify the source of gene fragments found in water samples. Further investigation of these gene loci from individual cyanobacterial cells allowed for improved analysis of the genetic diversity within microcystin producers as well as a method to predict microcystin variants for individuals. These efforts have also identified the source of the novel mcyA genotype previously termed Microcystis-like that is pervasive in the Laurentian Great Lakes and they predict the microcystin variant(s) that it produces.Microcystin-producing cyanobacteria are common nuisance organisms in harmful algal blooms in freshwaters around the world (4). This genetically diverse group (based on 16S rRNA, mcyA, mcyD, and mcyE gene sequences [6, 10, 15, 16, 22]) ranges in morphology from unicellular and colonial cocci to large filamentous strands. Many species can produce a variety of secondary metabolites that can act as hepatatoxins upon ingestion by animals (e.g., variants of microcystin) (4, 33). Microcystin production reduces the water quality in reservoirs used by human populations and fishery resources, and production of these toxins by this group of cyanobacteria makes them important organisms for continued observation and study (4, 33, 36). Much effort has been expended over the past 15 years to characterize the genomic and structural components of the microcystin (mcy) synthetase operon responsible for the production of microcystins. Several complete DNA sequences of the mcy synthetase operon are currently available in GenBank (3, 11, 29, 31).Although the mechanisms of microcystin production are now better understood, recent analyses of mcyA gene fragments from Lakes Erie and Ontario indicated a microcystin toxin producer of unknown phylogeny (7, 28). This discrepancy suggested a need for improved molecular characterization of naturally occurring microcystin producers, which spurred our research to identify the source of several unusual mcyA fragments from the cyanobacterial community (7, 28). It was apparent from initial sequence data that these mcyA gene fragments, termed Microcystis-like, were highly similar to those from Microcystis spp. (colonial or unicellular cocci). However, they contained a 6-nucleotide insert consistent with mcyA genes from filamentous cyanobacteria (e.g., Anabaena, Nostoc, and Planktothrix) (28). These preliminary findings suggested that these unusual mcyA fragments either came from (i) a novel species or strain, (ii) an ancestral Microcystis, (iii) the highly unlikely hybridization of colonial cocci and filamentous cyanobacteria, or (iv) a chimera of cocci and filamentous PCR products. To identify the source of these mcyA gene fragments from uncultured cyanobacteria, we used culture-independent methods to amplify and isolate long regions of the mcy synthetase operon for the simultaneous analysis of two genes, mcyA and mcyB, in one individual from a population. This approach ensures that both genes are contained on the same DNA molecule, thus allowing for more continuous sequence information to use in comparative phylogenetic analyses than previously described. We also envisioned that this mcy gene combination would provide an improved diagnostic tool for determining the genetic potential of naturally occurring cyanobacteria to produce specific microcystin variants by comparing the phylogenetic marker in mcyA to the predictor of amino acid incorporation (via an adenylation domain) in mcyB1.     

Myostatin基因突变激发骨骼肌发育机制研究进展

肌肉生长抑制素基因（myostatin,MSTN）是骨骼肌发育的负调节因子,在不同物种中具有高度保守性。自然突变或通过基因编辑技术对该基因进行操作,均可以获得肌肉异常发达的动物个体。研究表明,MSTN基因突变可以通过多种调控途径影响肌肉发育过程。因此,从成肌细胞增殖、分化、蛋白质合成分解代谢、组蛋白修饰以及巨噬细胞极化等5个方面对MSTN突变促进肌肉发育的机理进行综述,以期为农业动物育种新材料生产及重大恶病质的治疗提供借鉴。     

Identification of Specific Genes and Pathways by a Comparative Transcriptomic Study of Hypodermal and Body Muscle Development

Russian Journal of Developmental Biology - The development of diverse tissues may be determined by cell fate at the level of cell lineage, especially for hypodermal (Hyp) and body muscle (BM)...     

甘蓝型油菜ZF-HD基因家族的鉴定与系统进化分析

ZF-HD是一类植物特有的转录因子, 在植物生长发育及胁迫响应过程中发挥重要作用。利用生物信息学方法, 在甘蓝型油菜(Brassica napus)基因组中鉴定到62个ZF-HD基因, 其中83.9%的基因缺乏内含子, 93.5%的BnZF-HD等电点大于7, 预测定位于细胞核的蛋白大多由100个以上氨基酸组成。根据进化关系可将其分为6个亚群, 在每个亚群中, 甘蓝(B. oleracea)和白菜(B. rapa)的ZF-HD基因数量相等或近似相等, 而甘蓝型油菜的ZF-HD基因数量接近或等同于甘蓝和白菜的ZF-HD基因数量之和。同一亚群的motif数量和类型高度相似。共线性分析结果显示, 全基因组三倍体化使ZF-HD基因在二倍体祖先得到扩张, 而异源多倍体化又进一步使甘蓝型油菜ZF-HD基因家族扩张。K_a/K_s值说明大多数ZF-HD基因在进化过程中受到了纯化选择。所有BnZF-HD基因都具有光响应元件, 2/3的基因具有MeJA、ABA和厌氧诱导顺式作用元件, 推测这些基因可能参与相关生物学过程。研究结果为进一步挖掘该家族基因的生物学功能奠定基础, 同时为揭示多基因家族在异源多倍体中的进化式样提供借鉴。