Optimizing the creation of base populations for aquaculture breeding programs using phenotypic and genomic data and its consequences on genetic progress

The success of an aquaculture breeding program critically depends on the way in which the base population of breeders is constructed since all the genetic variability for the traits included originally in the breeding goal as well as those to be included in the future is contained in the initial founders. Traditionally, base populations were created from a number of wild strains by sampling equal numbers from each strain. However, for some aquaculture species improved strains are already available and, therefore, mean phenotypic values for economically important traits can be used as a criterion to optimize the sampling when creating base populations. Also, the increasing availability of genome-wide genotype information in aquaculture species could help to refine the estimation of relationships within and between candidate strains and, thus, to optimize the percentage of individuals to be sampled from each strain. This study explores the advantages of using phenotypic and genome-wide information when constructing base populations for aquaculture breeding programs in terms of initial and subsequent trait performance and genetic diversity level. Results show that a compromise solution between diversity and performance can be found when creating base populations. Up to 6% higher levels of phenotypic performance can be achieved at the same level of global diversity in the base population by optimizing the selection of breeders instead of sampling equal numbers from each strain. The higher performance observed in the base population persisted during 10 generations of phenotypic selection applied in the subsequent breeding program.

Direct Vlasov simulations of electron-attracting cylindrical Langmuir probes in flowing plasmas

Current collection by positively polarized cylindrical Langmuir probes immersed in flowing plasmas is analyzed using a non-stationary direct Vlasov-Poisson code. A detailed description of plasma density spatial structure as a function of the probe-to-plasma relative velocity U is presented. Within the considered parametric domain, the well-known electron density maximum close to the probe is weakly affected by U. However, in the probe wake side, the electron density minimum becomes deeper as U increases and a rarified plasma region appears. Sheath radius is larger at the wake than at the front side. Electron and ion distribution functions show specific features that are the signature of probe motion. In particular, the ion distribution function at the probe front side exhibits a filament with positive radial velocity. It corresponds to a population of rammed ions that were reflected by the electric field close to the positively biased probe. Numerical simulations reveal that two populations of trapped electrons exist: one orbiting around the probe and the other with trajectories confined at the probe front side. The latter helps to neutralize the reflected ions, thus explaining a paradox in past probe theory.

Transient excess of MYC activity can elicit genomic instability and tumorigenesis

Overexpression of the MYC protooncogene has been implicated in the genesis of diverse human tumors. Tumorigenesis induced by MYC has been attributed to sustained effects on proliferation and differentiation. Here we report that MYC may also contribute to tumorigenesis by destabilizing the cellular genome. A transient excess of MYC activity increased tumorigenicity of Rat1A cells by at least 50-fold. The increase persisted for >30 days after the return of MYC activity to normal levels. The brief surfeit of MYC activity was accompanied by evidence of genomic instability, including karyotypic abnormalities, gene amplification, and hypersensitivity to DNA-damaging agents. MYC also induced genomic destabilization in normal human fibroblasts, although these cells did not become tumorigenic. Stimulation of Rat1A cells with MYC accelerated their passage through G1/S. Moreover, MYC could force normal human fibroblasts to transit G1 and S after treatment with N-(phosphonoacetyl)-l-aspartate (PALA) at concentrations that normally lead to arrest in S phase by checkpoint mechanisms. Instead, the cells subsequently appeared to arrest in G2. We suggest that the accelerated passage through G1 was mutagenic but that the effect of MYC permitted a checkpoint response only after G2 had been reached. Thus, MYC may contribute to tumorigenesis through a dominant mutator effect.

Extent of genomic rearrangement after genome duplication in yeast

Whole-genome duplication approximately 108 years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70–100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150–200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Discrimination of the closely related A and D genomes of the hexaploid oat Avena sativa L.

A satellite DNA sequence, As120a, specific to the A-genome chromosomes in the hexaploid oat, Avena sativa L., was isolated by subcloning a fragment with internal tandem repeats from a plasmid, pAs120, that had been obtained from an Avena strigosa (As genome) genomic library. Southern and in situ hybridization showed that sequences with homology to sequences within pAs120 were dispersed throughout the genome of diploid (A and C genomes), tetraploid (AC genomes), and hexaploid (ACD genomes) Avena species. In contrast, sequences homologous to As120a were found in two A-genome species (A. strigosa and Avena longiglumis) and in the hexaploid A. sativa whereas this sequence was little amplified in the tetraploid Avena murphyi and was absent in the remaining A- and C-genome diploid species. In situ hybridization of pAs120a to hexaploid oat species revealed the distribution of elements of the As120a repeated family over both arms of 14 of 42 chromosomes of this species. By using double in situ hybridization with pAs120a and a C genome-specific probe, three sets of 14 chromosomes were revealed corresponding to the A, C, and D genomes of the hexaploid species. Simultaneous in situ hybridizations with pAs120a and ribosomal probes were used to assign the SAT chromosomes of hexaploid species to their correct genomes. This work reports a sequence able to distinguish between the closely related A and D genomes of hexaploid oats. This sequence offers new opportunities to analyze the relationships of Avena species and to explore the possible evolution of various polyploid oat species.

Human ATP-binding cassette transporter 1 (ABC1): Genomic organization and identification of the genetic defect in the original Tangier disease kindred

Tangier disease is characterized by low serum high density lipoproteins and a biochemical defect in the cellular efflux of lipids to high density lipoproteins. ABC1, a member of the ATP-binding cassette family, recently has been identified as the defective gene in Tangier disease. We report here the organization of the human ABC1 gene and the identification of a mutation in the ABC1 gene from the original Tangier disease kindred. The organization of the human ABC1 gene is similar to that of the mouse ABC1 gene and other related ABC genes. The ABC1 gene contains 49 exons that range in size from 33 to 249 bp and is over 70 kb in length. Sequence analysis of the ABC1 gene revealed that the proband for Tangier disease was homozygous for a deletion of nucleotides 3283 and 3284 (TC) in exon 22. The deletion results in a frameshift mutation and a premature stop codon starting at nucleotide 3375. The product is predicted to encode a nonfunctional protein of 1,084 aa, which is approximately half the size of the full-length ABC1 protein. The loss of a Mnl1 restriction site, which results from the deletion, was used to establish the genotype of the rest of the kindred. In summary, we report on the genomic organization of the human ABC1 gene and identify a frameshift mutation in the ABC1 gene of the index case of Tangier disease. These results will be useful in the future characterization of the structure and function of the ABC1 gene and the analysis of additional ABC1 mutations in patients with Tangier disease.

Probes bound to myosin Cys-707 rotate during length transients in contraction

It is widely conjectured that muscle shortens because portions of myosin molecules (the “cross-bridges”) impel the actin filament to which they transiently attach and that the impulses result from rotation of the cross-bridges. Crystallography indicates that a cross-bridge is articulated–consisting of a globular catalytic/actin-binding domain and a long lever arm that may rotate. Conveniently, a rhodamine probe with detectable attitude can be attached between the globular domain and the lever arm, enabling the observer to tell whether the anchoring region rotates. Well-established signature effects observed in shortening are tension changes resulting from the sudden release or quick stretch of active muscle fibers. In this investigation we found that closely correlated with such tension changes are changes in the attitude of the rhodamine probes. This correlation strongly supports the conjecture about how shortening is achieved.

Panhandle PCR strategy to amplify MLL genomic breakpoints in treatment-related leukemias

Panhandle PCR amplifies genomic DNA with known 5′ and unknown 3′ sequences from a template with an intrastrand loop schematically shaped like a pan with a handle. We used panhandle PCR to clone MLL genomic breakpoints in two pediatric treatment-related leukemias. The karyotype in a case of treatment-related acute lymphoblastic leukemia showed the t(4;11)(q21;q23). Panhandle PCR amplified the translocation breakpoint at position 2158 in intron 6 in the 5′ MLL breakpoint cluster region (bcr). The karyotype in a case of treatment-related acute myeloid leukemia was normal, but Southern blot analysis showed a single MLL gene rearrangement. Panhandle PCR amplified the breakpoint at position 1493 in MLL intron 6. Screening of somatic cell hybrid and radiation hybrid DNAs by PCR and reverse transcriptase-PCR analysis of the leukemic cells indicated that panhandle PCR identified a fusion of MLL intron 6 with a previously uncharacterized sequence in MLL intron 1, consistent with a partial duplication. In both cases, the breakpoints in the MLL bcr were in Alu repeats, and there were Alu repeats in proximity to the breakpoints in the partner DNAs, suggesting that Alu sequences were relevant to these rearrangements. This study shows that panhandle PCR is an effective method for cloning MLL genomic breakpoints in treatment-related leukemias. Analysis of additional pediatric cases will determine whether breakpoint distribution deviates from the predilection for 3′ distribution in the bcr that has been found in adult cases.

Fluorescence correlation analysis of probe diffusion simplifies quantitative pathogen detection by PCR

A sensitive, labor-saving, and easily automatable nonradioactive procedure named APEX-FCS (amplified probe extension detected by fluorescence correlation spectroscopy) has been established to detect specific in vitro amplification of pathogen genomic sequences. As an example, Mycobacterium tuberculosis genomic DNA was subjected to PCR amplification with the Stoffel fragment of Thermus aquaticus DNA polymerase in the presence of nanomolar concentrations of a rhodamine-labeled probe (third primer), binding to the target in between the micromolar amplification primers. The probe becomes extended only when specific amplification occurs. Its low concentration avoids false-positives due to unspecific hybridization under PCR conditions. With increasing portion of extended probe molecules, the probe’s average translational diffusion properties gradually change over the course of the reaction, reflecting amplification kinetics. Following PCR, this change from a stage of high to a stage of low mobility can directly be monitored during a 30-s measurement using a fluorescence correlation spectroscopy device. Quantitation down to 10 target molecules in a background of 2.5 μg unspecific DNA without post-PCR probe manipulations could be achieved with different primer/probe combinations. The assay holds the promise to concurrently perform amplification, probe hybridization, and specific detection without opening the reaction chamber, if sealable foils are used.

Genomic analysis of human and mouse TCL1 loci reveals a complex of tightly clustered genes

TCL1 and TCL1b genes on human chromosome 14q23.1 are activated in T cell leukemias by translocations and inversions at 14q32.1, juxtaposing them to regulatory elements of T cell receptor genes. In this report we present the cloning, mapping, and expression analysis of the human and murine TCL1/Tcl1 locus. In addition to TCL1 and TCL1b, the human locus contains two additional genes, TCL1-neighboring genes (TNG) 1 and 2, encoding proteins of 141 and 110 aa, respectively. Both genes show no homology to any known genes, but their expression profiles are very similar to those of TCL1 and TCL1b. TNG1 and TNG2 also are activated in T cell leukemias with rearrangements at 14q32.1. To aid in the development of a mouse model we also have characterized the murine Tcl1 locus and found five genes homologous to human TCL1b. Tcl1b1–Tcl1b5 proteins range from 117 to 123 aa and are 65–80% similar, but they show only a 30–40% similarity to human TCL1b. All five mouse Tcl1b and murine Tcl1 mRNAs are abundant in mouse oocytes and two-cell embryos but rare in various adult tissues and lymphoid cell lines. These data suggest a similar or complementary function of these proteins in early embryogenesis.

A model system to study genomic imprinting of human genes

Somatic-cell hybrids have been shown to maintain the correct epigenetic chromatin states to study developmental globin gene expression as well as gene expression on the active and inactive X chromosomes. This suggests the potential use of somatic-cell hybrids containing either a maternal or a paternal human chromosome as a model system to study known imprinted genes and to identify as-yet-unknown imprinted genes. Testing gene expression by using reverse transcription followed by PCR, we show that functional imprints are maintained at four previously characterized 15q11–q13 loci in hybrids containing a single human chromosome 15 and at two chromosome 11p15 loci in hybrids containing a single chromosome 11. In contrast, three γ-aminobutyric acid type A receptor subunit genes in 15q12–q13 are nonimprinted. Furthermore, we have found that differential DNA methylation imprints at the SNRPN promoter and at a CpG island in 11p15 are also maintained in somatic-cell hybrids. Somatic-cell hybrids therefore are a valid and powerful system for studying known imprinted genes as well as for rapidly identifying new imprinted genes.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5′ flanking sequence (−855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5′ flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.

The onset and extent of genomic instability in sporadic colorectal tumor progression

Cancer cell genomes contain alterations beyond known etiologic events, but their total number has been unknown at even the order of magnitude level. By sampling colorectal premalignant polyp and carcinoma cell genomes through use of the technique inter-(simple sequence repeat) PCR, we have found genomic alterations to be considerably more abundant than expected, with the mean number of genomic events per carcinoma cell totaling approximately 11,000. Colonic polyps early in the tumor progression pathway showed similar numbers of events. These results indicate that, as with certain hereditary cancer syndromes, genomic destabilization is an early step in sporadic tumor development. Together these results support the model of genomic instability being a cause rather than an effect of malignancy, facilitating vastly accelerated somatic cell evolution, with the observed orderly steps of the colon cancer progression pathway reflecting the consequences of natural selection.

Cloning of chlorophyllase, the key enzyme in chlorophyll degradation: Finding of a lipase motif and the induction by methyl jasmonate

Chlorophyllase (Chlase) is the first enzyme involved in chlorophyll (Chl) degradation and catalyzes the hydrolysis of ester bond to yield chlorophyllide and phytol. In the present study, we isolated the Chlase cDNA. We synthesized degenerate oligo DNA probes based on the internal amino acid sequences of purified Chlase from Chenopodium album, screened the C. album cDNA library, and cloned a cDNA (CaCLH, C. album chlorophyll-chlorophyllido hydrolase). The deduced amino acid sequence (347 aa residues) had a lipase motif overlapping with an ATP/GTP-binding motif (P-loop). CaCLH possibly was localized in the extraplastidic part of the cell, because a putative signal sequence for endoplasmic reticulum is at the N terminus. The amino acid sequence shared 37% identity with a function-unknown gene whose mRNA is inducible by coronatine and methyl jasmonate (MeJA) in Arabidopsis thaliana (AtCLH1). We expressed the gene products of AtCLH1 and of CaCLH in Escherichia coli, and they similarly exhibited Chlase activity. Moreover, we isolated another full-length cDNA based on an Arabidopsis genomic fragment and expressed it in E. coli, demonstrating the presence of the second Arabidopsis CLH gene (AtCLH2). No typical feature of signal sequence was identified in AtCLH1, whereas AtCLH2 had a typical signal sequence for chloroplast. AtCLH1 mRNA was induced rapidly by a treatment of MeJA, which is known to promote senescence and Chl degradation in plants, and a high mRNA level was maintained up to 9 h. AtCLH2, however, did not respond to MeJA.