Integrated cytogenetic BAC map of the genome of the gray, short-tailed opossum, Monodelphis domestica.

The generation of high-quality genome assemblies for numerous species is advancing at a rapid pace. As the number of genome assemblies increases, so does our ability to investigate genome relationships and their contributions to unraveling complex biological, evolutionary, and biomedical processes. A key process in the generation of a genome assembly is to determine and verify the precise physical location and order of the large sequence blocks (scaffolds) that result from the assembly. For organisms of relatively recent common ancestry this process may be achieved largely through comparative sequence alignment. However, as the evolutionary distance between species lengthens, the use of comparative sequence alignment becomes increasingly less reliable. Simultaneous cytogenetic mapping, using multicolor fluorescence in-situ hybridization (FISH) analysis, offers an alternative means to define the cytogenetic location and relative order of DNA sequences, thereby anchoring the genome sequence to the karyotype. In this article we report the molecular cytogenetic locations of 415 bacterial artificial chromosome (BAC) clones that served to anchor sequence scaffolds of the gray, short-tailed opossum (Monodelphis domestica) to its karyotype, which enabled accurate integration of these regions into the genome assembly.

Linkage mapping and physical localization of the major histocompatibility complex region of the marsupial Monodelphis domestica.

We used genetic linkage mapping and fluorescence in situ hybridization (FISH) to conduct the first analysis of genic organization and chromosome localization of the major histocompatibility complex (MHC) of a marsupial, the gray, short-tailed opossum Monodelphis domestica. Family based linkage analyses of two M. domestica MHC Class I genes (UA1, UG) and three MHC Class II genes (DAB, DMA, and DMB) revealed that these genes were tightly linked and positioned in the central region of linkage group 3 (LG3). This cluster of MHC genes was physically mapped to the centromeric region of chromosome 2q by FISH using a BAC clone containing the UA1 gene. An interesting finding from the linkage analyses is that sex-specific recombination rates were virtually identical within the MHC region. This stands in stark contrast to the genome-wide situation, wherein males exhibit approximately twice as much recombination as females, and could have evolutionary implications for maintaining equality between males and females in the ability to generate haplotype diversity in this region. These analyses also showed that three non-MHC genes that flank the MHC region on human chromosome 6, myelin oligodendrocyte glycoprotein (MOG), bone morphogenetic protein 6 (BMP6), and prolactin (PRL), are split among two separate linkage groups (chromosomes) in M. domestica. Comparative analysis with eight other vertebrate species suggests strong conservation of the BMP6-PRL synteny among birds and mammals, although the BMP6-PRL-MHC-ME1 synteny is not conserved.

Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral.

BACKGROUND: Genetic testing for hereditary cancer syndromes contributes to the medical management of patients who may be at increased risk of one or more cancers. BRCA1 and BRCA2 testing for hereditary breast and ovarian cancer is one such widely used test. However, clinical testing methods with high sensitivity for deleterious mutations in these genes also detect many unclassified variants, primarily missense substitutions. METHODS: We developed an extension of the Grantham difference, called A-GVGD, to score missense substitutions against the range of variation present at their position in a multiple sequence alignment. Combining two methods, co-occurrence of unclassified variants with clearly deleterious mutations and A-GVGD, we analysed most of the missense substitutions observed in BRCA1. RESULTS: A-GVGD was able to resolve known neutral and deleterious missense substitutions into distinct sets. Additionally, eight previously unclassified BRCA1 missense substitutions observed in trans with one or more deleterious mutations, and within the cross-species range of variation observed at their position in the protein, are now classified as neutral. DISCUSSION: The methods combined here can classify as neutral about 50% of missense substitutions that have been observed with two or more clearly deleterious mutations. Furthermore, odds ratios estimated for sets of substitutions grouped by A-GVGD scores are consistent with the hypothesis that most unclassified substitutions that are within the cross-species range of variation at their position in BRCA1 are also neutral. For most of these, clinical reclassification will require integrated application of other methods such as pooled family histories, segregation analysis, or validated functional assay.

Bone mineral density, carotid artery intimal medial thickness, and the vitamin D receptor BsmI polymorphism in Mexican American women.

Low bone mineral density (BMD) is a predictor of cardiovascular mortality, suggesting that osteoporosis and cardiovascular disease may share common risk factors. We assessed the relationship between BMD and intimal medial thickening (IMT) of the common carotid artery, a marker of sub-clinical atherosclerosis, in 471 women examined as part of the San Antonio Family Osteoporosis Study, a population-based study of osteoporosis risk conducted in Mexican American families. Because of the documented role of vitamin D metabolism in bone metabolism and its possible role in cardiovascular function, we further evaluated whether allelic variation at the vitamin D receptor locus (VDR) influenced joint variation in BMD and IMT. The association of BMD with IMT depended on age, with low BMD being correlated with high IMT in older women, but with low IMT in younger women [age by IMT interaction effects significant at the spine (P = 0.042), radius ultradistal (P = 0.010), and hip (P = 0.006)]. In all women, the VDR BsmI BB genotype was associated with significantly higher forearm BMD (P = 0.005 for both radius ultradistal and midpoint), higher IMT (P = 0.05), and higher spine BMD in older women (P = 0.06), but not with hip BMD. The association of the VDR genotype with IMT was independent of its association with BMD. Although a functional consequence of the BsmI polymorphism on vitamin D metabolism has not been established, these findings support a possible biological relationship among VDR, bone metabolism, and atherosclerosis. We conclude that VDR polymorphisms may be one of multiple factors influencing the joint risk of atherosclerosis and osteoporosis.

Genome-wide linkage analysis of blood pressure in Mexican Americans.

The genetic mechanisms that control variation in blood pressure level are largely unknown. One of the first steps in understanding those mechanisms is the localization of the genes that have a significant effect on blood pressure. We performed genome scans of systolic (SBP) and diastolic blood pressure (DBP) on a population-based sample of families in the San Antonio Family Heart Study. A likelihood-based Mendelian model incorporating genotype-specific effects of sex, age, age(2), BMI, and blood pressure (SBP or DBP, as appropriate) as covariates was used to perform two-point lodscore (Z) linkage on 399 polymorphic markers. Results showed that the genotype-specific covariate effects were highly significant for both SBP and DBP. Linkage results showed that a quantitative trait locus (QTL) influencing DBP was significantly linked to D2S1790 (Z = 3.92, theta = 0.00) and showed suggestive linkage to D8S373 (Z = 1.92, theta = 0.00). A QTL influencing SBP showed suggestive linkage to D21S1440 (Z = 2.82, theta = 0.00) and D18S844 (Z = 2.09, theta = 0.11). Without the genotype-specific effects in the model, the linkage to D2S1790 was not even suggestive (Z = 1.33, theta = 0.09); thus genotype-specific modeling was crucial in detecting this linkage. A comparison with linkage studies based in other populations showed that the significant linkage to D2S1790 has been replicated at the same marker in the Quebec Family Study. The replicated significant linkage at D2S1790 may begin to establish the locations of the genes that significantly affect blood pressure across several human ethnic groups.

Genome-wide linkage analysis of pulse pressure in Mexican Americans.

Pulse pressure, a measure of aortic stiffness, is a strong predictor of cardiovascular mortality. To locate genes that affect pulse pressure, we performed genetic analysis on randomly ascertained families in the San Antonio Family Heart Study. Pulse pressure was defined as the difference between systolic and diastolic blood pressures. Likelihood methods were used to construct a model that had both single-locus and polygenic components for 46 families (1308 individuals). The single-locus component included sex-specific and genotype-specific effects of both age and body mass index. Using this model, we then performed 2-point linkage analysis in 10 families (440 individuals) that were among the largest of the 46 families and that had been genotyped for 399 polymorphic markers. The model that contained only the polygenic component and simple effects of the covariates showed pulse pressure heritability of 0.21. When the single-locus component was added, the sex-specific and genotype-specific effects of age and body mass index were highly significant (P<0.002). The full model accounted for 73% of the total variation of pulse pressure. Linkage analysis using this model with each marker revealed 4 markers with lod scores >1.9, which is the Lander-Kruglyak suggestive linkage standard. D21S1440 had a lod score of 2.78 with a recombination fraction (theta) of 0.02. D7S1799 had a lod score of 2.04 (theta=0.01), D8S1100 had a lod score of 1.98 (theta=0.08), and D18S844 had a lod score of 1.95 (theta=0.11). These results are highly correlated with results involving systolic blood pressure, indicating that pulse pressure may not be genetically distinct from systolic blood pressure.

Linkage of essential hypertension to the angiotensinogen locus in Mexican Americans.

Essential hypertension has been linked to a highly polymorphic marker at the angiotensinogen locus, and association with a polymorphism in this locus has been found in some populations. We tested the hypothesis that these same polymorphic markers are linked to essential hypertension in Mexican Americans. The data comprised all the affected relative pairs in 46 extended families chosen at random from a low-income barrio in San Antonio. Specifically, we searched for linkage by testing for excessive marker alleles shared identical by descent (IBD) among hypertensive relative pairs. When women taking oral contraceptives or hormones were excluded, the affected relative pairs shared a significant excess of alleles IBD for the highly heterozygous GT repeat polymorphism (P=.038) and were marginally significant for the M235T variant (P=.079), which has a much lower heterozygosity (0.43 versus 0.85 for the GT repeat). We also assayed plasma levels of angiotensinogen and, using likelihood methods, found no significant association (P=.43) between plasma levels of angiotensinogen and M235T genotypes. These results support the linkage of essential hypertension to the angiotensinogen locus but do not indicate a specific role for the M235T variant.

Interspecific variation at the Y-linked RPS4Y locus in hominoids: implications for phylogeny.

Within- and between-species variation in restriction endonuclease recognition sites was examined at the Y-linked RPS4Y locus of six hominoid species: human (Homo sapiens), gorilla (Gorilla gorilla), chimpanzee (Pan troglodytes), bonobo (Pan paniscus), orangutan (Pongo pygmaeus), and gibbon (Hylobates lar). RPS4Y is an expressed gene that maps to the non-recombining region of the Y chromosome. An approximately 1,490 base pair fragment of the RPS4Y gene, including all of intron 3, was amplified by PCR from DNA extracted from each of the six species. Forty-seven restriction sites were identified on the six-species composite map derived from double-digest restriction analyses of the amplified fragment. As expected, maximum parsimony analysis indicated that chimpanzee and bonobo are the two most closely related living hominoids. The same analysis suggested that the closest living relative of Homo is Gorilla, not Pan, although support for this relationship was relatively weak. These results disagree with recently published phylogenies based on analyses of mtDNA sequences (Horai et al. [1995] Proc. Natl. Acad. Sci. U.S.A. 88:7401-7404) and the Y-linked ZFY locus (Dorit et al. [1995] Science 268:1183-1185). A combined data set derived from three distinct Y-linked loci-RPS4Y, SRY, and ZFY-was also analyzed. The maximum parsimony topology for the combined data provided only weak support for a shared common ancestor for Homo and Pan subsequent to divergence from the Gorilla lineage. Taken together, the data from the Y chromosome do not provide unequivocal support for any single, dichotomously branching species tree linking Homo, Pan, and Gorilla.

Effects of a major gene for apolipoprotein A-I concentration are thyroid hormone dependent in Mexican Americans.

Apolipoprotein A-I (apoA-I) is the principal protein component of HDL cholesterol. The thyroid hormone triiodothryonine (T3) is known to be a potent mediator of expression of the apoA-I structural gene (APOA1). Using complex segregation analysis, we detected a major gene influencing plasma concentration of apoA-I and examined its interaction with T3 serum level in Mexican Americans participating in the San Antonio Family Heart Study. Strong evidence for a major locus with two alleles (A and a) determining apoA-I level was obtained when interaction with T3 was allowed. The major gene appears not to be linked to the APOA1 structural locus. Genotypes differed significantly in their relationships to T3 level. The AA and Aa genotypes showed a positive relationship with T3 level, while the rarer aa homozygote showed a strong negative relationship with T3. The relative variance in apoA-I concentration due to this major gene varied from 56% to 18%, depending on T3 level. On average, the major gene accounts for 30% of apoA-I variation, and shared-household effects account for an additional 11%. These findings suggest that thyroid hormone has an important role in the genetic control of lipoprotein metabolism.

A new DNA marker, U15557, is linked to protease inhibitor and adenylate kinase-1 in the laboratory opossum, Monodelphis domestica.

A 323-bp DNA fragment (U15557) was isolated, cloned, and sequenced after polymerase chain reaction (PCR) amplification from Monodelphis domestica genomic DNA. A HindIII restriction fragment length polymorphism was identified in this species using the U15557 PCR, fragment as a hybridization probe. DNA samples exhibited either a 6.4 kb band, a 7.2 kb band, or both bands simultaneously. Behaviour of these two variants in family studies was consistent with codominant autosomal inheritance. Linkage between this marker and the loci encoding protease inhibitor (PI) and adenylate kinase 1 (AK1) was found in M. domestica.

Major gene with sex-specific effects influences fat mass in Mexican Americans.

Increased adiposity has repeatedly been identified as a major risk factor for a variety of chronic diseases. However, the question still remains whether the amount of adipose tissue itself is genetically mediated. To address this question, a segregation analysis, using maximum likelihood techniques as implemented in the computer program Pedigree Analysis Package (PAP), was performed on fat mass (kilograms of body fat) in a large sample of extended Mexican American families residing in San Antonio, TX. The only model not rejected was a Mendelian mixed model for fat mass, incorporating genotype x sex interaction. In males the major gene accounted for 37% of the total variance compared with 43% in females. In both sexes homozygous recessive individuals have a fat mass more than double that of individuals of the other two genotypes. It was possible to reject linkage of the anonymous major gene for fat mass with several candidate loci for obesity. However, tentative evidence of linkage was detected with markers on both chromosomes 2 and 11, thereby providing hypotheses for future testing.

Developmental progression of Gpd expression from the inactive X chromosome of the Virginia opossum.

Metatherian (marsupial) mammals possess a non-random form of X-chromosome inactivation in which the paternally-derived X is always the one inactivated. To examine the progression of X-linked gene expression during metatherian development, we compared relative levels of the maternally and paternally encoded Gpd gene products in heterozygous female Virginia opossums (Didelphis virginiana) across a major portion of the developmental period. Panels of tissues obtained from fetuses, newborns, and pouch young were examined via polyacrylamide gel electrophoresis of the G6PD protein. As in adults, G6PD phenotypes in these developmental stages were highly skewed in favor of the maternal allele product, but in some tissues there was a marked increase in paternal allele expression with advancing developmental age. However, even by 42 days of post-partum development, expression of the paternal Gpd allele had not attained the adult, tissue-specific activity pattern. Our findings indicate remarkable developmental changes in the activity of the paternal allele in several tissues/organs continuing well into mid pouch-life stages and beyond. Specifically we found that 1) a substantially repressed paternal Gpd gene is present in the cells of female stage 29 fetuses and later developmental stages, 2) the activity state of the paternal Gpd gene is not fixed during early embryonic development in this species, 3) major changes in paternal Gpd expression occur in advanced developmental stages and comprise a maturation of the gene expression pattern during ontogeny, and 4) alterations of paternal Gpd allele activity during development occur in a tissue-specific manner.

X-chromosome replication patterns in adult, newborn and prenatal opossums.

Somatic cells from the opossums Monodelphis domestica and Didelphis virginiana were labelled with 5-bromodeoxyuridine (BrdU), treated with colchicine, stained with acridine orange and examined using fluorescence microscopy. BrdU-incorporated metaphase spreads from females of M. domestica at developmental stages from late bilaminar blastocysts to adults showed replication asynchrony of the two (acrocentric) X chromosomes. The long arm of one X chromosome was the latest replicating region in the entire chromosome complement and is presumed to represent transcriptional inactivation and X dosage compensation. The minute short arm of the same X, which contains a nucleolar organizer region, was earlier replicating and synchronous with the short arm of its homologue and is thus assumed to escape inactivation. BrdU-incorporated spreads from cells of fetuses, neonates and adults of D. virginiana also showed a late replicating (submetacentric) X chromosome. The pattern was different from that of M. domestica because of the different morphology and the presence of large blocks of constitutive heterochromatin in both homologues. The timing and pattern of replication of the single X in males of both species resembled the earlier replicating X in females. The array of molecular techniques now available offers the best means for investigating X-chromosome replication and activity states of X-linked genes in the earliest stages of marsupial embryogenesis.

X-linked glucose-6-phosphate dehydrogenase (G6PD) and autosomal 6-phosphogluconate dehydrogenase (6PGD) polymorphisms in baboons.

Electrophoretic polymorphisms of glucose-6-phosphate dehydrogenase (G6PD) and 6-phosphogluconate dehydrogenase (6PGD) were examined in captive colonies of five subspecies of baboons (Papio hamadryas). Phenotype frequencies and family data verified the X-linked inheritance of the G6PD polymorphism. Insufficient family data were available to confirm autosomal inheritance of the 6PGD polymorphism, but the electrophoretic patterns of variant types (putative heterozygotes) suggested the codominant expression of alleles at an autosomal locus. Implications of the G6PD polymorphism are discussed with regard to its utility as a marker system for research on X-chromosome inactivation during baboon development and for studies of clonal cell proliferation and/or cell selection during the development of atherosclerotic lesions in the baboon model.

Biochemical characteristics and subcellular localizations of rat liver neuraminidase isozymes: a paradox resolved.

A striking discrepancy in the abilities of two analytical approaches (fluorometric and electrophoretic) to detect the effect of a gene, Neu-2, on rat liver neuraminidase phenotypes led us to examine the biochemical and physical properties of the liver isozymes NEU-1 and NEU-2 that might be responsible for this difference. Cell fractionation via Percoll gradient centrifugation revealed NEU-1 activity almost exclusively in the lysosomal cell fraction, while NEU-2 was strictly cytosolic in distribution. The two isozymes were also found to differ in pH activity curves and optima (optima: 4.6-4.8 and 5.4-5.8 for NEU-1 and NEU-2, respectively) and in solubility characteristics (NEU-2 highly soluble; NEU-1 relatively insoluble but solubilized by freezing/thawing). Both isozymes were found to be freeze-thaw stable in crude, whole-cell extracts, but NEU-1 was destabilized in the enriched (partially purified) lysosomal subcellular fraction. Consideration of these properties relative to those described previously for unidentified cytosolic and membrane bound (lysosomal) rat liver neuraminidases (Tulsiani, D. R. P., and Carubelli, R., J. Biol. Chem. 245:1821, 1970) leads us to believe that NEU-2 also is destabilized by partial purification and that NEU-1 and NEU-2 have very different relative abundances within the cell. The biochemical and physical differences between NEU-1 and NEU-2 can account for the discrepant abilities of the fluorometric and electrophoretic approaches to detect the effects of Neu-2. Ways to increase the sensitivity of the fluorometric approach for quantitative assays of specific NEU-1 and NEU-2 activity are discussed.

X-linked gene expression in metatherian fibroblasts: evidence from the Gpd and Pgk-A loci of the Virginia opossum and the red-necked wallaby.

Fibroblasts cultured from ear pinna biopsies of Virginia opossums (Didelphis virginiana) and red-necked wallabies (Macropus rufogriseus) were examined electrophoretically to determine the relative expression levels of the maternally and paternally derived alleles at X-linked, enzyme-coding loci. Only the maternally derived allele was expressed at the Pgk-A locus in fibroblasts of heterozygous D. virginiana (M. rufogriseus not examined), but fibroblasts of both species exhibited evidence of paternal allele expression at the Gpd locus. Furthermore, the heterozygous G6PD phenotypes in both species were skewed in favor of the maternal gene product, as expected if the paternal allele is only partially (incompletely) expressed. For M. rufogriseus this result is contrary to a previous finding which suggested equal expression of both Gpd alleles in cultured fibroblasts of this species. The present results suggest that X-linked genes in metatherian fibroblasts are subject to the same kind of determinate, paternal allele inactivation, incomplete at some loci, described previously for X-linked genes in adult tissues and that the pattern of paternal X-linked gene expression in these cells is independent of the patterns in the tissues from which the fibroblasts are derived.

Molecular studies of marsupial X chromosomes reveal limited sequence homology of mammalian X-linked genes.

To explore the extent to which the X chromosome has been conserved during mammalian evolution, we compared six loci that are X-linked in the human genome with the corresponding genes of the North American marsupial, the Virginia opossum (Didelphis virginiana). Our analysis shows that in the opossum genome there are sequences highly homologous to those of human cDNAs for housekeeping genes, glucose-6-phosphoribosyltransferase (HPRT), phosphoglycerate kinase A (PGK1), and alpha-galactosidase A (GLA). However, ornithine transcarbamylase and blood clotting Factor IX--tissue-specific genes that are X-linked in eutherians mammals--have no highly conserved homologs in the marsupial genome. By cloning opossum G6PD and HPRT, we found that these genes are X-linked in the opossum and that homologous sequences are limited to coding regions. As all genomic fragments hybridizing with the human GLA probe show dosage effects, it is likely that the opossum counterpart is X-linked. Finally, the pattern of hybridization suggests that the autosomal pseudogenes of HPRT and PGK1 in the opossum have remained highly homologous to the human X-linked genes.

X-linked gene expression in the Virginia opossum: differences between the paternally derived Gpd and Pgk-A loci.

Expression of X-linked glucose-6-phosphate dehydrogenase (G6PD) and phosphoglycerate kinase-A (PGK-A) in the Virginia opossum (Didelphis virginiana) was studied electrophoretically in animals from natural populations and those produced through controlled laboratory crosses. Blood from most of the wild animals exhibited a common single-banded phenotype for both enzymes. Rare variant animals, regardless of sex, exhibited single-banded phenotypes different in mobility from the common mobility class of the respective enzyme. The laboratory crosses confirmed the allelic basis for the common and rare phenotypes. Transmission of PGK-A phenotypes followed the pattern of determinate (nonrandom) inactivation of the paternally derived Pgk-A allele, and transmission of G6PD also was consistent with this pattern. A survey of tissue-specific expression of G6PD phenotypes of heterozygous females revealed, in almost all tissues, three-banded patterns skewed in favor of the allele that was expressed in blood cells. Three-banded patterns were never observed in males or in putatively homozygous females. These patterns suggest simultaneous, but unequal, expression of the maternally and paternally derived Gpd alleles within individual cells (i.e., partial paternal allele expression). The absence of such partial expression was noted in a parallel survey of females heterozygous at the Pgk-A locus. Thus, it appears that Gpd and Pgk-A are X-linked in D. virginiana and subject to preferential paternal allele inactivation, but that dosage compensation may not be complete for all paternally derived X-linked genes. The data establish the similarity between the American and Australian marsupial patterns of X-linked gene regulation and, thus, support the hypothesis that this form of dosage compensation was present in the early marsupial lineage that gave rise to these modern marsupial divisions. In addition, the data provide the first documentation of the differential expression of two X-linked genes in a single marsupial species. Because of its combination of X-linked variation, high fecundity, and short generation time, D. virginiana is a unique model for pursuing questions about marsupial gene regulation that have been difficult to approach through studies of Australian species.