Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ~2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2×10-9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3×10-12, and -0.16 SD per G allele, P = 1.2×10-15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10-9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10-6 and rs2707466: OR = 1.22, P = 7.2×10-6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16-/- mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5×10-13<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.

The MS risk allele of CD40 is associated with reduced cell-membrane bound expression in antigen presenting cells: Implications for gene function

Human genetic and animal studies have implicated the costimulatory molecule CD40 in the development of multiple sclerosis (MS). We investigated the cell specific gene and protein expression variation controlled by the CD40 genetic variant(s) associated with MS, i.e. the T-allele at rs1883832. Previously we had shown that the risk allele is expressed at a lower level in whole blood, especially in people with MS. Here, we have defined the immune cell subsets responsible for genotype and disease effects on CD40 expression at the mRNA and protein level. In cell subsets in which CD40 is most highly expressed, B lymphocytes and dendritic cells, the MS-associated risk variant is associated with reduced CD40 cell-surface protein expression. In monocytes and dendritic cells, the risk allele additionally reduces the ratio of expression of full-length versus truncated CD40 mRNA, the latter encoding secreted CD40. We additionally show that MS patients, regardless of genotype, express significantly lower levels of CD40 cell-surface protein compared to unaffected controls in B lymphocytes. Thus, both genotype-dependent and independent down-regulation of cell-surface CD40 is a feature of MS. Lower expression of a co-stimulator of T cell activation, CD40, is therefore associated with increased MS risk despite the same CD40 variant being associated with reduced risk of other inflammatory autoimmune diseases. Our results highlight the complexity and likely individuality of autoimmune pathogenesis, and could be consistent with antiviral and/or immunoregulatory functions of CD40 playing an important role in protection from MS. © 2015 Field et al.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Identification of mechanism, PECARN risk factors and injury patterns in severe paediatric cervical spine injuries in Queensland

Introduction Canadian C spine rule and NEXUS criteria have identified risk factors for cervical spine injury in adults but not for children. PECARN has developed an 8 variable model for cervical spine injury in children. We sought to identify the mechanism, prevalence of PECARN risk factors, injury patterns, and management of severe Paediatric cervical spine injuries presenting to the major children’s hospitals in Brisbane, Australia. Methods This a retrospective study of the children with cervical spine injuries who presented directly or were referred to the major children’s hospitals in Brisbane over 5 years. Results There were 38 patients with 18 male and 20 female.The mean age was 8.6 years. They were divided into two groups according to their age, (Group 1 < =8 years had 18 (47%) patients, while group 2 (9-15 years) had 20 (53%) patients. Motor vehicle related injuries were the most common (61%) in Group 1 while it was sporting injuries (50%) in group 2. All patients in group 1 had upper cervical injury (C0-C2) while subaxial injuries were most common in group 2 (66.6%). 82% of the patients had 2 or more PECARN risk factors. 18 children (47%) had normal neurological assessment at presentation, 6 (16%) had radicular symptoms, 11 (29%) could not be assessed as they had already been intubated due to the severity of the injury, 3 (8%) had incomplete cord injury. 29 (69%) patients had normal neurological assessment at final follow up and 2 children died from their injuries. Conclusion Our study confirms that younger children sustain upper cervical injuries most commonly secondary to motor vehicle accidents, while the older sustain subaxial injuries from sporting activities. The significant prevalence of the PECARN risk factors among this cohort of patients have led to them being incorporated into a protocol at these hospitals used to assess patients with suspected cervical spinal injury.

Identification of major loci controlling clinical manifestations of ankylosing spondylitis

Objective. To identify genomic regions linked with determinants of age at symptom onset, disease activity, and functional impairment in ankylosing spondylitis (AS). Methods. A whole genome linkage scan was performed in 188 affected sibling pair families with 454 affected individuals. Traits assessed were age at symptom onset, disease activity assessed by the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), and functional impairment assessed by the Bath Ankylosing Spondylitis Functional Index (BASFI). Parametric and nonparametric quantitative linkage analysis was performed using parameters defined in a previous segregation study. Results. Heritabilities of the traits studied in this data set were as follows: BASDAI 0.49 (P = 0.0001, 95% confidence interval [95% CI] 0.23-0.75), BASFI 0.76 (P = 10-7, 95% CI 0.49-1.0), and age at symptom onset 0.33 (P = 0.005, 95% CI 0.04-0.62). No linkage was observed between the major histocompatibility complex (MHC) and any of the traits studied (logarithm of odds [LOD] score <1.0). "Significant" linkage (LOD score 4.0) was observed between a region on chromosome 18p and the BASDAI. Age at symptom onset showed "suggestive" linkage to chromosome 11p (LOD score 3.3). Maximum linkage with the BASFI was seen at chromosome 2q (LOD score 2.9). Conclusion. In contrast to the genetic determinants of susceptibility to AS, clinical manifestations of the disease measured by the BASDAI, BASFI, and age at symptom onset are largely determined by a small number of genes not encoded within the MHC.

The association between ambient temperature and children's lung function in Baotou, China

The objective of this study is to examine the association between ambient temperature and children’s lung function in Baotou, China. We recruited 315 children (8–12 years) from Baotou, China in the spring of 2004, 2005, and 2006. They performed three successive forced expiratory measurements three times daily (morning, noon, and evening) for about 5 weeks. The highest peak expiratory flow (PEF) was recorded for each session. Daily data on ambient temperature, relative humidity, and air pollution were monitored during the same period. Mixed models with a distributed lag structure were used to examine the effects of temperature on lung function while adjusting for individual characteristics and environmental factors. Low temperatures were significantly associated with decreases in PEF. The effects lasted for lag 0–2 days. For all participants, the cumulative effect estimates (lag 0–2 days) were −1.44 (−1.93, −0.94) L/min, −1.39 (−1.92, −0.86) L/min, −1.40 (−1.97, −0.82) L/min, and −1.28 (−1.69, −0.88) L/min for morning, noon, evening, and daily mean PEF, respectively, associated with 1 °C decrease in daily mean temperature. Generally, the effects of temperature were slightly stronger in boys than in girls for noon, evening, and daily mean PEF, while the effects were stronger in girls for morning PEF. PM2.5 had joint effects with temperature on children’s PEF. Higher PM2.5 increased the impacts of low temperature. Low ambient temperatures are associated with lower lung function in children in Baotou, China. Preventive health policies will be required for protecting children from the cold weather.

Identification of bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei reveals distinct temperature-driven patterns

Change in temperature is often a major environmental factor in triggering waterborne disease outbreaks. Previous research has revealed temporal and spatial patterns of bacterial population in several aquatic ecosystems. To date, very little information is available on aquaculture environment. Here, we assessed environmental temperature effects on bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei (FASFL). Water samples were collected over a one-year period, and aquatic bacteria were characterized by polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) and 16S rDNA pyrosequencing. Resulting DGGE fingerprints revealed a specific and dynamic bacterial population structure with considerable variation over the seasonal change, suggesting that environmental temperature was a key driver of bacterial population in the FASFL. Pyrosequencing data further demonstrated substantial difference in bacterial community composition between the water at higher (WHT) and at lower (WLT) temperatures in the FASFL. Actinobacteria, Proteobacteria and Bacteroidetes were the highest abundant phyla in the FASFL, however, a large number of unclassified bacteria contributed the most to the observed variation in phylogenetic diversity. The WHT harbored remarkably higher diversity and richness in bacterial composition at genus and species levels when compared to the WLT. Some potential pathogenenic species were identified in both WHT and WLT, providing data in support of aquatic animal health management in the aquaculture industry.

The koala immunological toolkit: Sequence identification and comparison of key markers of the koala (Phascolarctos cinereus) immune response

The koala (Phascolarctos cinereus) is an Australian marsupial that continues to experience significant population declines. Infectious diseases caused by pathogens such as Chlamydia are proposed to have a major role. Very few species-specific immunological reagents are available, severely hindering our ability to respond to the threat of infectious diseases in the koala. In this study, we utilise data from the sequencing of the koala transcriptome to identify key immunological markers of the koala adaptive immune response and cytokines known to be important in the host response to chlamydial infection in other species. This report describes the identification and preliminary sequence analysis of (1) T lymphocyte glycoprotein markers (CD4, CD8); (2) IL-4, a marker for the Th2 response; (3) cytokines such as IL-6, IL-12 and IL-1β, that have been shown to have a role in chlamydial clearance and pathology in other hosts; and (4) the sequences for the koala immunoglobulins, IgA, IgG, IgE and IgM. These sequences will enable the development of a range of immunological reagents for understanding the koala’s innate and adaptive immune responses, while also providing a resource that will enable continued investigations into the origin and evolution of the marsupial immune system.

VIP Barcoding: Composition vector-based software for rapid species identification based on DNA barcoding

Species identification based on short sequences of DNA markers, that is, DNA barcoding, has emerged as an integral part of modern taxonomy. However, software for the analysis of large and multilocus barcoding data sets is scarce. The Basic Local Alignment Search Tool (BLAST) is currently the fastest tool capable of handling large databases (e.g. >5000 sequences), but its accuracy is a concern and has been criticized for its local optimization. However, current more accurate software requires sequence alignment or complex calculations, which are time-consuming when dealing with large data sets during data preprocessing or during the search stage. Therefore, it is imperative to develop a practical program for both accurate and scalable species identification for DNA barcoding. In this context, we present VIP Barcoding: a user-friendly software in graphical user interface for rapid DNA barcoding. It adopts a hybrid, two-stage algorithm. First, an alignment-free composition vector (CV) method is utilized to reduce searching space by screening a reference database. The alignment-based K2P distance nearest-neighbour method is then employed to analyse the smaller data set generated in the first stage. In comparison with other software, we demonstrate that VIP Barcoding has (i) higher accuracy than Blastn and several alignment-free methods and (ii) higher scalability than alignment-based distance methods and character-based methods. These results suggest that this platform is able to deal with both large-scale and multilocus barcoding data with accuracy and can contribute to DNA barcoding for modern taxonomy. VIP Barcoding is free and available at http://msl.sls.cuhk.edu.hk/vipbarcoding/.

Development and validation of a questionnaire to measure health professionals' attitudes toward identification of female victimes of domestic violence

Background Domestic violence against women is a major public health problem and violations of women’s human rights. Health professionals could play an important role in screening for the victims. From the evidence to date, it is unclear whether health professionals do play an active role in identification of the victims. Objectives To develop a reliable and valid instrument to measure health professionals’ attitude to identifying female victims of domestic violence. Methods A primary questionnaire was constructed in accordance with established guidelines using the Theory of Planned Behaviour Ajzen (1975) to develop an instrument to measure health professionals’ attitudes in identifying female victim of DV. An expert panel was used to establish content validity. Focus groups amongst a group of health professionals (N = 5) of the target population were performed to confirm face validity. A pilot study (N = 30 nurses and doctors) was undertaken to elicit the feasibility and reliability of the questionnaire. The questionnaire was also administered a second time after one week to check the stability of the tests. Results Feedbacks of the expert panel’s and group discussion confirmed that the questionnaire had the content and face validity. Cronbach’s alpha values for all the items were greater than 0.7. Strong correlations between the direct and indirect measures confirmed that the indirect measures were well constructed. High test-retest correlations confirmed that the measures were reliable in the sense of temporal stability. Significance This tool has the potential to be used by researchers in expanding the knowledge base in this important area.

The function of the RNA-binding protein TEL1 in moss reveals ancient regulatory mechanisms of shoot development

The shoot represents the basic body plan in land plants. It consists of a repeated structure composed of stems and leaves. Whereas vascular plants generate a shoot in their diploid phase, non-vascular plants such as mosses form a shoot (called the gametophore) in their haploid generation. The evolution of regulatory mechanisms or genetic networks used in the development of these two kinds of shoots is unclear. TERMINAL EAR1-like genes have been involved in diploid shoot development in vascular plants. Here, we show that disruption of PpTEL1 from the moss Physcomitrella patens, causes reduced protonema growth and gametophore initiation, as well as defects in gametophore development. Leafy shoots formed on ΔTEL1 mutants exhibit shorter stems with more leaves per shoot, suggesting an accelerated leaf initiation (shortened plastochron), a phenotype shared with the Poaceae vascular plants TE1 and PLA2/LHD2 mutants. Moreover, the positive correlation between plastochron length and leaf size observed in ΔTEL1 mutants suggests a conserved compensatory mechanism correlating leaf growth and leaf initiation rate that would minimize overall changes in plant biomass. The RNA-binding protein encoded by PpTEL1 contains two N-terminus RNA-recognition motifs, and a third C-terminus non-canonical RRM, specific to TEL proteins. Removal of the PpTEL1 C-terminus (including this third RRM) or only 16–18 amino acids within it seriously impairs PpTEL1 function, suggesting a critical role for this third RRM. These results show a conserved function of the RNA-binding PpTEL1 protein in the regulation of shoot development, from early ancestors to vascular plants, that depends on the third TEL-specific RRM.

Autosomal dominant spondylocostal dysostosis is caused by mutation in TBX6

In humans, congenital spinal defects occur with an incidence of 0.5-1 per 1000 live births. One of the most severe syndromes with such defects is spondylocostal dysostosis (SCD). Over the past decade, the genetic basis of several forms of autosomal recessive SCD cases has been solved with the identification of four causative genes (DLL3, MESP2, LFNG and HES7). Autosomal dominant forms of SCD have also been reported, but to date no genetic etiology has been described for these. Here, we have used exome capture and next-generation sequencing to identify a stoploss mutation in TBX6 that segregates with disease in two generations of one family. We show that this mutation has a deleterious effect on the transcriptional activation activity of the TBX6 protein, likely due to haploinsufficiency. In mouse, Tbx6 is essential for the patterning of the vertebral precursor tissues, somites; thus, mutation of TBX6 is likely to be causative of SCD in this family. This is the first identification of the genetic cause of an autosomal dominant form of SCD, and also demonstrates the potential of exome sequencing to identify genetic causes of dominant diseases even in small families with few affected individuals.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.