Asymptomatic bacteriuria Escherichia coli strain 83972 carries mutations in the foc locus and is unable to express F1C fimbriae

Escherichia coli is the most common organism associated with asymptomatic bacteriuria (ABU). In contrast to uropathogenic E. coli (UPEC), which causes symptomatic urinary tract infection (UTI), very little is known about the mechanisms by which these strains colonize the urinary tract. Bacterial adhesion conferred by specific surface-associated adhesins is normally considered as a prerequisite for colonization of the urinary tract. The prototype ABU E coli strain 83972 was originally isolated from a girl who had carried it asymptomatically for 3 years. This study characterized the molecular status of one of the primary adhesion factors known to be associated with UTI, namely F1C fimbriae, encoded by the foc gene cluster. F1C fimbriae recognize receptors present in the human kidney and bladder. Expression of the foc genes was found to be up-regulated in human urine. It was also shown that although strain 83972 contains a seemingly intact foc gene cluster, F1C fimbriae are not expressed. Sequencing and genetic complementation revealed that the focD gene, encoding a component of the F1C transport and assembly system, was non-functional, explaining the inability of strain 83972 to express this adhesin. The data imply that E. coli 83972 has lost its ability to express this important colonization factor as a result of host-driven evolution. The ancestor of the strain seems to have been a pyelonephritis strain of phylogenetic group B2. Strain 83972 therefore represents an example of bacterial adaptation from pathogenicity to commensalism through virulence factor loss.

A diffusion approach to approximating preservation probabilities for gene duplicates

Consider a haploid population and, within its genome, a gene whose presence is vital for the survival of any individual. Each copy of this gene is subject to mutations which destroy its function. Suppose one member of the population somehow acquires a duplicate copy of the gene, where the duplicate is fully linked to the original gene's locus. Preservation is said to occur if eventually the entire population consists of individuals descended from this one which initially carried the duplicate. The system is modelled by a finite state-space Markov process which in turn is approximated by a diffusion process, whence an explicit expression for the probability of preservation is derived. The event of preservation can be compared to the fixation of a selectively neutral gene variant initially present in a single individual, the probability of which is the reciprocal of the population size. For very weak mutation, this and the probability of preservation are equal, while as mutation becomes stronger, the preservation probability tends to double this reciprocal. This is in excellent agreement with simulation studies.

GPNN: Power studies and applications of a neural network method for detecting gene-gene interactions in studies of human disease

Background: The identification and characterization of genes that influence the risk of common, complex multifactorial disease primarily through interactions with other genes and environmental factors remains a statistical and computational challenge in genetic epidemiology. We have previously introduced a genetic programming optimized neural network (GPNN) as a method for optimizing the architecture of a neural network to improve the identification of gene combinations associated with disease risk. The goal of this study was to evaluate the power of GPNN for identifying high-order gene-gene interactions. We were also interested in applying GPNN to a real data analysis in Parkinson's disease. Results: We show that GPNN has high power to detect even relatively small genetic effects (2-3% heritability) in simulated data models involving two and three locus interactions. The limits of detection were reached under conditions with very small heritability (

Replication of association of IL1 gene complex members with ankylosing spondylitis in Taiwanese Chinese

Objective: To test the association of interleukin 1 ( IL1) gene family members with ankylosing spondylitis ( AS), previously reported in Europid subjects, in an ethnically remote population. Methods: 200 Taiwanese Chinese AS patients and 200 ethnically matched healthy controls were genotyped for five single nucleotide polymorphisms ( SNPs) and the IL1RN. VNTR, markers previously associated with AS. Allele, genotype, and haplotype frequencies were compared between cases and controls. Results: Association of alleles and genotypes of the markers IL1F10.3, IL1RN. 4, and IL1RN. VNTR was observed with AS ( p < 0.05). Haplotypes of pairs of these markers and of the markers IL1RN.6/1 and IL1RN.6/2 were also significantly associated with AS. The strongest associations observed were with the marker IL1RN. 4, and with the two-marker haplotype IL1RN.4-IL1RN.VNTR ( both p = 0.004). Strong linkage disequilibrium was observed between all marker pairs except those involving IL1B-511 ( D' 0.4 to 0.9, p < 0.01). Conclusions: The IL1 gene cluster is associated with AS in Taiwanese Chinese. This finding provides strong statistical support that the previously observed association of this gene cluster with AS is a true positive finding.

Genorne sequence of the thermostable Newcastle disease virus (strain I-2) reveals a possible phenotypic locus

The complete genome sequence of the Australian 1-2 heat-tolerant Newcastle disease virus (NDV) vaccine (master seed stocks) was determined and compared to the sequence of the parent virus from which it had been derived after exposure of the parent stock at 56 degrees C for 30 min. Nucleotide changes were observed at a number of positions with synonymous mutations being greater than those observed for non-synonymous mutations. Sequence data for the HN gene of a parental culture of V4 and two heat-tolerant variants of V4 were obtained. These were compared with the data for the 1-2 viruses and with published sequences for parental V4 and for a number of ND vaccine strains. Sequence analyses did not reveal the ARG 303 deletion in the HN protein, previously claimed to be responsible for the thermostable phenotype. No consistent changes were detected that would indicate involvement of the HN protein in heat resistance. The majority of alterations were observed in the L protein of the virus and it is proposed that these alterations were responsible for the heat-tolerant phenotype of the 1-2 NDV vaccine. (c) 2005 Elsevier B.V. All rights reserved.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Microvariation at the human D1S80 locus

Microvariant allelic polymorphisms have been known since 1966 when Harris, Hubby and Lewontin described the huge store of genetic variation detectable at the polypeptide level. Later Jeffreys used MVR (minisatellite variant repeat) analysis to describe the variation hidden within minisatellite VNTRs and to propose a mutational mechanism.^ The questions I have asked follow these traditions: (1) How much microvariant polymorphism exists at the discrete allele minisatellite D1S80 locus? (2) Do alleles or groups of alleles associate randomly with the flanking markers to form haplotypes? (3) What mechanisms might explain mutations at this locus? What are the phylogenetic relationships among the alleles?^ The minisatellite locus D1S80 (1p35-36), GenBank sequence (Accession # D28507), is a highly polymorphic Variable Number of Tandem Repeat (VNTR) based on a 16 base core. D1S80 alleles are electrophoretically separable into discontinuous sets of equivalent length alleles. Sequence variation or minor length variation within these classes was expected: I have sought to determine the nature of this microvariant heterogeneity by sequencing nominal and variant alleles.^ Alleles were analyzed by Single-Strand Conformation Polymorphism (SSCP) analysis. Sequences were determined to ascertain whether sequence variation or size variation is the major cause of altered electrophoretic migration of microvariant D1S80 alleles. Twenty three alleles from 14 previously typed individuals were sequenced. The individuals were from African American, Caucasian, or Hispanic databases.^ A Tsp509 I restriction site, previously reported as a Hinf I flanking polymorphism, and a 3$\sp\prime$ flanking region BsoF I restriction site polymorphism were identified. There appears to be a strong association of the 5$\sp\prime$ flanking region Hinf I(+) and Tsp509 I(-) site and the 3$\sp\prime$ flanking region BsoF I(-) site with the 18 allele, while the 24 tends to be associated with the Hinf I(-), Tsp509 I(+) and BsoF I(+) sites.^ The general conclusion for this locus is clearly the closer you look, the more you find. D1S80 allelic polymorphisms are primarily due to variation in the number of repeat units and to sequence variation among repeats. The sequenced based gene tree depicts two major classes of alleles which conform to the two most common alleles, reflecting either equivalent age or population size bottlenecks. ^

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes

Background - Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5–10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required. Methods - We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia. Results - We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study. Conclusions - A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.

Roles of CTCF and YY1 in T Cell Receptor Gene Rearrangement And T Cell Development

Diversity of T cell receptors (TCR) and immunoglobulins (Ig) is generated by V(D)J recombination of antigen receptor (AgR) loci. The Tcra-Tcrd locus is of particular interest because it displays a nested organization of Tcrd and Tcra gene segments and V(D)J recombination follows an intricate developmental program to assemble both TCRδ and TCRα repertoires. However, the mechanisms that dictate the developmental regulation of V(D)J recombination of the Tcra-Tcrd locus remain unclear.

We have previously shown that CCCTC-binding factor (CTCF) regulates Tcra gene transcription and rearrangement through organizing chromatin looping between CTCF- binding elements (CBEs). This study is one of many showing that CTCF functions as a chromatin organizer and transcriptional regulator genome-wide. However, detailed understanding of the impact of specific CBEs is needed to fully comprehend the biological function of CTCF and how CTCF influences the generation of the TCR repertoire during thymocyte development. Thus, we generated several mouse models with genetically modified CBEs to gain insight into the CTCF-dependent regulation of the Tcra-Tcrd locus. We revealed a CTCF-dependent chromatin interaction network at the Tcra-Tcrd locus in double-negative thymocytes. Disruption of a discrete chromatin loop encompassing Dδ, Jδ and Cδ gene segments allowed a single Vδ segment to frequently contact and rearrange to diversity and joining gene segments and dominate the adult TCRδ repertoire. Disruption of this loop also narrowed the TCRα repertoire, which, we believe, followed as a consequence of the restricted TCRδ repertoire. Hence, a single CTCF-mediated chromatin loop directly regulates TCRδ diversity and indirectly regulates TCRα diversity. In addition, we showed that insertion of an ectopic CBE can modify chromatin interactions and disrupt the rearrangement of particular Vδ gene segments. Finally, we investigated the role of YY1 in early T cell development by conditionally deleting YY1 in developing thymocytes. We found that early ablation of YY1 caused severe developmental defects in the DN compartment due to a dramatic increase in DN thymocyte apoptosis. Furthermore, late ablation of YY1 resulted in increased apoptosis of DP thymocytes and a restricted TCRα repertoire. Mechanistically, we showed that p53 was upregulated in both DN and DP YY1-deficient thymocytes. Eliminating p53 in YY1-deficient thymocytes rescued the survival and developmental defects, indicating that these YY1-dependent defects were p53-mediated. We conclude that YY1 is required to maintain cell viability during thymocyte development by thwarting the accumulation of p53.

Overall, this thesis work has shown that CTCF-dependent looping provides a central framework for lineage- and developmental stage-specific regulation of Tcra-Tcrd gene expression and rearrangements. In addition, we identified YY1 as a novel regulator of thymocyte viability.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk

Characterization of IGH locus breakpoints in multiple myeloma indicates a subset of translocations appear to occur in pregerminal center B cells.

Translocations in myeloma are thought to occur solely in mature B cells in the germinal center through class switch recombination (CSR). We used a targeted captured technique followed by massively parallel sequencing to determine the exact breakpoints in both the immunoglobulin heavy chain (IGH) locus and the partner chromosome in 61 presentation multiple myeloma samples. The majority of samples (62%) have a breakpoint within the switch regions upstream of the IGH constant genes and are generated through CSR in a mature B cell. However, the proportion of CSR translocations is not consistent between cytogenetic subgroups. We find that 100% of t(4;14) are CSR-mediated; however, 21% of t(11;14) and 25% of t(14;20) are generated through DH-JH recombination activation gene-mediated mechanisms, indicating they occur earlier in B-cell development at the pro-B-cell stage in the bone marrow. These 2 groups also generate translocations through receptor revision, as determined by the breakpoints and mutation status of the segments used in 10% and 50% of t(11;14) and t(14;20) samples, respectively. The study indicates that in a significant number of cases the translocation-based etiological events underlying myeloma may arise at the pro-B-cell hematological progenitor cell level, much earlier in B-cell development than was previously thought.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Characterization of Borrelia burgdorferi Gene Product BBD18 for Its Role(s) in Pathogen Biology and Infectivity

bbd18 is a differentially expressed Borrelia burgdorferi gene that is transcribed at almost undetectable levels in spirochetes grown in vitro but dramatically upregulated during tick infection. The gene also displays low yet detectable expression at various times in tissues of murine hosts. As the gene product bears no homology to known proteins, its biological significance remains enigmatic. To understand the gene function, we created isogenic bbd18-deletion mutants as well as genetically-complemented isolates from an infectious wild-type B. burgdorferi strain. Compared to parental isolates, bbd18 mutants - but not complemented spirochetes - displayed slower in vitro growth. The bbd18 mutants also reflect significantly reduced ability to persist or remain undetectable both in immunocompetent and SCID mice, yet were able to survive in ticks. This suggests BBD18 function is essential in mammalian hosts but redundant in the arthropod vector. Notably, although bbd18 expression and in vitro growth defects are restored in the complemented isolates, their phenotype is similar to the mutants - being unable to persist in mice but able to survive in ticks. Despite low expression in cultured wild-type B. burgdorferi, bbd18 deletion downregulated several genes. Interestingly, expression of some, including ospD and bbi39, could be complemented, while that of others could not be restored via bbd18 re-expression. Correspondingly, bbd18 mutants displayed altered production of several proteins, and similar to RNA levels, some were restored in the bbd18 complement and others not. To understand how bbd18 deletion results in apparently permanent and noncomplementable phenotypic defects, we sought to genetically disturb the DNA topology surrounding the bbd18 locus without deleting the gene. Spirochetes with an antibiotic cassette inserted downstream of the gene, between bbd17 and bbd18, were significantly attenuated in mice, while a similar upstream insertion, between bbd18 and bbd19, did not affect infectivity, suggesting that an unidentified cis element downstream of bbd18 may encode a virulence-associated factor critical for infection.

Cloning of the ω-secalin gene family in a wheat 1BL/1RS translocation line using BAC clone sequencing

Background: Wheat 1BL/1RS translocation lines are planted around the world for their disease resistance and high yield. Most of them are poor in bread making, which is partially caused by ω-secalins that are encoded by the ω-secalin gene family, which is located on the short arm of rye chromosome 1R (1RS). However, information on the structure and evolution of the ω-secalin gene family is still limited. Results: We first generated a physicalmap of the ω-secalin gene family covering 195 kb of the Sec-1 locus based on sequencing three bacterial artificial chromosome (BAC) clones of the 1BL/1RS translocation wheat cultivar Shimai 15. A BAC contig was constructed spanning 168 kb of the Sec-1 locus on 1RS. Twelve ω-secalin genes were arranged in a head-to-tail fashion, separated by 8.2–21.6 kb spacers on the contig, whereas six other ω-secalin genes were arranged head-to-tail, separated by 8.2–8.4 kb of spacers on clone BAC125. The 18 ω-secalin genes can be classified into six types among which eight ω-secalin genes were expressed during seed development. The ω-secalin genes with the 1074-bp open reading frame (ORF) represented the main population. Except for two pseudogenes, the N-terminal of the ω-secalin gene was conserved, whereas variations in the C-terminal led to a change in ORF length. The spacers can be sorted into two classes. Class-1 spacers contained conserved and non-conservative sequences. Conclusion: The ω-secalin gene family consisted of at least 18 members in the 1BL/1RS translocation line cv. Shimai 15. Eight ω-secalin genes were expressed during seed development. Eighteen members may originate from a progenitor with a 1,074-bp ORF. The spacers differed in length and sequence conservation.