IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.

Identification of genetic risk factors for cerebellar mutism in pediatric brain tumor patients

Following posterior fossa surgery for resection of childhood medulloblastoma and primitive neuroectodermal tumor (M/PNET), cerebellar mutism (CM) may develop. This is a condition of absent or diminished speech in a conscious patient with no evidence of oral apraxia, which can be accompanied by other symptoms of the posterior fossa syndrome complex, which includes ataxia and hypotonia. Little is known about the etiology. Therefore, we conducted a SNP, gene, and pathway-level analysis to assess the role of host genetic variation on the risk of CM in M/PNET subjects following treatment. Cases (n= 20) and controls (n= 53) were recruited from the Childhood Cancer Epidemiology and Prevention Center, in Houston, TX. DNA samples were genotyped using the Illumina Human 1M Quad SNP chip. Ten pathways were identified from logistic regression used to identify the marginal effect of each SNP on CM risk. The minP test was conducted to identify associations between SNPs categorized to genes and CM risk. Pathways were assessed to determine if there was a significant enrichment of genes in the pathway compared to all other pathways. There were 78 genes that reached the threshold of min P ≤0.05 in 948 genes. The Neurotoxicity pathway was the most significant pathway after adjusting for multiple comparisons (q=0.040 and q=0.005, using Fisher's exact test and a test of proportions, respectively). Most genes within the Neurotoxicity pathway that reached a threshold of minP ≤0.05 were known to have an apoptosis function, possibly inducing neuronal apoptosis in the dentatothalamocortical pathway, and may be important in CM etiology in this population. This is the first study to assess the potential role of genetic risk factors on CM. As an exploratory study, these results should be replicated in a larger sample. ^

Evidence of cryptic and pseudocryptic speciation in the Paracalanus parvus species complex (Crustacea, Copepoda, Calanoida), supplementary material

Introduction Many marine planktonic crustaceans such as copepods have been considered as widespread organisms. However, the growing evidence for cryptic and pseudo-cryptic speciation has emphasized the need of re-evaluating the status of copepod species complexes in molecular and morphological studies to get a clearer picture about pelagic marine species as evolutionary units and their distributions. This study analyses the molecular diversity of the ecologically important Paracalanus parvus species complex. Its seven currently recognized species are abundant and also often dominant in marine coastal regions worldwide from temperate to tropical oceans. Results COI and Cytochrome b sequences of 160 specimens of the Paracalanus parvus complex from all oceans were obtained. Furthermore, 42 COI sequences from GenBank were added for the genetic analyses. Thirteen distinct molecular operational taxonomic units (MOTU) and two single sequences were revealed with cladistic analyses (Maximum Likelihood, Bayesian Inference), of which seven were identical with results from species delimitation methods (barcode gaps, ABDG, GMYC, Rosenberg's P(AB)). In total, 10 to 12 putative species were detected and could be placed in three categories: (1) temperate geographically isolated, (2) warm-temperate to tropical wider spread and (3) circumglobal warm-water species. Conclusions The present study provides evidence of cryptic or pseudocryptic speciation in the Paracalanus parvus complex. One major insight is that the species Paracalanus parvus s.s. is not panmictic, but may be restricted in its distribution to the northeastern Atlantic.

Human midbrain precursors activate the expected developmental genetic program and differentiate long-term to functional A9 dopamine neurons in vitro. Enhancement by Bcl-XL

Understanding the molecular programs of the generation of human dopaminergic neurons (DAn) from their ventral mesencephalic (VM) precursors is of key importance for basic studies, progress in cell therapy, drug screening and pharmacology in the context of Parkinson's disease. The nature of human DAn precursors in vitro is poorly understood, their properties unstable, and their availability highly limited. Here we present positive evidence that human VM precursors retaining their genuine properties and long-term capacity to generate A9 type Substantia nigra human DAn (hVM1 model cell line) can be propagated in culture. During a one month differentiation, these cells activate all key genes needed to progress from pro-neural and prodopaminergic precursors to mature and functional DAn. For the first time, we demonstrate that gene cascades are correctly activated during differentiation, resulting in the generation of mature DAn. These DAn have morphological and functional properties undistinguishable from those generated by VM primary neuronal cultures. In addition, we have found that the forced expression of Bcl-XL induces an increase in the expression of key developmental genes (MSX1, NGN2), maintenance of PITX3 expression temporal profile, and also enhances genes involved in DAn long-term function, maintenance and survival (EN1, LMX1B, NURR1 and PITX3). As a result, Bcl-XL anticipates and enhances DAn generation.

Genetic characterization of chestnut (Castanea sativa Mill.) orchards and traditional nut varieties in El Bierzo, a glacial refuge and major cultivation site in Northwestern Spain

This paper presents a detailed genetic study of Castanea sativa in El Bierzo, a major nut production region with interesting features. It is located within a glacial refuge at one extreme of the distribution area (northwest Spain); it has a centenary tradition of chestnut management; and more importantly, it shows an unusual degree of genetic isolation. Seven nuclear microsatellite markers were selected to analyze the genetic variability and structure of 169 local trees grafted for nut production. We analyzed in the same manner 62 local nuts. The selected loci were highly discriminant for the genotypes studied, giving a combined probability of identity of 6.1 × 10−6. An unprecedented density of trees was sampled for this project over the entire region, and nuts were collected representing 18 cultivars marketed by local producers. Several instances of misclassification by local growers were detected. Fixation index estimates and analysis of molecular variance (AMOVA) data are supportive of an unexpectedly high level of genetic differentiation in El Bierzo, larger than that estimated in a previous study with broader geographical scope but based on limited local sampling (Pereira-Lorenzo et al., Tree Genet Genomes 6: 701–715, 2010a). Likewise, we have determined that clonality due to grafting had been previously overestimated. In line with these observations, no significant spatial structure was found using both a model-based Bayesian procedure and Mantel’s tests. Taken together, our results evidence the need for more fine-scale genetic studies if conservation strategies are to be efficiently improved.

Ecological and genetic determinants of Pepino mosaic virus emergence

Virus emergence is a complex phenomenon, which generally involves spread to a new host from a wild host, followed by adaptation to the new host. Although viruses account for the largest fraction of emerging crop pathogens, knowledge about their emergence is incomplete. We address here the question of whether Pepino mosaic virus (PepMV) emergence as a major tomato pathogen worldwide could have involved spread from wild to cultivated plant species and host adaptation. For this, we surveyed natural populations of wild tomatoes in southern Peru for PepMV infection. PepMV incidence, genetic variation, population structure, and accumulation in various hosts were analyzed. PepMV incidence in wild tomatoes was high, and a strain not yet reported in domestic tomato was characterized. This strain had a wide host range within the Solanaceae, multiplying efficiently in most assayed Solanum species and being adapted to wild tomato hosts. Conversely, PepMV isolates from tomato crops showed evidence of adaptation to domestic tomato, possibly traded against adaptation to wild tomatoes. Phylogenetic reconstructions indicated that the most probable ancestral sequence came from a wild Solanum species. A high incidence of PepMV in wild tomato relatives would favor virus spread to crops and its efficient multiplication in different Solanum species, including tomato, allowing its establishment as an epidemic pathogen. Later, adaptation to tomato, traded off against adaptation to other Solanum species, would isolate tomato populations from those in other hosts.

Identification of novel susceptibility loci for inflammatory bowel disease on chromosomes 1p, 3q, and 4q: Evidence for epistasis between 1p and IBD1

The idiopathic inflammatory bowel diseases, Crohn’s disease (CD) and ulcerative colitis (UC), are chronic, frequently disabling diseases of the intestines. Segregation analyses, twin concordance, and ethnic differences in familial risks have established that CD and UC are complex, non-Mendelian, related genetic disorders. We performed a genome-wide screen using 377 autosomal markers, on 297 CD, UC, or mixed relative pairs from 174 families, 37% Ashkenazim. We observed evidence for linkage at 3q for all families (multipoint logarithm of the odds score (MLod) = 2.29, P = 5.7 × 10−4), with greatest significance for non-Ashkenazim Caucasians (MLod = 3.39, P = 3.92 × 10−5), and at chromosome 1p (MLod = 2.65, P = 2.4 × 10−4) for all families. In a limited subset of mixed families (containing one member with CD and another with UC), evidence for linkage was observed on chromosome 4q (MLod = 2.76, P = 1.9 × 10−4), especially among Ashkenazim. There was confirmatory evidence for a CD locus, overlapping IBD1, in the pericentromeric region of chromosome 16 (MLod = 1.69, P = 2.6 × 10−3), particularly among Ashkenazim (MLod = 1.51, P = 7.8 × 10−3); however, positive MLod scores were observed over a very broad region of chromosome 16. Furthermore, evidence for epistasis between IBD1 and chromosome 1p was observed. Thirteen additional loci demonstrated nominal (MLod > 1.0, P < 0.016) evidence for linkage. This screen provides strong evidence that there are several major susceptibility loci contributing to the genetic risk for CD and UC.

Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid

The genetic basis of heterosis was investigated in an elite rice hybrid by using a molecular linkage map with 150 segregating loci covering the entire rice genome. Data for yield and three traits that were components of yield were collected over 2 years from replicated field trials of 250 F2:3 families. Genotypic variations explained from about 50% to more than 80% of the total variation. Interactions between genotypes and years were small compared with the main effects. A total of 32 quantitative trait loci (QTLs) were detected for the four traits; 12 were observed in both years and the remaining 20 were detected in only one year. Overdominance was observed for most of the QTLs for yield and also for a few QTLs for the component traits. Correlations between marker heterozygosity and trait expression were low, indicating that the overall heterozygosity made little contribution to heterosis. Digenic interactions, including additive by additive, additive by dominance, and dominance by dominance, were frequent and widespread in this population. The interactions involved large numbers of marker loci, most of which individually were not detectable on single-locus basis; many interactions among loci were detected in both years. The results provide strong evidence that epistasis plays a major role as the genetic basis of heterosis.

Sulfate reduction in higher plants: Molecular evidence for a novel 5′-adenylylsulfate reductase

Sulfate-assimilating organisms reduce inorganic sulfate for Cys biosynthesis. There are two leading hypotheses for the mechanism of sulfate reduction in higher plants. In one, adenosine 5′-phosphosulfate (APS) (5′-adenylylsulfate) sulfotransferase carries out reductive transfer of sulfate from APS to reduced glutathione. Alternatively, the mechanism may be similar to that in bacteria in which the enzyme, 3′-phosphoadenosine-5′-phosphosulfate (PAPS) reductase, catalyzes thioredoxin (Trx)-dependent reduction of PAPS. Three classes of cDNA were cloned from Arabidopsis thaliana termed APR1, -2, and -3, that functionally complement a cysH, PAPS reductase mutant strain of Escherichia coli. The coding sequence of the APR clones is homologous with PAPS reductases from microorganisms. In addition, a carboxyl-terminal domain is homologous with members of the Trx superfamily. Further genetic analysis showed that the APR clones can functionally complement a mutant strain of E. coli lacking Trx, and an APS kinase, cysC. mutant. These results suggest that the APR enzyme may be a Trx-independent APS reductase. Cell extracts of E. coli expressing APR showed Trx-independent sulfonucleotide reductase activity with a preference for APS over PAPS as a substrate. APR-mediated APS reduction is dependent on dithiothreitol, has a pH optimum of 8.5, is stimulated by high ionic strength, and is sensitive to inactivation by 5′-adenosinemonophosphate (5′-AMP). 2′-AMP, or 3′-phosphoadenosine-5′-phosphate (PAP), a competitive inhibitor of PAPS reductase, do not affect activity. The APR enzymes may be localized in different cellular compartments as evidenced by the presence of an amino-terminal transit peptide for plastid localization in APR1 and APR3 but not APR2. Southern blot analysis confirmed that the APR clones are members of a small gene family, possibly consisting of three members.

Evidence that Myb-related CDC5 proteins are required for pre-mRNA splicing

The conserved CDC5 family of Myb-related proteins performs an essential function in cell cycle control at G2/M. Although c-Myb and many Myb-related proteins act as transcription factors, herein, we implicate CDC5 proteins in pre-mRNA splicing. Mammalian CDC5 colocalizes with pre-mRNA splicing factors in the nuclei of mammalian cells, associates with core components of the splicing machinery in nuclear extracts, and interacts with the spliceosome throughout the splicing reaction in vitro. Furthermore, genetic depletion of the homolog of CDC5 in Saccharomyces cerevisiae, CEF1, blocks the first step of pre-mRNA processing in vivo. These data provide evidence that eukaryotic cells require CDC5 proteins for pre-mRNA splicing.

Genome scan of human systemic lupus erythematosus: Evidence for linkage on chromosome 1q in African-American pedigrees

Systemic lupus erythematosus (SLE) is an autoimmune disorder characterized by production of autoantibodies against intracellular antigens including DNA, ribosomal P, Ro (SS-A), La (SS-B), and the spliceosome. Etiology is suspected to involve genetic and environmental factors. Evidence of genetic involvement includes: associations with HLA-DR3, HLA-DR2, Fcγ receptors (FcγR) IIA and IIIA, and hereditary complement component deficiencies, as well as familial aggregation, monozygotic twin concordance >20%, λs > 10, purported linkage at 1q41–42, and inbred mouse strains that consistently develop lupus. We have completed a genome scan in 94 extended multiplex pedigrees by using model-based linkage analysis. Potential [log10 of the odds for linkage (lod) > 2.0] SLE loci have been identified at chromosomes 1q41, 1q23, and 11q14–23 in African-Americans; 14q11, 4p15, 11q25, 2q32, 19q13, 6q26–27, and 12p12–11 in European-Americans; and 1q23, 13q32, 20q13, and 1q31 in all pedigrees combined. An effect for the FcγRIIA candidate polymorphism) at 1q23 (lod = 3.37 in African-Americans) is syntenic with linkage in a murine model of lupus. Sib-pair and multipoint nonparametric analyses also support linkage (P < 0.05) at nine loci detected by using two-point lod score analysis (lod > 2.0). Our results are consistent with the presumed complexity of genetic susceptibility to SLE and illustrate racial origin is likely to influence the specific nature of these genetic effects.

Genome scan for teratogen-induced clefting susceptibility loci in the mouse: Evidence of both allelic and locus heterogeneity distinguishing cleft lip and cleft palate

Nonsyndromic clefting of the lip and palate in humans has a highly complex etiology, with both multiple genetic loci and exposure to teratogens influencing susceptibility. Previous studies using mouse models have examined only very small portions of the genome. Here we report the findings of a genome-wide search for susceptibility genes for teratogen-induced clefting in the AXB and BXA set of recombinant inbred mouse strains. We compare results obtained using phenytoin (which induces cleft lip) and 6-aminonicotinamide (which induces cleft palate). We use a new statistical approach based on logistic regression suitable for these categorical data to identify several chromosomal regions as possible locations of clefting susceptibility loci, and we review candidate genes located within each region. Because cleft lip and cleft palate do not frequently co-aggregate in human families and because these structures arise semi-independently during development, these disorders are usually considered to be distinct in etiology. Our data, however, implicate several of the same chromosomal regions for both forms of clefting when teratogen-induced. Furthermore, different parental strain alleles are usually associated with clefting of the lip versus that of the palate (i.e., allelic heterogeneity). Because several other chromosomal regions are associated with only one form of clefting, locus heterogeneity also appears to be involved. Our findings in this mouse model suggest several priority areas for evaluation in human epidemiological studies.

A Genetic Analysis of Interactions with Spc110p Reveals Distinct Functions of Spc97p and Spc98p, Components of the Yeast γ-Tubulin Complex

The spindle pole body (SPB) in Saccharomyces cerevisiae functions as the microtubule-organizing center. Spc110p is an essential structural component of the SPB and spans between the central and inner plaques of this multilamellar organelle. The amino terminus of Spc110p faces the inner plaque, the substructure from which spindle microtubules radiate. We have undertaken a synthetic lethal screen to identify mutations that enhance the phenotype of the temperature-sensitive spc110–221 allele, which encodes mutations in the amino terminus. The screen identified mutations in SPC97 and SPC98, two genes encoding components of the Tub4p complex in yeast. The spc98–63 allele is synthetic lethal only with spc110 alleles that encode mutations in the N terminus of Spc110p. In contrast, the spc97 alleles are synthetic lethal with spc110 alleles that encode mutations in either the N terminus or the C terminus. Using the two-hybrid assay, we show that the interactions of Spc110p with Spc97p and Spc98p are not equivalent. The N terminus of Spc110p displays a robust interaction with Spc98p in two different two-hybrid assays, while the interaction between Spc97p and Spc110p is not detectable in one strain and gives a weak signal in the other. Extra copies of SPC98 enhance the interaction between Spc97p and Spc110p, while extra copies of SPC97 interfere with the interaction between Spc98p and Spc110p. By testing the interactions between mutant proteins, we show that the lethal phenotype in spc98–63 spc110–221 cells is caused by the failure of Spc98–63p to interact with Spc110–221p. In contrast, the lethal phenotype in spc97–62 spc110–221 cells can be attributed to a decreased interaction between Spc97–62p and Spc98p. Together, these studies provide evidence that Spc110p directly links the Tub4p complex to the SPB. Moreover, an interaction between Spc98p and the amino-terminal region of Spc110p is a critical component of the linkage, whereas the interaction between Spc97p and Spc110p is dependent on Spc98p.

Chromosomal breakage-fusion-bridge events cause genetic intratumor heterogeneity

It has long been known that rearrangements of chromosomes through breakage-fusion-bridge (BFB) cycles may cause variability of phenotypic and genetic traits within a cell population. Because intercellular heterogeneity is often found in neoplastic tissues, we investigated the occurrence of BFB events in human solid tumors. Evidence of frequent BFB events was found in malignancies that showed unspecific chromosome aberrations, including ring chromosomes, dicentric chromosomes, and telomeric associations, as well as extensive intratumor heterogeneity in the pattern of structural changes but not in tumors with tumor-specific aberrations and low variability. Fluorescence in situ hybridization analysis demonstrated that chromosomes participating in anaphase bridge formation were involved in a significantly higher number of structural aberrations than other chromosomes. Tumors with BFB events showed a decreased elimination rate of unstable chromosome aberrations after irradiation compared with normal cells and other tumor cells. This result suggests that a combination of mitotically unstable chromosomes and an elevated tolerance to chromosomal damage leads to constant genomic reorganization in many malignancies, thereby providing a flexible genetic system for clonal evolution and progression.

Altered recognition mutants of the response regulator PhoB: A new genetic strategy for studying protein–protein interactions

Two-component regulatory systems require highly specific interactions between histidine kinase (transmitter) and response regulator (receiver) proteins. We have developed a novel genetic strategy that is based on tightly regulated synthesis of a given protein to identify domains and residues of an interacting protein that are critical for interactions between them. Using a reporter strain synthesizing the nonpartner kinase VanS under tight arabinose control and carrying a promoter-lacZ fusion activated by phospho-PhoB, we isolated altered recognition (AR) mutants of PhoB showing enhanced activation (phosphorylation) by VanS as arabinose-dependent Lac+ mutants. Changes in the PhoBAR mutants cluster in a “patch” near the proposed helix 4 of PhoB based on the CheY crystal structure (a homolog of the PhoB receiver domain) providing further evidence that helix 4 lies in the kinase-regulator interface. Based on the CheY structure, one mutant has an additional change in a region that may propagate a conformational change to helix 4. The overall genetic strategy described here may also be useful for studying interactions of other components of the vancomycin resistance and Pi signal transduction pathways, other two-component regulatory systems, and other interacting proteins. Conditionally replicative oriRR6Kγ attP “genome targeting” suicide plasmids carrying mutagenized phoB coding regions were integrated into the chromosome of a reporter strain to create mutant libraries; plasmids encoding mutant PhoB proteins were subsequently retrieved by P1-Int-Xis cloning. Finally, the use of similar genome targeting plasmids and P1-Int-Xis cloning should be generally useful for constructing genomic libraries from a wide array of organisms.