Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Diversity of Mycoplasma hyopneumoniae in pig farms revealed by direct molecular typing of clinical material

Mycoplasma hyopneumoniae is the etiological agent of enzootic pneumonia in swine. Various reports indicate that different strains are circulating in the swine population. We investigated the variety of M. hyopneumoniae strains by a newly developed genetic typing method based on the polyserine repeat motif of the LppS homolog P146. PCR amplification using M. hyopneumoniae specific, conserved primers flanking the region encoding the repeat motif, followed by sequencing and cluster analysis was carried out. The study included strains isolated from different geographic regions as well as lysates from lung swabs from a series of pig farms in Switzerland. High diversity of M. hyopneumoniae was observed but farms being in close geographic or operative contact generally seemed to be affected by the same strains. Moreover, analysis of multiple samples from single pig farms indicated that these harbored the same, farm-specific strain. The results indicate that multiple strains of M. hyopneumoniae are found in the swine population but that specific strains or clones are responsible for local outbreaks. The method presented is a highly reproducible epidemiologic tool allowing direct typing of M. hyopneumoniae from clinical material without prior isolation and cultivation of strains.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

Whole-body MRI of multiple myeloma: comparison of different MRI sequences in assessment of different growth patterns

PURPOSE: To determine sensitivity, specificity and inter-observer variability of different whole-body MRI (WB-MRI) sequences in patients with multiple myeloma (MM). METHODS AND MATERIALS: WB-MRI using a 1.5T MRI scanner was performed in 23 consecutive patients (13 males, 10 females; mean age 63+/-12 years) with histologically proven MM. All patients were clinically classified according to infiltration (low-grade, n=7; intermediate-grade, n=7; high-grade, n=9) and to the staging system of Durie and Salmon PLUS (stage I, n=12; stage II, n=4; stage III, n=7). The control group consisted of 36 individuals without malignancy (25 males, 11 females; mean age 57+/-13 years). Two observers independently evaluated the following WB-MRI sequences: T1w-TSE (T1), T2w-TIRM (T2), and the combination of both sequences, including a contrast-enhanced T1w-TSE with fat-saturation (T1+/-CE/T2). They had to determine growth patterns (focal and/or diffuse) and the MRI sequence that provided the highest confidence level in depicting the MM lesions. Results were calculated on a per-patient basis. RESULTS: Visual detection of MM was as follows: T1, 65% (sensitivity)/85% (specificity); T2, 76%/81%; T1+/-CE/T2, 67%/88%. Inter-observer variability was as follows: T1, 0.3; T2, 0.55; T1+/-CE/T2, 0.55. Sensitivity improved depending on infiltration grade (T1: 1=60%; 2=36%; 3=83%; T2: 1=70%; 2=71%; 3=89%; T1+/-CE/T2: 1=50%; 2=50%; 3=89%) and clinical stage (T1: 1=58%; 2=63%; 3=79%; T2: 1=58%; 2=88%; 3=100%; T1+/-CE/T2: 1=50%; 2=63%; 3=100%). T2w-TIRM sequences achieved the best reliability in depicting the MM lesions (65% in the mean of both readers). CONCLUSIONS: T2w-TIRM sequences achieved the highest level of sensitivity and best reliability, and thus might be valuable for initial assessment of MM. For an exact staging and grading the examination protocol should encompass unenhanced and enhanced T1w-MRI sequences, in addition to T2w-TIRM.

Functional properties of multiple isoforms of human divalent metal-ion transporter 1 (DMT1)

DMT1 (divalent metal-ion transporter 1) is a widely expressed metal-ion transporter that is vital for intestinal iron absorption and iron utilization by most cell types throughout the body, including erythroid precursors. Mutations in DMT1 cause severe microcytic anaemia in animal models. Four DMT1 isoforms that differ in their N- and C-termini arise from mRNA transcripts that vary both at their 5'-ends (starting in exon 1A or exon 1B) and at their 3'-ends giving rise to mRNAs containing (+) or lacking (-) the 3'-IRE (iron-responsive element) and resulting in altered C-terminal coding sequences. To determine whether these variations result in functional differences between isoforms, we explored the functional properties of each isoform using the voltage clamp and radiotracer assays in cRNA-injected Xenopus oocytes. 1A/IRE+-DMT1 mediated Fe2+-evoked currents that were saturable (K(0.5)(Fe) approximately 1-2 microM), temperature-dependent (Q10 approximately 2), H+-dependent (K(0.5)(H) approximately 1 muM) and voltage-dependent. 1A/IRE+-DMT1 exhibited the provisional substrate profile (ranked on currents) Cd2+, Co2+, Fe2+, Mn2+>Ni2+, V3+>>Pb2+. Zn2+ also evoked large currents; however, the zinc-evoked current was accounted for by H+ and Cl- conductances and was not associated with significant Zn2+ transport. 1B/IRE+-DMT1 exhibited the same substrate profile, Fe2+ affinity and dependence on the H+ electrochemical gradient. Each isoform mediated 55Fe2+ uptake and Fe2+-evoked currents at low extracellular pH. Whereas iron transport activity varied markedly between the four isoforms, the activity for each correlated with the density of anti-DMT1 immunostaining in the plasma membrane, and the turnover rate of the Fe2+ transport cycle did not differ between isoforms. Therefore all four isoforms of human DMT1 function as metal-ion transporters of equivalent efficiency. Our results reveal that the N- and C-terminal sequence variations among the DMT1 isoforms do not alter DMT1 functional properties. We therefore propose that these variations serve as tissue-specific signals or cues to direct DMT1 to the appropriate subcellular compartments (e.g. in erythroid cells) or the plasma membrane (e.g. in intestine).

Registration of Sub-Sequence and Multi-Camera Reconstructions for Camera Motion Estimation

This paper presents different application scenarios for which the registration of sub-sequence reconstructions or multi-camera reconstructions is essential for successful camera motion estimation and 3D reconstruction from video. The registration is achieved by merging unconnected feature point tracks between the reconstructions. One application is drift removal for sequential camera motion estimation of long sequences. The state-of-the-art in drift removal is to apply a RANSAC approach to find unconnected feature point tracks. In this paper an alternative spectral algorithm for pairwise matching of unconnected feature point tracks is used. It is then shown that the algorithms can be combined and applied to novel scenarios where independent camera motion estimations must be registered into a common global coordinate system. In the first scenario multiple moving cameras, which capture the same scene simultaneously, are registered. A second new scenario occurs in situations where the tracking of feature points during sequential camera motion estimation fails completely, e.g., due to large occluding objects in the foreground, and the unconnected tracks of the independent reconstructions must be merged. In the third scenario image sequences of the same scene, which are captured under different illuminations, are registered. Several experiments with challenging real video sequences demonstrate that the presented techniques work in practice.

Clustering subspecies of Aeromonas salmonicida using IS630 typing.

BACKGROUND The insertion element IS630 found in Aeromonas salmonicida belongs to the IS630-Tc1-mariner superfamily of transposons. It is present in multiple copies and represents approximately half of the IS present in the genome of A. salmonicida subsp. salmonicida A449. RESULTS By using High Copy Number IS630 Restriction Fragment Length Polymorphism (HCN-IS630-RFLP), strains of various subspecies of Aeromonas salmonicida showed conserved or clustering patterns, thus allowing their differentiation from each other. Fingerprints of A. salmonicida subsp. salmonicida showed the highest homogeneity while 'atypical' A. salmonicida strains were more heterogeneous. IS630 typing also differentiated A. salmonicida from other Aeromonas species. The copy number of IS630 in Aeromonas salmonicida ranges from 8 to 35 and is much lower in other Aeromonas species. CONCLUSIONS HCN-IS630-RFLP is a powerful tool for subtyping of A. salmonicida. The high stability of IS630 insertions in A. salmonicida subsp. salmonicida indicates that it might have played a role in pathoadaptation of A. salmonicida which has reached an optimal configuration in the highly virulent and specific fish pathogen A. salmonicida subsp. salmonicida.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Genetically controlled variation of "acid" beta-galactosidase detected in Rattus norvegicus by isoelectric focusing.

Two genetically variant forms of rat "acid" beta-galactosidase were found to differ in isoelectric point and pH dependence, but not in thermostability or sensitivity to inhibition by p-mercuribenzoate (PMB). The results of two backcrosses and an intercross indicated that the isoelectric focusing phenotypes are controlled by two codominant alleles at a single autosomal locus, for which we propose the name Glb-1. No significant linkage between Glb-1 and albino (LG I), brown (LG II), or hooded (LG VI) was observed. Strain-specific differences in total levels of kidney beta-galactosidase were detected, but it is not yet known whether the variation is controlled by genes linked to Glb-1. Experiments in which organ homogenates were incubated with neuraminidase indicated that the genetically variant forms do not result from differences in sialylation, though sialylation does appear to be largely responsible for the presence of multiple bands within each phenotype and for differences in the banding patterns of beta-galactosidases derived from different organs. The beta-galactosidase present in the bands used for Glb-1 typing resembles human GM1 gangliosidase (GLB1) with respect to pH optimum, substrate specificity, and susceptibility to inhibition by PMB. It also appears that Glb-1 is homologous with the Bgl-e locus of the mouse. In rats as in mice the genetically variant bands of beta-galactosidase are active at acid pH and have relatively high isoelectric points. In both species these bands are readily detectable in kidney homogenates, and can be revealed in homogenates of liver or spleen following treatment with neuraminidase. The presence of the same beta-galactosidase bands in homogenates of rat kidney and small intestine as well as in neuraminidase-treated homogenates of liver and spleen suggests that the Glb-1 variants differ by one or more point mutations in the structural gene for "acid" beta-galactosidase.

Multiple independent genetic factors at NOS1AP modulate the QT interval in a multi-ethnic population.

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6 x 10(-5)) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP x sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63 x 10(-8)), as well as the sex-interaction with rs16847548 (P = 8.68 x 10(-6)). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.

MULTIPLE POSTTRANSCRIPTIONAL REGULATORY FEATURES CONTROL EXPRESSION OF ETHANOLAMINE UTILIZATION GENES IN ENTEROCOCCUS FAECALIS

Enterococcus faecalis is a Gram-positive bacterium that lives as a commensal organism in the mammalian gastrointestinal tract, but can behave as an opportunistic pathogen. Our lab discovered that mutation of the eutK gene attenuates virulence of E. faecalis in the C. elegans model host. eutK is part of the ethanolamine metabolic pathway which was previously unknown in E. faecalis. I discovered the presence of two unique posttranscriptional regulatory features that control expression of eut locus genes. The first feature I found is an AdoCBL riboswitch, a cis-acting RNA regulatory element that acts as a positive regulator of gene expression. The second feature I discovered is a unique two-component system, EutVW. The EutV response regulator contains an ANTAR family domain, which binds RNA to trigger transcriptional antitermination. I determined that induction of expression of several genes in the eut locus is dependent on ethanolamine, AdoCBL and the two-component system. AdoCBL and ethanolamine are both required for induction of eut locus gene expression. Additionally, I discovered eutG is regulated by a unique mechanism of antitermination. Both the AdoCBL riboswitch and EutV response regulator control the expression of the downstream gene eutG. EutV potentially acts through a novel antitermination mechanism in which a dimer of EutV binds to a pair of mRNA stem loops forming an antitermination complex. My data show a unique mechanism by which two environmental signals are integrated by two different posttranscriptional regulators to regulate a single locus.

3D MPRAGE improves classification of cortical lesions in multiple sclerosis.

BACKGROUND: Gray matter lesions are known to be common in multiple sclerosis (MS) and are suspected to play an important role in disease progression and clinical disability. A combination of magnetic resonance imaging (MRI) techniques, double-inversion recovery (DIR), and phase-sensitive inversion recovery (PSIR), has been used for detection and classification of cortical lesions. This study shows that high-resolution three-dimensional (3D) magnetization-prepared rapid acquisition with gradient echo (MPRAGE) improves the classification of cortical lesions by allowing more accurate anatomic localization of lesion morphology. METHODS: 11 patients with MS with previously identified cortical lesions were scanned using DIR, PSIR, and 3D MPRAGE. Lesions were identified on DIR and PSIR and classified as purely intracortical or mixed. MPRAGE images were then examined, and lesions were re-classified based on the new information. RESULTS: The high signal-to-noise ratio, fine anatomic detail, and clear gray-white matter tissue contrast seen in the MPRAGE images provided superior delineation of lesion borders and surrounding gray-white matter junction, improving classification accuracy. 119 lesions were identified as either intracortical or mixed on DIR/PSIR. In 89 cases, MPRAGE confirmed the classification by DIR/PSIR. In 30 cases, MPRAGE overturned the original classification. CONCLUSION: Improved classification of cortical lesions was realized by inclusion of high-spatial resolution 3D MPRAGE. This sequence provides unique detail on lesion morphology that is necessary for accurate classification.

Bicarbonate enhances expression of the endocarditis and biofilm associated pilus locus, ebpR-ebpABC, in Enterococcus faecalis.

BACKGROUND: We previously identified ebpR, encoding a potential member of the AtxA/Mga transcriptional regulator family, and showed that it is important for transcriptional activation of the Enterococcus faecalis endocarditis and biofilm associated pilus operon, ebpABC. Although ebpR is not absolutely essential for ebpABC expression (100-fold reduction), its deletion led to phenotypes similar to those of an ebpABC mutant such as absence of pili at the cell surface and, consequently, reduced biofilm formation. A non-piliated ebpABC mutant has been shown to be attenuated in a rat model of endocarditis and in a murine urinary tract infection model, indicating an important participation of the ebpR-ebpABC locus in virulence. However, there is no report relating to the environmental conditions that affect expression of the ebpR-ebpABC locus. RESULTS: In this study, we examined the effect of CO2/HCO3(-), pH, and the Fsr system on the ebpR-ebpABC locus expression. The presence of 5% CO2/0.1 M HCO3(-) increased ebpR-ebpABC expression, while the Fsr system was confirmed to be a weak repressor of this locus. The mechanism by which the Fsr system repressed the ebpR-ebpABC locus expression appears independent of the effects of CO2(-) bicarbonate. Furthermore, by using an ebpA::lacZ fusion as a reporter, we showed that addition of 0.1 M sodium bicarbonate to TSBG (buffered at pH 7.5), but not the presence of 5% CO2, induced ebpA expression in TSBG broth. In addition, using microarray analysis, we found 73 genes affected by the presence of sodium bicarbonate (abs(fold) > 2, P < 0.05), the majority of which belong to the PTS system and ABC transporter families. Finally, pilus production correlated with ebpA mRNA levels under the conditions tested. CONCLUSIONS: This study reports that the ebp locus expression is enhanced by the presence of bicarbonate with a consequential increase in the number of cells producing pili. Although the molecular basis of the bicarbonate effect remains unclear, the pathway is independent of the Fsr system. In conclusion, E. faecalis joins the growing family of pathogens that regulates virulence gene expression in response to bicarbonate and/or CO2.

Genome sequence of Fusobacterium nucleatum subspecies polymorphum - a genetically tractable fusobacterium.

Fusobacterium nucleatum is a prominent member of the oral microbiota and is a common cause of human infection. F. nucleatum includes five subspecies: polymorphum, nucleatum, vincentii, fusiforme, and animalis. F. nucleatum subsp. polymorphum ATCC 10953 has been well characterized phenotypically and, in contrast to previously sequenced strains, is amenable to gene transfer. We sequenced and annotated the 2,429,698 bp genome of F. nucleatum subsp. polymorphum ATCC 10953. Plasmid pFN3 from the strain was also sequenced and analyzed. When compared to the other two available fusobacterial genomes (F. nucleatum subsp. nucleatum, and F. nucleatum subsp. vincentii) 627 open reading frames unique to F. nucleatum subsp. polymorphum ATCC 10953 were identified. A large percentage of these mapped within one of 28 regions or islands containing five or more genes. Seventeen percent of the clustered proteins that demonstrated similarity were most similar to proteins from the clostridia, with others being most similar to proteins from other gram-positive organisms such as Bacillus and Streptococcus. A ten kilobase region homologous to the Salmonella typhimurium propanediol utilization locus was identified, as was a prophage and integrated conjugal plasmid. The genome contains five composite ribozyme/transposons, similar to the CdISt IStrons described in Clostridium difficile. IStrons are not present in the other fusobacterial genomes. These findings indicate that F. nucleatum subsp. polymorphum is proficient at horizontal gene transfer and that exchange with the Firmicutes, particularly the Clostridia, is common.

Automatic quantitation of multiple sclerosis lesions on MR images

A two-pronged approach for the automatic quantitation of multiple sclerosis (MS) lesions on magnetic resonance (MR) images has been developed. This method includes the design and use of a pulse sequence for improved lesion-to-tissue contrast (LTC) and seeks to identify and minimize the sources of false lesion classifications in segmented images. The new pulse sequence, referred to as AFFIRMATIVE (Attenuation of Fluid by Fast Inversion Recovery with MAgnetization Transfer Imaging with Variable Echoes), improves the LTC, relative to spin-echo images, by combining Fluid-Attenuated Inversion Recovery (FLAIR) and Magnetization Transfer Contrast (MTC). In addition to acquiring fast FLAIR/MTC images, the AFFIRMATIVE sequence simultaneously acquires fast spin-echo (FSE) images for spatial registration of images, which is necessary for accurate lesion quantitation. Flow has been found to be a primary source of false lesion classifications. Therefore, an imaging protocol and reconstruction methods are developed to generate "flow images" which depict both coherent (vascular) and incoherent (CSF) flow. An automatic technique is designed for the removal of extra-meningeal tissues, since these are known to be sources of false lesion classifications. A retrospective, three-dimensional (3D) registration algorithm is implemented to correct for patient movement which may have occurred between AFFIRMATIVE and flow imaging scans. Following application of these pre-processing steps, images are segmented into white matter, gray matter, cerebrospinal fluid, and MS lesions based on AFFIRMATIVE and flow images using an automatic algorithm. All algorithms are seamlessly integrated into a single MR image analysis software package. Lesion quantitation has been performed on images from 15 patient volunteers. The total processing time is less than two hours per patient on a SPARCstation 20. The automated nature of this approach should provide an objective means of monitoring the progression, stabilization, and/or regression of MS lesions in large-scale, multi-center clinical trials. ^