The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes.

The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5'UTR, where exons 1A, 1A', 1C, 1D, and a novel non-coding exon 1A" were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5'UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3' splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5'UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.

Characterization and expression analysis of KIT and MITF-M genes in llamas and their relation to white coat color.

The llama (Lama glama) is a fiber-producing species that presents a wide range of coat colors, among which white is one of the most important for the textile industry. However, there is little information about the molecular mechanisms that control the white phenotype in this species. In domestic mammals, a white coat is usually produced by mutations in the KIT proto-oncogene receptor tyrosine kinase (KIT) and microphthalmia-associated transcription factor (MITF) genes. In this work we have sequenced and described the coding regions of KIT and MITF-M, the melanocyte-specific isoform, and the two transcriptional variants MITF-M(-) and MITF-M(+). Moreover, we studied the expression of these genes in the skin of white and colored llamas. Although no variants were revealed to be associated with white coat color, significant differences between phenotypes were observed in the expression levels of KIT and MITF-M. Interestingly, white llamas expressed less MITF-M(+) than did colored ones, which is consistent with a consequent reduction in the synthesis of melanin. Even though our results indicate that downregulation of KIT and MITF-M expression is involved in white phenotype production in llamas, the causative gene of white coat color remains unknown.

Molecular characterization of the llama FGF5 gene and identification of putative loss of function mutations.

Llama, the most numerous domestic camelid in Argentina, has good fiber-production ability. Although a few genes related to other productive traits have been characterized, the molecular genetic basis of fiber growth control in camelids is still poorly understood. Fibroblast growth factor 5 (FGF5) is a secreted signaling protein that controls hair growth in humans and other mammals. Mutations in the FGF5 gene have been associated with long-hair phenotypes in several species. Here, we sequenced the llama FGF5 gene, which consists of three exons encoding 813 bp. cDNA analysis from hair follicles revealed the expression of two FGF5 alternative spliced transcripts, in one of which exon 2 is absent. DNA variation analysis showed four polymorphisms in the coding region: a synonymous SNP (c.210A>G), a single base deletion (c.348delA), a 12-bp insertion (c.351_352insCATATAACATAG) and a non-sense mutation (c.499C>T). The deletion was always found together with the insertion forming a haplotype and producing a putative truncated protein of 123 amino acids. The c.499C>T mutation also leads to a premature stop codon at position 168. In both cases, critical functional domains of FGF5, including one heparin binding site, are lost. All animals analyzed were homozygous for one of the deleterious mutations or compound heterozygous for both (i.e. c.348delA, c.351_352insCATATAACATAG/c.499T). Sequencing of guanaco samples showed that the FGF5 gene encodes a full-length 270-amino acid protein. These results suggest that FGF5 is likely functional in short-haired wild species and non-functional in the domestic fiber-producing species, the llama.

Genetic diversity and conservation status of managed vicuña (Vicugna vicugna) populations in Argentina.

The vicuña (Vicugna vicugna) was indiscriminately hunted for more than 400 years and, by the end of 1960s, it was seriously endangered. At that time, a captive breeding program was initiated in Argentina by the National Institute of Agricultural Technology (INTA) with the aim of preserving the species. Nowadays, vicuñas are managed in captivity and in the wild to obtain their valuable fiber. The current genetic status of Argentinean vicuña populations is virtually unknown. Using mitochondrial DNA and microsatellite markers, we assessed levels of genetic diversity of vicuña populations managed in the wild and compared it with a captive population from INTA. Furthermore, we examined levels of genetic structure and evidence for historical bottlenecks. Overall, all populations revealed high genetic variability with no signs of inbreeding. Levels of genetic diversity between captive and wild populations were not significantly different, although the captive population showed the lowest estimates of allelic richness, number of mitochondrial haplotypes, and haplotype diversity. Significant genetic differentiation at microsatellite markers was found between free-living populations from Jujuy and Catamarca provinces. Moreover, microsatellite data also revealed genetic structure within the Catamarca management area. Genetic signatures of past bottlenecks were detected in wild populations by the Garza Williamson test. Results from this study are discussed in relation to the conservation and management of the species.

South American camelid illegal traffic detection by means of molecular markers.

South American camelids comprise the wild species guanaco and vicuña and their respective domestic relatives llama and alpaca. The aim of the present study was to determine by DNA analysis to which of these species belong a herd of camelids confiscated from a llama breeder but alleged to be alpacas by the prosecution, and to evaluate the usefulness of mitochondrial and autosomal DNA markers to solve judicial cases involving camelid taxa. Cytochrome b and cytochrome oxidase I mitochondrial genes and 7 STR were analyzed in 25 confiscated samples. Mitochondrial results were inconclusive because 18 of the sequestered samples presented haplotypes that corresponded to the guanaco haplogroup and the remaining seven belonged to a vicuña linage. Microsatellite data of casework samples and llama reference samples revealed different genetic profiles by the presence of private alleles at two microsatellites suggesting that the confiscated animals could be alpaca, or at least alpaca hybrids instead of pure llama.

Identification of camelid specific residues in mitochondrial ATP synthase subunits.

ATP synthase is an enzyme involved in oxidative phosphorylation from prokaryotic to eukaryotic cells. In mammals it comprises at least 16 subunits from which the mitochondrial encoded ATP6 and ATP8 are essential. Mitochondrial genes variations have been suggested to allow rapid human and animal adaptation to new climates and dietary conditions (Mishmar et al. 2003). Camelidae taxa are uniquely adapted to extremely hot and dry climates of African-Asian territories and to cold and hypoxic environments of the South American Andean region. We sequenced and analyzed ATP6 and ATP8 genes in all camelid species. Based on the available structural data and evolutionary conservation of the deduced proteins we identified features proper of the group. In Old World camels the ATP8, important in the assembly of the F0 complex, showed a number of positively charged residues higher than in the other aligned species. In ATP6 we found the camelid specific substitutions Q47H and I106V that occur in sites highly conserved in other species. We speculate that these changes may have functional importance.

2006 GEP-ISFG collaborative exercise on mtDNA: reflections about interpretation, artefacts, and DNA mixtures.

We report the results of the seventh edition of the GEP-ISFG mitochondrial DNA (mtDNA) collaborative exercise. The samples submitted to the participant laboratories were blood stains from a maternity case and simulated forensic samples, including a case of mixture. The success rate for the blood stains was moderate ( approximately 77%); even though four inexperienced laboratories concentrated about one-third of the total errors. A similar success was obtained for the analysis of mixed samples (78.8% for a hair-saliva mixture and 69.2% for a saliva-saliva mixture). Two laboratories also dissected the haplotypes contributing to the saliva-saliva mixture. Most of the errors were due to reading problems and misinterpretation of electropherograms, demonstrating once more that the lack of a solid devised experimental approach is the main cause of error in mtDNA testing.

Mutation rates at Y chromosome specific microsatellites.

A collaborative work was carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG) to estimate Y-STR mutation rates. Seventeen Y chromosome STR loci (DYS19, DYS385, DYS389I and II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS460, DYS461, DYS635 [GATA C4], GATA H4, and GATA A10) were analyzed in a sample of 3,026 father/son pairs. Among 27,029 allele transfers, 54 mutations were observed, with an overall mutation rate across the 17 loci of 1.998 x 10(-3) (95% CI, 1.501 x 10(-3) to 2.606 x 10(-3)). With just one exception, all of the mutations were single-step, and they were observed only once per gametogenesis. Repeat gains were more frequent than losses, longer alleles were found to be more mutable, and the mutation rate seemed to increase with the father's age. Hum Mutat 26(6), 520-528, 2005. (c) 2005 Wiley-Liss, Inc.

Genetic diversity and differentiation of guanaco populations from Argentina inferred from microsatellite data.

Genotype data from 14 microsatellite markers were used to assess the genetic diversity and differentiation of four guanaco populations from Argentine Patagonia. These animals were recently captured in the wild and maintained in semi-captivity for fibre production. Considerable genetic diversity in these populations was suggested by the finding of a total of 162 alleles, an average mean number of alleles per locus ranging from 6.50 to 8.19, and H(e) values ranging from 0.66 to 0.74. Assessment of population differentiation showed moderate but significant values of F(ST)=0.071 (P=0.000) and R(ST)=0.083 (P=0.000). An amova test showed that the genetic variation among populations was 5.6% while within populations it was 94.4%. A number of 6.6 migrants per generation may support these results. Unambiguous individual assignment to original populations was obtained for the Pilcaniyeu, Las Heras and La Esperanza populations. The erroneous assignment of 18.75% Rio Mayo individuals to the Las Heras population can be explained by the low genetic differentiation found between these two populations. Thirty-nine of 56 loci per population combinations were in Hardy--Weinberg disequilibrium because of guanaco heterozygote deficiency, which may be explained by population subdivision. The high level of genetic diversity of the guanacos analysed here indicates that the Patagonian guanaco constitutes an important genetic resource for conservation or economic utilization programmes.

Phylogenetic relationships among Robertsonian karyomorphs of Graomys griseoflavus (Rodentia, Muridae) by mitochondrial cytochrome b DNA sequencing.

Graomys griseoflavus (Waterhouse 1837) is a phyllotine murid rodent with a Robertsonian autosomal polymorphism, having been described 2n = 42, 41, 38, 37, 36, 35 and 34 karyomorphs, and proposed a chromosomal divergence pathway accounted by four sequential Robertsonian fusions. Sequences of a fragment (422 bp long) of the cytochrome b (cyt b) mitochondrial gene and its 5' flanking region (tRNA Glu) were obtained for 19 Graomys griseoflavus from different karyomorphs to infer phylogenetic relationships by using maximum parsimony. Outgroups considered for this analysis were the phyllotine rodents Phyllotis xanthopygus and Eligmodontia typus cyt b sequences. Three trees were produced showing the 2n = 38-34 karyomorphs grouped in a single clade while the 2n = 42-41 animals formed a different one. This is in agreement with a hypothesis of a single origin for 2n = 38-34 Robertsonian karyomorphs from the ancestral 2n = 42.

Variability of the F13B locus in South American populations.

A population study for the F13B locus was carried out in a total of 396 South American individuals. The analysis comprised new data from 5 Amerindian populations, existing data from 3 Amerindian populations, and I urban sample from La Plata, Argentina. In both pooled Amerindian and La Plata samples, 6 alleles were found. The individual Amerindian samples showed a lower number of alleles, changes in modal alleles, and restricted variability. Interpopulation comparisons revealed significant differences among samples from distinct geographical regions. Differences among the groups were also corroborated by the F(ST) statistic. Data support the hypothesis that genetic drift and gene flow influence Amerindian differentiation.

Allele frequencies of six STR loci in Argentine populations.

Allele frequencies of six short tandem repeat (STR) loci were determined in a Caucasian urban sample of La Plata city and three Amerindian sample populations of Argentina. Allele frequencies showed differences between urbans and Amerindians, and among Amerindians as well. The degree of genetic differentiation of subpopulations was mainly due to the Amerindian contribution. Mapuche, Mocovi, and pooled Amerindian populations showed little evidence of HW disequilibrium, and association of alleles. In the urban sample, there is no evidence of population substructuring. Forensic probabilities of exclusion and matching showed high differences between the population groups. Finally, La Plata sample did not show differences with Caucasians from other geographic regions.

Characterization of ancestral and derived Y-chromosome haplotypes of New World native populations.

We analyze the allelic polymorphisms in seven Y-specific microsatellite loci and a Y-specific alphoid system with 27 variants (alphah I-XXVII), in a total of 89 Y chromosomes carrying the DYS199T allele and belonging to populations representing Amerindian and Na-Dene linguistic groups. Since there are no indications of recurrence for the DYS199C-->T transition, it is assumed that all DYS199T haplotypes derive from a single individual in whom the C-->T mutation occurred for the first time. We identified both the ancestral founder haplotype, 0A, of the DYS199T lineage and seven derived haplogroups diverging from the ancestral one by one to seven mutational steps. The 0A haplotype (5.7% of Native American chromosomes) had the following constitution: DYS199T, alphah II, DYS19/13, DYS389a/10, DYS389b/27, DYS390/24, DYS391/10, DYS392/14, and DYS393/13 (microsatellite alleles are indicated as number of repeats). We analyzed the Y-specific microsatellite mutation rate in 1,743 father-son transmissions, and we pooled our data with data in the literature, to obtain an average mutation rate of.0012. We estimated that the 0A haplotype has an average age of 22,770 years (minimum 13,500 years, maximum 58,700 years). Since the DYS199T allele is found with high frequency in Native American chromosomes, we propose that 0A is one of the most prevalent founder paternal lineages of New World aborigines.

Molecular analysis of chromosomal polymorphism in the South American cricetid, Graomys griseoflavus.

Graomys griseoflavus is a South American phyllotine rodent widespread in Argentina that shows a high frequency of Robertsonian fusions (RFs). DNA restriction with EcoRI produced a 250-bp repeated family (EG250) specific for the genus. Southern hybridization and sequencing analysis indicate that the EG250 family is heterogeneous, comprising at least two subfamilies. In situ hybridized EG250 probe showed a centromere location in almost all chromosomes. In all karyomorphs C-banding was negative, but restriction enzyme banding (Re-banding) with Alul and Mbol showed centromeric blocks in the autosomes that will generate Robertsonian fusions. Thus, we found three groups of chromosomes: (a) EG250 and Re-banding negative; (b) EG250 positive and Re-banding negative; and (c) EG250 and Re-banding positive. We consider that group (b) is more the result of chromatin condensation state than that of the frequency of recognition sites for the enzymes used. Restriction enzyme blocks would appear in regions with heterochromatic EG250 subfamilies, while lack of banding would be due to decondensed EG250 subfamilies becoming an easier target for chromosomal restriction. It is suggested that heterochromatic EG250 DNA provides a favourable molecular environment for Robertsonian fusion occurrence.

An assessment of the relationships among species of Camelidae by satellite DNA comparisons.

Tandem satellite arrays and interspersed repetitive DNA components of the New World camelids guanaco, llama, alpaca, and vicuña and the Old World bactrian camel have been identified and compared. Southern hybridizations, using camel restriction fragments as probes, indicated that satellite DNAs in all camelids examined have been conserved since the last common ancestor about 5-10 MY ago. The hybridization profiles, however, varied from totally identical (MspI-sat) to highly differentiated (PstI-sat and EcoRI-sat) between Old and New World species. Repetitive DNA patterns specific of South American camelids were identified by most of the vicuña and guanaco probes and (a) llama and guanaco have undifferentiable patterns, supporting the view that the former is a domesticated form of the latter; (b) vicuña patterns were species-specific and in agreement with its position in a separate taxonomic unit; (c) the presence in alpaca of BamHI, TaqI and EcoRI patterns that are intermediate between those of the species above, suggested that the origin of the alpaca may be found in a cross-breed between the guanaco and vicuña.

Restriction site patterns in the ribosomal DNA of camelidae.

The restriction map of rDNA from South American camelids and the Bactrian camel was analyzed by digestion of high-molecular-weight DNA with endonucleases EcoRI,BamHI and the two combined followed by Southern blot hybridization with probes for the 18S and 28S rDNA sequences. We scored a total of 17 restriction sites, six of which were mapped conserved in all the species. The other eleven corresponded to spacer regions and revealed variations between these taxa. The study showed that the two groups differ in the length of the internal transcribed spacer. Also they showed the existence of two regions of fast evolution on the opposite termini of the external spacer. A restriction site present at low frequency in the non-transcribed spacer of guanaco and llama was the only difference encountered within the South American group.

Evolution of zinc finger-Y and zinc finger-X genes in oryzomyne-akodontine rodents (Cricetidae).

Zinc finger-Y (Zfy) and zinc finger-X (Zfx) genes were analyzed by Southern blotting in male and female specimens of 10 species belonging to the oryzomyne-akodontine stock of Cricetidae rodents. DNA fragments were used as characters to construct a parsimony tree of the genes. Zfx and Zfy trees in general coincide with the evolutionary history of the taxa. Both trees show Oryzomys longicaudatus genes as the outgroup whereas Akodon xanthorrhinus genes are also distant from those of the other species. Oxymycterus rufus and Bolomys obscurus share related sequences, while genes from the other six Akodon species form a group of their own. It was found that 9 out of the 10 species analyzed show Zfy amplification in a range varying from 2 to 24 copies and with a pattern that is clade specific. The estimation of the average changes per character strongly suggests that Zfy has evolved more rapidly than Zfx; our estimates of the rate of nucleotide substitution are 4.6 times higher for Zfy than for Zfx.

Effect of the metal chelating agent o-phenanthroline on the DNA and chromosome damage induced by bleomycin in Chinese hamster ovary cells.

A 30-min pulse treatment with bleomycin to Chinese hamster ovary cells in culture produces DNA degradation and chromosomal aberrations in a dose-dependent manner. Bleomycin also induces a long-lasting effect on the cell cycle producing a lengthening of two or more cycles after the treatment. The presence of o-phenanthroline, which chelates metal ions, totally inhibits DNA cleavage and the appearance of chromosome aberrations while partially correcting the lengthening of the cell cycle. These findings suggest that an important cellular target for bleomycin is the DNA. Chromosomal aberrations are a secondary effect resulting from DNA cleavage. On the other hand, the increase in the duration of the cell cycle is probably induced by DNA degradation and, perhaps, by damage to other cellular structures.

Ag-NOR staining and in situ hybridization of rDNA in the chromosomes of the South American camelids.

The location and frequency of Ag-stained NORs and sites of rDNA hybridization were studied in the chromosomes of the South American camelids. In the four camelids these regions occur distally on chromosomes 18, 21, and 27 and the smallest biarmed elements. Quantitative analysis of NOR distribution showed variations between both cells and species. In llama, guanaco and alpaca the NORs number averaged 6 per cell, this being higher than in vicuña where the average was 3. Relative frequencies of NOR-bearing chromosomes in the four camelids were similar. Yet, in vicuña virtual absence of NOR sites on one of the smallest biarmed pairs was observed. The rDNA sites assessed in llama and vicuña by in situ hybridization with cloned 18S DNA were coincident with the NOR locations and with the frequencies characteristics for each species. Moreover, varying the exposure time of the autoradiographs, labeling patterns specific for each camelid were observed. Grain counts on individual chromosomes indicated that under our conditions one month exposure is enough to demonstrate all the rDNA sites available in the complement of llama. Conversely, at least two months are necessary to show the total sites existing in vicuña. Most probably this finding reflects the presence of variations in the amount of copies of the ribosomal genes per chromosome.