The first mitochondrial genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997 (Diptera: Sciomyzoidea: Sepsidae), with mitochondrial genome phylogeny of cyclorrhapha

Sepsid flies (Diptera: Sepsidae) are important model insects for sexual selection research. In order to develop mitochondrial (mt) genome data for this significant group, we sequenced the first complete mt genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997. The circular 15,878 bp mt genome is typical of Diptera, containing all 37 genes usually present in bilaterian animals. We discovered inaccurate annotations of fly mt genomes previously deposited on GenBank and thus re-annotated all published mt genomes of Cyclorrhapha. These re-annotations were based on comparative analysis of homologous genes, and provide a statistical analysis of start and stop codon positions. We further detected two 18 bp of conserved intergenic sequences from tRNAGlu-tRNAPhe and ND1-tRNASer(UCN) across Cyclorrhapha, which are the mtTERM binding site motifs. Additionally, we compared automated annotation software MITOS with hand annotation method. Phylogenetic trees based on the mt genome data from Cyclorrhapha were inferred by Maximum-likelihood and Bayesian methods, strongly supported a close relationship between Sepsidae and the Tephritoidea.

Oldest pathology in a tetrapod bone illuminates the origin of terrestrial vertebrates

The origin of terrestrial tetrapods was a key event in vertebrate evolution, yet how and when it occurred remains obscure, due to scarce fossil evidence. Here, we show that the study of palaeopathologies, such as broken and healed bones, can help elucidate poorly understood behavioural transitions such as this. Using high-resolution finite element analysis, we demonstrate that the oldest known broken tetrapod bone, a radius of the primitive stem tetrapod Ossinodus pueri from the mid-Viséan (333 million years ago) of Australia, fractured under a high-force, impact-type loading scenario. The nature of the fracture suggests that it most plausibly occurred during a fall on land. Augmenting this are new osteological observations, including a preferred directionality to the trabecular architecture of cancellous bone. Together, these results suggest that Ossinodus, one of the first large (>2m length) tetrapods, spent a significant proportion of its life on land. Our findings have important implications for understanding the temporal, biogeographical and physiological contexts under which terrestriality in vertebrates evolved. They push the date for the origin of terrestrial tetrapods further back into the Carboniferous by at least two million years. Moreover, they raise the possibility that terrestriality in vertebrates first evolved in large tetrapods in Gondwana rather than in small European forms, warranting a re-evaluation of this important evolutionary event.

High connectivity in rastrelliger kanagurta: Influence of historical signatures and migratory behaviour inferred from mtDNA Cytochrome b

Phylogeographic patterns and population structure of the pelagic Indian mackerel, Rastrelliger kanagurta were examined in 23 populations collected from the Indonesian-Malaysian Archipelago (IMA) and the West Indian Ocean (WIO). Despite the vast expanse of the IMA and neighbouring seas, no evidence for geographical structure was evident. An indication that R. kanagurta populations across this region are essentially panmictic. This study also revealed that historical isolation was insufficient for R. kanagurta to attain migration drift equilibrium. Two distinct subpopulations were detected between the WIO and the IMA (and adjacent populations); interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the IMA and WIO regions suggest historical isolation as a result of fluctuations in sea levels during the late Pleistocene. This occurrence resulted in the evolution of a phylogeographic break for this species to the north of the Andaman Sea.

Comparative analysis and distribution of Omega-3 lcPUFA Biosynthesis genes in marine molluscs

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.

Resampling-based approaches to study variation in morphological modularity

Modularity has been suggested to be connected to evolvability because a higher degree of independence among parts allows them to evolve as separate units. Recently, the Escoufier RV coefficient has been proposed as a measure of the degree of integration between modules in multivariate morphometric datasets. However, it has been shown, using randomly simulated datasets, that the value of the RV coefficient depends on sample size. Also, so far there is no statistical test for the difference in the RV coefficient between a priori defined groups of observations. Here, we (1), using a rarefaction analysis, show that the value of the RV coefficient depends on sample size also in real geometric morphometric datasets; (2) propose a permutation procedure to test for the difference in the RV coefficient between a priori defined groups of observations; (3) show, through simulations, that such a permutation procedure has an appropriate Type I error; (4) suggest that a rarefaction procedure could be used to obtain sample-size-corrected values of the RV coefficient; and (5) propose a nearest-neighbor procedure that could be used when studying the variation of modularity in geographic space. The approaches outlined here, readily extendable to non-morphometric datasets, allow study of the variation in the degree of integration between a priori defined modules. A Java application – that will allow performance of the proposed test using a software with graphical user interface – has also been developed and is available at the Morphometrics at Stony Brook Web page (http://life.bio.sunysb.edu/morph/).

The gut microbial community of midas cichlid fish in repeatedly evolved limnetic-benthic species pairs

Gut bacterial communities are now known to influence a range of fitness related aspects of organisms. But how different the microbial community is in closely related species, and if these differences can be interpreted as adaptive is still unclear. In this study we compared microbial communities in two sets of closely related sympatric crater lake cichlid fish species pairs that show similar adaptations along the limnetic-benthic axis. The gut microbial community composition differs in the species pair inhabiting the older of two crater lakes. One major difference, relative to other fish, is that in these cichlids that live in hypersaline crater lakes, the microbial community is largely made up of Oceanospirillales (52.28%) which are halotolerant or halophilic bacteria. This analysis opens up further avenues to identify candidate symbiotic or co-evolved bacteria playing a role in adaptation to similar diets and life-styles or even have a role in speciation. Future functional and phylosymbiotic analyses might help to address these issues.

Bat species comparisons based on external morphology: A test of traditional versus geometric morphometric approaches

External morphology is commonly used to identify bats as well as to investigate flight and foraging behavior, typically relying on simple length and area measures or ratios. However, geometric morphometrics is increasingly used in the biological sciences to analyse variation in shape and discriminate among species and populations. Here we compare the ability of traditional versus geometric morphometric methods in discriminating between closely related bat species – in this case European horseshoe bats (Rhinolophidae, Chiroptera) – based on morphology of the wing, body and tail. In addition to comparing morphometric methods, we used geometric morphometrics to detect interspecies differences as shape changes. Geometric morphometrics yielded improved species discrimination relative to traditional methods. The predicted shape for the variation along the between group principal components revealed that the largest differences between species lay in the extent to which the wing reaches in the direction of the head. This strong trend in interspecific shape variation is associated with size, which we interpret as an evolutionary allometry pattern.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.66 x 10-10, odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10-4. OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6×10-5, OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10-5, OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.261025, OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1×10-4, OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10-5, OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9×10-4, OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

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.

N-ethyl-N-nitrosourea (ENU) induced mutations within the Klotho gene lead to ectopic calcification and reduced lifespan in mouse models

Ectopic calcification (EC), which is the pathological deposition of calcium and phosphate in extra-skeletal tissues, may be associated with hypercalcaemic and hyperphosphataemic disorders, or it may occur in the absence of metabolic abnormalities. In addition, EC may be inherited as part of several monogenic disorders and studies of these have provided valuable insights into the metabolic pathways regulating mineral metabolism. For example, studies of tumoural calcinosis, a disorder characterised by hyperphosphataemia and progressive EC, have revealed mutations of fibroblast growth factor 23 (FGF23), polypeptide N-acetyl galactosaminyltransferase 3 (GALNT3) and klotho (KL), which are all part of a phosphate-regulating pathway. However, such studies in humans are limited by the lack of available large families with EC, and to facilitate such studies we assessed the progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) for EC. This identified two mutants with autosomal recessive forms of EC, and reduced lifespan, designated Ecalc1 and Ecalc2. Genetic mapping localized the Ecalc1 and Ecalc2 loci to a 11.0 Mb region on chromosome 5 that contained the klotho gene (Kl), and DNA sequence analysis identified nonsense (Gln203Stop) and missense (Ile604Asn) Kl mutations in Ecalc1 and Ecalc2 mice, respectively. The Gln203Stop mutation, located in KL1 domain, was severely hypomorphic and led to a 17-fold reduction of renal Kl expression. The Ile604Asn mutation, located in KL2 domain, was predicted to impair klotho protein stability and in vitro expression studies in COS-7 cells revealed endoplasmic reticulum retention of the Ile604Asn mutant. Further phenotype studies undertaken in Ecalc1 (kl203X/203X) mice demonstrated elevations in plasma concentrations of phosphate, FGF23 and 1,25-dihydroxyvitamin D. Thus, two allelic variants of Kl that develop EC and represent mouse models for tumoural calcinosis have been established. © 2015 Esapa et al.

A mouse with an N-ethyl-N-nitrosourea (ENU) induced Trp589Arg Galnt3 mutation represents a model for hyperphosphataemic familial tumoural calcinosis

Mutations of UDP-N-acetyl-alpha-D-galactosamine polypeptide N-acetyl galactosaminyl transferase 3 (GALNT3) result in familial tumoural calcinosis (FTC) and the hyperostosis-hyperphosphataemia syndrome (HHS), which are autosomal recessive disorders characterised by soft-tissue calcification and hyperphosphataemia. To facilitate in vivo studies of these heritable disorders of phosphate homeostasis, we embarked on establishing a mouse model by assessing progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU), and identified a mutant mouse, TCAL, with autosomal recessive inheritance of ectopic calcification, which involved multiple tissues, and hyperphosphataemia; the phenotype was designated TCAL and the locus, Tcal. TCAL males were infertile with loss of Sertoli cells and spermatozoa, and increased testicular apoptosis. Genetic mapping localized Tcal to chromosome 2 (62.64-71.11 Mb) which contained the Galnt3. DNA sequence analysis identified a Galnt3 missense mutation (Trp589Arg) in TCAL mice. Transient transfection of wild-type and mutant Galnt3-enhanced green fluorescent protein (EGFP) constructs in COS-7 cells revealed endoplasmic reticulum retention of the Trp589Arg mutant and Western blot analysis of kidney homogenates demonstrated defective glycosylation of Galnt3 in Tcal/Tcal mice. Tcal/Tcal mice had normal plasma calcium and parathyroid hormone concentrations; decreased alkaline phosphatase activity and intact Fgf23 concentrations; and elevation of circulating 1,25-dihydroxyvitamin D. Quantitative reverse transcriptase-PCR (qRT-PCR) revealed that Tcal/Tcal mice had increased expression of Galnt3 and Fgf23 in bone, but that renal expression of Klotho, 25-hydroxyvitamin D-1α-hydroxylase (Cyp27b1), and the sodium-phosphate co-transporters type-IIa and -IIc was similar to that in wild-type mice. Thus, TCAL mice have the phenotypic features of FTC and HHS, and provide a model for these disorders of phosphate metabolism. © 2012 Esapa 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.

Polymorphisms in the receptor tyrosine kinase MERTK gene are associated with Multiple Sclerosis susceptibility

Multiple sclerosis (MS) is a debilitating, chronic demyelinating disease of the central nervous system affecting over 2 million people worldwide. The TAM family of receptor tyrosine kinases (TYRO3, AXL and MERTK) have been implicated as important players during demyelination in both animal models of MS and in the human disease. We therefore conducted an association study to identify single nucleotide polymorphisms (SNPs) within genes encoding the TAM receptors and their ligands associated with MS. Analysis of genotype data from a genome-wide association study which consisted of 1618 MS cases and 3413 healthy controls conducted by the Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) revealed several SNPs within the MERTK gene (Chromosome 2q14.1, Accession Number NG_011607.1) that showed suggestive association with MS. We therefore interrogated 28 SNPs in MERTK in an independent replication cohort of 1140 MS cases and 1140 healthy controls. We found 12 SNPs that replicated, with 7 SNPs showing p-values of less than 10-5 when the discovery and replication cohorts were combined. All 12 replicated SNPs were in strong linkage disequilibrium with each other. In combination, these data suggest the MERTK gene is a novel risk gene for MS susceptibility. © 2011 Ma et al.